cs_avg liion.cs_pce1 "Average intercalated lithium concentration"
theta_poor (cs_rich-cs_avg)/(cs_rich-cs_poor) "Volume fraction of lithium-poor particles"
E phis-phil "Electrode potential vs Li/Li+"
eta_rich E-Eeq_rich(cs_rich) "Overpotential, lithium-rich particles"
eta_poor E-Eeq_poor(cs_poor) "Overpotential, lithium-poor particles"
i0_poor 2*i0_ref*sqrt(cl/cl_0*cs_poor*(csmax-cs_poor)/csmax^2) "Exchange current density, lithium-poor particles"
i0_rich 2*i0_ref*sqrt(cl/cl_0*cs_rich*(csmax-cs_rich)/csmax^2) "Exchange current density, lithium-rich particles"
iloc_poor_expr i0_poor*(exp(0.5*F_const*eta_poor/(R_const*T))-exp(-0.5*F_const*eta_poor/(R_const*T))) "Local current density expression, lithium-poor particles"
loc_rich_expr i0_rich*(exp(0.5*F_const*eta_rich/(R_const*T))-exp(-0.5*F_const*eta_rich/(R_const*T))) "Local current density expression, lithium-rich particles"
iloc_poor "if(cs_avg<cs_avg_tp_max, iloc_poor_expr,0)" "Local current density, lithium-poor particles"
iloc_rich "if(cs_avg>cs_avg_tp_min, loc_rich_expr,0)" "Local current density, lithium-rich particles"
R_poor "if(cs_poor>cs_poor_max && eta_poor<0,0, -iloc_poor*Av/F_const)" "Lithium intercalation rate , lithium-poor particles"
R_rich "if(cs_rich<cs_rich_min && eta_rich>0,0, -iloc_rich*Av/F_const)" "Lithium intercalation rate, lithium-rich particles"
iloc theta_poor*iloc_poor+(1-theta_poor)*iloc_rich "Volume-averaged local current density"