MA Alexander Grunewald: Difference between revisions

From Arbeitsgruppe Kuiper
Jump to navigation Jump to search
m (→‎Setups: Lothars Plots)
m (→‎lin_T-Grad: init.c -> Code-Fragmente)
Line 14: Line 14:


==''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=

Revision as of 09:40, 3 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} \]

Auflösung 600 mit und ohne Dämpfung

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