Loading [MathJax]/jax/output/CommonHTML/jax.js

/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        
ρl0 Cl αl νl νvol_l
ρg0 γ P0 νg νvol_g

Definition

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: Psh=P0 (Real)

[Pa]
ρl0 Liquid reference density.

(Real)

[kgm3]
Cl Liquid bulk modulus.

(Real)

[Pa]
αl Initial mass fraction liquid proportion.
= 0.0
Gas
= 1.0
Liquid

(Real)

 
νl Shear kinematic viscosity ( =μ/ρl0 ).

(Real)

[m2s]
νvol_l Volumetric viscosity (kinematic) for liquid: =3λ+2μρ

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

(Real)

[m2s]
ρg0 Reference gas density.

(Real)

[kgm3]
γ Perfect gas constant.

(Real)

 
P0 Reference gas initial pressure.

(Real)

[Pa]
νg Shear kinematic viscosity = =μ/ρg0 .

(Real)

[m2s]
νvol_g Volumetric viscosity (kinematic) for gas: =3λ+2μρ

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

(Real)

[m2s]

Example

Comments

  1. For each cycle equilibrium is computed between liquid and gas: ΔPl=ΔPg
    Liquid EOS:(1)
    Pl=P0+Clμ

    with bulk modulus Cl=ρl0(cl0)2

    Where,
    ρl0
    Initial liquid density
    cl0
    Initial sound speed
    Gas EOS:(2)
    Pvγ=P0vγ0=constant
    Where,
    P
    Pressure
    v=VV0
    Specific volume
    Stress Tension(3)
    σ=P+τ
    (4)
    τij=λ˙εkkδij+2μ˙εij
    Where,
    σ
    Stress tension
    P
    Spherical stress tension
    τ
    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 Psh=P0 . This means that computed pressure is shifted so that at initial time output pressure value relative to P0 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=P0
    Expansion: P=ΔP<0 P<P0
    Compression: P=ΔP>0 P>P0
  6. Internal Energy is not computed with LAW37. Nevertheless, adiabatic condition dE=PdV is used through polytropic process equation:
    (5)
    PVγ=const.
  7. LAW37 is not compatible with elementary boundary conditions (/EBCS).