MA Alexander Grunewald/Code-Fragmente - Revision history

Lothar.brendel: +up

2026-03-03T08:43:52Z

+up

2026-03-03T08:42:43Z

+analytische Anfangsprofile

← Older revision		Revision as of 10:42, 3 March 2026
Line 102:		Line 102:
	#define USER_DEF_PARAMETERS 5		#define USER_DEF_PARAMETERS 5
	</syntaxhighlight>		</syntaxhighlight>

			= analytische Anfangsprofile =

			== 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_0pow(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.

2026-02-19T17:52:00Z

in init.c:UserDefBoundary(): tidy up

← Older revision		Revision as of 19:52, 19 February 2026
Line 6:		Line 6:
	{		{
	if(side == 0){ // is called before the other sides		if(side == 0){ // is called before the other sides
	~~int kl,jl,il;~~
	double Ymin = g_domBeg[JDIR],Ymax = g_domEnd[JDIR],		double Ymin = g_domBeg[JDIR],Ymax = g_domEnd[JDIR],
	dzs = g_inputParam[Dampingzones_size_cu],		dzs = g_inputParam[Dampingzones_size_cu],
	tau = g_inputParam[Damping_time_cu];		tau = g_inputParam[Damping_time_cu];
	~~RBox dom_box;~~
	if (dzs > 0 && g_intStage == 1) {		if (dzs > 0 && g_intStage == 1) {
			int kl,jl,il;
			RBox dom_box;
	DOM_LOOP(kl,jl,il) {		DOM_LOOP(kl,jl,il) {
	double damp = 0,y = grid->x[JDIR][jl];		double damp = 0,y = grid->x[JDIR][jl];

2026-02-19T17:44:37Z

in init.c:UserDefBoundary(): vdots

← Older revision		Revision as of 19:44, 19 February 2026
Line 4:		Line 4:

	<syntaxhighlight lang="C">		<syntaxhighlight lang="C">
			{
	if(side == 0){ // is called before the other sides		if(side == 0){ // is called before the other sides
	int kl,jl,il;		int kl,jl,il;
Line 27:		Line 27:
	ApplyRestrictions(d, grid);		ApplyRestrictions(d, grid);
	}		}
	else if(side == X2_END) { ~~...~~		else if(side == X2_END) {
	</syntaxhighlight>		</syntaxhighlight>
			⋮

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

2026-02-19T17:42:16Z

in user_defined_parameters.h: Werte

← Older revision		Revision as of 19:42, 19 February 2026
Line 32:		Line 32:
	== in user_defined_parameters.h ==		== in user_defined_parameters.h ==

			⋮
	<syntaxhighlight lang="C">		<syntaxhighlight lang="C">
	#define Dampingzones_size_cu 5		#define Dampingzones_size_cu 5
Line 42:		Line 43:
	</syntaxhighlight>		</syntaxhighlight>

	Für die beiden neuen Parameter sind natürlich in pluto.ini noch Werte einzusetzen.		Für die beiden neuen Parameter sind natürlich in pluto.ini noch Werte einzusetzen:

			* Sinnvoll sind \(0.02\) für die <code>Dampingzones_size_cu</code>, das ist die Breite einer jeden Zone, es bleibt also eine \(0.06\) breite, ungedämpfte Zone in der Mitte. Zum Innenrand einer Zone wird die Dämpfung linear schwächer.
			* Für die charakteristische Zeit <code>Damping_time_cu</code> gilt je kleiner, desto stärker die Dämpfung; \(5\cdot 10^{-7}\) ist schon recht stark.

	= Störung in der Beschleunigung =		= Störung in der Beschleunigung =

2026-02-19T16:53:08Z

Dämpfungszonen

← Older revision		Revision as of 18:53, 19 February 2026
Line 1:		Line 1:
			= Dämpfungzonen =

			== in init.c:UserDefBoundary() ==

			<syntaxhighlight lang="C">

			if(side == 0){ // is called before the other sides
			int kl,jl,il;
			double Ymin = g_domBeg[JDIR],Ymax = g_domEnd[JDIR],
			dzs = g_inputParam[Dampingzones_size_cu],
			tau = g_inputParam[Damping_time_cu];
			RBox dom_box;
			if (dzs > 0 && g_intStage == 1) {
			DOM_LOOP(kl,jl,il) {
			double damp = 0,y = grid->x[JDIR][jl];
			if (y < Ymin+dzs)
			damp = Ymin+dzs-y;
			else if (y > Ymax-dzs)
			damp = y-(Ymax-dzs);
			damp = 1/(dzstau);
			d->Vc[VX2][kl][jl][il] = MAX(1-dampg_dt,0.0);
			}
			RBoxDefine (IBEG,IEND,JBEG,JEND,KBEG,KEND,CENTER,&dom_box);
			PrimToCons3D(d->Vc, d->Uc, &dom_box);
			}

			ApplyRestrictions(d, grid);
			}
			else if(side == X2_END) { ...
			</syntaxhighlight>

			== in user_defined_parameters.h ==

			<syntaxhighlight lang="C">
			#define Dampingzones_size_cu 5
			#define Damping_time_cu 6

			//
			// Specify number of user defined parameters:
			//
			#define USER_DEF_PARAMETERS 7
			</syntaxhighlight>

			Für die beiden neuen Parameter sind natürlich in pluto.ini noch Werte einzusetzen.


	= Störung in der Beschleunigung =		= Störung in der Beschleunigung =

2026-02-12T11:46:14Z

2026-02-12T11:45:33Z

2026-02-12T11:42:18Z

in body_force.c: Einrückung

← Older revision		Revision as of 13:42, 12 February 2026
Line 22:		Line 22:

	double g_grav = G_GravityConstant * SolarMass / (SolarRadius*SolarRadius),		double g_grav = G_GravityConstant * SolarMass / (SolarRadius*SolarRadius),
	perturbation_mode = g_inputParam[IC_Perturbation_mode], // dimensionless and integer		perturbation_mode = g_inputParam[IC_Perturbation_mode], // dimensionless and integer
	perturbation_amplitude = g_inputParam[IC_Perturbation_Amplitude_g]*g_grav,		perturbation_amplitude = g_inputParam[IC_Perturbation_Amplitude_g]*g_grav,
	Xmin = grid->xl_glob[IDIR][grid->gbeg[IDIR]],		Xmin = grid->xl_glob[IDIR][grid->gbeg[IDIR]],
	Xmax = grid->xr_glob[IDIR][grid->gend[IDIR]],		Xmax = grid->xr_glob[IDIR][grid->gend[IDIR]],
	Ymin = grid->xl_glob[JDIR][grid->gbeg[JDIR]],		Ymin = grid->xl_glob[JDIR][grid->gbeg[JDIR]],
	Ymax = grid->xr_glob[JDIR][grid->gend[JDIR]],		Ymax = grid->xr_glob[JDIR][grid->gend[JDIR]],
	ampy = perturbation_amplitude * exp(- 10.0 * fabs(0.5 - (x2-Ymin) / (Ymax-Ymin)));		ampy = perturbation_amplitude * exp(- 10.0 * fabs(0.5 - (x2-Ymin) / (Ymax-Ymin)));

	g[JDIR] += ampy * sin(perturbation_mode * 2.0M_PI x1 / (Xmax-Xmin));		g[JDIR] += ampy * sin(perturbation_mode * 2.0M_PI x1 / (Xmax-Xmin));

2026-02-12T11:41:14Z

in init.c: soll

← Older revision		Revision as of 13:41, 12 February 2026
Line 7:		Line 7:
	<code>data->Vc[VX2][kl][jl][il] += perturbation_amplitude * sin(perturbation_mode * 2.0M_PI x1 / (Xmax-Xmin)) * ampy;</code>		<code>data->Vc[VX2][kl][jl][il] += perturbation_amplitude * sin(perturbation_mode * 2.0M_PI x1 / (Xmax-Xmin)) * ampy;</code>

	muss natürlich raus. (Und alles was dazugehört, ~~kann~~ auch ~~ruhig~~ raus.)		muss natürlich raus. (Und alles was dazugehört, auch raus.)

	== in body_force.c ==		== in body_force.c ==