/MAT/LAW37 (BIPHAS)

Block Format Keyword Describes the hydrodynamic bi-material liquid gas material. It is not recommended to use multi-material laws with Radioss single precision engine.

Format

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/MAT/LAW37/mat_ID/unit_ID or /MAT/BIPHAS/mat_ID/unit_ID
mat_title
Psh        
ρ 0 l C l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaWGSbaabeaaaaa@37DC@ α l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaadYgaaeqaaaaa@38B3@ ν l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyVd42aaS baaSqaaiaadYgaaeqaaaaa@38CC@ ν v o l _ l MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8akY=xipgYlh9vqqj=hEeei0xXdbb a9frFf0=yqFf0dbba91qpepeI8k8fiI+fsY=rqaqpepae9pg0Firpe pesP0xe9Fve9Fve9qapdbaGaaiGadiWaamaaceGaaqaacaqbaaGcba GaeqyVd42aaSbaaSqaaiaadAhacaWGVbGaamiBaiaaygW7caGGFbGa aGzaVlaadYgaaeqaaaaa@4382@
ρ 0 g MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdi3aa0 baaSqaaiaaicdaaeaacaWGNbaaaaaa@398A@ γ P 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiuamaaBa aaleaacaaIWaaabeaaaaa@37B2@ ν g MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyVd42aaS baaSqaaiaadYgaaeqaaaaa@38CC@ ν v o l _ g MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8akY=xipgYlh9vqqj=hEeei0xXdbb a9frFf0=yqFf0dbba91qpepeI8k8fiI+fsY=rqaqpepae9pg0Firpe pesP0xe9Fve9Fve9qapdbaGaaiGadiWaamaaceGaaqaacaqbaaGcba GaeqyVd42aaSbaaSqaaiaadAhacaWGVbGaamiBaiaaygW7caGGFbGa aGzaVlaadYgaaeqaaaaa@4382@

Definitions

Field Contents SI Unit Example
mat_ID Material identifier.

(Integer, maximum 10 digits)

 
unit_ID Unit Identifier

(Integer, maximum 10 digits)

 
mat_title Material title.

(Character, maximum 100 characters)

 
Psh Pressure shift. 5

Default: P s h = P 0 (Real)

[ Pa ]
ρ 0 l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdi3aa0 baaSqaaiaaicdaaeaacaWGSbaaaaaa@398F@ Liquid reference density.

(Real)

[ kg m 3 ]
C l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBa aaleaacaWGSbaabeaaaaa@37DC@ Liquid bulk modulus.

(Real)

[ Pa ]
α l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySde2aaS baaSqaaiaadYgaaeqaaaaa@38B3@ Initial mass fraction liquid proportion.
= 0.0
Gas
= 1.0
Liquid

(Real)

 
ν l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyVd42aaS baaSqaaiaadYgaaeqaaaaa@38CC@ Shear kinematic viscosity ( = μ / ρ 0 l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyypa0Jaeq iVd0Maai4laiabeg8aYnaaDaaaleaacaaIWaaabaGaamiBaaaaaaa@3CFE@ ).

(Real)

[ m 2 s ]
ν v o l _ l MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8akY=xipgYlh9vqqj=hEeei0xXdbb a9frFf0=yqFf0dbba91qpepeI8k8fiI+fsY=rqaqpepae9pg0Firpe pesP0xe9Fve9Fve9qapdbaGaaiGadiWaamaaceGaaqaacaqbaaGcba GaeqyVd42aaSbaaSqaaiaadAhacaWGVbGaamiBaiaaygW7caGGFbGa aGzaVlaadYgaaeqaaaaa@4382@ Volumetric viscosity (kinematic) for liquid: = 3 λ + 2 μ ρ

Default = 0: Stokes Hypothesis ( 3 λ + 2 μ = 0 )

(Real)

[ m 2 s ]
ρ 0 g MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdi3aa0 baaSqaaiaaicdaaeaacaWGNbaaaaaa@398A@ Reference gas density.

(Real)

[ kg m 3 ]
γ Perfect gas constant.

(Real)

 
P 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiuamaaBa aaleaacaaIWaaabeaaaaa@37B2@ Reference gas initial pressure.

(Real)

[ Pa ]
ν g MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyVd42aaS baaSqaaiaadYgaaeqaaaaa@38CC@ Shear kinematic viscosity = = μ / ρ 0 g MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyypa0Jaeq iVd0Maai4laiabeg8aYnaaDaaaleaacaaIWaaabaGaam4zaaaaaaa@3CF9@ .

(Real)

[ m 2 s ]
ν v o l _ g MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8akY=xipgYlh9vqqj=hEeei0xXdbb a9frFf0=yqFf0dbba91qpepeI8k8fiI+fsY=rqaqpepae9pg0Firpe pesP0xe9Fve9Fve9qapdbaGaaiGadiWaamaaceGaaqaacaqbaaGcba GaeqyVd42aaSbaaSqaaiaadAhacaWGVbGaamiBaiaaygW7caGGFbGa aGzaVlaadYgaaeqaaaaa@4382@ Volumetric viscosity (kinematic) for gas: = 3 λ + 2 μ ρ

Default = 0: Stokes Hypothesis ( 3 λ + 2 μ = 0 ).

(Real)

[ m 2 s ]

Example

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  1. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                   g                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW37/1/1
air (unit g_mm_ms)
#                Psh
                   0
#             RHO_l0                 C_l             ALPHA_l                NU_l            NU_VOL_l
                .001                2089                   0            1.787e-3                   0
#             RHO_G0               GAMMA                  P0                NU_g            NU_VOL_g
             1.22E-6                 1.4                  .1             1.48e-2                   
/ALE/MAT/1
#     Modif. factor.
                   0
/MAT/LAW37/2/1
water (unit g_mm_ms)
#                Psh
                   0
#             RHO_l0                 C_l             ALPHA_l                NU_l            NU_VOL_l
                .001                2089                   1            1.787e-3                   0
#             RHO_G0               GAMMA                  P0                NU_g            NU_VOL_g
             1.22E-6                 1.4                  .1             1.48e-2                   
/ALE/MAT/2
#     Modif. factor.
                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#enddata
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Comments

  1. For each cycle equilibrium is computed between liquid and gas: Δ P l = Δ P g
    Liquid EOS:(1)
    P l = P 0 + C l μ

    with bulk modulus C l = ρ 0 l ( c 0 l ) 2

    Where,
    ρ 0 l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdi3aa0 baaSqaaiaaicdaaeaacaWGSbaaaaaa@398F@
    Initial liquid density
    c 0 l
    Initial sound speed
    Gas EOS:(2)
    P v γ = P 0 v 0 γ = c o n s t a n t
    Where,
    P
    Pressure
    v = V V 0
    Specific volume
    Stress Tension(3)
    σ = P + τ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbc9v8 qqaqFr0xb9pg0xb9qqaqFn0dXdHiVcFbIOFHK8Feea0dXdar=Jb9hs 0dXdHuk9fr=xfr=xfrpeWZqaaiaaciWacmaadaGabiaaeaGaauaaaO qaaiaaho8acqGH9aqpcaWHqbGaey4kaSIaaCiXdaaa@3FA8@
    (4)
    τ i j = λ ε ˙ k k δ i j + 2 μ ε ˙ i j MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8akY=xipgYlh9vqqj=hEeei0xXdbb a9frFf0=yqFf0dbba91qpepeI8k8fiI+fsY=rqaqpepae9pg0Firpe pesP0xe9Fve9Fve9qapdbaGaaiGadiWaamaaceGaaqaacaqbaaGcba GaeqiXdq3aaSbaaSqaaiaadMgacaWGQbaabeaakiabg2da9iabeU7a Sjqbew7aLzaacaWaaSbaaSqaaiaadUgacaWGRbaabeaakiaayIW7cq aH0oazdaWgaaWcbaGaamyAaiaadQgaaeqaaOGaey4kaSIaaGOmaiab eY7aTjqbew7aLzaacaWaaSbaaSqaaiaadMgacaWGQbaabeaaaaa@5084@
    Where,
    σ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbc9v8 qqaqFr0xb9pg0xb9qqaqFn0dXdHiVcFbIOFHK8Feea0dXdar=Jb9hs 0dXdHuk9fr=xfr=xfrpeWZqaaiaaciWacmaadaGabiaaeaGaauaaaO qaaiaaho8aaaa@3B97@
    Stress tension
    P MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbc9v8 qqaqFr0xb9pg0xb9qqaqFn0dXdHiVcFbIOFHK8Feea0dXdar=Jb9hs 0dXdHuk9fr=xfr=xfrpeWZqaaiaaciWacmaadaGabiaaeaGaauaaaO qaaiaahcfaaaa@3B21@
    Spherical stress tension
    τ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbc9v8 qqaqFr0xb9pg0xb9qqaqFn0dXdHiVcFbIOFHK8Feea0dXdar=Jb9hs 0dXdHuk9fr=xfr=xfrpeWZqaaiaaciWacmaadaGabiaaeaGaauaaaO qaaiaahs8aaaa@3B98@
    Viscous stress tension
    λ
    Shear viscosity (kinematic)
    μ
    Dilatation viscosity (kinematic)
  2. Multi-material LAW51 (/MAT/MULTIMAT) can emulate LAW37.

    Anti-diffusive technique using MUSCL method is under implementation. It is already available for LAW51 (/ALE/MUSCL).

  3. If inlet/outlet are needed LAW51 must be used instead. LAW51 has multi- material inlet/outlet formulations.
  4. In Animation files:

    /ANIM/ELEM/DENS37 enables to output submaterial densities.

    /ANIM/BRIC/VFRAC allows to output volumetric fractions for each submaterial.

    Refer to Time Histories in the Radioss Starter Input Manual, or to /ANIM/Eltyp/Restype in the Radioss Engine Input Manual.

  5. By Default P s h = P 0 . This means that computed pressure is shifted so that at initial time output pressure value relative to P 0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiuamaaBa aaleaacaaIWaaabeaaaaa@37B2@ value is 0.0.
    To use total pressure, Psh must be nullify using a small value such as -1e-20 Pressure P could be exported in animation with /ANIM/ELEM/P.
      Relative Pressure

    Psh=-P0 (Default)

    Total Pressure

    Psh = -1e-20

    Initial Condition: P = Δ P = 0 P = P 0
    Expansion: P = Δ P < 0 P < P 0
    Compression: P = Δ P > 0 P > P 0
  6. Internal Energy is not computed with LAW37. Nevertheless, adiabatic condition d E = P d V is used through polytropic process equation:
    (5)
    P V γ = c o n s t .
  7. LAW37 is not compatible with elementary boundary conditions (/EBCS).