MA Alexander Grunewald
Programme & Anleitungen
- Handbuch zu PLUTO worauf belt basiert
- Handbuch zu Visit, welches die vtk-Dateien anzeigen kann
- Handbuch zu gnuplot, einem Plot-Klassiker
- Cheat-Sheet zu Emacs, einem Editor-Klassiker
- Cheat-Sheet zu tmux, einem Terminal-"Vervielfacher"
Setups
(belt installiert auf saturn und auf Laptop)
lin_T-Grad
init.c
relevante Zeilen für die initialen Druck- und Dichte-Profile:
if(jl <= 2) data->Vc[PRS][kl][jl-1][il] = P_0;
data->Vc[PRS][kl][jl][il] = (1.0 - xi2) / (1.0 + xi2) * data->Vc[PRS][kl][jl-1][il];
data->Vc[RHO][kl][jl][il] = data->Vc[PRS][kl][jl][il] / (R_UniversalGasConstant / MolarMass(data, kl, jl, il) * Tgas);
Lothar: Die erste davon kann entfallen, die dritte bleibt so, und in der zweiten soll der Druckwert direkt zugewiesen werden: data->Vc[PRS][kl][jl][il] = P_0*pow(1-w*(x2-Ymin),q/w); aber es müssen \(q\) und \(w\) zuvor noch deklariert werden (double q,w; weiter oben) und vor der DOM_LOOP aus den Parametern berechnet werden.
Notizen & Ergebnisse
- Start belt: 1. Esc+&, 2. ./belt
- g zum aktualisieren
\[ \mathrm{p(y)}=\mathrm{p_0}(1-\mathrm{wy)}^\frac{q}{w}, \mathrm{q}=\frac{\mathrm g\mathrm m}{\mathrm{k_B}\mathrm{T_0}}=7,7\cdot 10^{-8}\mathrm{m^{-1}}, \ \mathrm{Einheit} \ \mathrm{q}=\frac{\mathrm{mkg}}{\mathrm{s^2}\mathrm{J}}=\frac{1}{\mathrm{m}}, \mathrm{w}=\frac{\mathrm{u}}{\mathrm{T_0}}=1,29\cdot 10^{-8}\mathrm{m^{-1}}, \mathrm{Einheit} \ \mathrm{w}=\frac{1}{\mathrm{m}}=\mathrm{m^{-1}}, \frac{\mathrm{q}}{\mathrm{w}}=5,98 \] \[ \mathrm{\rho}=\frac{\mathrm{m}\mathrm{p_0}}{\mathrm{k_B}\mathrm{T_0}}(1-\mathrm{wy})^{\frac{\mathrm{q}}{\mathrm{w}}-1} \]