Add sponge as source term for compressible NS.

loganoz · Feb 28, 2024 · 02b6464 · 02b6464
1 parent 1261df3
commit 02b6464
Show file tree

Hide file tree

Showing 3 changed files with 490 additions and 3 deletions.
diff --git a/Solver/src/NavierStokesSolver/SpatialDiscretization.f90 b/Solver/src/NavierStokesSolver/SpatialDiscretization.f90
@@ -347,6 +347,7 @@ subroutine TimeDerivative_ComputeQDot( mesh , particles, t)
       use WallFunctionConnectivity
          use TripForceClass, only: randomTrip
          use ActuatorLine, only: farm
+         use SpongeClass, only: sponge
          implicit none
          type(HexMesh)              :: mesh
          type(Particles_t)          :: particles
@@ -545,6 +546,8 @@ subroutine TimeDerivative_ComputeQDot( mesh , particles, t)
                end associate
             end do
 !$omp end do
+            ! for the sponge, loops are in the internal subroutine as values are precalculated
+            call sponge % addSource(mesh)
 !
 !           Add Particles source
 !           ********************
@@ -1334,4 +1337,4 @@ SUBROUTINE computeBoundaryFlux(f, time, mesh)
 
       end subroutine computeBoundaryFlux
 
-end module SpatialDiscretization
+end module SpatialDiscretization
diff --git a/Solver/src/libs/sources/Makefile b/Solver/src/libs/sources/Makefile
@@ -10,7 +10,8 @@ MOD_FLAG?=-J
 
 LIBOBJS_NS = ./build_ns/ActuatorLine.o \
 			 ./build_ns/mkl_dfti.o \
-			 ./build_ns/TripForceClass.o 
+			 ./build_ns/TripForceClass.o \
+			 ./build_ns/SpongeClass.o 
 
 LIBOBJS_NSSA =  ./build_nssa/ActuatorLine.o \
 				./build_nssa/mkl_dfti.o \
@@ -157,4 +158,4 @@ clean: FORCE
 #############################################
 #	Dummy procedure to force a rule     #
 #############################################
-FORCE:
+FORCE: