/INIMAP2D

Block Format Keyword Map 2D velocity and thermodynamic values into a 3D axi-symmetric space for use with /ALE/EULER materials.

Format

(1)	(2)	(3)	(5)	(6)
/INIMAP2D/`form`/`inimap2d_ID`
`inimap2d_title`
`node_ID`₁	`node_ID`₂	`node_ID`₃
`grbric_ID`	`grquad_ID`	`grtria_ID`
`f2d_ID`_v	`Fscale`_v
`f2d_ID`_vf1	`f2d_ID`_rho1	`Fscale`_rho1	`f2d_ID`_{p_e1}	`Fscale`_{p_e1}
etc	etc	etc	etc	etc
`f2d_ID`_vfN	`f2d_ID`_rhoN	`Fscale`_rhoN	`f2d_ID`_{p_eN}	`Fscale`_{p_eN}

Definitions

Field	Contents	SI Unit Example
`form`	Initial state formulation. VE Mass density, specific internal energy and velocity are provided. VP Mass density, pressure and velocity are provided.
`inimap2d_ID`	Inimap2d block identifier. (Integer, maximum 10 digits)
`inimap2d_title`	Inimap2d block title. (Character, maximum 100 characters)
`node_ID`₁	Node 1 identifier. (Integer, maximum 10 digits)
`node_ID`₂	Node 2 identifier. (Integer, maximum 10 digits)
`node_ID`₃	Node 3 identifier. (Integer, maximum 10 digits)
`grbric_ID`	Brick group on which initialization is performed. (Integer)
`grquad_ID`	Quad group on which initialization is performed. (Integer)
`grtria_ID`	Tria group on which initialization is performed. (Integer)
`f2d_ID`_v	Function identifier (/FUNC_2D) for velocity initialization. (Integer)
`Fscale`_v	Velocity scale factor. 4 (Integer)	$[\frac{m}{s}]$
Repeat the following lines for each /MAT/LAW151 or /MAT/LAW51 submaterial 1 through N
`f2d_ID`_vf1	2D function identifier (/FUNC_2D) for volume fraction $f_{v f} (t)$ . = 0 $Volume Fraction = F s c a l e_{v f}$ > 0 $Volume Fraction = F s c a l e_{v f} \cdot f_{v f} (t)$ (Integer)
`f2d_ID`_rhoi	2D function identifier (/FUNC_2D) for density $f_{r h o} (t)$ . = 0 $ρ = F s c a l e_{r h o}$ > 0 $ρ = F s c a l e_{r h o} \cdot f_{r h o} (t)$ (Integer)
`Fscale`_rhoi	Density scale factor. 4 (Real)	$[\frac{kg}{m^{3}}]$
`f2d_ID`_{p_ei}	If `form`=VE, 2D function identifier (/FUNC_2D) for energy $f_{p_e} (t)$ . If `form`=VP, 2D function identifier (/FUNC_2D) for pressure $f_{p_e} (t)$ . = 0 $P or E = F s c a l e_{p_e}$ > 0 $P or E = F s c a l e_{p_e} \cdot f_{p_e} (t)$ (Integer)
`Fscale`_{p_e1}	If `form`=VE, Energy scale factor. 4 If `form`=VP, Pressure scale factor. (Real)	$[Pa]$ or $[\frac{J}{m^{3}}]$

Example

/INIMAP2D/VE/1
INIMAP2D1
# node_ID1  node_ID2  node_ID3
    130602    131262     96789
#grbric_ID grquad_ID grsh3n_ID
         1         0         0
#    fct_u
         3
#  fct_alp   fct_rho                      fct_eint
         0         1                             0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNC_2D/1
for density initialization
#      dim
         1
#                  X                   Y                  Z1
#            X-coord             Y-Coord             density 
          0.00000000          0.00000000        109.00000000
          0.00195695          0.00000000        114.00000000
          0.00391389          0.00000000        111.00000000
          0.00587084          0.00000000        109.00000000
          0.00782779          0.00000000        114.00000000
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNC_2D/3
for velocity initialization
#      dim
         2
#                  X                   Y                  Z1                  Z2
#            X-coord             Y-Coord          velocity x          velocity y
          0.00097656          0.00097656         -0.00026680         -0.02174114
          0.00097656          0.00292969         -0.00080037         -0.02173459
          0.00097656          0.00488281         -0.00133382         -0.02172149
          0.00097656          0.00683594         -0.00186707         -0.02170185
          0.00097656          0.00878906         -0.00240003         -0.02167567
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata

Comments

node_ID₁, node_ID₂ and node_ID₃ define the local coordinate system used for mapping.
Where, node_ID₁ amd node_ID₂ define $X'$
node_ID₁ and node_ID₃ define $Y "$
- $Z' = X' Λ Y "$
- $Y' = Z^{'} Λ X^{'}$
Figure 1. Moving System with Dir=X
$Y'$ is the axis of symmetry for the 3D axi-symmetric mapped data.
The coordinates of 2D data that will be mapped are defined in the local $X'$ and $Y'$ system using the /FUNC_2D X and Y input.
When mapping scalar quantities such as density, pressure and energy, use /FUNC_2D with dim=1 and define the scalar value using Z1. For vector quantities such as velocity, use /FUNC_2D with dim=2 and define the vector using Z1 and Z2 (Figure 3 and Figure 4).
If any inputs for veloctiy, volume fraction, density, energy, or pressure are constant, then the function can be left undefined and the constant values can be entered in the scale factor fields, Fscale_v, Fscale_rhoi, or Fscale_{p_e1}. If functions are defined, then the default value for the scale factors are 1.0.
When form=VE, the model is initialized from density, specific internal energy and velocity. This formulation is available for all mono-material ALE / EULER material laws and with material /MAT/LAW151 (MULTIFLUID) and /MAT/LAW51 (MULTIMAT).
When form=VP, the model is initialized from density, pressure and velocity. This formulation is available for all mono-material ALE / EULER material laws whose equation of state is provided through the /EOS card and with material /MAT/LAW151 and /MAT/LAW51 only if I_form= 12 is used.
If explosive material LAW5 (JWL) is used as submaterial, then mapping will be proceeded considering all explosive was burnt (burn fraction = 1.0). The mapping will be done with detonation products only with a burn fraction = 1.0.

Figure 2. 2D Scalar Input

Figure 3. Resulting Mapping onto a 50 x 50 x 50 Cubic Mesh

Figure 4. 2D Scalar Velocity Profile

Figure 5. Resulting mapping on a 50 x 50 x 50 cubic mesh (bottom)