Merge pull request #2425 from vkrm1612/issue_2342

pcarruscag · web-flow · commit 23286b078500 · 2025-01-24T19:41:27.000-08:00
CTurbSASolver.cpp changed as mentioned
diff --git a/SU2_CFD/src/solvers/CTurbSASolver.cpp b/SU2_CFD/src/solvers/CTurbSASolver.cpp
@@ -500,7 +500,7 @@ void CTurbSASolver::BC_HeatFlux_Wall(CGeometry *geometry, CSolver **solver_conta
          Res_Wall = coeff*RoughWallBC*Area;
          LinSysRes.SubtractBlock(iPoint, &Res_Wall);
 
-         su2double Jacobian_i = (laminar_viscosity*Area)/(0.03*Roughness_Height*sigma);
+         su2double Jacobian_i = (laminar_viscosity /density *Area)/(0.03*Roughness_Height*sigma);
          Jacobian_i += 2.0*RoughWallBC*Area/sigma;
          if (implicit) Jacobian.AddVal2Diag(iPoint, -Jacobian_i);
       }

Original file line number	Diff line number	Diff line change
`@@ -500,7 +500,7 @@ void CTurbSASolver::BC_HeatFlux_Wall(CGeometry geometry, CSolver *solver_conta`
`500`	`500`	`Res_Wall = coeffRoughWallBCArea;`
`501`	`501`	`LinSysRes.SubtractBlock(iPoint, &Res_Wall);`
`502`	`502`
`503`		`- su2double Jacobian_i = (laminar_viscosityArea)/(0.03Roughness_Height*sigma);`
	`503`	`+ su2double Jacobian_i = (laminar_viscosity /density Area)/(0.03Roughness_Height*sigma);`
`504`	`504`	`Jacobian_i += 2.0RoughWallBCArea/sigma;`
`505`	`505`	`if (implicit) Jacobian.AddVal2Diag(iPoint, -Jacobian_i);`
`506`	`506`	`}`