Axisymmetric viscous source term bug in energy equation

[mskim0401]

Dear all,

In my opinion, there is a bug in axisymmetric viscous source term of energy equation.

From the equation D-55 of reference in PR#1192 , the energy equation should be

...
+v*TWO3*(2*PrimVar_Grad_i[2][1]-PrimVar_Grad_i[1][0]-v*yinv))
...

below is current code in void CSourceAxisymmetric_Flow::ResidualDiffusion()

residual[3] -= Volume*(yinv*(total_viscosity_i*(u*(PrimVar_Grad_i[2][0]+PrimVar_Grad_i[1][1])
                                                 +v*TWO3*(2*PrimVar_Grad_i[1][1]-PrimVar_Grad_i[1][0]
                                                 -v*yinv+U_i[0]*turb_ke_i))
                                                 +total_conductivity_i*PrimVar_Grad_i[0][1])
                                                 -TWO3*(AuxVar_Grad_i[1][1]+AuxVar_Grad_i[2][0]));

I appreciate if someone let me know correct thing if I miss something.
Thank you.

[FlorianDm]

Oh thank you very much for pointing this out. I don't know how it happened but I will address it right away

[FlorianDm]

fixed #1366