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).
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
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