From 66e4f84da6d7168858e963f2eb4550694a0e93ef Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Thu, 17 Oct 2024 07:40:20 -0400 Subject: [PATCH 1/2] fix center for flame_wave in 1D --- Exec/science/flame_wave/problem_initialize.H | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/Exec/science/flame_wave/problem_initialize.H b/Exec/science/flame_wave/problem_initialize.H index 61eb43808c..8c0d1fe5f8 100644 --- a/Exec/science/flame_wave/problem_initialize.H +++ b/Exec/science/flame_wave/problem_initialize.H @@ -190,6 +190,10 @@ void problem_initialize () problem::center[d] = 0.5_rt * (problo[d] + probhi[d]); } +#if AMREX_SPACEDIM == 1 + problem::center[0] = 0.0_rt; +#endif + #if AMREX_SPACEDIM == 2 // for axisymmetry, put the x-center on the x-axis // and the y-center at 0, so the height computation is okay From 29aed12a4ddc5bf023c110839cf6e38b7758c989 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Thu, 17 Oct 2024 08:21:31 -0400 Subject: [PATCH 2/2] a 1D inputs to test HSE --- Exec/science/flame_wave/inputs.1d | 156 ++++++++++++++++++++++++++++++ 1 file changed, 156 insertions(+) create mode 100644 Exec/science/flame_wave/inputs.1d diff --git a/Exec/science/flame_wave/inputs.1d b/Exec/science/flame_wave/inputs.1d new file mode 100644 index 0000000000..63e5a27caa --- /dev/null +++ b/Exec/science/flame_wave/inputs.1d @@ -0,0 +1,156 @@ +# ------------------ INPUTS TO MAIN PROGRAM ------------------- +max_step = 250000000 +stop_time = 3.0 + +# PROBLEM SIZE & GEOMETRY +geometry.is_periodic = 0 0 0 +geometry.prob_lo = 0 0 0 + +geometry.coord_sys = 0 # 0 => cart, 1 => RZ 2=>spherical + + geometry.prob_hi = 2.4576e4 + amr.n_cell = 192 + + # >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< + # 0 = Interior 3 = Symmetry + # 1 = Inflow 4 = SlipWall + # 2 = Outflow 5 = NoSlipWall + # >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< + castro.lo_bc = 3 + castro.hi_bc = 2 + + castro.fill_ambient_bc = 1 + castro.ambient_fill_dir = 0 + castro.ambient_outflow_vel = 1 + +#castro.allow_non_unit_aspect_zones = 1 + +castro.domain_is_plane_parallel = 1 + +# WHICH PHYSICS +castro.do_hydro = 1 +castro.do_react = 1 +castro.do_rotation = 0 +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 +castro.use_pslope = 1 +castro.pslope_cutoff_density = 1.e4 + +gravity.gravity_type = ConstantGrav +gravity.const_grav = -1.5e14 + +castro.rotational_period = 0.0 +castro.rotation_include_centrifugal = 0 + +castro.diffuse_temp = 0 +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 = 5.e6 + +castro.react_T_min = 6.e7 + +castro.sponge_upper_density = 1.e2 +castro.sponge_lower_density = 1.e0 +castro.sponge_timescale = 1.e-7 + +# GPU options +castro.hydro_memory_footprint_ratio = 3 + +# 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.use_retry = 1 +castro.max_subcycles = 16 + +castro.retry_small_density_cutoff = 10.0 + +castro.abundance_failure_tolerance = 0.1 +castro.abundance_failure_rho_cutoff = 1.0 + +# DIAGNOSTICS & VERBOSITY +castro.sum_interval = 10 # 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 = 0 # static grids ftw +amr.blocking_factor = 32 # 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_1000Hz_25cm_chk # root name of checkpoint file +#if AMREX_SPACEDIM == 3 +amr.check_int = 25 # number of timesteps between checkpoints +#elif AMREX_SPACEDIM == 2 +amr.check_int = 250 # number of timesteps between checkpoints +#endif + +amr.checkpoint_files_output = 1 + +# PLOTFILES +amr.plot_file = flame_wave_1000Hz_25cm_plt # root name of plotfile +amr.plot_per = 5.e-3 # number of seconds between plotfiles +amr.derive_plot_vars = ALL +amr.plot_files_output = 1 + +amr.small_plot_file = flame_wave_1000Hz_25cm_smallplt # root name of plotfile +amr.small_plot_per = 5.e-4 # number of seconds between plotfiles +amr.small_plot_vars = density Temp + +#if AMREX_SPACEDIM == 3 + amr.derive_small_plot_vars = abar magvel enuc +#elif AMREX_SPACEDIM == 2 + amr.derive_small_plot_vars = abar x_velocity y_velocity z_velocity enuc +#endif + +fab.format = NATIVE_32 + + +# problem initialization + +problem.dtemp = 1.2e9 +problem.x_half_max = 2.56e4 +problem.x_half_width = 2048.e0 + +problem.dens_base = 3.43e6 + +problem.T_star = 2.e8 +problem.T_hi = 2.e8 +problem.T_lo = 8.e6 + +problem.H_star = 2000.e0 +problem.atm_delta = 50.0 + +problem.fuel1_name = "helium-4" +problem.fuel1_frac = 1.0 + +problem.ash1_name = "nickel-56" +problem.ash1_frac = 1.0 + +problem.low_density_cutoff = 1.e-4 + +problem.tag_by_density = 0 +problem.refine_height = 3600.0 +problem.max_base_tagging_level = 3 + +# Microphysics + +integrator.rtol_spec = 1.e-6 +integrator.atol_spec = 1.e-6 + +network.use_tables = 1