E_aSi_Li a*E_Li+(1-a)*E_aSi "Young's modulus of LixSi"
nu_aSi_Li a*nu_Li+(1-a)*nu_aSi "Poisson's ratio of LixSi"
a xmax*c_bar/(xmax*c_bar+1) "Atomic fraction of Li atoms in a-Si"
sY "sYsat + (sY0-sYsat)*exp(-c_bar/cs)" "Yield strength"
c_bar c/cmax "Normalized concentration"
Vol0 pi*((Ri+W_Si)^2-Ri^2)*L0 "Initial volume of anode"
A0 2*pi*(Ri+W_Si)*L0 "Initial area of anode"
M D0*Omega0/(R_const*T0)*ts.c*(1-c_bar) "Mobility coefficient (spatial)"
Psi_ch "R_const*T0*cmax*(c_bar*log(c_bar)+(1-c_bar)*log(1-c_bar))+cmax*(a2*c_bar^2 +a3*c_bar^3+a4*c_bar^4+a5*c_bar^5+a6*c_bar^6+a7*c_bar^7)" "Chemical free energy"
SOC intV(c_bar)/intV(1) State-of-charge
Jsw 1+Omega0*(laginterp(1,c)-c_init) "Swelling volume ratio"
I -F_const*J_Li "Applied current"
I0 F_const*k0*sqrt(1-c_bar)*sqrt(c_bar) "Concentration-dependent exchange current"
V Veq+eta Voltage
Veq V0-ts.mu_c/F_const "Equilibrium potential"
eta (2*R_const*T0/F_const)*asinh(-0.5*I/I0) Overpotential
J -Vol0/A0*Crate/3600[s]*cmax "Lithium flux"
J_Li fluxDirection*J "Lithium flux with sign"
mu_def d(solid.Ws,c) "Deformation contribution to the chemical potential"