added example case

Author: wilfonba

examples/2D_hardcodied_ic/case.py

@@ -0,0 +1,88 @@
+#!/usr/bin/env python3
+
+import json
+
+eps = 1e-9
+
+# Configuring case dictionary
+print(json.dumps({
+    # Logistics ================================================
+    'run_time_info'                : 'T',
+    # ==========================================================
+
+    # Computational Domain Parameters ==========================
+    'x_domain%beg'                 : 0.E+00,
+    'x_domain%end'                 : 4.E+00,
+    'stretch_x'                    : 'T',
+    'a_x'                          : 7,
+    'x_a'                          : -2,
+    'x_b'                          : 2,
+    'y_domain%beg'                 : 0.E+00,
+    'y_domain%end'                 : 4.E+00,
+    'stretch_y'                    : 'T',
+    'a_y'                          : 7,
+    'y_a'                          : -2,
+    'y_b'                          : 2,
+    'm'                            : 199,
+    'n'                            : 199,
+    'p'                            : 0,
+    'dt'                           : 5.E-07,
+    't_step_start'                 : 0,
+    't_step_stop'                  : 100000,
+    't_step_save'                  : 1000,
+    # ==========================================================
+
+    # Simulation Algorithm Parameters ==========================
+    'num_patches'                  : 1,
+    'model_eqns'                   : 3,
+    'alt_soundspeed'               : 'F',
+    'num_fluids'                   : 2,
+    'adv_alphan'                   : 'T',
+    'mpp_lim'                      : 'T',
+    'mixture_err'                  : 'T',
+    'time_stepper'                 : 3,
+    'weno_order'                   : 5,
+    'weno_eps'                     : 1.E-16,
+    'mapped_weno'                  : 'T',
+    'null_weights'                 : 'F',
+    'mp_weno'                      : 'F',
+    'riemann_solver'               : 2,
+    'wave_speeds'                  : 1,
+    'avg_state'                    : 2,
+    'bc_x%beg'                     : -2,
+    'bc_x%end'                     : -2,
+    'bc_y%beg'                     : -6,
+    'bc_y%end'                     : -6,
+    # ==========================================================
+
+    # Formatted Database Files Structure Parameters ============
+    'format'                       : 1,
+    'precision'                    : 2,
+    'prim_vars_wrt'                :'T',
+    'parallel_io'                  :'T',
+    # ==========================================================
+                                                            
+    # Patch 1: Base ============================================
+    'patch_icpp(1)%geometry'       : 7,
+    'patch_icpp(1)%hcid'           : 200,
+    'patch_icpp(1)%x_centroid'     : 4.,
+    'patch_icpp(1)%y_centroid'     : 4.,
+    'patch_icpp(1)%length_x'       : 8.,
+    'patch_icpp(1)%length_y'       : 8.,
+    'patch_icpp(1)%vel(1)'         : 0.,
+    'patch_icpp(1)%vel(2)'         : 0.,
+    'patch_icpp(1)%pres'           : 100,
+    'patch_icpp(1)%alpha_rho(1)'   : (1-eps)*1000,
+    'patch_icpp(1)%alpha_rho(2)'   : eps*1,
+    'patch_icpp(1)%alpha(1)'       : 1-eps,
+    'patch_icpp(1)%alpha(2)'       : eps,
+    # ===========================================================
+
+    # Fluids Physical Parameters ===============================
+    'fluid_pp(1)%gamma'            : 1.E+00/(2.35E+00-1.E+00),
+    'fluid_pp(1)%pi_inf'           : 2.35E+00*1.E+09/(2.35E+00-1.E+00),
+    'fluid_pp(2)%gamma'            : 1.E+00/(1.4E+00-1.E+00),
+    'fluid_pp(2)%pi_inf'           : 0.E+00,
+    # ==========================================================
+}))
+# ==============================================================================

src/pre_process/include/1dHardcodedIC.fpp

@@ -1,3 +1,9 @@
+#:def Hardcoded1DVariables()
+    ! Place any declaration of intermediate variables here
+    
+
+#:enddef
+
 #:def Hardcoded1D()
 
     if (patch_icpp(patch_id)%hcid == 100) then

src/pre_process/include/2dHardcodedIC.fpp

@@ -1,6 +1,25 @@
+#:def Hardcoded2DVariables()
+
+    real(kind(0d0)) :: eps
+
+    eps = 1e-9
+
+#:enddef
+
 #:def Hardcoded2D()
 
-    if (patch_icpp(patch_id)%hcid == 200) then
+    if (patch_icpp(patch_id)%hcid == 200) then ! 2D_hardcoded_ic example case
+
+        if (y_cc(j) <= (-x_cc(i)**3 + 1)**(1d0/3d0)) then
+            ! Volume Fractions
+            q_prim_vf(advxb)%sf(i, j, 0) = eps
+            q_prim_vf(advxe)%sf(i, j, 0) = 1d0-eps
+            ! Denssities
+            q_prim_vf(contxb)%sf(i, j, 0) = eps*1000d0
+            q_prim_vf(contxe)%sf(i, j, 0) = (1d0-eps)*1d0
+            ! Pressure
+            q_prim_vf(E_idx)%sf(i, j, 0) = 1d0
+        end if
 
     end if
 

src/pre_process/include/3dHardcodedIC.fpp

@@ -1,3 +1,9 @@
+#:def Hardcoded3DVariables()
+    ! Place any declaration of intermediate variables here
+    
+
+#:enddef
+
 #:def Hardcoded3D()
 
     if (patch_icpp(patch_id)%hcid == 300) then

src/pre_process/m_patches.fpp

@@ -785,6 +785,8 @@ contains
         ! Generic loop iterators
         integer :: i, j, k
 
+        @:Hardcoded1DVariables()
+
         pi_inf = fluid_pp(1)%pi_inf
         gamma = fluid_pp(1)%gamma
         lit_gamma = (1d0 + gamma)/gamma
@@ -890,6 +892,8 @@ contains
 
         integer :: i, j, k !< generic loop iterators
 
+        @:Hardcoded2DVariables()
+
         pi_inf = fluid_pp(1)%pi_inf
         gamma = fluid_pp(1)%gamma
         lit_gamma = (1d0 + gamma)/gamma
@@ -949,6 +953,8 @@ contains
 
         integer :: i, j, k !< generic loop iterators
 
+        @:Hardcoded3DVariables()
+
         pi_inf = fluid_pp(1)%pi_inf
         gamma = fluid_pp(1)%gamma
         lit_gamma = (1d0 + gamma)/gamma