PLUTO Test Problems  4.4-patch2
Functions
init.c File Reference

Planar shock wave problem with particles. More...

Functions

void InitDomain (Data *d, Grid *grid)
 
void Analysis (const Data *d, Grid *grid)
 
void BackgroundField (double x1, double x2, double x3, double *B0)
 
void UserDefBoundary (const Data *d, RBox *box, int side, Grid *grid)
 

Detailed Description

Set the initial condition for a Planar shocks wave problem. Ambient unshocked density and pressure are constant everywhere and equal to rho0 and Pg. While rho0 is an input parameter, Pg=1.e-4 (for MHD runs) and Pg = 1.0e-8 (for RMHD runs). Also for RMHD runs, a transverse velocity of 0.01 is set initially in the Ambient medium.

The shocked gas (x < 1.0) is allowed to move from left to right and the fluid quantities in this shocked region are given as input Parameters. The fluid quantities in the shocked region are estimated by solving the Rankine-Hugoniot conditions for specified shock velocities.

The input parameters read from pluto.ini are labeled as:

The available configuration refer to: -#01 Cartesian (2.5D, MHD, Ideal, THETA_B = 30.0, BETA = 100.) -#02 Cartesian (2.5D, RMHD, Taub, THETA_B = 3.0, BETA = 0.01) -#03 Cartesian (2.5D, RMHD, Taub, THETA_B = 83.0, BETA = 0.01) -#04 Same as #02, but with a perturbation in the upstream.

Author
B. Vaidya (bvaid.nosp@m.ya@i.nosp@m.iti.a.nosp@m.c.in) A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it)
Date
Sep 19, 2019

References:

  • "A Particle Module for the PLUTO Code. II. Hybrid Framework for Modeling Nonthermal Emission from Relativistic Magnetized Flows", Vaidya et al. (2018) ApJS [ Sec. 4.1 - 4.2 ]

Function Documentation

◆ Analysis()

void Analysis ( const Data *  d,
Grid *  grid 
)

Perform runtime data analysis.

Parameters
[in]dthe PLUTO Data structure
[in]gridpointer to array of Grid structures

◆ BackgroundField()

void BackgroundField ( double  x1,
double  x2,
double  x3,
double *  B0 
)

Define the component of a static, curl-free background magnetic field.

Parameters
[in]x1position in the 1st coordinate direction $x_1$
[in]x2position in the 2nd coordinate direction $x_2$
[in]x3position in the 3rd coordinate direction $x_3$
[out]B0array containing the vector componens of the background magnetic field

◆ InitDomain()

void InitDomain ( Data *  d,
Grid *  grid 
)

Assign initial condition by looping over the computational domain. Called after the usual Init() function to assign initial conditions on primitive variables. Value assigned here will overwrite those prescribed during Init().

◆ UserDefBoundary()

void UserDefBoundary ( const Data *  d,
RBox *  box,
int  side,
Grid *  grid 
)

Assign user-defined boundary conditions.

Parameters
[in,out]dpointer to the PLUTO data structure containing cell-centered primitive quantities (d->Vc) and staggered magnetic fields (d->Vs, when used) to be filled.
[in]boxpointer to a RBox structure containing the lower and upper indices of the ghost zone-centers/nodes or edges at which data values should be assigned.
[in]sidespecifies the boundary side where ghost zones need to be filled. It can assume the following pre-definite values: X1_BEG, X1_END, X2_BEG, X2_END, X3_BEG, X3_END. The special value side == 0 is used to control a region inside the computational domain.
[in]gridpointer to an array of Grid structures.