regenerate test data, gpu update

Author: conraddelgado

examples/2D_ibm_circle_periodic/case.py

@@ -17,13 +17,13 @@
 M = 1.2
 
 # reynolds number
-Re = 1500.0
+Re = 400.0
 
 # pressure
 P = 101325
 
 # density
-rho = 1.225
+rho = 1.0
 
 # fluid velocity
 v1 = M * np.sqrt(gam_a * P / rho)
@@ -32,8 +32,8 @@
 mu = rho * v1 * D / Re
 
 dt = 1.0e-06
-Nt = 1000  # int(10 * L / v1 / dt)
-t_save = 10  # Nt // 300
+Nt = 5000 
+t_save = 20 
 
 Nx = 511  # to fully resolve requires ~ 40-60 cells across sphere diameter
 Ny = Nx

examples/3D_ibm_sphere_periodic/case.py

@@ -8,7 +8,7 @@
 D = 0.1
 
 # sphere density
-rho_s = 200  # water 1000 kg/m^3
+rho_s = 200 
 
 # domain length
 L = 10 * D
@@ -17,13 +17,13 @@
 M = 1.2
 
 # reynolds number
-Re = 1500.0
+Re = 400.0
 
 # pressure
 P = 101325
 
 # density
-rho = 1.225
+rho = 1.0
 
 # fluid velocity
 v1 = M * np.sqrt(gam_a * P / rho)
@@ -32,8 +32,8 @@
 mu = rho * v1 * D / Re
 
 dt = 2.0e-06
-Nt = int(4 * L / v1 / dt)
-t_save = Nt // 300
+Nt = 2000
+t_save = Nt // 250
 
 Nx = 127  # to fully resolve requires ~ 40-60 cells across sphere diameter
 Ny = Nx
@@ -46,7 +46,7 @@
         f"patch_ib({1})%geometry": 8,
         f"patch_ib({1})%x_centroid": 0.5 * L,
         f"patch_ib({1})%y_centroid": 0.5 * L,
-        f"patch_ib({1})%z_centroid": 0.0 * L,
+        f"patch_ib({1})%z_centroid": 0.5 * L,
         f"patch_ib({1})%radius": D / 2,
         f"patch_ib({1})%slip": "F",
         f"patch_ib({1})%moving_ibm": 0,

src/common/m_mpi_common.fpp

@@ -1919,6 +1919,7 @@ contains
             end if
         end if
 #endif
+        $:GPU_UPDATE(device='[domain_glb]')
     end subroutine s_mpi_global_domain_bounds
 #endif