18.06.2013 Views

Boundary Lyer Theory

Boundary Lyer Theory

Boundary Lyer Theory

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

FRO XXII. Tho incornprcs~iblo turbulent boundary layer 11. Tho calculaLion of two-climcnsional turbulent boundary lnycrs (iX I<br />

tribution U(x). Since the valuc of the shape factor, N(x) = 1, has already been<br />

assigned, the only quantity that we need to calculate is the local Reynolds number<br />

determined by eqn. (22.11a). Since m = const, we can contract the two terms on<br />

the Icft-hand side and solve the problcm by performing two integrations, one each<br />

for Rz(x) and R3(x). In ~ont~racted form these are<br />

The relations and numerical values to be used for the exponents i, e and n, as<br />

well M for the modified lrinematie viscosity v' are listed in Table 22.2, separately<br />

for the momentum and for the energy method. The constant of integration is<br />

DifTercntiat.ing eqn. (22.17) wit11 rc~pcctr to x and taking into account thc contmctions<br />

dcfincd in Tablrs 22.1 antl 22.2 wc can tlcmonstrat~ eonsist~ency with eqn. (22.11a).<br />

In the case of the moment,um m&hod, eqn. (22.17) becomes identical with eqn. (22.16)<br />

if wc put i = n, R = Rz antl IZ = Ez.<br />

Comparing the nurnerical data of Table 22.2, we find far-reaching agreement.<br />

In spite of considerable difl'erences in the assumptions for the shape factor and for<br />

the shearing stress at the wall we discover that thc two explicit cquat.ions for the<br />

calculation of the momentum thickness are equivalent. The following specific numeri-<br />

cal values can be recommcndctl:<br />

The cnergy mcthod is discussed below.<br />

Analogy with self-similar eolutions (ener~y method): Self-similar solutions in<br />

boundary-layer thcory are generally described as equilibrium /lows when they occur<br />

in turbulent motion. They are charnctcrized by the fact that the velocity profiles<br />

u/U nt varying positions x become similar for certain velocity distributions U(x) of<br />

the ontcr flow. 'J'his means that, the shape factor H(x) remains constant with z, that<br />

is that dlfldx = 0. Figure 22.7 implics that all quantities which depend on x in<br />

gcnrml n~ust bcrome constant for such equilibrium boundary layers.<br />

We now suI)st,itutc in cqn. (22.10) the expression for cn from cqn. (22.813) and<br />

note that the integration with respect to x can be performed in closed form with b =<br />

const and p' = const. The Reynolds nnmbcr formcd with the energy thickness is<br />

thus given by<br />

. .<br />

I hc numerical vn.lucs for 6, r, n, and v' are to Ip sclcctcd in accordance with the<br />

relations in Table 22.2 and l'ig. 22.7s.<br />

I<br />

In the spceial rase of separation-prone flows for which Hs < I1 < Ilk, we find,<br />

for example, that. 1 4- b$ = 1.004 and I + b, = 1.152. Tliesc two values differ by<br />

about. 5%. Such a disrrcpancy can t)c tlisregardcd in view of the i~ncert~ainties in-<br />

Iirrent in s~wh appro xi mat,^ mctzhods. In othcr words, this signifies that it is possiblc<br />

to perform calculations using numerical values based on t,hc flat-plat,c ani~logy. 'rhc<br />

quantity P' = PIPQ) that appears in eqn. (22.19) also depends only wcaltly on the shape<br />

factor; by way of approximation, we let it be P' = I. 'rhns, the calculabion of ~(IIc<br />

cncrgy t,hickness with the aid of eqn. (22.10) can be bascci on tdlc followitig numerical<br />

valucs :<br />

R += R3: h = 0.152; e == 2.3 ; n = 3.3; r' = 78 v; P' =- 1 . (22.20)<br />

With thrse assumption, eqn. (22.19) transforms into cqn. (22.17) bearing in mint1 that<br />

i = b, R = R3, and /C = "3, as cxpcctcd.<br />

Refcrence [114] diows t,hat eqn. (22.1 la) suffices by it,sclf to solvc the prnbleln<br />

when the energy mct,llocI is used. 13y contrast, when t.hc nionlent.urn mctliotl is usctl,<br />

t.he coupling between cqns. (22.1 la) and (22.11b) cannot hc disregardcdt. 'Hie Iatkr<br />

Icads us t,o the trivial result that Rz = IIz3R3 in view of thc dcfnit.ion Ifzn -- 02/fi3.<br />

At. t.hr icvcl of approximation - -<br />

considcrcd so far, thc momcntum mctliotl t.urns out<br />

t,o be itlcntical witJr t.11~ cncrgy mct.liot1. Ncvcrthclcsrr, tho two proc:rtl~~~.c,s tlilli,~,<br />

csscntially from one another in that tlic momentum mcthod employs tho t,wo basic<br />

cquat,ions (22.11~~) and (22.1 lb), whcrras t,hc cncrgy mct,hotl gct,s by with cqn. (22.1 la)<br />

nlonc. As far as tlic devclop~nent of furthcr approximat.ions in the for111 of simple<br />

inkgrals is concerned, we have cxhaust,ctl t.hc pot,ent.inl inlierent in thc rnomrnt.um<br />

mct,horl. In the energy mcthod, cqn. (22.1 1b) is used to tlcrivc n forninla for t,hc shape<br />

fnclor by closed-form int.rgrnl.ion, as wc arc about, t.o how.<br />

Integration method due to E. Truckenbrodt: 13. Trucltcnbrodt [I 11, 1141 worlrcd<br />

out an npproximnt,~ mc:l,hotl for 1.11~ cxplicil inttcyrnl.iol~ of t,llc cqi~n.Iio~is of t.111.1~111v1it~<br />

boundary layers which serves to obtain the boundary-layer t,l~iclrness (cnergy tliickness)<br />

as well as the (modified) shape fact,or. The first vcrsion of the method [Ill] has<br />

proved to be practicable for calculat,ions in engineering applications. It was, therefore,<br />

tliougl~t useful to modify it in the light of more rccent discoveries. Employing<br />

the classification introduced above, we tlescribc this as an energy method. Tlrc method<br />

can be used for two-dimensional as wellasforaxiallysymmct~ricflows,c/. Sco. XXlX cl.<br />

Tlic method is based on cqn. (22.11a) which is used to calculate the Iteynolds<br />

number formed with the energy thickness. Hence we put R = R3, m = 2 and (IS =<br />

@a = (21~) /? Rab in accordance with Table 22.1. If we assume 1,hat 6 = const and<br />

p(r) is known, we can integrate eqn. (22.11a) with rcspcct to x and obtain<br />

1:<br />

1 B4x)<br />

{R3(4I1+* = -;; [C1(2)1C where E3(x) = E3(x1) + JP' Un dx. (22.21)<br />

The numerical values in cqn. (22.20) are valid up to P'(.T) = P(:c)/Pw. Ilowcvcr, an<br />

inspection of Fig. 22.7b shows that /I' docs not clcviatc much from the valuc 1.0,<br />

and we may calculate with /I' = 1.0 by way of approximation. In this case, eqn.<br />

(22.21) transforms into cqn. (22.19). Thus, if no great dernands of accuracy are made<br />

on the values of the Reynolds number, we obtain<br />

--<br />

z<br />

2%<br />

\ lllltbl<br />

t By way of amplificstion, we mcnt,ion that the result, derived in Tal~le 6 of [I141 in also valid<br />

when a -= const nnd b = ronut,.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!