diff --git a/Exec/science/flame_wave/inputs.H_He b/Exec/science/flame_wave/inputs.H_He new file mode 100644 index 0000000000..4a9b0f8da3 --- /dev/null +++ b/Exec/science/flame_wave/inputs.H_He @@ -0,0 +1,94 @@ +# ------------------ INPUTS TO MAIN PROGRAM ------------------- +max_step = 9900000 +stop_time = 3.0 + +# PROBLEM SIZE & GEOMETRY +geometry.is_periodic = 0 0 +geometry.coord_sys = 1 # 0 => cart, 1 => RZ 2=>spherical +geometry.prob_lo = 0 0 +geometry.prob_hi = 1.2288e5 3.072e4 +amr.n_cell = 768 192 + +# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< +# 0 = Interior 3 = Symmetry +# 1 = Inflow 4 = SlipWall +# 2 = Outflow 5 = NoSlipWall +# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< +castro.lo_bc = 3 1 +castro.hi_bc = 2 1 + +castro.yl_ext_bc_type = "hse" +castro.hse_interp_temp = 1 +castro.hse_reflect_vels = 1 + +castro.yr_ext_bc_type = "interp" + +# WHICH PHYSICS +castro.do_hydro = 1 +castro.do_react = 1 +castro.do_rotation = 1 +castro.do_grav = 1 +castro.do_sponge = 1 + +castro.small_temp = 1.e6 +castro.small_dens = 1.e-5 + +castro.ppm_type = 1 +castro.grav_source_type = 2 + +gravity.gravity_type = ConstantGrav +gravity.const_grav = -1.5e14 + +castro.rotational_period = 0.0005 +castro.rotation_include_centrifugal = 0 + +castro.diffuse_temp = 1 +castro.diffuse_cutoff_density_hi = 5.e4 +castro.diffuse_cutoff_density = 2.e4 + +castro.diffuse_cond_scale_fac = 1.0 + +castro.react_rho_min = 1.e2 +castro.react_rho_max = 1.5e7 + +castro.react_T_min = 6.e7 + + +# TIME STEP CONTROL +castro.cfl = 0.8 # cfl number for hyperbolic system +castro.init_shrink = 0.1 # scale back initial timestep +castro.change_max = 1.1 # max time step growth +castro.dt_cutoff = 5.e-20 # level 0 timestep below which we halt +castro.use_retry = 1 +castro.hard_cfl_limit = 0 +castro.max_subcycles = 16 + +# DIAGNOSTICS & VERBOSITY +castro.sum_interval = 1 # timesteps between computing mass +castro.v = 1 # verbosity in Castro.cpp +amr.v = 1 # verbosity in Amr.cpp + +# REFINEMENT / REGRIDDING +amr.max_level = 2 # maximum level number allowed +amr.ref_ratio = 4 2 2 2 # refinement ratio +amr.regrid_int = 2 2 2 2 # how often to regrid +amr.blocking_factor = 16 # block factor in grid generation +amr.max_grid_size = 128 +amr.n_error_buf = 2 2 2 2 # number of buffer cells in error est + +# CHECKPOINT FILES +amr.check_file = flame_wave_chk # root name of checkpoint file +amr.check_int = 250 # number of timesteps between checkpoints + +# PLOTFILES +amr.plot_file = flame_wave_plt # root name of plotfile +amr.plot_per = 2.e-3 # number of seconds between plotfiles +amr.derive_plot_vars = ALL + +amr.small_plot_file = flame_wave_smallplt # root name of plotfile +amr.small_plot_per = 5.e-5 # number of seconds between plotfiles +amr.small_plot_vars = density Temp enuc +amr.derive_small_plot_vars = abar x_velocity y_velocity z_velocity + +#PROBIN FILENAME +amr.probin_file = probin.H_He diff --git a/Exec/science/flame_wave/probin.H_He b/Exec/science/flame_wave/probin.H_He new file mode 100644 index 0000000000..b73a05bb29 --- /dev/null +++ b/Exec/science/flame_wave/probin.H_He @@ -0,0 +1,53 @@ +&fortin + + dtemp = 1.2e9 + x_half_max = 2.048d4 + x_half_width = 2048.d0 + + dx_model = 10.d0 + + dens_base = 1.e6 + + T_star = 1.d8 + T_hi = 2.d8 + T_lo = 8.d6 + + H_star = 2000.d0 + atm_delta = 50.0 + + fuel1_name = "hydrogen-1" + fuel1_frac = 0.75d0 + + fuel2_name = "helium-4" + fuel2_frac = 0.25d0 + + ash1_name = "nickel-56" + ash1_frac = 1.0d0 + + low_density_cutoff = 1.d-4 + + cutoff_density = 2.5e4 + max_hse_tagging_level = 3 + max_base_tagging_level = 2 + + X_min = 1.e-2 + + x_refine_distance = 6.144e4 +/ + +&sponge + + sponge_upper_density = 1.d2 + sponge_lower_density = 1.d0 + sponge_timescale = 1.d-7 + +/ + + +&extern + + rtol_spec = 1.d-6 + atol_spec = 1.d-6 + + use_tables = T +/ \ No newline at end of file