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

Hydrodynamic jet propagation in 2D cylindrical coordinates. More...

Functions

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

Detailed Description

This problem considers the propagation of a hydrodynamic jet into a static uniform medium with constant density and pressure. The ambient density, in units of the jet density, is prescribed to be $ \rho_a = \eta $ where $\eta$ is the ambient/jet density ratio. The ambient pressure is $ p = 1/\Gamma $ when an IDEAL EoS is used or it is computed from temperature as $ p = p(\rho_a, T_a) $ for the PVTE_LAW EoS (here Ta is the ambient temperature). These values are set in Init() while the jet inflow is set through a user-defined boundary condition at the lower z-boundary. A simple top-hat injection nozzle is used.

The configuration is defined in terms of the following parameters:

defined in pluto.ini. The reference density and length are given by the jet density and radius while the reference velocity is 1 Km/s. The actual numerical values are needed only when using the PVTE_LAW EoS.

hd_jet.jpg
Pressure (left) and density (right) maps for configuration #01 at t=15
Author
A. Mignone (migno.nosp@m.ne@t.nosp@m.o.inf.nosp@m.n.it)
Date
Feb 25, 2019

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

◆ UserDefBoundary()

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

Assign user-defined boundary conditions in the lower boundary ghost zones. The profile is top-hat:

\[ V_{ij} = \left\{\begin{array}{ll} V_{\rm jet} & \quad\mathrm{for}\quad r_i < 1 \\ \noalign{\medskip} \mathrm{Reflect}(V) & \quad\mathrm{otherwise} \end{array}\right. \]

where $ V_{\rm jet} = (\rho,v,p)_{\rm jet} = (1,M,1/\Gamma)$ and M is the flow Mach number (the unit velocity is the jet sound speed, so $ v = M$).