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

MHD Shock Cloud interaction. More...

Functions

void InitDomain (Data *d, Grid *grid)
 

Detailed Description

Set the initial condition for the 2D or 3D shock-cloud problem. The shock propagates in the x-direction and it is initially located at x=0.6 with post- and -pre shock value equal to

\[ \left(\begin{array}{l} \rho \\ \noalign{\medskip} v_x \\ \noalign{\medskip} p \\ \noalign{\medskip} B_y \\ \noalign{\medskip} B_z \end{array}\right) = \left(\begin{array}{l} 3.86859 \\ \noalign{\medskip} 0 \\ \noalign{\medskip} 167.345 \\ \noalign{\medskip} B_{\rm post} \\ \noalign{\medskip} -B_{\rm post} \\ \noalign{\medskip} \end{array}\right) \quad\mathrm{for}\quad x < 0.6\,,\qquad\qquad \left(\begin{array}{l} \rho \\ \noalign{\medskip} v_x \\ \noalign{\medskip} p \\ \noalign{\medskip} B_y \\ \noalign{\medskip} B_z \end{array}\right) = \left(\begin{array}{l} 1 \\ \noalign{\medskip} -11.2536 \\ \noalign{\medskip} 1 \\ \noalign{\medskip} B_{\rm pre} \\ \noalign{\medskip} -B_{\rm pre} \\ \noalign{\medskip} \end{array}\right) \quad\mathrm{for}\quad x > 0.6\,,\qquad \]

while the remaining vector components are 0. The cloud has radius $R$ centered at $r_0 = (x_0,y_0,z_0) = (0.8, 0.5, 0.5)$ and larger density $\rho = 10$.

The interaction may be divided into two phases: 1) the collapse stage where the front of the cloud is strongly compressed and two fast shocks are generated and 2) the reexpansion phase which begins when the transmitted fast shock overtake the back of the cloud.

The runtime parameters that are read from pluto.ini are

A list of the available configurations is given in the following table:

Conf.GEOMETRY DIMT. STEPPINGRECONSTRUCTIONdivBAMR
#01 CARTESIAN 2 RK2 LINEAR CT NO
#02 CARTESIAN 2 ChTr LINEAR CT NO
#03 CARTESIAN 2 HANCOCK LINEAR 8W NO
#04 CARTESIAN 3 ChTr LINEAR CT NO
#05 CARTESIAN 3 RK3 LINEAR CT NO
#06 CARTESIAN 3 ChTr LINEAR GLMNO
#07 CARTESIAN 3 RK3 WENO3_FD GLMNO
#08 CARTESIAN 3 HANCOK LINEAR GLMYES
#09 CARTESIAN 2 RK2 LINEAR CT NO
mhd_shock_cloud.02.jpg
Density map with overplotted field lines for configuration #02
Author
A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)
Date
Sept 17, 2014

References:

Function Documentation

◆ 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().