diff --git a/Exec/science/nova/README.md b/Exec/science/nova/README.md new file mode 100644 index 0000000000..d12152006f --- /dev/null +++ b/Exec/science/nova/README.md @@ -0,0 +1,16 @@ +# `nova` + +This sets up a classical nova envelope on the surface of a carbon/oxygen +white dwarf using initial models provided by Ami Glasner. + +It is recommended to begin with the `inputs_nova_t7` inputs file. + + +## Publications + +This problem setup was used in: + +* *Multidimensional Nova Simulations with an Extended Buffer and Lower + Initial Mixing Temperatures*, A. Smith Clark & M. Zingale, 2025, + Open Journal of Astrophysics, volume 8, + https://doi.org/10.33232/001c.136890 diff --git a/Exec/science/nova/inputs_nova b/Exec/science/nova/inputs_nova deleted file mode 100644 index b18afa3259..0000000000 --- a/Exec/science/nova/inputs_nova +++ /dev/null @@ -1,108 +0,0 @@ -# ------------------ INPUTS TO MAIN PROGRAM ------------------- -max_step = 100000 - -# PROBLEM SIZE & GEOMETRY -geometry.is_periodic = 1 0 -geometry.coord_sys = 0 # 0 => cart, 1 => RZ 2=>spherical -geometry.prob_lo = 0 0 -geometry.prob_hi = 8.192e7 4.096e7 -amr.n_cell = 256 128 - -# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< -# 0 = Interior 3 = Symmetry -# 1 = Inflow 4 = SlipWall -# 2 = Outflow 5 = NoSlipWall -# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< -castro.lo_bc = 0 1 -castro.hi_bc = 0 2 - -castro.yl_ext_bc_type = 1 -castro.hse_interp_temp = 1 -castro.hse_reflect_vels = 1 - -castro.fill_ambient_bc = 1 -castro.ambient_fill_dir = 1 -castro.ambient_outflow_vel = 1 - -# WHICH PHYSICS -castro.do_hydro = 1 #hydrodynamics -castro.do_react = 1 #microphysics -castro.add_ext_src = 0 #external sources -castro.do_grav = 1 #gravity -castro.do_sponge = 0 #using sponge -castro.ppm_type = 1 -castro.riemann_solver = 0 -castro.use_flattening = 1 - -castro.small_temp = 1.0e6 -gravity.gravity_type = ConstantGrav -gravity.const_grav = -7.34e8 -castro.react_rho_min = 1.0e2 -castro.react_T_min = 5.0e6 - -# TIME STEP CONTROL -castro.cfl = 0.5 # cfl number for hyperbolic system -#castro.init_shrink = 0.1 # scale back initial timestep -castro.change_max = 1.1 # max time step growth - -# 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 = 2 2 2 2 # refinement ratio -amr.regrid_int = 2 2 2 2 # how often to regrid -amr.blocking_factor = 64 # block factor in grid generation -amr.max_grid_size = 256 -amr.n_error_buf = 2 2 2 2 # number of buffer cells in error est - -# CHECKPOINT FILES -amr.check_file = chk # root name of checkpoint file -amr.check_int = 200 # number of timesteps between checkpoints - -# PLOTFILES -amr.plot_file = plt # root name of plotfile -amr.plot_int = 1000 # number of timesteps between plotfiles -amr.derive_plot_vars = ALL - -# DATA LOG -amr.data_log = log - -# problem initialization - -problem.model_name = "glasner_T51600.hse" - -problem.apply_vel_field = 1 -problem.velpert_scale = 2.0e5 -problem.velpert_amplitude = 1.0e6 -problem.num_vortices = 4 - -# Microphysics - -integrator.jacobian = 2 -integrator.rtol_spec = 1.e-6 -integrator.atol_spec = 1.e-6 -integrator.rtol_enuc = 1.e-6 - -# refinement - -amr.refinement_indicators = dengrad tempgrad - -amr.refine.dengrad.gradient = 30.0 -amr.refine.dengrad.field_name = density -amr.refine.dengrad.max_level = 2 - -#amr.refine.denerr.value_greater = 0.1 -#amr.refine.denerr.field_name = density -#amr.refine.denerr.max_level = 2 - -amr.refine.tempgrad.gradient = 20000.0 -amr.refine.tempgrad.field_name = Temp -amr.refine.tempgrad.max_level = 2 - -#amr.refine.temperr.value_greater = 1.e6 -#amr.refine.temperr.field_name = Temp -#amr.refine.temperr.max_level = 2 diff --git a/Exec/science/nova/inputs_nova_en b/Exec/science/nova/inputs_nova_en deleted file mode 100644 index 747d789f2b..0000000000 --- a/Exec/science/nova/inputs_nova_en +++ /dev/null @@ -1,110 +0,0 @@ -# ------------------ INPUTS TO MAIN PROGRAM ------------------- -max_step = 100000 - -# PROBLEM SIZE & GEOMETRY -geometry.is_periodic = 1 0 -geometry.coord_sys = 0 # 0 => cart, 1 => RZ 2=>spherical -geometry.prob_lo = 0 0 -geometry.prob_hi = 8.192e7 4.096e7 -amr.n_cell = 256 128 - -# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< -# 0 = Interior 3 = Symmetry -# 1 = Inflow 4 = SlipWall -# 2 = Outflow 5 = NoSlipWall -# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< -castro.lo_bc = 0 1 -castro.hi_bc = 0 2 - -castro.yl_ext_bc_type = 1 -castro.hse_interp_temp = 1 -castro.hse_reflect_vels = 1 - -castro.fill_ambient_bc = 1 -castro.ambient_fill_dir = 1 -castro.ambient_outflow_vel = 1 - -# WHICH PHYSICS -castro.do_hydro = 1 #hydrodynamics -castro.do_react = 1 #microphysics -castro.add_ext_src = 0 #external sources -castro.do_grav = 1 #gravity -castro.do_sponge = 0 #using sponge -castro.ppm_type = 1 -castro.riemann_solver = 0 -castro.use_flattening = 1 - -castro.small_temp = 1.0e5 -castro.small_dens = 1.0e-6 -castro.retry_small_density_cutoff = 1.0 - -gravity.gravity_type = ConstantGrav -gravity.const_grav = -8.84e8 -castro.react_rho_min = 1.0e2 -castro.react_T_min = 1.0e6 - -# TIME STEP CONTROL -castro.cfl = 0.5 # cfl number for hyperbolic system -#castro.init_shrink = 0.1 # scale back initial timestep -castro.change_max = 1.1 # max time step growth - -# 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 = 1 # maximum level number allowed -amr.ref_ratio = 2 2 2 2 # refinement ratio -amr.regrid_int = 2 2 2 2 # how often to regrid -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 = chk # root name of checkpoint file -amr.check_int = 200 # number of timesteps between checkpoints - -# PLOTFILES -amr.plot_file = plt # root name of plotfile -amr.plot_int = 500 # number of timesteps between plotfiles -amr.derive_plot_vars = ALL - -# DATA LOG -amr.data_log = log - -# problem initialization - -problem.model_name = "glasner_enriched1600.hse" - -problem.apply_vel_field = 0 -#problem.velpert_scale = 5.0e5 -#problem.velpert_scale_ratio = 7.0 -#problem.velpert_amplitude = 7.0e5 -#problem.num_vortices = 8 - -problem.apply_temp_field = 1 -problem.width = 1.0e6 -problem.amplitude = 0.5 -problem.num_vortices = 10 -problem.width_ratio = 1.8 - -# Microphysics - -integrator.jacobian = 2 -integrator.rtol_spec = 1.e-6 -integrator.atol_spec = 1.e-6 -integrator.rtol_enuc = 1.e-6 - -# refinement - -amr.refinement_indicators = dens temp - -amr.refine.dens.max_level = 1 -amr.refine.dens.value_greater = 1.0e2 -amr.refine.dens.field_name = density - -amr.refine.temp.max_level = 1 -amr.refine.temp.value_greater = 1.0e7 -amr.refine.temp.field_name = Temp \ No newline at end of file