Implement the initialization of VOF using a parameter incflo.vof_init…

… defined in input file

test_2d/inputs.dropflight

@@ -0,0 +1,105 @@
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#            SIMULATION STOP            #
+#.......................................#
+#stop_time               = 2.2426211256
+stop_time               =   100e-6         # Max (simulated) time to evolve
+#stop_time =1.66481717925811447992
+#max_step                = 10          # Max number of time steps
+steady_state            =   0           # Steady-state solver? 
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#         TIME STEP COMPUTATION         #
+#.......................................#
+#incflo.fixed_dt        =   1e-7         # Use this constant dt if > 0
+incflo.cfl              =   0.1        # CFL factor
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#            INPUT AND OUTPUT           #
+#.......................................#
+amr.plot_int            =  400     # Steps between plot files
+amr.check_int           =   1000        # Steps between checkpoint files
+amr.restart             =   ""          # Checkpoint to restart from 
+amr.plt_p   = 1
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#               PHYSICS                 #
+#.......................................#
+incflo.gravity          =   0.  0. 0.   # Gravitational force (3D)
+incflo.ro_0             =   1.0e-3     # Reference density 
+incflo.ntrac            =   1
+incflo.fluid_model      =   "newtonian" # Fluid model (rheology)
+incflo.mu               =   1.8e-4         # Dynamic viscosity coefficient
+incflo.vof_advect_tracer=   1
+#incflo.mu_s             = 0.0182574185835055
+#incflo.ro_s             =  1.0
+incflo.mu_s             =   1e-2 
+incflo.ro_s             =   1.
+#incflo.mu_s             =   0.0182574185835055
+#incflo.ro_s             =   .1
+
+incflo.sigma            =   72.
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#        ADAPTIVE MESH REFINEMENT       #
+#.......................................#
+amr.n_cell              =  32 64   # Grid cells at coarsest AMRlevel
+amr.max_level           =   0           # Max AMR level in hierarchy 
+amr.max_grid_size       =  64 128  
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#              GEOMETRY                 #
+#.......................................#
+geometry.prob_lo        =   -60.e-4 0.           # Lo corner coordinates
+geometry.prob_hi        =    60.e-4 240.e-4      # Hi corner coordinates
+geometry.is_periodic    =   0 0       # Periodicity x y z (0/1)
+geometry.coord_sys   =  0
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#           INITIAL CONDITIONS          #
+#.......................................#
+incflo.probtype         =   1109        #
+incflo.ic_u             =   0.
+incflo.ic_v             =   100.
+incflo.ic_w             =   0.
+
+#initialize the VOF field with AMREX EB functions
+incflo.vof_init_with_eb = 1
+#it is a spherical droplet
+incflo.vof_init          ="x*x+(y-35.e-4)*(y-35.e-4)-20.e-4*20.e-4"
+
+# Boundary conditions
+ylo.type                =   "mi"
+#ylo.tracer              = 0
+#ylo.pressure            =   0
+ylo.velocity           = 0. .1 0.  
+yhi.type                =   "po"
+#yhi.velocity           =0. .1 0.
+yhi.pressure            =   0.  
+#yhi.tracer              =  0.
+xlo.type                =   "nsw"
+#xlo.pressure            =   0.
+xhi.type                =   "nsw"
+#xhi.pressure            =   0.0
+
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#          NUMERICAL PARAMETERS         #
+#.......................................#
+incflo.steady_state_tol = 1.e-5         # Tolerance for steady-state
+amrex.fpe_trap_invalid  =   1           # Trap NaNs
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#              VERBOSITY                #
+#.......................................#
+incflo.verbose          =   2           # incflo_level
+amrex.fpe_trap_overflow=1
+incflo.use_cc_proj      =   1
+incflo.advect_momentum =  0
+incflo.number_of_averaging = 0
+incflo.diffusion_type = 2
+incflo.output_drop_frequence = 10
+incflo.godunov_include_diff_in_forcing = 1
+incflo.use_mac_phi_in_godunov = 0
+incflo.godunov_use_forces_in_trans = 0
+incflo.initial_iterations =0
+incflo.do_initial_proj = 0
+amrex.fpe_trap_zero=1

test_2d/inputs.droplet

@@ -0,0 +1,95 @@
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#            SIMULATION STOP            #
+#.......................................#
+stop_time               =   8         # Max (simulated) time to evolve
+max_step                =   10000          # Max number of time steps
+steady_state            =   0           # Steady-state solver? 
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#         TIME STEP COMPUTATION         #
+#.......................................#
+#incflo.fixed_dt        =   .1         # Use this constant dt if > 0
+incflo.cfl              =   0.1        # CFL factor
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#            INPUT AND OUTPUT           #
+#.......................................#
+amr.plot_int            =  40           # Steps between plot files
+amr.check_int           =   1000        # Steps between checkpoint files
+amr.restart             =   ""          # Checkpoint to restart from 
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#               PHYSICS                 #
+#.......................................#
+incflo.gravity          =   0.  0. 0.   # Gravitational force (3D)
+incflo.ro_0             =   1.0         # Reference density 
+incflo.ntrac            =   1
+incflo.fluid_model      =   "newtonian" # Fluid model (rheology)
+incflo.mu               =   1.e-3         # Dynamic viscosity coefficient
+incflo.vof_advect_tracer=   1
+incflo.mu_s             =   1.0 
+incflo.ro_s             =   1000.0
+incflo.sigma            =   10.
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#        ADAPTIVE MESH REFINEMENT       #
+#.......................................#
+amr.n_cell              =   32 32      # Grid cells at coarsest AMRlevel
+amr.max_level           =   0           # Max AMR level in hierarchy 
+amr.max_grid_size       =   32 32 
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#              GEOMETRY                 #
+#.......................................#
+geometry.prob_lo        =   0.  0.    # Lo corner coordinates
+geometry.prob_hi        =   1.  1.    # Hi corner coordinates
+geometry.is_periodic    =   1 1       # Periodicity x y z (0/1)
+geometry.coord_sys   =  0
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#           INITIAL CONDITIONS          #
+#.......................................#
+incflo.probtype         =   1109        #
+incflo.ic_u             =   0.
+incflo.ic_v             =   0.
+incflo.ic_w             =   0.
+
+#initialize the VOF field with AMREX EB functions
+incflo.vof_init_with_eb = 1
+#the interface is a cosin function
+incflo.vof_init    ="if (x>0.3 and x<0.7 and y>0.3 and y<0.7,-1.,1.)"
+
+
+# Boundary conditions
+#ylo.type                =   "po"
+#ylo.pressure            =   0.  
+#yhi.type                =   "po"
+#yhi.pressure            =   0.  
+#xlo.type                =   "po"
+#xlo.pressure            =   0.
+#xhi.type                =   "po"
+#xhi.pressure            =   0.0
+
+
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#          NUMERICAL PARAMETERS         #
+#.......................................#
+incflo.steady_state_tol = 1.e-5         # Tolerance for steady-state
+amrex.fpe_trap_invalid  =   1           # Trap NaNs
+
+#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
+#              VERBOSITY                #
+#.......................................#
+incflo.verbose          =   2           # incflo_level
+amrex.fpe_trap_overflow=1
+incflo.use_cc_proj      =   1
+incflo.advect_momentum =  0
+incflo.number_of_averaging = 1
+incflo.diffusion_type = 2
+incflo.godunov_include_diff_in_forcing = 1
+incflo.use_mac_phi_in_godunov = 0
+incflo.godunov_use_forces_in_trans = 0
+incflo.initial_iterations =0
+incflo.do_initial_proj = 0
+amrex.fpe_trap_zero=1