Struct PolynomialEos 

pub struct PolynomialEos {
    pub rho0: f64,
    pub a: [f64; 3],
    pub b: [f64; 3],
    pub t: [f64; 2],
}

Polynomial EOS commonly used in hydrocode simulations.

p = a₁μ + a₂μ² + a₃μ³ + (b₀ + b₁μ + b₂μ²)ρe (compressed) p = t₁μ + t₂μ² (expanded, e<0) where μ = ρ/ρ₀ - 1.

Fields

rho0: f64

Reference density ρ₀ (kg/m³).

a: [f64; 3]

Compression coefficients a₁, a₂, a₃ (Pa).

b: [f64; 3]

Energy coupling b₀, b₁, b₂.

t: [f64; 2]

Tension coefficients t₁, t₂ (Pa).

Implementations

impl PolynomialEos

pub fn new(rho0: f64, a: [f64; 3], b: [f64; 3], t: [f64; 2]) -> Self

Create a polynomial EOS.

pub fn mu(&self, density: f64) -> f64

Compression μ = ρ/ρ₀ - 1.

pub fn pressure_energy(&self, density: f64, energy: f64) -> f64

Pressure from density and energy.

Trait Implementations

impl Clone for PolynomialEos

fn clone(&self) -> PolynomialEos

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for PolynomialEos

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl EquationOfState for PolynomialEos

fn pressure(&self, density: f64) -> f64

Compute pressure from density (Pa).

fn sound_speed(&self, density: f64) -> f64

Compute sound speed from density (m/s).

fn density_from_pressure(&self, pressure: f64) -> f64

Compute density from pressure (inverse of pressure, if invertible).

impl Copy for PolynomialEos