add shock burning ref to docs (#3165)

this adds Ryan's RNAAS to the section on shock burning

Author: zingale

Docs/source/reactions.rst

@@ -93,7 +93,10 @@ compiled with::
 
 in the ``GNUmakefile``.  This will allocate storage for a shock flag in the conserved
 state array.  This flag is computed via a multidimensional shock detection algorithm
-described in :cite:`doubledet2024`.  A zone is tagged as a shock if the following
+described in :cite:`doubledet2024`.  A test problem demonstrating the effect of
+this flag is also shown in :cite:`brady:2025`.
+
+A zone is tagged as a shock if the following
 conditions are true:
 
 .. math::
@@ -355,5 +358,3 @@ In ``Castro_react.cpp``, the flow is:
     * $U(\mu_p) = \mathtt{burn\_state.mu\_p}$
 
     * $U(\mu_n) = \mathtt{burn\_state.mu\_n}$
-
-

Docs/source/refs.bib

@@ -1150,3 +1150,38 @@ @ARTICLE{amrex-ecp
       adsnote = {Provided by the SAO/NASA Astrophysics Data System}
 }
 
+
+@article{brady:2025,
+doi = {10.3847/2515-5172/add686},
+url = {https://doi.org/10.3847/2515-5172/add686},
+year = {2025},
+month = {may},
+publisher = {The American Astronomical Society},
+volume = {9},
+number = {5},
+pages = {113},
+author = {Brady, Ryan and Zingale, Michael},
+title = {Numerical Treatment of Shock-induced Nuclear Burning in Double Detonation Type Ia Supernovae},
+journal = {Research Notes of the AAS},
+abstract = {We present a benchmark problem to assess the treatment of
+                  shock-induced nuclear burning in the context of
+                  double detonation Type Ia supernovae. In a
+                  stratified white dwarf model, we implement a
+                  shock-detection criterion that suppresses burning in
+                  zones characterized by compression and significant
+                  pressure gradients, controlled by a tunable
+                  parameter, fshock. One-dimensional simulations,
+                  using the open-source Castro suite, were conducted
+                  across three treatments—burning fully enabled, and
+                  burning suppressed with fshock = 2/3 and
+                  fshock = 1—across three spatial resolutions (5.0,
+                  2.5, and 0.3125 km). At the finest resolution, the
+                  burning-enabled and fshock = 1 models converge,
+                  while the fshock = 2/3 front continues to show
+                  slight offset behavior. Since most simulations are
+                  carried out at much lower resolutions, our tests
+                  support the idea that burning in shocks should
+                  always be disabled in practice. We also observe that
+                  the behavior of lower-resolution simulations remains
+                  extremely sensitive to the choice of fshock.}
+}