conductive_flame¶
This test suite case models a conductively-propagated deflagration wave (“flame”) in a high-density, degenerate carbon-oxygen mixture. It is similar to a calculation presented by Timmes & Woosley (1992).
Warning
Careful modeling of the flame speed requires a much larger nuclear network than the 21 isotope network used in this test. Timmes & Woosley (1992) use a 130 isotope network, Chamulak et al. (2007) use a 430 isotope network, and Schwab et al. (2020) use a 495 isotope network.
Unlike most MESA calculations, this models a small sphere of material
(instead of an entire star). The initial model is built in
run_star_extras
using the other_build_initial_model
hook. A
spatially uniform model with a given density, temperature, and
composition is constructed. A small hot spot is then added at the
center of the model. The properties of this initial model can be
controlled from the inlist.
use_other_build_initial_model = .true.
x_integer_ctrl(1) = 300 ! number of points
x_ctrl(1) = 10 ! mass (in g)
x_ctrl(2) = 1d10 ! initial density (cgs)
x_ctrl(3) = 2d8 ! initial temperature (cgs)
x_ctrl(4) = 0.0001 ! size (in q) of region to heat
x_ctrl(5) = 5d9 ! temperature of that region
x_ctrl(6) = 0.50 ! mass fraction c12
x_ctrl(7) = 0.50 ! mass fraction o16
The inner boundary is at r = 0. The outer boundary has a fixed
temperature and a fixed pressure equal to the initial pressure of the
material. This is achieved via the use_other_surface_PT
hook.
After an initial transient, the entire flame structure, approximately isobaric, propagates into the upstream fuel with a unique speed and width. The test succeeds if the flame successfully propagates through half of the domain.
Last-Updated: 2019-11-16 (mesa r12307) by Josiah Schwab