/INIMAP1D

Block Format Keyword Map 1D velocities and thermodynamic values into 2D or 3D space using spherical, cylindrical, or planar mapping methods. The mapped results are used with /ALE/EULER materials.

Format

(1)	(2)	(3)	(5)	(6)
/INIMAP1D/`form`/`inimap1d_ID`
`inimap1d_title`
`type`
`node_ID`₁	`node_ID`₂
`grbric_ID`	`grquad_ID`	`grtria_ID`
`fct_ID`_v	`Fscale`_v
`fct_ID`_vf1	`fct_ID`_rho1	`Fscale`_rho1	`fct_ID`_{p_e1}	`Fscale`_{p_e1}
etc	etc	etc	etc	etc
`fct_ID`_vfN	`fct_ID`_rhoN	`Fscale`_rhoN	`fct_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.
`inimap1d_ID`	Inimap1d block identifier. (Integer, maximum 10 digits)
`inimap1d_title`	Inimap1d block title. (Character, maximum 100 characters)
`type`	Initial mapping type. 1 Planar mapping 2 Cylindrical mapping 3 Spherical mapping (Integer)
`node_ID`₁	Node 1 identifier. (Integer, maximum 10 digits)
`node_ID`₂	Node 2 identifier, only if `type`=1 or 2. (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)
`fct_ID`_v	Function identifier for velocity initialization. (Integer)
`Fscale`_v	Velocity scale factor. 2 (Integer)	$[\frac{m}{s}]$
Repeat the following lines for each /MAT/LAW151 or /MAT/LAW51 submaterial 1 through N
`fct_ID`_vf1	Function $f_{v f} (t)$ identifier for volume fraction. = 0 $V o l u m e F r a c t i o n = F s c a l e_{v f}$ > 0 $V o l u m e f r a c t i o n = F s c a l e_{v f} \cdot f_{v f} (t)$ (Integer)
`fct_ID`_rhoi	Function $f_{r h o} (t)$ identifier for density. = 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. (Real)	$[\frac{kg}{m^{3}}]$
`fct_ID`_{p_ei}	If `form`=VE, Energy function identifier. If `form`=VP, Pressure function identifier. = 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. 2 If `form`=VP, Pressure scale factor. (Real)	$[Pa]$ or $[\frac{J}{m^{3}}]$

Example

#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/INIMAP1D/VE/1
INIMAP1D for cylindrical mapping
#     Type
         2
# Node_ID1  Node_ID2
         1         4
#grbric_ID grquad_ID grsh3n_ID
         1         0         0
#   func_u
         3
# func_alp  func_rho    func_e
         0         1                             2          
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/INIMAP1D/VE/2
INIMAP1D for planar mapping
#     Type
         1
# Node_ID1  Node_ID2
      2602    393299
#grbric_ID grquad_ID grsh3n_ID
         2         0         0
#   func_u
         3
# func_alp  func_rho    func_e
         0         1                             2          
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/INIMAP1D2/VE/3
INIMAP1D for spherical mapping
#     Type
         3
# Node_ID1		 
      5203
#grbric_ID grquad_ID grsh3n_ID
         3         0         0
#   func_u
         3
# func_alp  func_rho    func_e
         0         1                             2          
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/1
RHO
          0.00000000          1.00000000
          0.25000000          1.00000000
          0.50000000          0.25000000
          0.65000000          0.30000000
          0.70000000          0.20000000
          0.80000000          0.20000000
          0.81000000          0.12500000
          1.00000000          0.12500000
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/2
EINT
          0.00000000          2.50000000
          0.25000000          2.50000000
          0.50000000          1.40000000
          0.65000000          1.60000000
          0.70000000          2.50000000
          0.80000000          2.50000000
          0.81000000          2.00000000
          1.00000000          2.00000000
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/3
VELOCITY
          0.00000000          0.00000000
          0.25000000          0.00000000
          0.50000000          1.16000000
          0.65000000          0.80000000
          0.70000000          0.70000000
          0.80000000          0.60000000
          0.81000000          0.00000000
          1.00000000          0.00000000
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata

Comments

The option can be used to initialize a 2D or 3D simulation from the results of a 1D simulation potentially obtained from an external solver.
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.
Three different initialization types are available depending on the type parameter:
- type = 1: a planar initialization is performed where the 1D results are mapped onto a plane which is normal to the vector defined from node_ID₁ to node_ID₂. The abscissa values in an input function are the distance along the plane’s normal vector. The ordinate values from the input function are a constant value for the entire plane at a given distance along the normal. See Figure 1 and Figure 2.
- type = 2: a cylindrical initialization is performed where the local z cylindrical axis is defined from node_ID₁ to node_ID₂. The abscissa values in an input function are the radial distance from the local z cylindrical axis. The ordinate values from the input function are a constant value for the entire cylindrical surface at given radial distance. See Figure 3 and Figure 4.
- type = 3: a spherical initialization is performed, using the origin defined by node_ID₁. The abscissa values in an input function are the radial distance from origin. The ordinate values from the input function are a constant value for the entire spherical surface at given radial distance. See Figure 5 and Figure 6.
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 1. Mass Density (with input function) Planar Mapping on a 50 x 50 x 50 Cubic Mesh

Figure 2. Velocity (with input function) Planar Mapping on a 50 x 50 x 50 Cubic Mesh

Figure 3. Mass Density (with input function) Cylindrical Mapping on a 50 x 50 x 50 Cubic Mesh

Figure 4. Velocity (with input function) Cylindrical Mapping on a 50 x 50 x 50 Cubic Mesh

Figure 5. Mass Density (with input function) Spherical Mapping on a 50 x 50 x 50 Cubic Mesh

Figure 6. Velocity (with input function) Spherical Mapping on a 50 x 50 x 50 Cubic Mesh