init.c File Reference

Hydrodynamics blast wave problem. More...

#include "pluto.h"

Include dependency graph for init.c:

Go to the source code of this file.

Functions
void	Init (double *v, double x1, double x2, double x3)
void	Analysis (const Data *d, Grid *grid)
void	UserDefBoundary (const Data *d, RBox *box, int side, Grid *grid)

Detailed Description

Hydrodynamics blast wave problem.

Set up an isothermal blast wave in a non-uniform, randomnly perturbed density field. The background perturbation is interpolated read from an input file generated externally.

Authors

A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)
P. Tzeferacos (petro.nosp@m.s.tz.nosp@m.efera.nosp@m.cos@.nosp@m.ph.un.nosp@m.ito..nosp@m.it)

Date

May 06, 2014

Definition in file init.c.

Function Documentation

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

Perform runtime data analysis.

Parameters

[in]	d	the PLUTO Data structure
[in]	grid	pointer to array of Grid structures

Definition at line 71 of file init.c.

{

}

void Init	(	double *	v,
		double	x1,
		double	x2,
		double	x3
	)

The Init() function can be used to assign initial conditions as as a function of spatial position.

Parameters

[out]	v	a pointer to a vector of primitive variables
[in]	x1	coordinate point in the 1st dimension
[in]	x2	coordinate point in the 2nd dimension
[in]	x3	coordinate point in the 3rdt dimension

The meaning of x1, x2 and x3 depends on the geometry:

Variable names are accessed by means of an index v[nv], where nv = RHO is density, nv = PRS is pressure, nv = (VX1, VX2, VX3) are the three components of velocity, and so forth.

Definition at line 19 of file init.c.

{
  static int first_call = 1;
  double r;

  #if GEOMETRY == CARTESIAN
   r = D_EXPAND(x1*x1, + x2*x2, + x3*x3);
   r = sqrt(r);
  #elif GEOMETRY == POLAR
   r = x1;
  #endif

  v[RHO] = 1.0;

  #if ADD_TURBULENCE == YES
  if (first_call){
    int k, input_var[256];
    
    input_var[0] = RHO;
    input_var[1] = -1;
    InputDataSet ("./grid0.out",input_var);
    InputDataRead("./rho0.dbl"," ");
    first_call = 0;
  }
  InputDataInterpolate(v, x1, x2, x3);  /* -- interpolate density from 
                                              input data file -- */
  #endif

  #if EOS == IDEAL
   v[PRS] = 1.0/g_gamma;  /* ambient pressure */
  #elif EOS == PVTE_LAW
   v[PRS] = 1.0;
  #endif

/* -- set a high pressure region inside r = 1 -- */

  if (r < 0.1) {
    v[RHO] = g_inputParam[RHO_IN];
    #if HAVE_ENERGY
     v[PRS] = g_inputParam[PRS_IN];
    #endif
  }

  v[VX1] = 0.0;
  v[VX2] = 0.0;
  v[VX3] = 0.0;

}

Here is the call graph for this function:

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

Assign user-defined boundary conditions.

Parameters

[in,out]	d	pointer to the PLUTO data structure containing cell-centered primitive quantities (d->Vc) and staggered magnetic fields (d->Vs, when used) to be filled.
[in]	box	pointer 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]	side	specifies 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]	grid	pointer to an array of Grid structures.

Definition at line 89 of file init.c.

{
  int   i, j, k, nv;
  real  *x1, *x2, *x3;

  x1 = grid[IDIR].x;
  x2 = grid[JDIR].x;
  x3 = grid[KDIR].x;

  if (side == 0) {    /* -- check solution inside domain -- */
    DOM_LOOP(k,j,i){};
  }

  if (side == X1_BEG){  /* -- X1_BEG boundary -- */
    X1_BEG_LOOP(k,j,i){}
  }

  if (side == X1_END){  /* -- X1_END boundary -- */
    X1_END_LOOP(k,j,i){}
  }

  if (side == X2_BEG){  /* -- X2_BEG boundary -- */
    X2_BEG_LOOP(k,j,i){}
  }

  if (side == X2_END){  /* -- X2_END boundary -- */
    X2_END_LOOP(k,j,i){}
  }

  if (side == X3_BEG){  /* -- X3_BEG boundary -- */
    X3_BEG_LOOP(k,j,i){}
  }
 
  if (side == X3_END) {  /* -- X3_END boundary -- */
    X3_END_LOOP(k,j,i){}
  }
}

Here is the call graph for this function: