Update SU2_CFD/src/solvers/CTurbSSTSolver.cpp

Shubham shukla · pcarruscag · web-flow · commit 3e8292c2de06 · 2025-12-27T23:54:07.000+05:30
Co-authored-by: Pedro Gomes &lt;38071223+pcarruscag@users.noreply.github.com&gt;
diff --git a/SU2_CFD/src/solvers/CTurbSSTSolver.cpp b/SU2_CFD/src/solvers/CTurbSSTSolver.cpp
@@ -272,7 +272,7 @@ void CTurbSSTSolver::Postprocessing(CGeometry *geometry, CSolver **solver_contai
         }
         shearStress = sqrt(shearStress);
 
-        const su2double FrictionVelocity = max(sqrt(shearStress/flowNodes->GetDensity(iPoint)), 1e-10);
+        const su2double FrictionVelocity = max(sqrt(shearStress/flowNodes->GetDensity(iPoint)), EPS);
         const su2double wall_dist = geometry->vertex[iMarker][iVertex]->GetNearestNeighborDistance();
 
         const su2double Derivative = flowNodes->GetLaminarViscosity(jPoint) * pow(nodes->GetSolution(jPoint, 0), 0.673) / max(wall_dist, EPS);

Original file line number	Diff line number	Diff line change
`@@ -272,7 +272,7 @@ void CTurbSSTSolver::Postprocessing(CGeometry geometry, CSolver *solver_contai`
`272`	`272`	`}`
`273`	`273`	`shearStress = sqrt(shearStress);`
`274`	`274`
`275`		`- const su2double FrictionVelocity = max(sqrt(shearStress/flowNodes->GetDensity(iPoint)), 1e-10);`
	`275`	`+ const su2double FrictionVelocity = max(sqrt(shearStress/flowNodes->GetDensity(iPoint)), EPS);`
`276`	`276`	`const su2double wall_dist = geometry->vertex[iMarker][iVertex]->GetNearestNeighborDistance();`
`277`	`277`
`278`	`278`	`const su2double Derivative = flowNodes->GetLaminarViscosity(jPoint) * pow(nodes->GetSolution(jPoint, 0), 0.673) / max(wall_dist, EPS);`