MA Alexander Grunewald: Difference between revisions

From Arbeitsgruppe Kuiper
Jump to navigation Jump to search
(graph)
No edit summary
 
(27 intermediate revisions by 2 users not shown)
Line 1: Line 1:
=Programme & Anleitungen=  
=Programme & Anleitungen=  


Line 10: Line 11:


(''belt'' installiert auf saturn und auf Laptop)
(''belt'' installiert auf saturn und auf Laptop)
[[MA Alexander Grunewald/Code-Fragmente|Code-Fragmente]] & [[MA Alexander Grunewald/Lothars Plots|Plots]] (Lothar)


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


=Notizen & Ergebnisse=
=Notizen & Ergebnisse=
Line 50: Line 42:
\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>
\]</nowiki>
==Geschwindigkeitsprofile==
==Geschwindigkeitsprofile==
[[File:MAAG_Numerisches und exaktes Profil.png|thumb|center|numerische und exakte Berechnung der 1D-Navier-Stokes-Gleichung]]
[[File:MAAG Vergleich Auflösung num und ex.png|thumb|center|Vergleich der numerischen und exakten Berechnung bei verschiedenen Auflösungen]]


[[File:MAAG_Auflösungstest.png|thumb|center|Auflösungstest für 300 und 600 Zellen ohne Störung]]
[[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 Vergleich Auflösung num und ex.png|thumb|Vergleich der numerischen und exakten Berechnung bei verschiedenen Auflösungen]][[File:MAAG Vergrößerung des hinteren Teils der verschiedenen Auflösungen.png|thumb|Vergrößerung des hinteren Teils der verschiedenen Auflösungen]]
[[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.001.png|thumb|center|Temperatur 1e+6-1e+5, Amplitude 1-0.001, s=3,87e-4 R_Sonne]]
[[File:MAAG 300 T6-5 Ampl1-0.001(zoom).png|thumb|center|Temperatur 1e+6-1e+5, Amplitude 1-0.001, s=3,87e-4 R_Sonne]]
[[File:MAAG 300 T7-4 Ampl1-0.001.png|thumb|center|Temperatur 1e+7-1e+4, Amplitude 1-0.001, s=3,88e-4 R_Sonne]]
[[File:MAAG 300 T7-4 Ampl1-0.001(zoom).png|thumb|center|Temperatur 1e+7-1e+4, Amplitude 1-0.001, s=3,88e-4 R_Sonne]]
[[File:MAAG 300 T6-5 Ampl0.01-0.0001.png|thumb|center|Temperatur 1e+6-1e+5, Amplitude 0.01-0.0001, s=3,88e-4 R_Sonne]]
[[File:MAAG 300 T7-4 Ampl0.01-0.0001.png|thumb|center|Temperatur 1e+7-1e+4, Amplitude 0.01-0.0001, s=3,88e-4 R_Sonne]]

Latest revision as of 14:10, 25 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.001, s=3,87e-4 R_Sonne
Temperatur 1e+6-1e+5, Amplitude 1-0.001, s=3,87e-4 R_Sonne
Temperatur 1e+7-1e+4, Amplitude 1-0.001, s=3,88e-4 R_Sonne
Temperatur 1e+7-1e+4, Amplitude 1-0.001, s=3,88e-4 R_Sonne
Temperatur 1e+6-1e+5, Amplitude 0.01-0.0001, s=3,88e-4 R_Sonne
Temperatur 1e+7-1e+4, Amplitude 0.01-0.0001, s=3,88e-4 R_Sonne