/EOS/PUFF

Format

(1)	(3)	(5)	(7)
/EOS/PUFF/`mat_ID`/`unit_ID`
`eos_title`
`C`₁	`C`₂	`C`₃	$γ_{0}$
`T`₁	`T`₂	`E`_s
`H`	`E`₀	$ρ_{0}$

Definitions

Field	Contents	SI Unit Example
`mat_ID`	Material identifier. (Integer, maximum 10 digits)
`unit_ID`	Unit Identifier. (Integer, maximum 10 digits)
`eos_title`	EOS title. (Character, maximum 100 characters)
`C`₁	`C`₁ coefficient. 1 (Real)	$[Pa]$
`C`₂	`C`₂ coefficient. (Real)	$[Pa]$
`C`₃	`C`₃ coefficient. (Real)	$[Pa]$
$γ_{0}$	Gruneisen gamma. (Real)
`T`₁	`T`₁ coefficient. 1 Default= `C`₁ (Real)	$[Pa]$
`T`₂	`T`₂ coefficient. (Real)	$[Pa]$
`E`_s	Sublimation energy per unit initial volume. (Real)	$[\frac{J}{m^{3}}]$
`H`	`H` coefficient. 1 (Real)
`E`₀	Initial energy per unit initial volume. (Real)	$[\frac{J}{m^{3}}]$
$ρ_{0}$	Reference density. Default = material density (Real)	$[\frac{kg}{m^{3}}]$

Example (Copper)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/UNIT/1
unit for mat
                   g                  cm                 mus
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/HYD_JCOOK/1/1
Copper
#              RHO_I               RHO_0
                 8.9                   0
#                  E                  nu
                 1.3                 .34
#                  A                   B                   n              epsmax              sigmax
                9E-4              .00292                 .31                   0                   0
#               Pmin
                   0
#                  C           EPS_DOT_0                   M               Tmelt                Tmax
                .025                1E-6                1.09                1356                1E30
#              RHOCP                                                         T_r
            3.432E-5                                                           0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/EOS/PUFF/1/1
Copper
#                 C1                  C2                  C3                  G0
              1.3597              2.7466              2.4479                1.97
#                 T1                  T2                  ES
                   0                   0             0.12282
#                  H                  E0               RHO_0
                   1             1.25675                   0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA

Comments

Let $μ = \frac{ρ}{ρ_{0}} - 1$ , $η = 1 + μ$ , $x = 1 - \frac{ρ_{0}}{ρ}$ and E, the internal energy per unit initial volume, the pressure is defined by:
Region 1: $μ \geq 0$ (1)
$P = (C_{1} μ + C_{2} μ^{2} + C_{3} μ^{3}) (1 - \frac{γ μ}{2}) + γ (1 + μ) E$

with $γ ρ = γ_{0} ρ_{0}$

Region 2: $μ < 0$ and $E \geq E_{s}$ (2)
$P = (T_{1} μ + T_{2} μ^{2}) (1 - \frac{γ μ}{2}) + γ (1 + μ) E$

Region 3: $μ < 0$ and $E \geq E_{s}$ (3)
$P = η [H + (γ_{0} - H) \sqrt{η}] [E - E_{s} (1 - \exp (\frac{N (η - 1)}{η^{2}}))]$

with $N = \frac{C_{1} η}{γ_{0} E_{s}}$
Equations of state are used by Radioss to compute the hydrodynamic pressure and are compatible with the material laws:
- /MAT/LAW3 (HYDPLA)
- /MAT/LAW4 (HYD_JCOOK)
- /MAT/LAW6 (HYDRO or HYD_VISC)
- /MAT/LAW10 (DPRAG1)
- /MAT/LAW12 (3D_COMP)
- /MAT/LAW49 (STEINB)
- /MAT/LAW102 (DPRAG2)
- /MAT/LAW103 (HENSEL-SPITTEL)