xA cet.x_Li/(cet.x_Li+cet.x_K) "Li2CO3 salt molar fraction"
D_Li_CO3 (5.73e-7-3.77e-7*xA)*exp(-4.5e4[J/mol]/(R_const*T))[m^2/s] "Maxwell-Stefan diffusivity"
D_K_CO3 (1.25e-8+4.86e-7*xA)*exp(-4.5e4[J/mol]/(R_const*T))[m^2/s] "Maxwell-Stefan diffusivity"
D_Li_K (3.13e-7-1e-7*xA)*exp(-4.49e4[J/mol]/(R_const*T))[m^2/s] "Maxwell-Stefan diffusivity"
i0_neg i0_neg_ref*xH2_neg^(0.5*beta_neg)*xH2O_neg^(0.5*(1-beta_neg))*xCO2_neg^(0.5*(1-beta_neg)) "Exchange current density"
phis_neg 0 "Electrode phase potential"
Eeq_neg Eeq_neg_ref-R_const*T/(2*F_const)*log(xH2_neg/(xH2O_neg*xCO2_neg)) "Equilibrium potential"
eta_neg phis_neg-cet.phil-Eeq_neg Overpotential
iloc_neg i0_neg*(exp((1-beta_neg)*F_const/(R_const*T)*eta_neg)-exp(-(1+beta_neg)*F_const/(R_const*T)*eta_neg)) "Local current density"
iv_neg iloc_neg*Av_neg "Volumetric current density"
i0_pos i0_pos_ref*xO2_pos^(0.5-0.25*beta_pos)*xCO2_pos^(-(1+0.5*beta_pos)) "Exchange current density"
phis_pos E_cell "Electrode phase potential"
Eeq_pos Eeq_pos_ref-R_const*T/(2*F_const)*log(1/(xCO2_pos*xO2_pos^0.5)) "Equilibrium potential"
eta_pos phis_pos-cet.phil-Eeq_pos Overpotential
iloc_pos i0_pos*(exp((1+beta_pos)*F_const/(R_const*T)*eta_pos)-exp(-(1-beta_pos)*F_const/(R_const*T)*eta_pos)) "Local current density"
iv_pos iloc_pos*Av_pos "Volumetric current density"