nearest-neighbor distance for wallfunction

Author: nijso beishuizen

SU2_CFD/src/solvers/CIncNSSolver.cpp

@@ -684,12 +684,9 @@ void CIncNSSolver::SetTau_Wall_WF(CGeometry *geometry, CSolver **solver_containe
 
       const su2double VelTangMod = GeometryToolbox::Norm(int(MAXNDIM), VelTang);
 
-      /*--- Compute normal distance of the interior point from the wall ---*/
+      /*--- Compute nearest-neighbor distance of the interior point from the wall ---*/
 
-      su2double WallDist[MAXNDIM] = {0.0};
-      GeometryToolbox::Distance(nDim, Coord, Coord_Normal, WallDist);
-
-      su2double WallDistMod = GeometryToolbox::Norm(int(MAXNDIM), WallDist);
+      su2double WallDistMod = geometry->vertex[iMarker][iVertex]->GetNearestNeighborDistance();
 
       su2double Density_Wall = nodes->GetDensity(iPoint);
 

SU2_CFD/src/solvers/CNSSolver.cpp

@@ -863,12 +863,9 @@ void CNSSolver::SetTau_Wall_WF(CGeometry *geometry, CSolver **solver_container,
 
       const su2double VelTangMod = GeometryToolbox::Norm(int(MAXNDIM), VelTang);
 
-      /*--- Compute normal distance of the interior point from the wall ---*/
+      /*--- Compute nearest-neighbor distance of the interior point from the wall ---*/
 
-      su2double WallDist[MAXNDIM] = {0.0};
-      GeometryToolbox::Distance(nDim, Coord, Coord_Normal, WallDist);
-
-      const su2double WallDistMod = GeometryToolbox::Norm(int(MAXNDIM), WallDist);
+      const su2double WallDistMod = geometry->vertex[iMarker][iVertex]->GetNearestNeighborDistance();
 
       su2double T_Wall = nodes->GetTemperature(iPoint);
       const su2double Conductivity_Wall = nodes->GetThermalConductivity(iPoint);