Add new iNS NACA test

Author: AbbBallout

.github/workflows/CI_sequential 3.yml

@@ -499,3 +499,19 @@ jobs:
           ./horses3d.ins Actuator.control
         if: '!cancelled()'
 
+# 8) INS NACA0012
+# -------------------------------
+
+      - name: Build iNS-NACA0012
+        working-directory: ./Solver/test/IncompressibleNS/NACA0012/SETUP
+        run: |
+          source /opt/intel/oneapi/setvars.sh || true
+          make MODE=${{matrix.mode}} COMPILER=${{matrix.compiler}} COMM=${{matrix.comm}} ENABLE_THREADS=${{matrix.enable_threads}} WITH_MKL=${{matrix.mkl}}
+        if: '!cancelled()'
+
+      - name: Run iNS-NACA0012
+        working-directory: ./Solver/test/IncompressibleNS/NACA0012
+        run: |
+          source /opt/intel/oneapi/setvars.sh || true
+          ./horses3d.ins NACA0012.control
+        if: '!cancelled()'

.gitignore

@@ -9,6 +9,7 @@ bin/HORSES3D
 make.inc
 *.hmesh
 *.bmesh
+*.pmesh
 *.hsol
 RESULTS
 

Solver/Makefile.in

@@ -41,7 +41,7 @@ GNU_HPC_FLAGS= -cpp -fallow-argument-mismatch -std=legacy -ffree-line-length-0 -
 		   -ftree-vectorizer-verbose=0 -fbackslash -D_has_Quad -march=znver3 -funroll-loops
 GNU_RELEASE_FLAGS= -cpp -ffree-line-length-0 -O3 -ftree-vectorize  -ftree-vectorizer-verbose=0 \
 		   -fbackslash -D_has_Quad
-GNU_DEBUG_FLAGS=  -ffree-line-length-0 -cpp -O0 -g -fcheck=all -Wno-unused-variable \
+GNU_DEBUG_FLAGS=  -ffree-line-length-0 -cpp -O0 -frecursive -g -fcheck=all -Wno-unused-variable \
 		  -fbacktrace -fbounds-check -finit-real=snan -Wall \
 		  -ffpe-trap=invalid,zero,overflow -fcheck-array-temporaries \
 		  -fbackslash -fcheck=pointer -D_has_Quad

Solver/configure

@@ -92,6 +92,7 @@ TEST_CASES="./Euler/BoxAroundCircle \
             ./IncompressibleNS/Cylinder \
             ./IncompressibleNS/CylinderVreman \
             ./IncompressibleNS/ActuatorLineInterpolation \
+            ./IncompressibleNS/NACA0012 \
             ./CahnHilliard/TJokisaari \
             ./Multiphase/Convergence \
             ./Multiphase/Convergence_cs \

Solver/test/IncompressibleNS/NACA0012/NACA0012_ins.control

@@ -0,0 +1,66 @@
+/* Control file for the implicit Navier-Stokes solver */
+Flow equations        = "iNS"
+mesh file name        = "../../TestMeshes/naca0012_p3.msh"
+Polynomial order      = 2 ! T
+Number of time steps  = 100
+Output Interval       = 10 
+autosave mode         = "iteration"
+autosave interval     = 1000
+
+Convergence tolerance = 1.d-10
+cfl                   = 0.1d0
+dcfl                  = 0.1d0
+Density (kg/m^3) = 1.0
+Viscosity (Pa.s)  = 1.0
+reference velocity (m/s)  = 100.0
+Artificial compressibility factor = 1000.0
+
+discretization nodes  = Gauss-Lobatto
+
+solution file name    = "RESULTS/naca.hsol"
+save gradients with solution = .true.
+restart               = .false.
+restart file name     = "RESULTS/restart.hsol"
+riemann solver         = "central"
+viscous discretization = "BR1"
+
+Riemann solver = central  
+Lambda stabilization = 1.0
+!Inviscid discretization = Split-form
+Averaging = Skew-symmetric 2 ! Pirozzoli/Morinishi/Kennedy-Gruber/Entropy conserving/Entropy and
+LES Model = Vreman
+                               
+simulation type       = time-accurate
+final time            = 40.d0
+!simulation type       = steady-state
+time integration        = explicit
+partitioning = SFC
+!-------------------------------------------------------------
+! Inflow
+!-------------------------------------------------------------
+#define boundary BACK
+   type = FreeSlipWall
+#end
+
+#define boundary FRONT
+   type = FreeSlipWall
+#end
+
+#define boundary AIRFOILT
+  type = NoSlipWall
+#end
+
+#define boundary AIRFOILB
+  type = NoSlipWall
+#end
+
+#define boundary INFLOW
+    type = Inflow
+#end
+
+#define boundary OUTFLOW
+    type = Outflow
+#end
+!-------------------------------------------------------------
+
+