Add moving imersed boundaries

.github/workflows/frontier/submit.sh

@@ -33,7 +33,7 @@ sbatch <<EOT
 #SBATCH -A CFD154                  # charge account
 #SBATCH -N 1                       # Number of nodes required
 $sbatch_device_opts
-#SBATCH -t 04:59:00                # Duration of the job (Ex: 15 mins)
+#SBATCH -t 05:59:00                # Duration of the job (Ex: 15 mins)
 #SBATCH -o$job_slug.out            # Combined output and error messages file
 #SBATCH -p extended                # Extended partition for shorter queues
 #SBATCH -W                         # Do not exit until the submitted job terminates.

.typos.toml

@@ -19,6 +19,7 @@ strang = "strang"
 Strang = "Strang"
 TKE = "TKE"
 HSA = "HSA"
+infp = "infp"
 
 [files]
 extend-exclude = ["docs/documentation/references*", "tests/", "toolchain/cce_simulation_workgroup_256.sh"]

.vscode/settings.json

@@ -49,7 +49,7 @@
     // Enable ONLY fortls language server
     "fortran.enableLanguageServer": true,
     "fortran.languageServer": "fortls",
-    "fortran.fortls.disabled": false,
+    "fortran.fortls.disabled": true,
     "fortran.fortls.path": "fortls",
 
     // Try to disable any built-in language features

CMakeLists.txt

@@ -328,6 +328,13 @@ macro(HANDLE_SOURCES target useCommon)
                                                "${CMAKE_BINARY_DIR}/modules/${target}/*.fpp")
     if (${useCommon})
         file(GLOB common_FPPs CONFIGURE_DEPENDS "${common_DIR}/*.fpp")
+
+        # If we're building post_process, exclude m_compute_levelset.fpp
+        if("${target}" STREQUAL "post_process")
+            list(FILTER common_FPPs EXCLUDE REGEX ".*/m_compute_levelset\.fpp$")
+            list(FILTER common_FPPs EXCLUDE REGEX ".*/m_ib_patches\.fpp$")
+        endif()
+
         list(APPEND ${target}_FPPs ${common_FPPs})
     endif()
 
@@ -443,7 +450,7 @@ function(MFC_SETUP_TARGET)
 
         if (ARGS_SILO)
             find_package(SILO REQUIRED)
-            target_link_libraries(${a_target} PRIVATE SILO::SILO)
+            target_link_libraries(${a_target} PRIVATE SILO::SILO stdc++)
         endif()
 
         if (ARGS_HDF5)

docs/documentation/gpuParallelization.md

@@ -31,7 +31,7 @@ Note: Ordering is not guaranteed or stable, so use key-value pairing when using
 
 - Data on the GPU has a reference counter
 - When data is referred to being allocated, it means that GPU memory is allocated if it is not already present in GPU memory. If a variable is already present, the reference counter is just incremented.
-- When data is referred to being dellocated, it means that the reference counter is decremented. If the reference counter is zero, then the data is actually deallocated from GPU memory
+- When data is referred to being deallocated, it means that the reference counter is decremented. If the reference counter is zero, then the data is actually deallocated from GPU memory
 - When data is referred to being attached, it means that the device pointer attaches to target if it not already attached. If pointer is already attached, then the attachment counter is just incremented
 - When data is referred to being detached, it means that the attachment counter is decremented. If attachment counter is zero, then actually detached
 

examples/1D_shuosher_analytical/case.py

@@ -57,7 +57,7 @@
             "patch_icpp(1)%pres": 10.3333,
             "patch_icpp(1)%alpha_rho(1)": 3.957143,
             "patch_icpp(1)%alpha(1)": 1.0,
-            # One anlytic patch to take care of -4 < x < 5
+            # One analytic patch to take care of -4 < x < 5
             # Patch 2 Analytic
             "patch_icpp(2)%geometry": 1,
             "patch_icpp(2)%x_centroid": 0.5,

examples/1D_shuosher_old/case.py

@@ -57,7 +57,7 @@
             "patch_icpp(1)%pres": 10.3333,
             "patch_icpp(1)%alpha_rho(1)": 3.957143,
             "patch_icpp(1)%alpha(1)": 1.0,
-            # One anlytic patch to take care of -4 < x < 5
+            # One analytic patch to take care of -4 < x < 5
             # Patch 2 Analytic
             "patch_icpp(2)%geometry": 1,
             "patch_icpp(2)%x_centroid": 0.5,

examples/1D_shuosher_teno5/case.py

@@ -59,7 +59,7 @@
             "patch_icpp(1)%pres": 10.333,
             "patch_icpp(1)%alpha_rho(1)": 3.857,
             "patch_icpp(1)%alpha(1)": 1.0,
-            # One anlytic patch to take care of 1 < x < 10
+            # One analytic patch to take care of 1 < x < 10
             # Patch 2 Analytic
             "patch_icpp(2)%geometry": 1,
             "patch_icpp(2)%x_centroid": 5.5,

examples/1D_shuosher_teno7/case.py

@@ -59,7 +59,7 @@
             "patch_icpp(1)%pres": 10.333,
             "patch_icpp(1)%alpha_rho(1)": 3.857,
             "patch_icpp(1)%alpha(1)": 1.0,
-            # One anlytic patch to take care of 1 < x < 10
+            # One analytic patch to take care of 1 < x < 10
             # Patch 2 Analytic
             "patch_icpp(2)%geometry": 1,
             "patch_icpp(2)%x_centroid": 5.5,

examples/1D_shuosher_wenojs5/case.py

@@ -58,7 +58,7 @@
             "patch_icpp(1)%pres": 10.333,
             "patch_icpp(1)%alpha_rho(1)": 3.857,
             "patch_icpp(1)%alpha(1)": 1.0,
-            # One anlytic patch to take care of 1 < x < 10
+            # One analytic patch to take care of 1 < x < 10
             # Patch 2 Analytic
             "patch_icpp(2)%geometry": 1,
             "patch_icpp(2)%x_centroid": 5.5,

examples/1D_shuosher_wenom5/case.py

@@ -58,7 +58,7 @@
             "patch_icpp(1)%pres": 10.333,
             "patch_icpp(1)%alpha_rho(1)": 3.857,
             "patch_icpp(1)%alpha(1)": 1.0,
-            # One anlytic patch to take care of 1 < x < 10
+            # One analytic patch to take care of 1 < x < 10
             # Patch 2 Analytic
             "patch_icpp(2)%geometry": 1,
             "patch_icpp(2)%x_centroid": 5.5,

examples/1D_shuosher_wenoz5/case.py

@@ -58,7 +58,7 @@
             "patch_icpp(1)%pres": 10.333,
             "patch_icpp(1)%alpha_rho(1)": 3.857,
             "patch_icpp(1)%alpha(1)": 1.0,
-            # One anlytic patch to take care of 1 < x < 10
+            # One analytic patch to take care of 1 < x < 10
             # Patch 2 Analytic
             "patch_icpp(2)%geometry": 1,
             "patch_icpp(2)%x_centroid": 5.5,

src/common/include/case.fpp

@@ -4,4 +4,5 @@
 
 ! For pre-process.
 #:def analytical()
+
 #:enddef

src/common/m_boundary_common.fpp

@@ -3,7 +3,7 @@
 !! @brief Contains module m_boundary_conditions_common
 
 !> @brief The purpose of the module is to apply noncharacteristic and processor
-!! boundary condiitons
+!! boundary conditions
 
 #:include 'macros.fpp'
 

src/pre_process/m_compute_levelset.fpp → src/common/m_compute_levelset.fpp

@@ -300,6 +300,7 @@ contains
                         end if
                     end do
 
+                    ! TODO :: This entire if/else tree can be simplified by just calling arrays by idx index
                     if (idx == 1) then
                         levelset%sf(i, j, 0, ib_patch_id) = side_dists(1)
                         if (f_approx_equal(side_dists(1), 0._wp)) then

src/common/m_derived_types.fpp

@@ -318,6 +318,11 @@ module m_derived_types
 
         real(wp) :: model_threshold !<
         !! Threshold to turn on smoothen STL patch.
+
+        !! Patch conditions for moving imersed boundaries
+        integer :: moving_ibm ! 0 for no moving, 1 for moving, 2 for moving on forced path
+        real(wp), dimension(1:3) :: vel
+
     end type ib_patch_parameters
 
     !> Derived type annexing the physical parameters (PP) of the fluids. These