Feladat: folytonos kezdeti feltétel diszkretizálása

Next: Feladat: FTCS módszer a Up: Diszkretizáció Previous: Diszkretizáció Tartalomjegyzék

Feladat: folytonos kezdeti feltétel diszkretizálása

Ebben a feladatban egy a

$\begin{displaymath} \rho(x)=\left\{ \begin{array}{ll} 1+0.3\sin [0.1 \pi (x-10... ...0\le x\le20$}\\ 1 & \mbox{egy\'ebk\'ent} \end{array}\right. \end{displaymath}$

(7)

képlettel definiált szinuszhullám alakú kezdeti feltételt generálunk az

tartományban. Indítsa el a VAC grafikus interfészét a

  weblink &

utasítással. Kattintson a RUN VACINI gombra, majd írja be a felsõ SET PROBLEM mögötti mezõbe, hogy

  hullam

és kattintson a SET PROBLEM gombra. Ezzel a kiválasztotta a ,,hullám'' problémához tartozó állományokat, így a kezdeti feltétel paramétereit tartalmazó inipar/hullam11 állományt. (A 11 arra utal, hogy a probléma 1 dimenziós és a sebességnek is csak 1 komponense különbözik 0-tól).

Mielõtt elindítaná a kezdeti feltételt létrehozó VACINI programnak a kontinuitási egyenletre konfigurált és lefordított vacini_rho11 változatát, tekintse meg a futtatás paramétereit definiáló inipar/hullam11 paraméterállományt, azaz kattintson az EDIT gombra, és nézze át a beállításokat. A rács méretét és a rácstávolságot az alábbi paraméterek határozzák meg:

grid
50                      -racspontok szama
1                       -minimum koordinata
50                      -maximum koordinata

Mint látható, 50 rácspontot használunk és $\Delta x=1$ . A sûrûségprofilt a

uniform
1.                      -konstans suruseg
perturbvar
1  0.3  10 20           -rho_, drho, jmin, jmax
0.31415 0.              -k=2*pi/lamba, fazis

utasításokkal állítottuk be, míg a

sebességet, és a diffúziós együtthatót az

eqpar
2.  0.                  -v sebesseg, D diffuzios allando

utasítás adja meg.

Miután áttekintette a paramétereket, kattintson a EXIT EDITOR (kilépés a szövegszerkeszõbõl) gombra, és a RUN (futtatás) gombra való kattintással futassa le a VACINI programot.

Válassza a RESULTS (eredmények) menüpontot a bal oldalon. Írja be, hogy 65 a DATA HEAD (adat eleje) mögötti mezõbe, kattintson a DATA HEAD gombra, és nézze meg a diszkrét koordináták és változók számszerû értékét! Indítsa el az IDL vizualizációs programot a START IDL gombra való kattintással. Az IDL program számára a most megjelent APPLICATION címû ablakban adhat utasításokat. Írja ide be, hogy

  .r getpict
  .r plotfunc

majd a függvény nevének írja be, hogy rho. A megjelenõ IDL 0 feliratú ablakban láthatja a kezdeti feltételt. Lépjen ki az IDL programból az

  exit

utasítás beírásával.

Next: Feladat: FTCS módszer a Up: Diszkretizáció Previous: Diszkretizáció Tartalomjegyzék

Gabor Toth 2000-09-04