New Golden Files

examples/3D_ibm/case.py

@@ -4,7 +4,7 @@
 Mu = 1.84E-05
 gam_a = 1.4
 
-D = 0.6E-03
+D = 0.1
 
 # Configuring case dictionary
 print(json.dumps({
@@ -16,23 +16,23 @@
     # For these computations, the cylinder is placed at the (0,0,0)
     # domain origin. 
     # x direction
-    'x_domain%beg'                 : -5.7*D,
-    'x_domain%end'                 : 15.0*D,
+    'x_domain%beg'                 : -5*D,
+    'x_domain%end'                 : 5.0*D,
     # y direction
-    'y_domain%beg'                 : -7.5*D,
-    'y_domain%end'                 : 7.5*D,
+    'y_domain%beg'                 : -2.5*D,
+    'y_domain%end'                 : 2.5*D,
     # z direction
-    'z_domain%beg'                 : -7.5*D,
-    'z_domain%end'                 : 7.5*D,
+    'z_domain%beg'                 : -2.5*D,
+    'z_domain%end'                 : 2.5*D,
 
     'cyl_coord'                    : 'F',
-    'm'                            : 100,
-    'n'                            : 50,
-    'p'                            : 50,
-    'dt'                           : 0.5E-6,
+    'm'                            : 512,
+    'n'                            : 256,
+    'p'                            : 256,
+    'dt'                           : 1.0E-6,
     't_step_start'                 : 0,
-    't_step_stop'                  : 4,  #3000
-    't_step_save'                  : 2,  #10
+    't_step_stop'                  : 3000,  #3000
+    't_step_save'                  : 300,  #10
     # ==========================================================================
 
     # Simulation Algorithm Parameters ==========================================
@@ -52,7 +52,7 @@
     # Correct errors when computing speed of sound
     'mixture_err'                  : 'T',
     # Use TVD RK3 for time marching
-    'time_stepper'                 : 1,
+    'time_stepper'                 : 2,
     # Reconstruct the primitive variables to minimize spurious
     # Use WENO5
     'weno_order'                   : 5,
@@ -93,19 +93,19 @@
     # Patch: Constant Tube filled with air =====================================
     # Specify the cylindrical air tube grid geometry
     'patch_icpp(1)%geometry'       : 9,
-    'patch_icpp(1)%x_centroid'     : 4.65*D,
+    'patch_icpp(1)%x_centroid'     : 0.0,
     # Uniform medium density, centroid is at the center of the domain
     'patch_icpp(1)%y_centroid'     : 0.0,
     'patch_icpp(1)%z_centroid'     : 0.0,
-    'patch_icpp(1)%length_x'       : 20.7*D,
-    'patch_icpp(1)%length_y'       : 15.0*D,
-    'patch_icpp(1)%length_z'       : 15.0*D,
+    'patch_icpp(1)%length_x'       : 10*D,
+    'patch_icpp(1)%length_y'       : 5*D,
+    'patch_icpp(1)%length_z'       : 5*D,
     # Specify the patch primitive variables 
-    'patch_icpp(1)%vel(1)'         : 1.0E+00,
+    'patch_icpp(1)%vel(1)'         : 527.2E+00,
     'patch_icpp(1)%vel(2)'         : 0.0E+00,
     'patch_icpp(1)%vel(3)'         : 0.0E+00,
-    'patch_icpp(1)%pres'           : 1.E+00,
-    'patch_icpp(1)%alpha_rho(1)'   : 1.E+00,
+    'patch_icpp(1)%pres'           : 10918.2549,
+    'patch_icpp(1)%alpha_rho(1)'   : 0.2199,
     'patch_icpp(1)%alpha(1)'       : 1.E+00,
     # # ========================================================================
 
@@ -141,11 +141,12 @@
     # Patch: Cylinder Immersed Boundary ========================================
     # Specify the cylindrical air tube grid geometry
     'patch_ib(1)%geometry'       : 8,
-    'patch_ib(1)%x_centroid'     : 0.0,
+    'patch_ib(1)%x_centroid'     : -3.0E-3,
     # Uniform medium density, centroid is at the center of the domain
     'patch_ib(1)%y_centroid'     : 0.0,
     'patch_ib(1)%z_centroid'     : 0.0,
     'patch_ib(1)%radius'         : D/2,
+    'patch_ib(1)%slip'           : 'T',
     # 'patch_ib(1)%length_x'         : 0.2E-03,
     # 'patch_ib(1)%length_y'         : 0.4E-03,
     # # ==========================================================================
@@ -154,7 +155,7 @@
     # Use the same stiffness as the air bubble
     'fluid_pp(1)%gamma'            : 1.E+00/(gam_a-1.E+00),  # 2.50(Not 1.40)
     'fluid_pp(1)%pi_inf'           : 0,
-    'fluid_pp(1)%Re(1)'            : 3700*2/D,
+    'fluid_pp(1)%Re(1)'            : 7535533.2,
     # 'fluid_pp(1)%Re(2)'            : 1.0e+6,
     # ==========================================================================
 }))

src/simulation/m_start_up.fpp

@@ -438,16 +438,6 @@ contains
 
         end if
 
-
-
-
-        ! if (any(ib_markers == 1)) then
-        !     print *, "Yay!", proc_rank
-        ! else
-        !     print *, "No :(", proc_rank
-        ! end if
-        ! ==================================================================
-
     end subroutine s_read_serial_data_files ! -------------------------------------
 
         !! @param q_cons_vf Conservative variables