Struct EnthalpyMethod 

pub struct EnthalpyMethod {
    pub rho: f64,
    pub cp: f64,
    pub latent_heat: f64,
    pub t_solidus: f64,
    pub t_liquidus: f64,
    pub temperature: f64,
    pub enthalpy: f64,
}

Enthalpy method for melting and solidification problems.

The total enthalpy H at temperature T includes sensible and latent heat:

H(T) = ρ·cp·T + ρ·L·f_l(T)

where f_l is the liquid fraction. The temperature-update is performed by inverting H(T) given the updated enthalpy.

Fields

rho: f64

Density [kg/m³].

cp: f64

Specific heat capacity [J/(kg·K)].

latent_heat: f64

Latent heat of fusion L [J/kg].

t_solidus: f64

Solidus temperature T_s [K].

t_liquidus: f64

Liquidus temperature T_l [K].

temperature: f64

Current temperature [K].

enthalpy: f64

Current enthalpy [J/m³].

Implementations

impl EnthalpyMethod

pub fn new( rho: f64, cp: f64, latent_heat: f64, t_solidus: f64, t_liquidus: f64, ) -> Self

Create a new enthalpy method model.

pub fn liquid_fraction(&self, temp: f64) -> f64

Liquid fraction as a linear function of temperature in mushy zone.

pub fn enthalpy_at(&self, temp: f64) -> f64

Total enthalpy H(T) [J/m³].

pub fn temperature_from_enthalpy(&self, h: f64) -> f64

Invert enthalpy to temperature (Newton-Raphson, max 50 iterations).

pub fn effective_cp(&self, temp: f64) -> f64

Effective (apparent) specific heat cp_eff [J/(kg·K)] including latent heat.

pub fn apply_heat_flux(&mut self, q: f64, dt: f64)

Apply a heat flux q [W/m³] over time step dt [s].

pub fn is_mushy(&self) -> bool

Check if material is currently in the mushy zone.

Trait Implementations

impl Clone for EnthalpyMethod

fn clone(&self) -> EnthalpyMethod

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for EnthalpyMethod

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

Formats the value using the given formatter.