Lothars Scratch - Revision history

Lothar.brendel: /* mögliche natürliche Einheiten */ Umrechnungsfaktor

2026-07-31T12:49:16Z

mögliche natürliche Einheiten: Umrechnungsfaktor

← Older revision		Revision as of 14:49, 31 July 2026
Line 17:		Line 17:

	In diesen Einheiten wäre also \(P_\infty = 1/\gamma\) und \(GM=2\).		In diesen Einheiten wäre also \(P_\infty = 1/\gamma\) und \(GM=2\).

			Einheit Sonnenmassen/Jahr in analysis.out: ''M''<sub>⊙</sub> = 1.9891·10<sup>33</sup>g, Year = 31557600s ⇒ Umrechungsfaktor = 6.303·10<sup>25</sup>

	= Alexanders Setups =		= Alexanders Setups =

2026-07-30T15:35:51Z

Notizen Elias: Korrektur nat. E.

← Older revision		Revision as of 17:35, 30 July 2026
Line 4:		Line 4:
	* barotrop: \(\rho=\rho(p)\), gültig fürs ideale Gas bei isothermen (\(\rho\propto p\)) und isentropen (\(\rho\propto p^{1/\gamma}\)) Zustandsänderungen		* barotrop: \(\rho=\rho(p)\), gültig fürs ideale Gas bei isothermen (\(\rho\propto p\)) und isentropen (\(\rho\propto p^{1/\gamma}\)) Zustandsänderungen

	~~mögliche natürliche Einheiten~~		== Sims ==

			=== mögliche natürliche Einheiten ===

	\[		\[
	\begin{align}		\begin{align}
	\text{Länge } r_\mathrm{B} &= \frac{GM}{2c^2_{\text{s},\infty}}\\		\text{Länge } r_\mathrm{B} &= \frac{GM}{2c^2_{\text{s},\infty}}\\
	\text{Geschwindigkeit} &= c_{\text{s},\infty} = \sqrt{\gamma\frac{RT}{M_\text{mol}}}\\		\text{Geschwindigkeit} &= c_{\text{s},\infty} = \sqrt{\gamma \frac{P_\infty}{\rho_\infty}}=\sqrt{\gamma\frac{RT_\infty}{M_\text{mol}}}\\
	\text{Dichte} &= \rho_\infty		\text{Dichte} &= \rho_\infty
	\end{align}		\end{align}
	\]		\]

			In diesen Einheiten wäre also \(P_\infty = 1/\gamma\) und \(GM=2\).

	= Alexanders Setups =		= Alexanders Setups =

2026-07-30T09:14:39Z

Notizen Elias: nat. E.

@@ Line 3: / Line 3: @@
 * \(S\) = ''spezifische'' Entropie, d.h. Entropie pro Masse
 * barotrop: \(\rho=\rho(p)\), gültig fürs ideale Gas bei isothermen (\(\rho\propto p\)) und isentropen (\(\rho\propto p^{1/\gamma}\)) Zustandsänderungen
 = Alexanders Setups =

2026-06-24T14:59:50Z

+Elias

← Older revision		Revision as of 16:59, 24 June 2026
Line 1:		Line 1:
			= Notizen Elias =

			* \(S\) = ''spezifische'' Entropie, d.h. Entropie pro Masse
			* barotrop: \(\rho=\rho(p)\), gültig fürs ideale Gas bei isothermen (\(\rho\propto p\)) und isentropen (\(\rho\propto p^{1/\gamma}\)) Zustandsänderungen

	= Alexanders Setups =		= Alexanders Setups =

2026-02-05T00:12:50Z

Alexanders Setups: +quadratischer und kubischer T-Verlauf

← Older revision		Revision as of 02:12, 5 February 2026
Line 75:		Line 75:

	(<code>mtot</code> in g/mol)		(<code>mtot</code> in g/mol)

			== Alternativen zum linearen Verlauf ==

			=== quadratisch ===

			\[
			T(y) = b\Delta T\left(\frac{y}{\Delta Y}-\frac{y^2}{\Delta Y^2}\right)+\frac{\Delta T}{\Delta Y}\,y+T_0
			\]
			mit
			\[
			T'(0)=(1+b)\frac{\Delta T}{\Delta Y}
			\quad,\quad
			T'(\Delta y/2)=\frac{\Delta T}{\Delta Y}
			\quad,\quad
			T'(\Delta y)=(1-b)\frac{\Delta T}{\Delta Y}
			\]

			=== kubisch ===

			\[
			T(y) = b\Delta T\left(2\frac{y^3}{\Delta Y^3}-3\frac{y^2}{\Delta Y^2}+\frac{y}{\Delta Y}\right)+\frac{\Delta T}{\Delta Y}\,y+T_0
			\]
			mit
			\[
			T'(0)=T'(\Delta Y)=(1+b)\frac{\Delta T}{\Delta Y}
			\quad,\quad
			T'(\Delta Y/2)=(1-b/2)\frac{\Delta T}{\Delta Y}
			\]
			und
			\[
			T(\Delta Y/2)=T_0+\Delta T/2
			\]

2026-01-15T14:47:54Z

analytisch: α≈6

← Older revision		Revision as of 16:47, 15 January 2026
Line 28:		Line 28:
	\alpha = \frac{mg}{k_\text{B}T'}		\alpha = \frac{mg}{k_\text{B}T'}
	\]		\]
	Für eine "Inversionslage" muss \(\alpha<1\) gelten.		Für eine "Inversionslage" muss \(\alpha<1\) gelten. Für die ursprünglichen Werte (\(T'=9\cdot 10^6 K/R_\odot, g=GM_\odot/R_\odot^2\)) gilt aber \(\alpha\approx 6\).

			<small>(/ (* 2.353e-3 6.674e-11 1.989e30) (* 8.314 9e6 6.963e8)) ≈ 5.995</small>

	=== pluto.ini ===		=== pluto.ini ===

2026-01-15T10:34:05Z

Alexanders Setups: +Molecules.conf

← Older revision		Revision as of 12:34, 15 January 2026
Line 65:		Line 65:
	IC_Perturbation_mode		IC_Perturbation_mode
	IC_Perturbation_Amplitude_cs		IC_Perturbation_Amplitude_cs


			=== Molecules.conf ===

			mtot = 2.353
			Ntot = 1
			f_DegreesOfFreedom = 3

			(<code>mtot</code> in g/mol)

2026-01-15T06:04:31Z

pluto.ini

← Older revision		Revision as of 08:04, 15 January 2026
Line 33:		Line 33:
	=== pluto.ini ===		=== pluto.ini ===

	* Grid = [0.0,0.1] x [1.0,1.1]		* Grid = [0.0,0.1] x [1.0,1.1] (LE = Sonnenradius)
	* Boundary = [periodic,periodic] x [reflective,userdef]		* Boundary = [periodic,periodic] x [reflective,userdef]
	* Anfangswerte:		* Anfangswerte:
Line 40:		Line 40:
	IC_Perturbation_mode +5.0		IC_Perturbation_mode +5.0
	IC_Perturbation_Amplitude_cs +0.3		IC_Perturbation_Amplitude_cs +0.3


	=== bodyforce.c ===		=== bodyforce.c ===

2026-01-14T22:46:40Z

lin_T-Grad: +analytisch

← Older revision		Revision as of 00:46, 15 January 2026
Line 18:		Line 18:
	* BC_Outflow_NoInflow für \(v_y\) und \(\rho \Rightarrow\) BC_ZeroGradient für \(P\)		* BC_Outflow_NoInflow für \(v_y\) und \(\rho \Rightarrow\) BC_ZeroGradient für \(P\)
	* BC_ZeroGradient für \(v_x\) (und \(v_z\))		* BC_ZeroGradient für \(v_x\) (und \(v_z\))


			==== analytisch ====

			Für einen linearen Verlauf \(T=T_0-T'\, y\) mit \(T'=(T_0-T_1)/L_y\) können Druck- und Dichteprofil analytisch bestimmt werden, Letzteres lautet
			\[
			\rho(y)=\frac{P_0\,m}{k_\text{B}T_0}\left(1-\frac{T'}{T_0}y\right)^{\alpha-1}
			\qquad,\qquad
			\alpha = \frac{mg}{k_\text{B}T'}
			\]
			Für eine "Inversionslage" muss \(\alpha<1\) gelten.


	=== pluto.ini ===		=== pluto.ini ===

	* Grid = [0.0,0.1] x [1.0,1.1]		* Grid = [0.0,0.1] x [1.0,1.1]
	* Boundary = [periodic,periodic] x [reflective,userdef]		* Boundary = [periodic,periodic] x [reflective,userdef]
Line 28:		Line 40:
	IC_Perturbation_mode +5.0		IC_Perturbation_mode +5.0
	IC_Perturbation_Amplitude_cs +0.3		IC_Perturbation_Amplitude_cs +0.3



Line 37:		Line 48:


			=== makefile_user_simulationdefs ===

	~~=== makefile_user_simulationdefs ===~~
	DIMENSIONS = 2		DIMENSIONS = 2
	GEOMETRY = CARTESIAN		GEOMETRY = CARTESIAN
Line 46:		Line 57:

	=== user_defined_parameters.h ===		=== user_defined_parameters.h ===

	UNIT_DENSITY (1e-0)		UNIT_DENSITY (1e-0)
	UNIT_LENGTH (CONST_Rsun)		UNIT_LENGTH (CONST_Rsun)

2026-01-14T18:16:18Z

InitDomain(): fix Druck-DGL

← Older revision		Revision as of 20:16, 14 January 2026
Line 6:		Line 6:

	==== <code>InitDomain()</code> ====		==== <code>InitDomain()</code> ====

	* lineares \(T\)-Profil gemäß <code>IC_Temperature_X2BEG_K</code>><code> IC_Temperature_X2END_K</code>		* lineares \(T\)-Profil gemäß <code>IC_Temperature_X2BEG_K</code>><code> IC_Temperature_X2END_K</code>
	* ideales Gas: isotherme Kompressibilität = \(1/P \Rightarrow \tilde c^2=P/\rho=k_BT/m\)		* ideales Gas: isotherme Kompressibilität = \(1/P \Rightarrow \tilde c^2=P/\rho=k_BT/m\)
	* Druck-DGL \(-\partial_y p=-pg/\tilde c^2(T) \Rightarrow\) ~~Druckprofil (~~numerisch) mit RB \(P(y_\text{min})=10^9\)bar		* Druck-DGL \(\partial_y P=-Pg/\tilde c^2(T) \Rightarrow P(y)\) numerisch mit RB \(P(y_\text{min})=10^9\)bar
	* Dichteprofil aus \(\rho=P/\tilde c^2(T)\), beide Profile ggf. von unten gedeckelt durch <code>RHO_min</code> aus ''Restrictions.conf''		* Dichteprofil aus \(\rho=P/\tilde c^2(T)\), beide Profile ggf. von unten gedeckelt durch <code>RHO_min</code> aus ''Restrictions.conf''
	* \(v_y\) gestört durch Sinus mit Wellenlänge = Systembreite/<code>IC_Perturbation_mode</code> und Amplitude <code>IC_Perturbation_Amplitude_cs</code> × \(\exp(-\vert y_\text{rel}-1/2\vert)\)		* \(v_y\) gestört durch Sinus mit Wellenlänge = Systembreite/<code>IC_Perturbation_mode</code> und Amplitude <code>IC_Perturbation_Amplitude_cs</code> × \(\exp(-\vert y_\text{rel}-1/2\vert)\)