Feladat: lökéshullám modellezése

Indítsa el a VAC grafikus interfészét weblink & utasítással, kattintson a RUN VACINI gombra, írja be a felsõ SET PROBLEM mögötti mezõbe, hogy

  lokes

és kattintson a SET PROBLEM gombra. Tekintse meg a inipar/lokes11 paraméterállományt. A négyszöghullám alakú nyomásprofilt a

grid
200             -racspontok szama
  1.            -minimum koordinata
200.            -maximum koordinata
uniform
1. 0. 1.        -rho v p
setvar
61 140          -jmin, jmax
3 100. T        -p_,   ertek, utolso

utasításokkal állítottuk be. Mint láthatja 200 rácspontot használunk. Miután áttekintette a paramétereket futassa le a vacini_hd11 programot.

Kattintson a RUN VAC menüpontra, és tekintse meg a par/lokes11 paraméterállományt. A hanghullámoknál is bevált TVD-MUSCL diszkretizációt használjuk az MC (más néven Woodward) meredekség limiterrel. Mindkét oldalon nyitott határfeltételeket adtunk meg. Újdonság a

 &paramlist
        slowsteps=10
        courantpar=0.9
 /

beállítás. A slowsteps=10 azt jelenti, hogy az elsõ 10 lépésben az idõlépést fokozatosan növeljük a végsõ $C=0.9$ maximális Courant számmal dinamikusan megadott értékhez. Erre azért van szükség, mert a kezdeti lépésekben a megoldás rendkívül gyorsan változik. Futassa le a vac_hd11 programot ezekkel a paraméterekkel, és nézze meg az eredményt az IDL-lel. Milyen gyorsan megy a lökéshullám? Mekkora a sûrûség kompressziója? Mennyit változik az entrópia a lökéshullámon illetve a kontakt diszkontinuitáson keresztül? Hány rácspont ábrázolja a lökéshullámot és a kontakt diszkontinuitást? Visszaverõdik-e a lökéshullám a zéró gradiensû határról?

Az érdekesség kedvéért próbálja meg a feladatot megoldani az elsõrendû upwind (typefull='tvdmu1')és a MacCormack (typefull='mc') módszerekkel is.