MA Alexander Grunewald: Difference between revisions
Jump to navigation
Jump to search
(→Geschwindigkeitsprofile: Präfix MAAG_) |
Alexander.G (talk | contribs) (Auflösungstest) |
||
| Line 1: | Line 1: | ||
[[File:Auflösungstest.png|thumb|Auflösungstest für 300 und 600 Zellen ohne Störung]] | |||
=Programme & Anleitungen= | |||
= Setups = | *[https://plutocode.ph.unito.it/userguide.pdf Handbuch] zu [https://plutocode.ph.unito.it PLUTO] worauf ''belt'' basiert | ||
*[https://visit-sphinx-github-user-manual.readthedocs.io/en/develop/index.html Handbuch] zu [https://visit-dav.github.io/visit-website/index.html Visit], welches die vtk-Dateien anzeigen kann | |||
*[http://www.gnuplot.info/docs_6.0/gnuplot6.html Handbuch] zu [http://gnuplot.info gnuplot], einem Plot-Klassiker | |||
*[https://www.gnu.org/software/emacs/refcards/pdf/refcard.pdf Cheat-Sheet] zu [https://de.wikipedia.org/wiki/Emacs Emacs], einem Editor-Klassiker | |||
*[https://tmux.info/docs/cheatsheet Cheat-Sheet] zu [https://github.com/tmux/tmux/wiki tmux], einem Terminal-"Vervielfacher" | |||
=Setups = | |||
(''belt'' installiert auf saturn und auf Laptop) | (''belt'' installiert auf saturn und auf Laptop) | ||
== ''lin_T-Grad'' == | ==''lin_T-Grad''== | ||
=== init.c === | ===init.c=== | ||
relevante Zeilen für die initialen Druck- und Dichte-Profile: | relevante Zeilen für die initialen Druck- und Dichte-Profile: | ||
| Line 24: | Line 26: | ||
Lothar: Die erste davon kann entfallen, die dritte bleibt so, und in der zweiten soll der Druckwert direkt zugewiesen werden: <code>data->Vc[PRS][kl][jl][il] = P_0*pow(1-w*(x2-Ymin),q/w); </code> aber es müssen \(q\) und \(w\) zuvor noch [https://de.wikibooks.org/wiki/C-Programmierung:_Variablen_und_Konstanten deklariert] werden (<code>double q,w;</code> weiter oben) und vor der <code>DOM_LOOP</code> aus den Parametern ''berechnet'' werden. | Lothar: Die erste davon kann entfallen, die dritte bleibt so, und in der zweiten soll der Druckwert direkt zugewiesen werden: <code>data->Vc[PRS][kl][jl][il] = P_0*pow(1-w*(x2-Ymin),q/w); </code> aber es müssen \(q\) und \(w\) zuvor noch [https://de.wikibooks.org/wiki/C-Programmierung:_Variablen_und_Konstanten deklariert] werden (<code>double q,w;</code> weiter oben) und vor der <code>DOM_LOOP</code> aus den Parametern ''berechnet'' werden. | ||
= Notizen & Ergebnisse = | =Notizen & Ergebnisse= | ||
* Start ''belt'': 1. Esc+&, 2. ./belt | *Start ''belt'': 1. Esc+&, 2. ./belt | ||
* g zum aktualisieren | *g zum aktualisieren | ||
== Analytische Lösung für 1D-Navier-Stokes-Gleichung == | ==Analytische Lösung für 1D-Navier-Stokes-Gleichung== | ||
\begin{align} | <nowiki>\begin{align} | ||
p(y)=p_{\mathrm{0}}(1-wy)^{\frac{q}{w}}, | p(y)=p_{\mathrm{0}}(1-wy)^{\frac{q}{w}}, | ||
\end{align} | \end{align} | ||
| Line 45: | Line 47: | ||
\begin{align} | \begin{align} | ||
\frac{q}{w}=5,98 | \frac{q}{w}=5,98 | ||
\end{align} | \end{align}</nowiki> | ||
\[ | <nowiki>\[ | ||
\rho=\frac{mp_0}{k_{\mathrm{B}}T_0}(1-wy)^{\frac{q}{w}-1} | \rho=\frac{mp_0}{k_{\mathrm{B}}T_0}(1-wy)^{\frac{q}{w}-1} | ||
\] | \]</nowiki> | ||
== Geschwindigkeitsprofile == | ==Geschwindigkeitsprofile== | ||
[[File:MAAG_Numerisches und exaktes Profil.png|thumb|center|numerische und exakte Berechnung der 1D-Navier-Stokes-Gleichung]] | [[File:MAAG_Numerisches und exaktes Profil.png|thumb|center|numerische und exakte Berechnung der 1D-Navier-Stokes-Gleichung]] | ||
Revision as of 09:06, 10 February 2026
File:Auflösungstest.png
Auflösungstest für 300 und 600 Zellen ohne Störung
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
\begin{align} p(y)=p_{\mathrm{0}}(1-wy)^{\frac{q}{w}}, \end{align} \begin{align} q&=\frac{gm}{k_{\mathrm{B}}T_{0}}=7,7\cdot 10^{-8}\mathrm{m^{-1}},\\ \mathrm{Einheit} \ q&=\frac{\mathrm{mkg}}{\mathrm{s^2}\mathrm{J}}=\frac{1}{\mathrm{m}}, \\ u&=\frac{T_{0}-T_1}{Y_1-Y_0}\\ \mathrm{Einheit} \ u&=\frac{\mathrm{K}}{\mathrm{m}}\\ w&=\frac{u}{T_0}=1,29\cdot 10^{-8}\mathrm{m^{-1}}\\ \mathrm{Einheit}\ w&=\frac{1}{\mathrm{m}}=\mathrm{m^{-1}} \end{align} \begin{align} \frac{q}{w}=5,98 \end{align}
\[ \rho=\frac{mp_0}{k_{\mathrm{B}}T_0}(1-wy)^{\frac{q}{w}-1} \]
