From 75a4c63af920686f7b0f12341ef3c808d9cf4076 Mon Sep 17 00:00:00 2001 From: Hua Tan Date: Sat, 7 Dec 2024 23:20:08 -0800 Subject: [PATCH] Implement the initialization of VOF using a parameter incflo.vof_init defined in input file --- test_2d/inputs.dropflight | 105 ++++++++++++++++++++++++++++++++++++++ test_2d/inputs.droplet | 95 ++++++++++++++++++++++++++++++++++ 2 files changed, 200 insertions(+) create mode 100644 test_2d/inputs.dropflight create mode 100644 test_2d/inputs.droplet diff --git a/test_2d/inputs.dropflight b/test_2d/inputs.dropflight new file mode 100644 index 00000000..9186f753 --- /dev/null +++ b/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 diff --git a/test_2d/inputs.droplet b/test_2d/inputs.droplet new file mode 100644 index 00000000..9dcaaa40 --- /dev/null +++ b/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