HL update

HL
su2code · Jan 14, 2025 · a583d7c · a583d7c
1 parent 448f5b9
commit a583d7c
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 8 deletions.
diff --git a/SU2_CFD/include/numerics/turbulent/transition/trans_sources.hpp b/SU2_CFD/include/numerics/turbulent/transition/trans_sources.hpp
@@ -422,7 +422,7 @@ class CSourcePieceWise_TransAFT final : public CNumerics {
         su2double tt = 0;
       }
 
-      const su2double HL = dist_i * dist_i / Laminar_Viscosity_i * HLGradTerm;
+      const su2double HL = dist_i * dist_i * Density_i / Laminar_Viscosity_i  * HLGradTerm;
       /*--- Cal H12, dNdRet, Ret0, D_H12, l_H12, m_H12, kv ---*/
       const su2double H12 = TransCorrelations.H12_Correlations(HL);
       const su2double dNdRet = TransCorrelations.dNdRet_Correlations(H12);
@@ -463,7 +463,7 @@ class CSourcePieceWise_TransAFT final : public CNumerics {
       Jacobian_i[0][0] = 0.0;
       Jacobian_i[0][1] = 0.0;
       Jacobian_i[1][0] = 0.0;
-      Jacobian_i[1][1] = ( -c_1 * StrainMag_i * F_onset * exp(TransVar_i[1]) -c_2 * VorticityMag * F_turb * c_3 * exp(TransVar_i[1]) ) * Volume;
+      Jacobian_i[1][1] = (-c_1 * StrainMag_i * F_onset * exp(TransVar_i[1]) - c_2 * VorticityMag * F_turb * c_3 * exp(TransVar_i[1])) * Volume;
     }
 
     AD::SetPreaccOut(Residual, nVar);

diff --git a/SU2_CFD/src/solvers/CTransAFTSolver.cpp b/SU2_CFD/src/solvers/CTransAFTSolver.cpp
@@ -202,14 +202,14 @@ void CTransAFTSolver::Postprocessing(CGeometry *geometry, CSolver **solver_conta
     const su2double c_3 = 50.0;
     su2double Temp1 = nodes->GetAuxVarGradient(iPoint, 0, 0);
     su2double Temp2 = nodes->GetAuxVarGradient(iPoint, 0, 1);
-    su2double Temp3 = flowNodes->GetDensity(iPoint) * flowNodes->GetVelocity(iPoint, 0) * nodes->GetAuxVarGradient(iPoint, 0, 0);
-    su2double Temp4 = flowNodes->GetDensity(iPoint) * flowNodes->GetVelocity(iPoint, 1) * nodes->GetAuxVarGradient(iPoint, 0, 1);
+    su2double Temp3 = flowNodes->GetVelocity(iPoint, 0) * nodes->GetAuxVarGradient(iPoint, 0, 0);
+    su2double Temp4 = flowNodes->GetVelocity(iPoint, 1) * nodes->GetAuxVarGradient(iPoint, 0, 1);
 
     if(nDim == 2) {
       nodes->SetAuxVar(iPoint, 1, Temp3 + Temp4);
     }
     else {
-      nodes->SetAuxVar(iPoint, 1, Temp3 + Temp4 + flowNodes->GetDensity(iPoint) * flowNodes->GetVelocity(iPoint, 2) * nodes->GetAuxVarGradient(iPoint, 0, 2));
+      nodes->SetAuxVar(iPoint, 1, Temp3 + Temp4 + flowNodes->GetVelocity(iPoint, 2) * nodes->GetAuxVarGradient(iPoint, 0, 2));
     }
 
     su2double VorticityMag = GeometryToolbox::Norm(3, flowNodes->GetVorticity(iPoint));
@@ -221,7 +221,7 @@ void CTransAFTSolver::Postprocessing(CGeometry *geometry, CSolver **solver_conta
     }
 
     /*--- Cal H12, Hk, dNdRet, Ret0 ---*/
-    const su2double HL = dist_i * dist_i / Laminar_Viscosity_i * HLGradTerm;
+    const su2double HL = dist_i * dist_i * Density_i / Laminar_Viscosity_i * HLGradTerm;
     const su2double H12 = TransCorrelations.H12_Correlations(HL);
     const su2double dNdRet = TransCorrelations.dNdRet_Correlations(H12);
     const su2double Ret0 = TransCorrelations.Ret0_Correlations(H12);
@@ -337,8 +337,8 @@ void CTransAFTSolver::Source_Residual(CGeometry *geometry, CSolver **solver_cont
     numerics->SetCoord(geometry->nodes->GetCoord(iPoint), nullptr);
     numerics->SetAuxVarGrad(nodes->GetAuxVarGradient(iPoint), nullptr);
 
-    su2double Temp1 = flowNodes->GetDensity(iPoint) * flowNodes->GetVelocity(iPoint, 0) * nodes->GetAuxVarGradient(iPoint, 0, 0);
-    su2double Temp2 = flowNodes->GetDensity(iPoint) * flowNodes->GetVelocity(iPoint, 1) * nodes->GetAuxVarGradient(iPoint, 0, 1);
+    su2double Temp1 = flowNodes->GetVelocity(iPoint, 0) * nodes->GetAuxVarGradient(iPoint, 0, 0);
+    su2double Temp2 = flowNodes->GetVelocity(iPoint, 1) * nodes->GetAuxVarGradient(iPoint, 0, 1);
     su2double Temp3 = nodes->GetAuxVarGradient(iPoint, 0, 0);
     su2double Temp4 = nodes->GetAuxVarGradient(iPoint, 0, 1);
     su2double Temp5 = nodes->GetAuxVarGradient(iPoint, 1, 0);