Struct PhaseChangeModel 

pub struct PhaseChangeModel {
    pub t_liquidus: f64,
    pub t_solidus: f64,
    pub latent_heat: f64,
    pub cp_solid: f64,
    pub cp_liquid: f64,
    pub density: f64,
}

Phase change (melting/solidification) model with latent heat.

Implements the enthalpy method for solving phase change problems. The apparent heat capacity is increased in the mushy zone to account for the latent heat release/absorption.

Fields

t_liquidus: f64

Liquidus temperature (melting point or end of mushy zone) [K].

t_solidus: f64

Solidus temperature (start of mushy zone) [K].

latent_heat: f64

Latent heat of fusion L [J/kg].

cp_solid: f64

Solid-phase specific heat [J/(kg·K)].

cp_liquid: f64

Liquid-phase specific heat [J/(kg·K)].

density: f64

Material density [kg/m³].

Implementations

impl PhaseChangeModel

pub fn new( t_solidus: f64, t_liquidus: f64, latent_heat: f64, cp_solid: f64, cp_liquid: f64, density: f64, ) -> Self

Create a new phase change model.

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

Liquid fraction at temperature T (linear in mushy zone).

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

Apparent (effective) specific heat in the enthalpy method [J/(kg·K)].

In the mushy zone, cp_eff = cp + L/(T_liq - T_sol).

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

Specific enthalpy h(T) relative to T = 0.

pub fn conductivity(&self, k_solid: f64, k_liquid: f64, temp: f64) -> f64

Thermal conductivity (interpolated between solid and liquid).