Copy the ib flag to the device (#1937)

Author: sbryngelson

src/simulation/m_global_parameters.fpp

@@ -991,7 +991,7 @@ contains
 
         $:GPU_UPDATE(device='[alt_soundspeed, acoustic_source, num_source]')
         $:GPU_UPDATE(device='[dt, sys_size, buff_size, eqn_idx, mpp_lim, bubbles_euler, hypoelasticity, alt_soundspeed, &
-                     & avg_state, model_eqns, mixture_err, grid_geometry, cyl_coord, mp_weno, weno_eps, teno_CT, low_Mach]')
+                     & avg_state, model_eqns, mixture_err, grid_geometry, cyl_coord, mp_weno, weno_eps, teno_CT, low_Mach, ib]')
         $:GPU_UPDATE(device='[riemann_hypo_ADC, ADC_kappa, hll_u_interface, hypo_hll_interface_rhs, hypo_nc_mode]')
 
         $:GPU_UPDATE(device='[Bx0]')

tests/09FDDC81/golden-metadata.txt

@@ -844,7 +844,7 @@ def alter_num_fluids(dimInfo):
             alter_low_Mach_correction()
             if num_fluids == 1:
                 alter_eos()
-            alter_ib(dimInfo)
+            alter_ib(dimInfo, num_fluids=num_fluids)
             if len(dimInfo[0]) > 1:
                 alter_igr()
 
@@ -1194,7 +1194,7 @@ def alter_ppn(dimInfo):
             if ARG("rdma_mpi"):
                 cases.append(define_case_d(stack, "2 MPI Ranks -> RDMA MPI", {"rdma_mpi": "T"}, ppn=2))
 
-    def alter_ib(dimInfo, six_eqn_model=False, viscous=False):
+    def alter_ib(dimInfo, six_eqn_model=False, viscous=False, num_fluids=1):
         for slip in [True, False]:
             stack.push(
                 "IBM",
@@ -1328,6 +1328,37 @@ def alter_ib(dimInfo, six_eqn_model=False, viscous=False):
                 )
             )
 
+        if len(dimInfo[0]) == 2 and not viscous and num_fluids == 1:
+            # IB-interior cells must not set the adaptive dt. A body moving through gas at rest carries its velocity into
+            # its interior, so there |u| + c exceeds anywhere in the fluid and dt shrinks ~8% if those cells are counted.
+            # The edges sit on cell faces and the body moves 0.2 dx, so no cell crosses the surface.
+            cases.append(
+                define_case_d(
+                    stack,
+                    "IBM -> Adaptive dt",
+                    {
+                        "ib": "T",
+                        "num_ibs": 1,
+                        "fd_order": 2,
+                        "n": 49,
+                        "cfl_adap_dt": "T",
+                        "cfl_target": 0.5,
+                        "n_start": 0,
+                        "t_save": 0.04,
+                        "t_stop": 0.04,
+                        **{f"patch_icpp({i})%{k}": v for i in range(1, 4) for k, v in [("pres", 1.0), ("alpha_rho(1)", 1.0), ("vel(1)", 0.0)]},
+                        "patch_ib(1)%geometry": 3,
+                        "patch_ib(1)%x_centroid": 0.5,
+                        "patch_ib(1)%y_centroid": 0.5,
+                        "patch_ib(1)%length_x": 0.2,
+                        "patch_ib(1)%length_y": 0.2,
+                        "patch_ib(1)%slip": "F",
+                        "patch_ib(1)%moving_ibm": 2,
+                        "patch_ib(1)%vel(1)": 0.1,
+                    },
+                )
+            )
+
         if len(dimInfo[0]) == 2 and not viscous:
             cases.append(
                 define_case_d(