Struct FuelCellStack 

pub struct FuelCellStack {
    pub e_rev: f64,
    pub i0: f64,
    pub tafel_slope: f64,
    pub r_ohmic: f64,
    pub i_lim: f64,
    pub m_conc: f64,
    pub n_conc: f64,
}

Hydrogen PEM fuel cell polarisation curve model.

Total overpotential = activation loss + ohmic loss + concentration loss.

Fields

e_rev: f64

Reversible (Nernst) cell voltage [V]

i0: f64

Exchange current density i₀ [A/m²]

tafel_slope: f64

Tafel slope for activation loss [V/decade] (positive value)

r_ohmic: f64

Ohmic area-specific resistance [Ω·m²]

i_lim: f64

Limiting current density i_L [A/m²]

m_conc: f64

Concentration loss empirical coefficient m [V]

n_conc: f64

Concentration loss empirical coefficient n [m²/A]

Implementations

impl FuelCellStack

pub fn new( e_rev: f64, i0: f64, tafel_slope: f64, r_ohmic: f64, i_lim: f64, m_conc: f64, n_conc: f64, ) -> Self

Create a new fuel cell stack model.

pub fn activation_loss(&self, j: f64) -> f64

Activation overpotential [V] at current density j [A/m²].

η_act = tafel_slope · log10(j / i0)

pub fn ohmic_loss(&self, j: f64) -> f64

Ohmic overpotential [V] at current density j [A/m²].

η_ohm = R_ohmic · j

pub fn concentration_loss(&self, j: f64) -> f64

Concentration overpotential [V] at current density j [A/m²].

η_conc = m · exp(n · j)

pub fn cell_voltage(&self, j: f64) -> f64

Cell terminal voltage [V] at current density j [A/m²].

V = E_rev − η_act − η_ohm − η_conc

Returns 0.0 if j ≥ i_lim (cell has stalled).

pub fn power_density(&self, j: f64) -> f64

Power density [W/m²] at current density j [A/m²].

pub fn efficiency(&self, j: f64) -> f64

Efficiency relative to reversible voltage.

η = V_cell / E_rev

Trait Implementations

impl Clone for FuelCellStack

fn clone(&self) -> FuelCellStack

Returns a duplicate of the value.

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

Performs copy-assignment from source.

impl Debug for FuelCellStack

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

Formats the value using the given formatter.