- Fix Normal Distance bug

Author: rois1995

Common/include/geometry/CGeometry.hpp

@@ -1843,6 +1843,13 @@ class CGeometry {
    */
   virtual void SetWallDistance(su2double val) {}
 
+  /*!
+   * \brief Set wall distances a specific value and at a specific point
+   *  \param[in] iPoint - i-th point to set the wall distance.
+   *  \param[in] val - new value for the wall distance at the input point.
+   */
+  virtual void SetWallDistance(unsigned long iPoint, su2double val) {}
+
   /*!
    * \brief Compute the distances to the closest vertex on viscous walls over the entire domain
    * \param[in] config_container - Definition of the particular problem.

Common/include/geometry/CPhysicalGeometry.hpp

@@ -811,6 +811,13 @@ class CPhysicalGeometry final : public CGeometry {
     }
   }
 
+  /*!
+   * \brief Set wall distance of i-th point to a specific value
+   */
+  void SetWallDistance(unsigned long iPoint, su2double val) override {
+    nodes->SetWall_Distance(iPoint, val);
+  }
+
   /*!
    * \brief For streamwise periodicity, find & store a unique reference node on the designated periodic inlet.
    * \param[in] config - Definition of the particular problem.

Common/src/geometry/CGeometry.cpp

@@ -4119,7 +4119,20 @@ void CGeometry::ComputeWallDistance(const CConfig* const* config_container, CGeo
        * This is necessary, because before a computed distance is set, it will be checked
        * whether the new distance is smaller than the currently stored one. ---*/
       CGeometry* geometry = geometry_container[iZone][iInst][MESH_0];
-      if (wallDistanceNeeded[iZone]) geometry->SetWallDistance(numeric_limits<su2double>::max());
+      if (wallDistanceNeeded[iZone]){
+        geometry->SetWallDistance(numeric_limits<su2double>::max());
+
+        const auto nMarker = geometry->GetnMarker();
+        for (auto iMarker = 0; iMarker < config_container[iZone]->GetnMarker_All(); iMarker++) {
+          if (!config_container[iZone]->GetViscous_Wall(iMarker)) continue;
+
+          SU2_OMP_FOR_STAT(OMP_MIN_SIZE)
+          for (auto iVertex = 0u; iVertex < geometry->nVertex[iMarker]; iVertex++) {
+            const auto iPoint = geometry->vertex[iMarker][iVertex]->GetNode();
+            geometry->SetWallDistance(iPoint, 0.0);
+          }
+        }
+      } 
     }
 
     /*--- Loop over all zones and compute the ADT based on the viscous walls in that zone ---*/

Common/src/geometry/CPhysicalGeometry.cpp

@@ -10559,7 +10559,9 @@ void CPhysicalGeometry::SetWallDistance(CADTElemClass* WallADT, const CConfig* c
 
         WallADT->DetermineNearestElement(nodes->GetCoord(iPoint), dist, markerID, elemID, rankID);
 
-        if (dist < nodes->GetWall_Distance(iPoint)) {
+        if (nodes->GetWall_Distance(iPoint) == 0) dist = 0;
+
+        if (dist < nodes->GetWall_Distance(iPoint) || dist == 0) {
           nodes->SetWall_Distance(iPoint, dist, rankID, iZone, markerID, elemID);
         }
       }