/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) (4) (5) (6) (7) (8) (9) (10)
/INIMAP2D/form/inimap2d_ID
inimap2d_title
node_ID1 node_ID2 node_ID3              
grbric_ID grquad_ID grtria_ID              
f2d_IDv Fscalev              
f2d_IDvf1 f2d_IDrho1 Fscalerho1 f2d_IDp_e1 Fscalep_e1      
etc etc etc etc etc      
f2d_IDvfN f2d_IDrhoN FscalerhoN f2d_IDp_eN Fscalep_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_ID1 Node 1 identifier.

(Integer, maximum 10 digits)

 
node_ID2 Node 2 identifier.

(Integer, maximum 10 digits)

 
node_ID3 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_IDv Function identifier (/FUNC_2D) for velocity initialization.

(Integer)

 
Fscalev Velocity scale factor. 4

(Integer)

[ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@
Repeat the following lines for each /MAT/LAW151 or /MAT/LAW51 submaterial 1 through N
f2d_IDvf1 2D function identifier (/FUNC_2D) for volume fraction f v f ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciOzamaaBa aaleaacaWG2bGaamOzaaqabaGcdaqadaqaaiaadshaaiaawIcacaGL Paaaaaa@3B81@ .
= 0
Volume Fraction = F s c a l e v f MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeOvaiaab+ gacaqGSbGaaeyDaiaab2gacaqGLbGaaeiiaiaabAeacaqGYbGaaeyy aiaabogacaqG0bGaaeyAaiaab+gacaqGUbGaeyypa0JaamOraiaado hacaWGJbGaamyyaiaadYgacaWGLbWaaSbaaSqaaiaadAhacaWGMbaa beaaaaa@4BF6@
> 0
Volume Fraction = F s c a l e v f f v f ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeOvaiaab+ gacaqGSbGaaeyDaiaab2gacaqGLbGaaeiiaiaabAeacaqGYbGaaeyy aiaabogacaqG0bGaaeyAaiaab+gacaqGUbGaeyypa0JaamOraiaado hacaWGJbGaamyyaiaadYgacaWGLbWaaSbaaSqaaiaadAhacaWGMbaa beaakiabgwSixlGacAgadaWgaaWcbaGaamODaiaadAgaaeqaaOWaae WaaeaacaWG0baacaGLOaGaayzkaaaaaa@53D3@

(Integer)

 
f2d_IDrhoi 2D function identifier (/FUNC_2D) for density f r h o ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciOzamaaBa aaleaacaWGYbGaamiAaiaad+gaaeqaaOWaaeWaaeaacaWG0baacaGL OaGaayzkaaaaaa@3C73@ .
= 0
ρ = F s c a l e r h o MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdiNaey ypa0JaamOraiaadohacaWGJbGaamyyaiaadYgacaWGLbWaaSbaaSqa aiaadkhacaWGObGaam4Baaqabaaaaa@412D@
> 0
ρ = F s c a l e r h o f r h o ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdiNaey ypa0JaamOraiaadohacaWGJbGaamyyaiaadYgacaWGLbWaaSbaaSqa aiaadkhacaWGObGaam4BaaqabaGccqGHflY1ciGGMbWaaSbaaSqaai aadkhacaWGObGaam4BaaqabaGcdaqadaqaaiaadshaaiaawIcacaGL Paaaaaa@49FD@

(Integer)

 
Fscalerhoi Density scale factor. 4

(Real)

[ kg m 3 ]
f2d_IDp_ei If form=VE, 2D function identifier (/FUNC_2D) for energy f p _ e ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciOzamaaBa aaleaacaWGWbGaai4xaiaadwgaaeqaaOWaaeWaaeaacaWG0baacaGL OaGaayzkaaaaaa@3C5D@ .
If form=VP, 2D function identifier (/FUNC_2D) for pressure f p _ e ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciOzamaaBa aaleaacaWGWbGaai4xaiaadwgaaeqaaOWaaeWaaeaacaWG0baacaGL OaGaayzkaaaaaa@3C5D@ .
= 0
P  or  E = F s c a l e p _ e MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiuaiaabc cacaqGVbGaaeOCaiaabccacaWGfbGaeyypa0JaamOraiaadohacaWG JbGaamyyaiaadYgacaWGLbWaaSbaaSqaaiaadchacaGGFbGaamyzaa qabaaaaa@4423@
> 0
P  or  E = F s c a l e p _ e f p _ e ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiuaiaabc cacaqGVbGaaeOCaiaabccacaWGfbGaeyypa0JaamOraiaadohacaWG JbGaamyyaiaadYgacaWGLbWaaSbaaSqaaiaadchacaGGFbGaamyzaa qabaGccqGHflY1ciGGMbWaaSbaaSqaaiaadchacaGGFbGaamyzaaqa baGcdaqadaqaaiaadshaaiaawIcacaGLPaaaaaa@4CDD@

(Integer)

 
Fscalep_e1 If form=VE, Energy scale factor. 4

If form=VP, Pressure scale factor.

(Real)

[ Pa ] or [ 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

  1. node_ID1, node_ID2 and node_ID3 define the local coordinate system used for mapping.

    Where, node_ID1 amd node_ID2 define X ' MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCiwaiGacE caaaa@3784@

    node_ID1 and node_ID3 define Y " MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbc9v8 qqaqFr0xb9pg0xb9qqaqFn0dXdHiVcFbIOFHK8Feea0dXdar=Jb9hs 0dXdHuk9fr=xfr=xfrpeWZqaaiaaciWacmaadaGabiaaeaGaauaaaO qaaiaahMfacaGGIaaaaa@3BD0@
    • Z ' = X ' Λ Y " MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbc9v8 qqaqFr0xb9pg0xb9qqaqFn0dXdHiVcFbIOFHK8Feea0dXdar=Jb9hs 0dXdHuk9fr=xfr=xfrpeWZqaaiaaciWacmaadaGabiaaeaGaauaaaO qaaiaahQfacaGGNaGaeyypa0JaaCiwaiaacEcacqqHBoatcaWHzbGa aiOiaaaa@4165@
    • Y ' = Z Λ X MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=HhbHc9v8 qqaqFr0xb9pg0xb9qqaqFn0dXdHiVcFbIOFHK8Feea0dXdar=Jb9hs 0dXdHuk9fr=xfr=xfrpeWZqaaiaaciWacmaadaGabiaaeaGaauaaaO qaaiaahMfacaGGNaGaeyypa0JabCOwayaafaGaeu4MdWKabCiwayaa faaaaa@4034@


    Figure 1. Moving System with Dir=X

    Y ' MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCiwaiGacE caaaa@3784@ is the axis of symmetry for the 3D axi-symmetric mapped data.

  2. The coordinates of 2D data that will be mapped are defined in the local X ' MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCiwaiGacE caaaa@3784@ and Y ' MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCiwaiGacE caaaa@3784@ system using the /FUNC_2D X and Y input.
  3. 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).
  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, Fscalev, Fscalerhoi, or Fscalep_e1. If functions are defined, then the default value for the scale factors are 1.0.
  5. 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).
  6. 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 Iform= 12 is used.
  7. 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)