Disable normals vectors for non-analytic velocity fields

ViennaTools · Feb 13, 2025 · 024810b · 024810b
1 parent 1af374b
commit 024810b
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/include/viennaps/psProcess.hpp b/include/viennaps/psProcess.hpp
@@ -291,6 +291,9 @@ template <typename NumericType, int D> class Process {
     advectionKernel.setSaveAdvectionVelocities(advectionParams.velocityOutput);
     advectionKernel.setDissipationAlpha(advectionParams.dissipationAlpha);
     advectionKernel.setIgnoreVoids(advectionParams.ignoreVoids);
+    // normals vectors are only necessary for analytical velocity fields
+    if (model->getVelocityField()->getTranslationFieldOptions() > 0)
+      advectionKernel.setCalculateNormalVectors(false);
 
     for (auto dom : domain->getLevelSets()) {
       meshConverter.insertNextLevelSet(dom);