15. Metoda konÄnih elementov in reÅ¡evanje BDE in PDE v Matlabu
15. Metoda konÄnih elementov in reÅ¡evanje BDE in PDE v Matlabu
15. Metoda konÄnih elementov in reÅ¡evanje BDE in PDE v Matlabu
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Zgled 3 za bvp4c<br />
Naslednje enačbe opisujejo pretok v dolgem navpičnem kanalu, kjer je tekoč<strong>in</strong>a vbrizgana skozi eno izmed<br />
stranic:<br />
f ′′′ [<br />
− R (f ′ ) 2 − ff ′′] + RA = 0,<br />
h ′′ + Rfh ′ + 1 = 0,<br />
θ ′′ + P fθ ′ = 0,<br />
pri čemer je R Raynoldsovo število <strong>in</strong> P = 0.7R. Ker parameter A spada med neznanke, imamo 8 robnih<br />
pogojev:<br />
f(0) = f ′ (0) = 0, f(1) = 1, f ′ (1) = 0,<br />
h(0) = h(1) = 0,<br />
θ(0) = 0, θ(1) = 1.<br />
Problem bi radi rešili pri R = 10000, vendar nimamo nobenih dobrih začetnih približkov. Zato uporabimo<br />
metodo zveznega nadaljevanja.<br />
približek za R = 1000, tega pa za R = 10000.<br />
Najprej problem rešimo za R = 100, potem to uporabimo kot začetni<br />
Bor Plestenjak - Numerična analiza 2005/06