BA Elias Saarmann - Revision history

Lothar.brendel: /* pluto.ini-Syntax */ x[23]-grid

2026-07-31T12:57:09Z

pluto.ini-Syntax: x[23]-grid

← Older revision		Revision as of 14:57, 31 July 2026
Line 207:		Line 207:

	* Alle Parameter müssen angegeben werden, auch wenn gar kein Zugriff mittels <code>g_inputParam[]</code> erfolgt.		* Alle Parameter müssen angegeben werden, auch wenn gar kein Zugriff mittels <code>g_inputParam[]</code> erfolgt.
	* Die Reihenfolge der Parameter '''~~muss~~''' ~~der Nummerierung~~ in user_defined_parameters.h folgen.		* Die Reihenfolge der Parameter muss der '''Nummerierung''' in user_defined_parameters.h folgen.
			* Auch für <code>DIMENSIONS = 1</code> sind die Intervalle in <code>x2-grid</code> und <code>x3-grid</code> nicht irrelevant, für <code>GEOMETRY = SPHERICAL</code> daher stets <code>0 1 u 3.141592653589793</code> bzw. <code>0 1 u 6.283185307179586</code>.

Elias.S: /* Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren */

2026-07-30T09:59:52Z

Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren

← Older revision		Revision as of 11:59, 30 July 2026
Line 165:		Line 165:
	* mehr Grid Zellen, optische Dichte \(\tau = \kappa \rho \Delta x\) der einzelnen Zellen sollte klein sein		* mehr Grid Zellen, optische Dichte \(\tau = \kappa \rho \Delta x\) der einzelnen Zellen sollte klein sein
	** Parameter wie oben, 300 Zellen -> Oszilationen werden schwächer, power law erkennbar, immernoch Probleme mit Zero-Gradient-Randbedingung		** Parameter wie oben, 300 Zellen -> Oszilationen werden schwächer, power law erkennbar, immernoch Probleme mit Zero-Gradient-Randbedingung
	** Parameter wie oben, 1000 Zellen -> nur noch sehr sehr schwache Oszillationen, power law klar ~~eknnebra~~, Zero-Gradient scheinbar efüllt		** Parameter wie oben, 1000 Zellen -> nur noch sehr sehr schwache Oszillationen, power law klar erkennbar, Zero-Gradient scheinbar efüllt
	Allerdings: Massenflux steigt linear, kein Steady State erreicht		Allerdings: Massenflux steigt linear, kein Steady State erreicht

Elias.S: /* Simulationsprobleme */

2026-07-29T13:50:05Z

Simulationsprobleme

← Older revision		Revision as of 15:50, 29 July 2026
Line 162:		Line 162:

	Läuft, allerings wieder viele negative Drücke in der Log file und Ergebnisse mit Oszilation und Verletzung der Zero Gradient Bedingung bei Druck, Temperatur und Strahlungsenergiedichte.		Läuft, allerings wieder viele negative Drücke in der Log file und Ergebnisse mit Oszilation und Verletzung der Zero Gradient Bedingung bei Druck, Temperatur und Strahlungsenergiedichte.

			* mehr Grid Zellen, optische Dichte \(\tau = \kappa \rho \Delta x\) der einzelnen Zellen sollte klein sein
			** Parameter wie oben, 300 Zellen -> Oszilationen werden schwächer, power law erkennbar, immernoch Probleme mit Zero-Gradient-Randbedingung
			** Parameter wie oben, 1000 Zellen -> nur noch sehr sehr schwache Oszillationen, power law klar eknnebra, Zero-Gradient scheinbar efüllt
			Allerdings: Massenflux steigt linear, kein Steady State erreicht

	=== Versuch \(\tau_\text{B} = 10^{-1}\) und \(\beta = 10^{-3}\) zu untersuchen ===		=== Versuch \(\tau_\text{B} = 10^{-1}\) und \(\beta = 10^{-3}\) zu untersuchen ===

Lothar.brendel: /* Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren */ 1e-9

2026-07-28T11:11:18Z

Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren: 1e-9

← Older revision		Revision as of 13:11, 28 July 2026
Line 152:		Line 152:
	Lösungsansätze:		Lösungsansätze:
	* grid cells auf 150 erhöhen		* grid cells auf 150 erhöhen
	Simulationsschritte werden schnell sehr klein dt = e-9 und Simulation kommt kaum voran		Simulationsschritte werden schnell sehr klein dt = '''1'''e-9 = 10<sup>-9</sup> und Simulation kommt kaum voran

	* Parameter varrieren, so dass Bondi Radius kleiner wird		* Parameter varrieren, so dass Bondi Radius kleiner wird

Elias.S: /* Simulationsprobleme */

2026-07-28T02:13:07Z

Simulationsprobleme

← Older revision		Revision as of 04:13, 28 July 2026
Line 143:		Line 143:


	Problem:		*Problem: Starke Schwankungen nahe Innenrand für <code> prs </code>, <code> Tgas </code> und <code> erad </code>, Zero-Gradient Randbedingungen werden nicht eingehalten. Negative Drücke in Log File:
	Starke Schwankungen nahe Innenrand für <code> prs </code>, <code> Tgas </code> und <code> erad </code>, Zero-Gradient Randbedingungen werden nicht eingehalten. Negative Drücke in Log File:		**Simulationsparameter: \[\rho_\infty = 3 \cdot 10^{-14} g/cm^3\] \[T_\infty = 2.7 K\] \[\kappa = 0.1 cm^2/g\] \[m = 1 m_\odot\] \[M = 2.01588 g/Mol\] \[r_\text{s} = 3 au\] \[r_\text{o} = 206265.0 au\] grid cells = 114
			**theoretische Parameter:
	* Simulationsparameter:
	\[\rho_\infty = 3 \cdot 10^{-14} g/cm^3\]
	\[T_\infty = 2.7 K\]
	\[\kappa = 0.1 cm^2/g\]
	\[m = 1 m_\odot\]
	\[M = 2.01588 g/Mol\]
	\[r_\text{s} = 3 au\]
	\[r_\text{o} = 206265.0 au\]
	grid cells = 114

	* theoretische Parameter:
	\[\tau_\text{B} = 1.28 \cdot 10^{3}\]		\[\tau_\text{B} = 1.28 \cdot 10^{3}\]
	\[\beta = 1.69 \cdot 10^{-3}\]		\[\beta = 1.69 \cdot 10^{-3}\]
	\[r_\text{B} = 28450.69 au\]		\[r_\text{B} = 28450.69 au\]


	Lösungsansätze:		Lösungsansätze:
Line 167:		Line 155:

	* Parameter varrieren, so dass Bondi Radius kleiner wird		* Parameter varrieren, so dass Bondi Radius kleiner wird
			**Simulationsparameter: \[\rho_\infty = 3 \cdot 10^{-12} g/cm^3\] \[T_\infty = 10 K\] \[\kappa = 0.01 cm^2/g\] \[m = 1 m_\odot\] \[M = 2.01588 g/Mol\] \[r_\text{s} = 3 au\] \[r_\text{o} = 206265.0 au\] grid cells = 150
	* Simulationsparameter:		**theoretische Parameter:
	\[\rho_\infty = 3 \cdot 10^{-12} g/cm^3\]
	\[T_\infty = 10 K\]
	\[\kappa = 0.01 cm^2/g\]
	\[m = 1 m_\odot\]
	\[M = 2.01588 g/Mol\]
	\[r_\text{s} = 3 au\]
	\[r_\text{o} = 206265.0 au\]
	grid cells = 150

	* theoretische Parameter:
	\[\tau_\text{B} = 3.45 \cdot 10^{3}\]		\[\tau_\text{B} = 3.45 \cdot 10^{3}\]
	\[\beta = 2.38 \cdot 10^{-3}\]		\[\beta = 2.38 \cdot 10^{-3}\]
Line 202:		Line 180:
	\[r_\text{B} = 28450.69 au\]		\[r_\text{B} = 28450.69 au\]

	* Problem:		* Problem: Zeitschritte werden sehr klein (dt = 1e-9) und die Simulation braucht ewig: nach fast 3e6 Simulationsschritten immer noch bei t = 0.23.
	Zeitschritte werden sehr klein (dt = 1e-9) und die Simulation braucht ewig: nach fast 3e6 Simulationsschritten immer noch bei t = 0.23.		**Simulationsparameter: \[\rho_\infty = 3 \cdot 10^{-18} g/cm^3\] \[T_\infty = 2.7 K\] \[\kappa = 0.1 cm^2/g\] \[m = 1 m_\odot\] \[M = 2.01588 g/Mol\] \[r_\text{s} = 3 au\] \[r_\text{o} = 206265.0 au\] grid cells = 114
			**theoretische Parameter:
			\[\tau_\text{B} = 1.28 \cdot 10^{-1}\]
			\[\beta = 1.69 \cdot 10^{-3}\]
			\[r_\text{B} = 28450.69 au\]

	== auf Elias' Laptop (Ubuntu 26) ==		== auf Elias' Laptop (Ubuntu 26) ==

Elias.S: /* Simulationsprobleme */

2026-07-28T02:08:46Z

Simulationsprobleme

← Older revision		Revision as of 04:08, 28 July 2026
Line 132:		Line 132:
	=== Versuch \(\tau_\text{B} = 10^{-3}\) und \(\beta = 10^{-3}\) zu reproduzieren ===		=== Versuch \(\tau_\text{B} = 10^{-3}\) und \(\beta = 10^{-3}\) zu reproduzieren ===

	* Simulationsparameter:		*Problem: Code bricht ab, dt wird zu klein
	\[\rho_\infty = 3 \cdot 10^{-20} g/cm^3\]		**Simulationsparameter: \[\rho_\infty = 3 \cdot 10^{-20} g/cm^3\] \[T_\infty = 2.7 K\] \[\kappa = 0.1 cm^2/g\] \[m = 1 m_\odot\] \[M = 2.01588 g/Mol\] \[r_\text{s} = 3 au\] \[r_\text{o} = 206265.0 au\] grid cells = 114
	\[T_\infty = 2.7 K\]		**theoretische Parameter:
	\[\kappa = 0.1 cm^2/g\]
	\[m = 1 m_\odot\]
	\[M = 2.01588 g/Mol\]
	\[r_\text{s} = 3 au\]
	\[r_\text{o} = 206265.0 au\]
	grid cells = 114

	* theoretische Parameter:
	\[\tau_\text{B} = 1.28 \cdot 10^{-3}\]		\[\tau_\text{B} = 1.28 \cdot 10^{-3}\]
	\[\beta = 1.69 \cdot 10^{-3}\]		\[\beta = 1.69 \cdot 10^{-3}\]
	\[r_\text{B} = 28450.69 au\]		\[r_\text{B} = 28450.69 au\]

	* Problem:
	~~Code bricht ab, dt wird zu klein~~

	=== Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren ===		=== Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren ===

Elias.S: Fortschritt bei Problemlösung

2026-07-28T02:05:11Z

Fortschritt bei Problemlösung

@@ Line 151: / Line 151: @@
 === Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\)  zu reproduzieren ===
 * Simulationsparameter:
@@ Line 167: / Line 172: @@
 \[r_\text{B} = 28450.69 au\]
-* Problem:
-Starke Schwankungen nahe Innenrand für <code> prs </code>, <code> Tgas </code> und <code> erad </code>, Zero-Gradient Randbedingungen werden nicht eingehalten.
+Lösungsansätze:
 === Versuch \(\tau_\text{B} = 10^{-1}\) und \(\beta = 10^{-3}\)  zu untersuchen ===

Elias.S: /* Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren */

2026-07-28T01:45:04Z

Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren

← Older revision		Revision as of 03:45, 28 July 2026
Line 153:		Line 153:

	* Simulationsparameter:		* Simulationsparameter:
	\[\rho_\infty = 3 \cdot 10^{-16} g/cm^3\]		\[\rho_\infty = 3 \cdot 10^{-14} g/cm^3\]
	\[T_\infty = 2.7 K\]		\[T_\infty = 2.7 K\]
	\[\kappa = 0.1 cm^2/g\]		\[\kappa = 0.1 cm^2/g\]

Elias.S: /* Simulationsprobleme */

2026-07-28T01:19:04Z

Simulationsprobleme

← Older revision		Revision as of 03:19, 28 July 2026
Line 132:		Line 132:
	=== Versuch \(\tau_\text{B} = 10^{-3}\) und \(\beta = 10^{-3}\) zu reproduzieren ===		=== Versuch \(\tau_\text{B} = 10^{-3}\) und \(\beta = 10^{-3}\) zu reproduzieren ===

	Simulationsparameter:		* Simulationsparameter:

	\[\rho_\infty = 3 \cdot 10^{-20} g/cm^3\]		\[\rho_\infty = 3 \cdot 10^{-20} g/cm^3\]
	\[T_\infty = 2.7 K\]		\[T_\infty = 2.7 K\]
Line 143:		Line 142:
	grid cells = 114		grid cells = 114

	theoretische Parameter:		* theoretische Parameter:


	\[\tau_\text{B} = 1.28 \cdot 10^{-3}\]		\[\tau_\text{B} = 1.28 \cdot 10^{-3}\]
	\[\beta = 1.69 \cdot 10^{-3}\]		\[\beta = 1.69 \cdot 10^{-3}\]
	\[r_\text{B} = 28450.69 au\]		\[r_\text{B} = 28450.69 au\]

	Problem:		* Problem:
	Code bricht ab, dt wird zu klein		Code bricht ab, dt wird zu klein

	=== Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren ===		=== Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren ===

	Simulationsparameter:		* Simulationsparameter:


	\[\rho_\infty = 3 \cdot 10^{-16} g/cm^3\]		\[\rho_\infty = 3 \cdot 10^{-16} g/cm^3\]
	\[T_\infty = 2.7 K\]		\[T_\infty = 2.7 K\]
Line 167:		Line 162:
	grid cells = 114		grid cells = 114

	theoretische Parameter:		* theoretische Parameter:


	\[\tau_\text{B} = 1.28 \cdot 10^{3}\]		\[\tau_\text{B} = 1.28 \cdot 10^{3}\]
	\[\beta = 1.69 \cdot 10^{-3}\]		\[\beta = 1.69 \cdot 10^{-3}\]
	\[r_\text{B} = 28450.69 au\]		\[r_\text{B} = 28450.69 au\]

	Problem:		* Problem:
	Starke Schwankungen nahe Innenrand für <code> prs </code>, <code> Tgas </code> und <code> erad </code>, Zero-Gradient Randbedingungen werden nicht eingehalten.		Starke Schwankungen nahe Innenrand für <code> prs </code>, <code> Tgas </code> und <code> erad </code>, Zero-Gradient Randbedingungen werden nicht eingehalten.

	=== Versuch \(\tau_\text{B} = 10^{-1}\) und \(\beta = 10^{-3}\) zu untersuchen ===		=== Versuch \(\tau_\text{B} = 10^{-1}\) und \(\beta = 10^{-3}\) zu untersuchen ===

	Simulationsparameter:		* Simulationsparameter:


	\[\rho_\infty = 3 \cdot 10^{-18} g/cm^3\]		\[\rho_\infty = 3 \cdot 10^{-18} g/cm^3\]
	\[T_\infty = 2.7 K\]		\[T_\infty = 2.7 K\]
Line 191:		Line 182:
	grid cells = 114		grid cells = 114

	theoretische Parameter:		* theoretische Parameter:


	\[\tau_\text{B} = 1.28 \cdot 10^{-1}\]		\[\tau_\text{B} = 1.28 \cdot 10^{-1}\]
	\[\beta = 1.69 \cdot 10^{-3}\]		\[\beta = 1.69 \cdot 10^{-3}\]
	\[r_\text{B} = 28450.69 au\]		\[r_\text{B} = 28450.69 au\]

	Problem:		* Problem:
	Zeitschritte werden sehr klein (dt = 1e-9) und die Simulation braucht ewig: nach fast 3e6 Simulationsschritten immer noch bei t = 0.23.		Zeitschritte werden sehr klein (dt = 1e-9) und die Simulation braucht ewig: nach fast 3e6 Simulationsschritten immer noch bei t = 0.23.

Elias.S: Aktuelle Simulationsprobleme näher beschrieben

2026-07-28T01:17:35Z

Aktuelle Simulationsprobleme näher beschrieben

← Older revision		Revision as of 03:17, 28 July 2026
Line 130:		Line 130:
	== Simulationsprobleme ==		== Simulationsprobleme ==

	=== Versuch ~~Wertebereich aus Fig. 5~~ zu reproduzieren ===		=== Versuch \(\tau_\text{B} = 10^{-3}\) und \(\beta = 10^{-3}\) zu reproduzieren ===

	Simulationsparameter:		Simulationsparameter:

			\[\rho_\infty = 3 \cdot 10^{-20} g/cm^3\]
			\[T_\infty = 2.7 K\]
			\[\kappa = 0.1 cm^2/g\]
			\[m = 1 m_\odot\]
			\[M = 2.01588 g/Mol\]
			\[r_\text{s} = 3 au\]
			\[r_\text{o} = 206265.0 au\]
			grid cells = 114

	theoretische Parameter:		theoretische Parameter:


			\[\tau_\text{B} = 1.28 \cdot 10^{-3}\]
			\[\beta = 1.69 \cdot 10^{-3}\]
			\[r_\text{B} = 28450.69 au\]

	Problem:		Problem:
	Code bricht ab, dt wird zu klein		Code bricht ab, dt wird zu klein

	=== Versuch ~~Wertebereich aus Fig. 6~~ zu reproduzieren ===		=== Versuch \(\tau_\text{B} = 10^{3}\) und \(\beta = 10^{-3}\) zu reproduzieren ===

	Simulationsparameter:		Simulationsparameter:


			\[\rho_\infty = 3 \cdot 10^{-16} g/cm^3\]
			\[T_\infty = 2.7 K\]
			\[\kappa = 0.1 cm^2/g\]
			\[m = 1 m_\odot\]
			\[M = 2.01588 g/Mol\]
			\[r_\text{s} = 3 au\]
			\[r_\text{o} = 206265.0 au\]
			grid cells = 114

	theoretische Parameter:		theoretische Parameter:


			\[\tau_\text{B} = 1.28 \cdot 10^{3}\]
			\[\beta = 1.69 \cdot 10^{-3}\]
			\[r_\text{B} = 28450.69 au\]

	Problem:		Problem:
	Starke Schwankungen nahe Innenrand für <code> prs </code>, <code> Tgas </code> und <code> erad </code>, Zero-Gradient Randbedingungen werden nicht eingehalten.		Starke Schwankungen nahe Innenrand für <code> prs </code>, <code> Tgas </code> und <code> erad </code>, Zero-Gradient Randbedingungen werden nicht eingehalten.

	=== Versuch ~~mittleren Wertebereich~~ zu untersuchen ===		=== Versuch \(\tau_\text{B} = 10^{-1}\) und \(\beta = 10^{-3}\) zu untersuchen ===

	Simulationsparameter:		Simulationsparameter:


			\[\rho_\infty = 3 \cdot 10^{-18} g/cm^3\]
			\[T_\infty = 2.7 K\]
			\[\kappa = 0.1 cm^2/g\]
			\[m = 1 m_\odot\]
			\[M = 2.01588 g/Mol\]
			\[r_\text{s} = 3 au\]
			\[r_\text{o} = 206265.0 au\]
			grid cells = 114

	theoretische Parameter:		theoretische Parameter:


			\[\tau_\text{B} = 1.28 \cdot 10^{-1}\]
			\[\beta = 1.69 \cdot 10^{-3}\]
			\[r_\text{B} = 28450.69 au\]

	Problem:		Problem: