step_variable (max(cMg,eps)*max(cOH,eps)^2/ksp)-1 "Step variable for precipitation reaction"
R_MgOH2 k*(max(cMg,eps)*(max(cOH,eps))^2-ksp)*step1(step_variable) "Precipitation reaction rate"
R_Mg R_MgOH2 "Reaction rate of Mg ions"
R_OH 2*R_MgOH2 "Reaction rate of OH ions"
vn_corprod R_MgOH2*M_MgOH2/rho_MgOH2 "Corrosion product velocity"
u_corprod vn_corprod*ls.intnormx*step1(0.6-phils) "Corrosion product velocity, x-component"
v_corprod max(vn_corprod*ls.intnormy+genext1(cd.vn*cd.ny),eps)*step1(0.6-phils) "Corrosion product velocity, y-component"