On solutions of the integral form of the second-order boundary-layer equations

On solutions of the integral form of the second-order boundary-layer equations

Int. J. Engn~ Sci. Vol. 8. pp. 15 1-l 56. Pergamon Press 1970. Printed in Great Britain ON SOLUTIONS OF THE INTEGRAL FORM OF THE SECOND-ORDER ...

351KB Sizes 1 Downloads 38 Views

Int. J. Engn~

Sci. Vol. 8. pp. 15 1-l 56.

Pergamon

Press

1970.

Printed

in Great

Britain

ON SOLUTIONS OF THE INTEGRAL FORM OF THE SECOND-ORDER BOUNDARY-LAYER EQUATIONS MICHAEL Virginia

Polytechnic

Institute,

J. WERLEt Engineering

and

R. THOMAS

Mechanics

Department,

DAVIS Blacksburg,

Va. 24061,

U.S.A

Abstract- Approximate solutions to the second-order boundary-layer equations for incompressible flow are considered. A simplified form of the integral momentum equation for two-dimensional and axisymmetric flow is presented. Solutions to the problem of shear flow over a flat plate are obtained for five different polynomial approximations to the second-order boundary-layer velocity profile. It is found that the results obtained depend rather erratically on the polynomial degree employed.

INTRODUCTION THE

PURPOSE of this note is to bring to light results of an investigation concerned with the approximate solution of the integral form of the second-order boundary-layer equations for laminar incompressible flow. In particular, attention will be centered on the sensitivity of the method to the order of the polynomial used to represent the firstand second-order boundary-layer velocity profiles. The insensitivity of the first-order solutions to the polynomial order (see Schlichting, [ 11) certainly gives impetus to the present approach because it indicates that it is more important to satisfy the integral momentum equation than to be concerned with the details of the velocity profile employed. It is not immediately apparent that this latter fact is directly applicable to the second-order solutions and it is felt that an investigation of this point offers one means of testing the validity of the polynomial approximations.

GOVERNING

EQUATIONS

The most general form of the integrated boundary-layer equations for two-dimensional or axisymmetric incompressible how was presented in reference [2]. Subsequent study (see reference [3]) has shown that considerable simplification is achieved through a simple redefinition of terms involved. For the sake of completeness, the general problem is briefly restated here. Thus, the first-order (or classical) integrated boundarylayer equation may be stated as+:

*Former address: Applied Aerodynamics Division, U.S. Naval Ordnance Laboratory, Silver Spring. Md. 20910, U.S.A. Sunless, otherwise noted, the nomenclature used here is identical to that of reference [2] with one added simplification, i.e. the inviscid flow properties evaluated at the surface will be written without their arguments. Thus, for example, U, is in reality, C/, (s, 0). 151

and. as was shown

in reference

[3]. the second-order

equation

can be written

as

In equation (2). the second-order contributions due to vorticity interaction can be identified by the coefficient B’(O), the displacement speed effect by iJz, the longitudinal curvature by K, and the transverse curvature by u& The origins of the new boundarylayer thicknesses, & and 6,. introduced in equation (2) are best understood in terms of the classical concepts of the displacement and momentum thicknesses. Simply stated. the displacement thickness is that amount by which a body must be thickened in order to correctly calculate the inviscid flow field about that shape. The appropriate asymptotic expansion for large Reynolds numbers is then given as; 6-

c6, +E’&+.

. ..

(3)

where 1 E-S

(4)

P

and R, is the flow’s characteristic Reynolds number. The related to the second-order displacement thickness i& through

new thickness the relation

8, is then

(5) These

two thicknesses

are defined

by the integrals

(6)

(7)

In an analogous

manner,

one can define a momentum

o- Eo,+&&+. where 0, is related

to H, through

thickness, ...

0 such that: (8)

the relation

(9) and

(10)

Second-order

boundary-layer

equations

153

with

TEST

PROBLEM

The test problem considered here is that of a semi-infinite flat plate placed in a stream which has a uniform (but symmetric) shear. The exact solution for such was first given by Murray [4] and was later discussed by Van Dyke [S]. Devan and Oberai [6] found that they could match these exact solutions within a few percent through use of fourth-order polynomial approximations to both the first- and second-order velocity profiles in conjunction with the appropriate integral momentum equations. In the present work the integral equations will be solved for six different polynomial approximations to the second-order velocity profile in order to test the validity of the concept being applied. For a free-stream velocity distribution of the form U(--3c,_v) = I $_oJ\>I

(12)

where o is constant. it has been shown by Van Dyke[5] that the solution of the firstorder boundary-layer equations is identical to that forw = 0, and that there is no secondorder correction to the inviscid flow. The vorticity interaction produces the only second-order viscous effect and the integral equation reduces to simply de, -72,.-w, ds where 6, is given by (see equation 8x-J;:

(13)

( I I ))

(1~,--wN)(I-2rr,)dN-t2~J;;

(I-u,)NdN,

( 14)

and

For the first-order profile, the Karman so that withr) = N/l, = N/C,.s”‘.

Pohlhausen

fourth-order

profile

is employed

11, = 27) - 2$ - n4. Solution of equation ( I) using equations momentum equation now becomes

Before

d

m

-[I ds

0

(16)

(6) and ( 16) gives C, = 5.84. The integral

(l-++4$+2r)4)(~~X-~N)dN

considering

solutions

of equation

(I 7) via polynomial

approximations

to

153

hl. J. WERLE

Us. it is fruitful to first consider be employed here are

and R. I’. DAVIS

an appropriate

Np(S,

set of boundary

conditions.

The ones to

0)

=

0,

0)

=

-+y

12)

=

cfJLJ.

(20)

&Js,

1,)

=

w.

(21)

L4,,v,,(s,

12)

=

0,

(22)

u,,v,v(s,

u,(s,

( IX)

( 19)

and

where I, represents the value of N at which the second-order quantities take on theii limit values (note C, is a constant resulting from the first-order solutions, C, = 1.752). With respect to these conditions. equation ( 18) represents the no-slip condition, equation (19) is required to satisfy the boundary-layer equation at the wall, equation (20) represents the outer matching condition and equations (21) and (22) are required to satisfy the boundary-layer equations at N = I,?. While admittedly, one could generate an infinite number of boundary conditions by successive differentiation of the secondorder boundary-layer equation: those presented above represent a natural set. That second-order profile which satisfies all of the conditions of equations (I 8-22), plus the integral momentum equation, will be the second-order counterpart of the first-order Karman Pohlhausen solutions. Since both 14, and uq appear within the integral of equation (I 7). it will be necessary to investigate the relationship between I, and i2. The possibility that 1, > I, has already been eliminated and it remains to consider the two distinct cases of I, = I, or I, < I,. This investigation begins by considering I, = I2 and employs polynomial profiles in 5 = N/l, = 7, taking the general second-order velocity profile to be of the form

(23) With equation

(23), the integral

momentum

equation

(equation

(I 7)) reduces

to

(24) where I = I,~f’(~)[l-4~+4~3-2~“] The second-order

d&

(25)

wall shear is given by

tThese conditions are based on the assumption that the first-order velocity goes to its limit value at or before the second-order velocity does so. This assumption appears a logical one since the first-order quantities in fact appear as coefficients in the second-order equations and therefore serve to drive the second-order velocity to its limit value.

Second-order

boundary-layer

Equation (24) has been solved for five separate in equation (25) and the results are given in Table conditions which have been used.

Boundary

Shear stress - ~~,,,/w

(IX). (20) (181, (20). (21)

-i)

l-5)’ Ic,GY+t)+t515(1-5)’

(18-21) (18-22)

[C,(31+l)+f515(1-i)” ioriginally

The exact solution

conditions

(18)

Cl5 c‘,i(l C,C(

polynomial representations off(c) 1, along with a list of the boundary

I

Table Profile-J‘(c)

due to Devan

and Oberia,

was given by Murray, T

1.55

equations

1.337 2.209 3.510 3.022t 3.734

(6).

141, as

pUz= 3.126 w,

(27)

which was predicted quite accurately using the fourth-order profile above. The rather erratic dependence of the solutions on the profile order certainly does little to build one’s confidence in this approach. As a final test of the polynomial approach a solution has also been obtained taking I, > I,, and a velocity distribution of the form (note that 77no longer equals 5)

Fitting a fourth-order (18-22)) yields

polynomial

f(5)

in 5 to all of the boundary

=85(1-L)“.

conditions

(equations

(29)

and the shear stress at the wall as rq(. --

2.414~.

(30)

Certainly this is a poor approximation to the exact result (equation (27)). While the above results serve as a test only of the vorticity interaction solution using the polynomial approximations, the conclusions reached above can be extended to some extent with little difficulty. As discussed in reference [7], the linearity of the second-order equations produces a breakdown of the general problem into four separate contributions: the displacement speed effect. the vorticity interaction effect, the effects due to longitudinal curvature, and those effects due to the presence of transverse curvature. As was shown in reference [3], the latter three of these problems all have velocity profiles of the same general form. Since the polynomial approach failed for the vorticity problem, and both the longitudinal and transverse curvature profiles are of the same form, (but with different slopes as N + m) it is reasonable to assume equally poor results will be obtained for the latter as well.

I”/,

kl..l.

and

wL-Kl.t~

I<. I.

I).\\

IS

It hxwld appear then that the polynomial approach needs little further investigation. the rather accurate results given by Devan and Oberai [6] must. unfortunately. he viewed as being rather t’ortuitou% in light of the present study. ;tnd

K El’ EK t-:N(‘t:S

[ I ] 11.S<‘HLICHTING. 121 M. J. WER1.E 131M. J. WERLE,

Borrrrc/tr,?-Ltrv~,~ 7‘lr~or~. p. 205. Pergamon Pre\\ and R. T. DAVIS. /nt../. E/Is/?~.SC,~. 4. 423 (I 9hh). U.S. Ntr~:n/Ord. Lab.R~,J. No. TR 6%‘; I 196X).

[4] J. D. b1URRAY.J. 151 M. VAN DYKE,J. [6] [7]

(I’)551

FlrritlM~~~h. II. 309( I‘961). FluidMcvl~. 14.481 ( IYhZ).

I.. DEVAN and M. M. ORERAI.A/AA.//r/2. 1X38( IYh4). hl.V.AN I)YKF../.FI/A;~M~‘~./~.I~. Ihl (IYh2).

H&um6Les auteurs Ctudient des solutions approchkea des iquations du second ordre de la couche limite d’un Ccoulement de fluide incompressible et proposent une forme simplifiee de I’inttgrale de la quantitC de mouvement applicable aux Ccoulements :I deux dimensions ou de rCvolution. Ils obtiennent des solutions du problkme de la couche de choc sur une plaque plane pour cinq approximations diffkrentes de la rCpartition des vitesses dans la couche limite effectuees g I’aide d’un polynbme et trouvent que les rCsultats obtenuh dtpendent assez irrCgul&ement du degrC du polynBme employ6.

Zusammenfassung-

Ungefahre Liiungen fiir die Grenzschichtgleichungen zweiter Ordnung fiir inkompressible StrGmung werden untersucht. Eine vereinfachte Form der integralen Momentgleichung fiir zweidimensionale und achsensymmetrische Striimung wird vorgelegt. Liisungen fiir das Problem der ScherstrC mung iiber eine flache Platte werden fiir fiinf verschiedene polynomische Anngherungen an das Grenzschichtgeschwindigkeitsprofil zweiter Ordnung erhalten. Es wird gefunden, dass die erhaltenen Resultate errdtisch vom polynomischen Grad, der verwendet wird, abhgngen. Summarioconsiderano le aoluzioni approssimate delle equazioni di strata limite di seconda grandezra per un flusso incomprimibile. Si presenta una forma semplificata dell’equazione di moment0 integrato per un flusso bidimensionale ed assisimmetrico. Si ottengono soluzioni al problema di flusso su un piatto liscio per cinque approssimazioni polinomiali differenti al profilo di velocith dello strato limite di second’ordine. Si scoprc the i risultati ottenuti dipendono in manieraalquanto erratica dal grado polinomiale impiegato. A6CTpPKT-PPaCCMaTpuBaroTCrr IlOp!$LIKa n0

KOflWECTBJ

‘GlRaVH .IJfR

IIJISI HPA2KHMat?MOTO

n0

CKOPOCTA

OIlIM60~HOfi

MnororneHa.

i-PaHH’IHOrO

MHHMYto

PelIIeHMs

npenCTaBiWHa

LlBkGKKeHFiHLIJIR ABj’Xpa3MepHOrO

TtYEHHEO 06TeKafOLWMy

KOHTYPa

CKOpee

npu6nEixeHHble TeSCHBII.

ypaBHeHH8

yIlpOIIZHH%l

11 OCeCAMMeTPH’lHOrO

IIJlOCKyIO IIJIrlTy rIpEi IllTA CJlOIl

3aBNCF,MOCTb

“0

BTOpOrO

IIOpJlAKa.

nOJ,yYeHHblX

TplHHSHOMY

+OpMa

PaSJIA’IHbIX

Te’ieHPiB.

OT

YTO

CTeIleHA

BTOpOrO

ypaBH’ZHEiSl

nOJIj”faIoTCR

MHO,-O’WeHHblX

YCTaHaBJlHBaeTCSl, pe3yJIbTaTOB

CJIOH)

AHTerpaJIbHOrO

pWleHrtH

npFf6mmemnx CJIeAyeT

C’IATaTb

yrIOTpe6neHHOrO