/INIMAP1D

Format

(1)	(2)	(3)
/INIMAP1D/`form`/`inimap1d_ID`
`inimap1d_title`
`type`
`node_ID`₁	`node_ID`₂
`grbric_ID`	`grquad_ID`	`grtria_ID`
`fct_ID`₁	`fct_ID`₂	`fct_ID`₃

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. 2 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`₁	Function identifier for density initialization. (Integer)
`fct_ID`₂	Function identifier for specific internal energy initialization (if `form`=VE). Pressure initialization (if `form`=VP). (Integer)
`fct_ID`₃	Function identifier for velocity initialization. (Integer)

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 grtria_ID
         1         0         0
#  fct_ID1   fct_ID2   fct_ID3
         1         2         3  
#---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 grtria_ID
         2         0         0
#  fct_ID1   fct_ID2   fct_ID3
         1         2         3  
#---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 grtria_ID
         3         0         0
#  fct_ID1   fct_ID2   fct_ID3
         1         2         3               
#---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.
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) when only one sub-material is present.
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 when only one sub-material is present.

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