BA Fynn Wawrzyniak

Einströmgeschwindigkeit konstant

Damit die Auflösung des Grids für verschiedene Machzahlen nicht geändert werden muss, wird die Einströmgeschwindikeit als konstant gewählt 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}

Des Weiteren gilt für den Druck folgender Zusammenhang mit der Temperatur:

\begin{align} p_{\infty} = \frac{\rho_{\infty} T_{\infty} R}{\mu} \end{align}

Hier ist \( R \) die Gaskonstante und \( \mu \) die molare Masse. Aus diesem beiden Gleichungen folgt diese Proportionalität, wenn für verschiedene Mach Zahlen die Einströmgeschwindigkeit konstant bleiben soll:

\begin{align} T \propto \frac{1}{\mathcal{M}^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.