add figures and references

docs/references.bib

@@ -263,3 +263,20 @@ @ARTICLE{Jin_1996
       adsnote = {Provided by the SAO/NASA Astrophysics Data System}
 }
 
+@article{Krumholz_2007,
+  title = {Equations and {{Algorithms}} for {{Mixed-frame Flux-limited Diffusion Radiation Hydrodynamics}}},
+  author = {{Krumholz}, Mark R. and {Klein}, Richard I. and {McKee}, Christopher F. and {Bolstad}, John},
+  year = {2007},
+  month = sep,
+  journal = {The Astrophysical Journal},
+  volume = {667},
+  pages = {626--643},
+  issn = {0004-637X},
+  doi = {10.1086/520791},
+  url = {https://ui.adsabs.harvard.edu/abs/2007ApJ...667..626K},
+  urldate = {2023-10-18},
+  abstract = {We analyze the mixed-frame equations of radiation hydrodynamics under the approximations of flux-limited diffusion and a thermal radiation field and derive the minimal set of evolution equations that includes all terms that are of leading order in any regime of nonrelativistic radiation hydrodynamics. Our equations are accurate to first order in v/c in the static diffusion regime. In contrast, we show that previous lower order derivations of these equations omit leading terms in at least some regimes. In comparison to comoving-frame formulations of radiation hydrodynamics, our equations have the advantage that they manifestly conserve total energy, making them very well suited to numerical simulations, particularly with adaptive meshes. For systems in the static diffusion regime, our analysis also suggests an algorithm that is both simpler and faster than earlier comoving-frame methods. We implement this algorithm in the Orion adaptive mesh refinement code and show that it performs well in a range of test problems.},
+  keywords = {Astrophysics,Hydrodynamics,Methods: Numerical,Radiative Transfer},
+  annotation = {ADS Bibcode: 2007ApJ...667..626K},
+  file = {/Users/cche/Googledrive2/Academic/Zotero/Code/QUOKKA/Krumholz_plus_2007_Equations_and_Algorithms_for_Mixed-frame2.pdf}
+}

docs/tests/index.rst

@@ -10,4 +10,6 @@ Listed here are the test problems that are included with Quokka. *This page is s
    shu_osher
    sms
    energy_exchange
+   radhydro_uniform_adv
+   radhydro_pulse
 

docs/tests/radhydro_pulse.rst

@@ -9,7 +9,7 @@ the RHD moment equations, in a fully-coupled RHD problem. The problems
 involve the advection of the a pulse of radiation energy in an optically
 thick (:math:`\tau \gg 1`) gas in both static (:math:`\beta \tau \gg 1`)
 and dynamic (:math:`\beta \tau \ll 1`) diffusion regimes, with a uniform
-background flow velocity Krumholz+07.
+background flow velocity :cite:`Krumholz_2007`.
 
 Parameters
 ----------
@@ -51,33 +51,33 @@ Results
 Static diffusion regime:
 
 .. figure:: attach/radhydro_pulse_temperature-1.png
-   :alt: radhydro_pulse_temperature-1
+   :alt: radhydro_pulse_temperature-static-diffusion
 
-   radhydro_pulse_temperature-1
+   radhydro_pulse_temperature-static-diffusion
 
 .. figure:: attach/radhydro_pulse_density-1.png
-   :alt: radhydro_pulse_density-1
+   :alt: radhydro_pulse_density-static-diffusion
 
-   radhydro_pulse_density-1
+   radhydro_pulse_density-static-diffusion
 
 .. figure:: attach/radhydro_pulse_velocity-1.png
-   :alt: radhydro_pulse_velocity-1
+   :alt: radhydro_pulse_velocity-static-diffusion
 
-   radhydro_pulse_velocity-1
+   radhydro_pulse_velocity-static-diffusion
 
 Dynamic diffusion regime:
 
 .. figure:: attach/radhydro_pulse_temperature.png
-   :alt: radhydro_pulse_temperature
+   :alt: radhydro_pulse_temperature-dynamic-diffusion
 
-   radhydro_pulse_temperature
+   radhydro_pulse_temperature-dynamic-diffusion
 
 .. figure:: attach/radhydro_pulse_density.png
-   :alt: radhydro_pulse_density
+   :alt: radhydro_pulse_density-dynamic-diffusion
 
-   radhydro_pulse_density
+   radhydro_pulse_density-dynamic-diffusion
 
 .. figure:: attach/radhydro_pulse_velocity.png
-   :alt: radhydro_pulse_velocity
+   :alt: radhydro_pulse_velocity-dynamic-diffusion
 
-   radhydro_pulse_velocity
+   radhydro_pulse_velocity-dynamic-diffusion

docs/tests/radhydro_uniform_adv.rst

@@ -28,13 +28,28 @@ Results
 
 With :math:`O(\beta \tau)` terms:
 
-|radhydro_uniform_advecting_temperature|
-|radhydro_uniform_advecting_velocity|
+.. figure:: attach/radhydro_uniform_advecting_temperature.png
+    :alt: A figure showing the radiation temperature and material temperature as a function of time.
+
+    The radiation temperature and matter temperatures, along with the exact solution.
+
+.. figure:: attach/radhydro_uniform_advecting_velocity.png
+    :alt: A figure showing the radiation velocity and material velocity as a function of time.
+
+    The matter velocity, along with the exact solution.
 
 Without :math:`O(\beta \tau)` terms:
 
-|radhydro_uniform_advecting_temperature-nobeta|
-|radhydro_uniform_advecting_velocity-nobeta|
+.. figure:: attach/radhydro_uniform_advecting_temperature-nobeta.png
+    :alt: A figure showing the radiation temperature and material temperature as a function of time.
+
+    The radiation temperature and matter temperatures, along with the exact solution.
+
+.. figure:: attach/radhydro_uniform_advecting_velocity-nobeta.png
+    :alt: A figure showing the radiation velocity and material velocity as a function of time.
+
+    The matter velocity, along with the exact solution.
+
 
 Physics
 -------