BA Fynn Wawrzyniak: Difference between revisions
(new) |
m (edit) |
||
| Line 1: | Line 1: | ||
= Einströmgeschwindigkeit Konstant = | = Einströmgeschwindigkeit Konstant = | ||
Damit die Auflösung des Grid's für verschiedene Machzahlen nicht geändert werden muss, wird Einströmgeschwindikeit als Konstant angesehen und die Schallgeschwindigkeit variiert. Für die Schallgeschwindigkeit gilt: | 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} | \begin{align} | ||
| Line 11: | Line 11: | ||
\begin{align} | \begin{align} | ||
p_{\infty} = \frac{\rho_{\infty}}{\gamma} \biggl( \frac{v_{\infty}}{\mathcal{M}} \biggr)^2 | 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} | \end{align} | ||
Revision as of 21:46, 25 April 2024
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.