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.

ccnt,rred nt, (m -1- 1, n). 1'110 two t,xptwsions are t.hertwpon cornt)ined in such a way<br />

thnt ternls of ortlcr Aq2 are elitnina1.rtl. The corrcsponclir~g difference qnot,icnks can<br />

I,c given the, form (index we 1 I omitted) :<br />

a',, 1<br />

{ I I I 2 T I I 1 7 I] I 0 ( A . 1 ) (9 69)<br />

il,l 2 /I 11,'<br />

where<br />

p1 - 1<br />

l2 (I -+- K), rPz - L ~,<br />

.-- 2 r1 J'~, I>,, -= I .<br />

Eq~tn.tions (9.69) nncl (9.70) rcrlr~cc t.o t,hc st.nntlartl form for cent,ral difl'crenccs when<br />

K = I.<br />

For the (-tlcrivnt,ivrs in cqnnt.inn (9.64) a simple bnclrwnrtl tlini.rrner formule is<br />

used<br />

y - Fsr 11. n -- Fnr, n<br />

E -<br />

-I om. (9.7 I )<br />

At --<br />

The 1nrgc.r I,rtlnc.at,ion error which appcn.t,s here is balancctl by t,hc it.crat.ivc scltclne<br />

proposc(1 for solving thr tlilkrcncc c-qt~at,ion. 'l'llc non-lincn.r t,ernis in r(lnnt,ion (9.04)<br />

Imvc to be rcplaced by lincarizccl diflkrencc quoticnt.s. Tho tcrlns fFIl and FFg may<br />

serve as exntnples and thcy are writ,t,cn as<br />

l'hc lincnrizrtl Iini1,r-clifl'ercnw qnotionts given nhove are su11st~it.ut~erl into the<br />

tlifl'crentinl ty~tnt.ion (9.G4) nnd Lhe result is multiplied througlt by A E to give n<br />

tlilYcrcnce equa.t,iotl. 'J'his is writden ns follows<br />

i. The n~ctl~otl of finite dilTwcnces 191<br />

111 cqun1,ions (9.75), 6 and 0 nre cvnlunlctl nt (111. 1- I), mtl ot~ly the vn.rir~l,lcx wit.11<br />

sttp(wcripll i ntt 11~1cInlc~1 through st~cccssivt: it,crntions. To s~~xxl-III) t,Itv il~r~~:t,l,ion<br />

proces.s t,he tcrrns (/S)t can be licJ)t constnrlt. (equnl to t . 1 ~ &luc nt t.hc prrvious<br />

shtion) unt,il initial convergence is nchicved.<br />

Method of nolution: Equations (9.74) rcprc~ont~ n ~ oof t N-1 si~ntlll.nncwus r~.lgc:-<br />

I~rnic equntionrr for the unlrnown k;ntl,n (n = 2, 3, . . ., N). At, cnch levcl IL t.l~rcc<br />

unknown quantities nppenr, namely Fnajl, .-I, Fm.kl,n and Fniit, ,,+I, but sincc<br />

F,+I,~ and Fm+l,~ nre known from the bonndnry conilitions, 1.11~ totnl nurnbcr of<br />

cquntions equals the nunlber of unlrnowns. The set of nlgclmic cqunt,ions rnn be<br />

writt.cn in so-callod tllrcc-tlingonnl matrix form. MnLriccs of Illis ttypcwhcro oK-tlir~go~tnl<br />

elements vnnish outdc the three-tlingonal band can bc inverted bg n sirnplc: and<br />

direct nlothod well suit.cd for digital con~putcrs.l'o end tlriseqna.tion (0.74) is rcwrit,toll<br />

in "stantlard form" (subscripts (m -1 I) ornittcd)<br />

Thc botlnclary conditions arc<br />

F1 = 0 nntl PN = 1, (9.70)<br />

wllere IL = I tlenotcs tbc wall and ?t = N thc edge of the bountlary Inycr. J1, is asst~tn-<br />

etl now t,llat a solution existst in tllc form<br />

The boundary condit,ion F1 = 0 nnd t,he rcyuirernent t,llnt rquation (9.77) sl~ould<br />

rcrnnin valid indcpcntlcnt,ly of the sl,cp size /Iq leads to<br />

A direct, colisrquencc of rqnntion (9.77) is that<br />

When the preceding expressions are substitutml int20 oqn. (!).741,), 1.11~ following<br />

relohion is obt.aincd<br />

By tncnns of equalion (9.81) and the condition (9.78), it becomes possil~lc t.o cotnpot,~

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

Saved successfully!

Ooh no, something went wrong!