diff --git a/docs/source_docs/conf.py b/docs/source_docs/conf.py
index 82aa9ba742cb265a3331558791c7690ba7944610..e4bbc5f011be06102d1e68748ad1790406511989 100644
--- a/docs/source_docs/conf.py
+++ b/docs/source_docs/conf.py
@@ -86,6 +86,9 @@ todo_include_todos = False
numfig = True
math_numfig = True
+numfig_secnum_depth = 1
+
+math_eqref_format = "Eq.{number}"
bibtex_bibfiles = ['refs.bib']
diff --git a/docs/source_docs/index.rst b/docs/source_docs/index.rst
index 939c4e9599d2af1f93f94b1293acdffbaeffbe96..b5be3101c893d972cd888964053df35b85be8ec7 100644
--- a/docs/source_docs/index.rst
+++ b/docs/source_docs/index.rst
@@ -81,6 +81,7 @@ and performance:
test_benchmarks/CITests
test_benchmarks/NightlyTests
test_benchmarks/qualitative_bencharks/index
+ test_benchmarks/verification/index
Notice
------
diff --git a/docs/source_docs/test_benchmarks/verification/chemistry/eulerian01.rst b/docs/source_docs/test_benchmarks/verification/chemistry/eulerian01.rst
new file mode 100644
index 0000000000000000000000000000000000000000..0ceb3a1d0aab0e90d649081d9f9f6d39d6b02d4f
--- /dev/null
+++ b/docs/source_docs/test_benchmarks/verification/chemistry/eulerian01.rst
@@ -0,0 +1,93 @@
+.. include:: ../../../definitions.rst
+
+Eulerian01: athermal, constant density
+--------------------------------------
+
+Description
+~~~~~~~~~~~
+
+Consider a system with a fluid comprised of two chemical species, :math:`A` and :math:`B`, with the transfer of species given by the reversible homogeneous reaction,
+
+.. math::
+ :label: chem_eulerian01_eq01
+
+ A(g) \leftrightarrow B(g)
+
+If both species have the same molecular weight, specific heat, and enthalpy of formation, then the mixture molecular weight, system temperature, and density are constant.
+The evolution of species mass in the absence of convection and diffusion is given by
+
+.. math::
+ :label: chem_eulerian01_eq02
+
+ \rho\frac{\partial X_{A}}{\partial t} = \dot{\omega}_{A}
+
+where :math:`\rho` is density;
+:math:`X_A` is the mass fraction of species :math:`A`; and
+:math:`\dot{\omega}_{A}` is the net rate of production (consumption if negative) of species :math:`A`.
+
+A manufactured solution for species :math:`A` is
+
+.. math::
+ :label: chem_eulerian01_eq03
+
+ X_A(t) = X_A^0 - \sin(6\pi t)/4
+
+where :math:`X_A^0 = X_A(t=0)` is the initial mass fraction of species :math:`A`. Taking the derivative with respect to time, and substituting the result into :eq:`chem_eulerian01_eq02` gives
+
+.. math::
+ :label: chem_eulerian01_eq04
+
+ \dot{\omega}_A = -\rho 3\pi \cos(6\pi t)/2
+
+which has units of kg\ |cdot|\ m\ :math:`^{-3}`\ |cdot| sec\ :math:`^{-1}`. Two modifications are required to implement this Eulerian reaction rate into MFIX-Exa:
+
+- Eq. :eq:`chem_eulerian01_eq04` needs to be divided by the molecular weight of species A because Eulerian reaction rates in MFIX-Exa are provided in units of mol\ |cdot|\ m\ :math:`^{-3}`\ |cdot| sec\ :math:`^{-1}`.
+- The negative sign is omitted because species :math:`A` is a reactant. Specifically, by inspection of chemical equation :eq:`chem_eulerian01_eq01`, species :math:`A` is a reactant when the rate is positive (forward reaction), and a product when the rate is negative (backward reaction).
+
+The resulting Eulerian reaction rate is
+
+.. math::
+ :label: chem_eulerian01_eq54
+
+ R_A = \frac{\rho 3\pi \cos(6\pi t)}{2 M_A}
+
+Setup
+~~~~~
+
+The domain is a single cuboid cell with 0.0025 m edge length and convection and diffusion terms are set to zero for all governing equations. The fluid has a constant temperature of 300 K at 60,795 Pa, and the initial composition is :math:`X_A=X_B=0.5`. The density is chosen to be 2 kg\ |cdot|\ m\ :math:`^{-3}`; therefore, the component and mixture molecular weights are computed from the ideal gas law :math:`M=\rho RT / p_{therm} \approx 0.0821` kg\ |cdot|\ mol\ :math:`^{-1}`.
+
+
+Results
+~~~~~~~
+
+The test is run at five fixed time-step sizes with a refinement factor of 2: 0.2, 0.1, 0.05, 0.025, 0.0125 milliseconds. The :math:`L_1`\ -norm is computed by dividing the sum of the absolute difference between the analytical and numerical solution over all time steps by the total number of time steps.
+
+.. math::
+
+ L_1 = \frac{\sum_n \| X_A(t_n) - X_A^{\mathrm{\texttt{MFIX-Exa}}}\|}{n}
+
+
+Forward Euler
+'''''''''''''
+
+The analytical solution and absolute error with dt = 0.2 msec are shown in :numref:`fig_chem_eulerian01_ForwardEuler_sol`. The :math:`L_1` norm is used to compute the observed order of accuracy shown in :numref:`fig_chem_eulerian01_ForwardEuler_ooa`.
+
+
+.. _fig_chem_eulerian01_ForwardEuler_sol:
+
+.. figure:: ./images/eulerian01_ForwardEuler_solution.png
+ :width: 85%
+ :align: center
+ :alt: Eulerian01 test numerical solution and absolute error
+
+ Numerical solution for species mass fraction and absolute error.
+
+
+.. _fig_chem_eulerian01_ForwardEuler_ooa:
+
+.. figure:: ./images/eulerian01_ForwardEuler_accuracy.png
+ :width: 45%
+ :align: center
+ :alt: Eulerian01 test observed order
+
+ Observed order of species mass fraction.
diff --git a/docs/source_docs/test_benchmarks/verification/chemistry/images/eulerian01_ForwardEuler_accuracy.png b/docs/source_docs/test_benchmarks/verification/chemistry/images/eulerian01_ForwardEuler_accuracy.png
new file mode 100644
index 0000000000000000000000000000000000000000..030b20f09a875e48c41589c9cc3beeda93f56eda
Binary files /dev/null and b/docs/source_docs/test_benchmarks/verification/chemistry/images/eulerian01_ForwardEuler_accuracy.png differ
diff --git a/docs/source_docs/test_benchmarks/verification/chemistry/images/eulerian01_ForwardEuler_solution.png b/docs/source_docs/test_benchmarks/verification/chemistry/images/eulerian01_ForwardEuler_solution.png
new file mode 100644
index 0000000000000000000000000000000000000000..463c3bc37674932ac38c8d22b410d304b8c4110d
Binary files /dev/null and b/docs/source_docs/test_benchmarks/verification/chemistry/images/eulerian01_ForwardEuler_solution.png differ
diff --git a/docs/source_docs/test_benchmarks/verification/index.rst b/docs/source_docs/test_benchmarks/verification/index.rst
new file mode 100644
index 0000000000000000000000000000000000000000..6d96dde0d626fb6e96dbd031c09c3a82f1f08bd1
--- /dev/null
+++ b/docs/source_docs/test_benchmarks/verification/index.rst
@@ -0,0 +1,27 @@
+====================
+Chemistry code tests
+====================
+
+The test cases presented in this chapter for the MFIX fluid solver are summarized below.
+
+.. "Eulerian01", "C", "A", " ", " ", ":math:`{\checkmark}`"
+
+
+.. |check| replace:: :math:`{\checkmark}`
+
++-------------------+---------+---------+--------+-----------+---------+---------+---------+
+| Test | Freq\ †| Data\ ‡ | Phases | Momentum | Density | Species | Energy |
++===================+=========+=========+========+===========+=========+=========+=========+
+| Eulerian01 | C | A | Fluid | | |check| | | |
++-------------------+---------+---------+--------+-----------+---------+---------+---------+
+
+†C-Incorporated into the continuous integration server; M-Monthly; Q-Quarterly; X-Manual; D-Disabled
+
+‡ A-Analytical solution; P-Published benchmark data; R-Regression data from previous code versions
+
+.. toctree::
+ :numbered:
+ :caption: Validation and Verification
+ :maxdepth: 1
+
+ chemistry/eulerian01