L 0.02[m] "Cell length"
H_ch 1e-3[m] "Channel height"
W_ch 0.7874e-3[m] "Channel width"
W_rib 0.90932e-3[m] "Rib width"
H_gdl 380e-6[m] "GDL width"
H_electrode 50e-6[m] "Porous electrode thickness"
H_membrane 100e-6[m] "Membrane thickness"
eps_gdl 0.4 "GDL porosity"
kappa_gdl 1.18e-11[m^2] "GDL permeability"
sigma_gdl 222[S/m] "GDL electric conductivity"
mu_anode 1.19e-5[Pa*s] "Anode viscosity"
mu_cathode 2.46e-5[Pa*s] "Cathode viscosity"
MH2 0.002[kg/mol] "Hydrogen molar mass"
MN2 0.028[kg/mol] "Nitrogen molar mass"
MH2O 0.018[kg/mol] "Water molar mass"
MO2 0.032[kg/mol] "Oxygen molar mass"
D_H2_H2O 9.15e-5*(T/307.1[K])^1.75[m^2/s] "H2-H2O Binary diffusion coefficient"
D_N2_H2O 2.56e-5*(T/307.15[K])^1.75[m^2/s] "N2-H2O Binary diffusion coefficient"
D_O2_N2 2.2e-5*(T/293.2[K])^1.75[m^2/s] "O2-N2 binary diffusion coefficient"
D_O2_H2O 2.82e-5*(T/308.1[K])^1.75[m^2/s] "O2-H2O binary diffusion coefficient"
T 180+273.15[K] "Cell temperature"
p_ref 1[atm] "Reference pressure"
V_cell 0.95 "Cell voltage"
eps_l 0.3 "Electrolyte phase volume fraction"
eps_cl 1-eps_l-eps_gdl "Open volume fraction for gas diffusion in porous electrodes"
kappa_cl kappa_gdl/5 "Permeability (porous electrode)"
sigma_m 9.825[S/m] "Membrane conductivity"
alpha_c 1 "Transfer coefficient, cathode"
Eeq_ref_H2 -(0-(T-25[degC])*131[J/(K*mol)])/(2*F_const) "Reference equilibrium potential, anode"
Eeq_ref_O2 -(-2*229[kJ/mol]-(T-25[degC])*(2*189[J/(K*mol)]-205[J/(K*mol)]))/(4*F_const) "Reference equilibrium potential, cathode"
c_ref p_ref/(R_const*T) "Reference concentrations"
T_hum 28[degC] "Humidification temperature"
p_vapor 0.61121*exp((18.678-(T_hum-0[degC])/234.5[K])*((T_hum-0[degC])/(257.14+T_hum-0[degC])))[kPa] "Water partial pressure"
x_H2O_in p_vapor/p_ref "Inlet water molar fraction"
x_H2_in 1-x_H2O_in "Inlet hydrogen molar fraction"
x_O2_in 0.21*(1-x_H2O_in) "Inlet oxygen molar fraction"
I_stoich 1[A/cm^2]*L*(W_rib+W_ch) "Operating current for gas flow calculations"
lambda_anode 1.2 "Anode stoichiometry"
lambda_cathode 2 "Cathode stoichiometry"
U_in_anode lambda_anode*I_stoich/(2*F_const)*(R_const*T)/p_ref/x_H2_in/(H_ch*W_ch) "Inlet velocity, anode"
U_in_cathode lambda_cathode*I_stoich/(4*F_const)*(R_const*T)/p_ref/x_O2_in/(H_ch*W_ch) "Inlet velocity, cathode"
i0_ref_c 1e-3[A/m^2] "Reference exchange current density, cathode"
i0_ref_a 1e2[A/m^2] "Reference exchange current density, anode"
Av 1e7[m^2/m^3] "Specific surface area"