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

Stellar wind test problem. More...

Functions

void InitDomain (Data *d, Grid *grid)
 

Detailed Description

Sets initial condition for a spherically symmetric radial wind blowing from the origin of coordinates, see [Mig14]. The initial condition consists of a constant-density ambient medium with

\[ \rho = \rho_a \,,\qquad p = p_a\,,\qquad \vec{v} = v_{csm}\hvec{j} \]

where v_csm is the velocity of the star with respect to the background. The wind is injeted using the INTERNAL_BOUNDARY where flow quantities are kept constant in time and equal to

\[ r^2v_r\rho = \rho_0 V_0^2r_0^2 \,,\qquad v_r = V_0\hvec{r} \,,\qquad p = \frac{c_s^2}{\Gamma}\rho^\Gamma \]

These value are defined through the UserDefBoundary() function when side is equal to 0.

Dimensions are chosen so that the spherical wind shell has radius 1, density 1 and velocity 1 ( $ r_0 = V_0 = \rho_0 = 1$).

The input parameters that control the problem dynamics are

Configurations #01-05 and #07-08 work in 2D cylindrical axisymmetric coordinates while conf. #06 is 3D Cartesian. An AMR setup is available with configuration #04.

hd_stellar_wind.08.jpg
Density map at the end of computation for configuration #8
Author
A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it)
Date
July 12, 2019

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