/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.
(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
- 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=constantWhere,- P
- Pressure
- v=VV0
- Specific volume
Stress Tension(3) σ=P+τ(4) τij=λ˙εkk δij+2μ˙εijWhere,- σ
- Stress tension
- P
- Spherical stress tension
- τ
- Viscous stress tension
- λ
- Shear viscosity (kinematic)
- μ
- Dilatation viscosity (kinematic)
- 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).
- If inlet/outlet are needed LAW51 must be used instead. LAW51 has multi- material inlet/outlet formulations.
- 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.
- 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 - Internal Energy is not computed with LAW37. Nevertheless, adiabatic
condition
dE=−P⋅dV
is used through polytropic process
equation:
(5) PVγ=const. - LAW37 is not compatible with elementary boundary conditions (/EBCS).