MA Alexander Grunewald: Difference between revisions

From Arbeitsgruppe Kuiper
Jump to navigation Jump to search
(→‎Geschwindigkeitsprofile: Präfix MAAG_, zentriert)
 
(41 intermediate revisions by 2 users not shown)
Line 1: Line 1:
= Programme & Anleitungen =


* [https://plutocode.ph.unito.it/userguide.pdf Handbuch] zu [https://plutocode.ph.unito.it PLUTO] worauf ''belt'' basiert
=Programme & Anleitungen=
* [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 =
*[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'' ==
[[MA Alexander Grunewald/Code-Fragmente|Code-Fragmente]] & [[MA Alexander Grunewald/Lothars Plots|Plots]] (Lothar)
 
=== init.c ===
 
relevante Zeilen für die initialen Druck- und Dichte-Profile:
<syntaxhighlight lang="C">
        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);
</syntaxhighlight>


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.
==''lin_T-Grad''==


= 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 37:
\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 ==
 
[[File:MAAG_Numerisches Profil ohne störung.png|thumb|center|Geschwindigkeitsprofil der numerischen Lösung]]
 
==Geschwindigkeitsprofile==
[[File:MAAG Vergleich Auflösung num und ex.png|thumb|center|Vergleich der numerischen und exakten Berechnung bei verschiedenen Auflösungen]]


[[File:MAAG_Exaktes Profil ohne Störung.png|thumb|center|Geschwindigkeitsprofil der exakten Lösung]]
[[File:MAAG Vergrößerung des hinteren Teils der verschiedenen Auflösungen.png|thumb|center|Vergrößerung des hinteren Teils der verschiedenen Auflösungen]][[File:MAAG verschiedene Auflösungen mit und ohne Dämpfung.png|thumb|center|verschiedene Auflösungen mit und ohne Dämpfung]]
[[File:MAAG zoom verschiedene Auflösungen mit und ohne Dämpfung.png|thumb|center|Vergrößerung verschiedenen Auflösungen mit und ohne Dämpfung]]
[[File:MAAG exakt 600 mit und ohne dämpfung.png|thumb|center|Auflösung 600 mit und ohne Dämpfung]]
[[File:MAAG exakt 600 mit und ohne dämpfung2.png|thumb|center|Auflösung 600 mit und ohne Dämpfung (Dämpfungszone mit 0,02 R_Sonne, Dämpfungszeit 4,9359*10^-7 a, Strecke 46,6*10^8m=6,7 R_Sonne)]]
[[File:MAAG Auflösungen ohne stoerung.png|thumb|center|parabolisch mit Dämpfung, ohne störung]]
[[File:MAAG Auflösungen mit stoerung ampl1.png|thumb|center|parabolisch, mit Dämpfung, mit Störung (Mode=1, Amplitude=1, Zeit=1*10^-6)]]
[[File:MAAG Auflösungen mit stoerung ampl10.png|thumb|center|parabolisch, mit Dämpfung, mit Störung (Mode=1, Amplitude=1, Zeit=1*10^-6)]]
[[File:MAAG 300 ampl1 dauerstörung1e-5 zelle75 t6-5.png|thumb|center|Temperatur 1e+6-1e+5, Amplitude 1, Dauer Störung 1e-5, (s_Anfang=7,33e-6,s_Mitte=3,87e-4, s_Ende=7,7e-4) R_Sonne]]
[[File:MAAG 300 ampl1 dauerstörung1e-5 zelle75 t7-4.png|thumb|center|Temperatur 1e+7-1e+4, Amplitude 1, Dauer Störung 1e-5, (s_Anfang=7,73e-6,s_Mitte=3,88e-4, s_Ende=7,69e-4) R_Sonne]]
[[File:MAAG 300 T6-5 Ampl1-0.0001konst2.png|thumb|center|Temperatur 1e+6-1e+5, Amplitude 1-0.0001, s=3,87e-4 R_Sonne]]
[[File:MAAG 300 T7-4 Ampl1-0.0001konst2.png|thumb|center|Temperatur 1e+7-1e+4, Amplitude 1-0.0001, s=3,87e-4 R_Sonne]]

Latest revision as of 12:14, 28 March 2026

Programme & Anleitungen

Setups

(belt installiert auf saturn und auf Laptop)

Code-Fragmente & Plots (Lothar)

lin_T-Grad

Notizen & Ergebnisse

  • Start belt: 1. Esc+&, 2. ./belt
  • g zum aktualisieren


Analytische Lösung für 1D-Navier-Stokes-Gleichung

\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} \]


Geschwindigkeitsprofile

Vergleich der numerischen und exakten Berechnung bei verschiedenen Auflösungen
Vergrößerung des hinteren Teils der verschiedenen Auflösungen
verschiedene Auflösungen mit und ohne Dämpfung
Vergrößerung verschiedenen Auflösungen mit und ohne Dämpfung
Auflösung 600 mit und ohne Dämpfung
Auflösung 600 mit und ohne Dämpfung (Dämpfungszone mit 0,02 R_Sonne, Dämpfungszeit 4,9359*10^-7 a, Strecke 46,6*10^8m=6,7 R_Sonne)
parabolisch mit Dämpfung, ohne störung
parabolisch, mit Dämpfung, mit Störung (Mode=1, Amplitude=1, Zeit=1*10^-6)
parabolisch, mit Dämpfung, mit Störung (Mode=1, Amplitude=1, Zeit=1*10^-6)
Temperatur 1e+6-1e+5, Amplitude 1, Dauer Störung 1e-5, (s_Anfang=7,33e-6,s_Mitte=3,87e-4, s_Ende=7,7e-4) R_Sonne
Temperatur 1e+7-1e+4, Amplitude 1, Dauer Störung 1e-5, (s_Anfang=7,73e-6,s_Mitte=3,88e-4, s_Ende=7,69e-4) R_Sonne
Temperatur 1e+6-1e+5, Amplitude 1-0.0001, s=3,87e-4 R_Sonne
Temperatur 1e+7-1e+4, Amplitude 1-0.0001, s=3,87e-4 R_Sonne
Retrieved from "https://wiki.uni-due.de/agk/index.php?title=MA_Alexander_Grunewald&oldid=1172"