e11 ux Strain rate 11 component e12 0.5*(uy+vx) Strain rate 12 component e22 vy Strain rate 22 component s11 -p+2*mu_s*e11+T11 Total stress 11 component s12 2*mu_s*e12+T12 Total stress 12 component s22 -p+2*mu_s*e22+T22 Total stress 22 component f11es 2*mu_p*e11-T11-Wi*(u*T11x+v*T11y-2*ux*T11-2*uy*T12) Source term for T11 f12es 2*mu_p*e12-T12-Wi*(u*T12x+v*T12y-vx*T11-uy*T22) Source term for T12 f22es 2*mu_p*e22-T22-Wi*(u*T22x+v*T22y-2*vx*T12-2*vy*T22) Source term for T22 esweak tau*2*mu_p*(es11test*es11+2*es12test*es12+es22test*es22) Weak stabilization term for extra stress normU sqrt(u*u+v*v) Magnitude of velocity normgradU sqrt(ux*ux+uy*uy+vx*vx+vy*vy) Magnitude of velocity gradients tau nojac(1/sqrt(1+(2*Wi*normU/h)^2+(Wi*normgradU)^2)) Stabilization parameter for extra stress es11 T11+Wi*(u*T11x+v*T11y-2*ux*T11-2*uy*T12)-2*mu_p*e11 Stabilization term for T11 es11test test(T11)+Wi*(u*test(T11x)+v*test(T11y)-2*ux*test(T11)-2*uy*test(T12))-2*mu_p*test(e11) Stabilization test term for T11 es12 T12+Wi*(u*T12x+v*T12y-vx*T11-uy*T22)-2*mu_p*e12 Stabilization term for T12 es12test test(T12)+Wi*(u*test(T12x)+v*test(T12y)-vx*test(T11)-uy*test(T22))-2*mu_p*test(e12) Stabilization test term for T12 es22 T22+Wi*(u*T22x+v*T22y-2*vx*T12-2*vy*T22)-2*mu_p*e22 Stabilization term for T22 es22test test(T22)+Wi*(u*test(T22x)+v*test(T22y)-2*vx*test(T12)-2*vy*test(T22))-2*mu_p*test(e22) Stabilization test term for T22