Skip to content

Commit b66271c

Browse files
committed
fix SA walldist
1 parent 2cb1465 commit b66271c

File tree

1 file changed

+3
-3
lines changed

1 file changed

+3
-3
lines changed

SU2_CFD/src/solvers/CTurbSASolver.cpp

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -275,7 +275,7 @@ void CTurbSASolver::Postprocessing(CGeometry *geometry, CSolver **solver_contain
275275
FrictionVelocity = sqrt(flowNodes->GetLaminarViscosity(iPoint)*VorticityMag);
276276
}
277277

278-
const su2double wall_dist = geometry->nodes->GetWall_Distance(jPoint);
278+
const su2double wall_dist = geometry->vertex[iMarker][iVertex]->GetNearestNeighborDistance();
279279
const su2double Derivative = nodes->GetSolution(jPoint, 0) / wall_dist;
280280
const su2double turbulence_index = Derivative / (FrictionVelocity * 0.41);
281281

@@ -309,7 +309,7 @@ void CTurbSASolver::Source_Residual(CGeometry *geometry, CSolver **solver_contai
309309
CNumerics **numerics_container, CConfig *config, unsigned short iMesh) {
310310

311311
bool axisymmetric = config->GetAxisymmetric();
312-
312+
313313
const bool implicit = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
314314
const bool harmonic_balance = (config->GetTime_Marching() == TIME_MARCHING::HARMONIC_BALANCE);
315315
const bool transition_BC = config->GetSAParsedOptions().bc;
@@ -389,7 +389,7 @@ void CTurbSASolver::Source_Residual(CGeometry *geometry, CSolver **solver_contai
389389
/*--- Set y coordinate ---*/
390390
numerics->SetCoord(geometry->nodes->GetCoord(iPoint), geometry->nodes->GetCoord(iPoint));
391391
}
392-
392+
393393
/*--- Compute the source term ---*/
394394

395395
auto residual = numerics->ComputeResidual(config);

0 commit comments

Comments
 (0)