BA Fynn Wawrzyniak

Einströmgeschwindigkeit Konstant

Damit die Auflösung des Grid's für verschiedene Machzahlen nicht geändert werden muss, wird die Einströmgeschwindikeit als Konstant angesehen und die Schallgeschwindigkeit variiert. Für die Schallgeschwindigkeit gilt:

\begin{align} c_{\infty} = \frac{v_{\infty}}{\mathcal{M}}=\sqrt{\gamma \cdot \frac{p_{\infty}}{\rho_{\infty}}} \end{align} \( \) Wenn die Einströmgeschwindigkeit \(v_{\infty}\) und die Dichte \( \rho_{\infty} \) als Konstanten gegeben sind und die Machzahl \( \mathcal{M} \) sowie der adiabatische Index \( \gamma \) als unabhängige Parameter einstellbar sind, lässt sich der Druck in folgender Form ausdrücken:

\begin{align} p_{\infty} = \frac{\rho_{\infty}}{\gamma} \biggl( \frac{v_{\infty}}{\mathcal{M}} \biggr)^2 \end{align}

Oder auch:

\begin{align} p_{\infty} = \frac{\rho_{\infty}}{\gamma} c_{\infty}^2 \end{align}

Probleme

Für Machzahlen \( \mathcal{M} > 10 \) tritt kurz nach Beginn der Simulation folgender Fehler auf:

-> ! ConsToPrim(): p(E) < 0 (-1.01e-07), @step = 552 (stage = 1); [i,j = 2, 7], [x1,x2 = 0.000001, 0.172788]

Ein Beispiel "RunFolder", in welchem dieser Fehler auftritt, ist in Nextcloud hochgeladen. In diesem sollte die "init.c" mit meiner Rechnung, die Log-Files und "pluto.ini" zu finden sein.