EOS_GRUNEISEN
Material properties
*EOS_GRUNEISEN
eosid, $S_1$, $\Gamma$, $L$, $p_{cut}$
$S_2$, $S_3$, $a$
eosid, $S_1$, $\Gamma$, $L$, $p_{cut}$
$S_2$, $S_3$, $a$
Parameter definition
Variable
Description
eosid
Unique EOS identification number
$S_1$
Linear Hugoniot slope coefficient
$\Gamma$
Gruneisen gamma
$L$
Characteristic element size in shock viscosity calculation
$p_{cut}$
Pressure cut-off
$S_2$
Quadratic coefficient ($\gamma$SPH only)
$S_3$
Cubic coefficient ($\gamma$SPH only)
$a$
Gruneisen coefficient ($\gamma$SPH only)
Description
This is the Mie-Gruneisen equation-of-state. Note that the linear bulk modulus, $K$, is determined from the elastic properties in the material command.
$p = \displaystyle{ \frac{K \eta}{(1-S_1 \eta)^2} \cdot \left( 1 - \frac{\Gamma \eta}{2} \right) + \Gamma \rho_0 e }$
where $\rho$ is the current density, $\rho_0$ is the initial density and $e$ is the specific internal energy. $\eta$ is a measure of the volumetric compression.
$\displaystyle{ \eta = 1 - \frac{\rho_0}{\rho}}$
Only $\gamma$SPH:
Pressure for compressed material:
$\displaystyle{ p = \frac{\rho_0 c^2 \mu \left[ 1 + (1 - \frac{\Gamma}{2})\mu - \frac{a}{2}\mu^2 \right]} {1 - (S_1 - 1)\mu - \frac{S_2 \mu^2}{\mu + 1} - \frac{S_3 \mu^3}{(\mu + 1)^2}} + (\Gamma + a\mu)E}$
Pressure for expanded material:
$\displaystyle{ p = \rho_0 c^2 \mu + (\Gamma + a\mu) E}$
where $\rho$ is the current density, $\rho_0$ is the initial density and $E$ is the internal energy. $\mu$ is a measure of the volumetric compression.
$\displaystyle{ \mu = \frac{\rho}{\rho_0} - 1}$
