diff --git a/src/lagrangian/intermediate/submodels/MPPIC/CorrectionLimitingMethods/absolute/absolute.H b/src/lagrangian/intermediate/submodels/MPPIC/CorrectionLimitingMethods/absolute/absolute.H
index 85fc18054e..637d99ec44 100644
--- a/src/lagrangian/intermediate/submodels/MPPIC/CorrectionLimitingMethods/absolute/absolute.H
+++ b/src/lagrangian/intermediate/submodels/MPPIC/CorrectionLimitingMethods/absolute/absolute.H
@@ -25,9 +25,11 @@ Class
     Foam::CorrectionLimitingMethods::absolute
 
 Description
-    Correction limiting method that limits the velocity correction to that of a
-    rebound with a coefficient of restitution $e$. The absolute velocity of the
-    particle is used when calculating the magnitude of the limited correction.
+    Correction limiting method based on the absolute particle velocity.
+
+    This method that limits the velocity correction to that of a rebound with a
+    coefficient of restitution \f$e\f$. The absolute velocity of the particle
+    is used when calculating the magnitude of the limited correction.
     The direction is calculated using the relative velocity.
 
 SourceFiles
diff --git a/src/lagrangian/intermediate/submodels/MPPIC/CorrectionLimitingMethods/relative/relative.H b/src/lagrangian/intermediate/submodels/MPPIC/CorrectionLimitingMethods/relative/relative.H
index 73c772f723..0d84e4eb39 100644
--- a/src/lagrangian/intermediate/submodels/MPPIC/CorrectionLimitingMethods/relative/relative.H
+++ b/src/lagrangian/intermediate/submodels/MPPIC/CorrectionLimitingMethods/relative/relative.H
@@ -25,10 +25,12 @@ Class
     Foam::CorrectionLimitingMethods::relative
 
 Description
-    Correction limiting method that limits the velocity correction to that of a
-    rebound with a coefficient of restitution $e$. The relative velocity of the
-    particle with respect to the mean value is used to calculate the direction
-    and magnitude of the limited velocity.
+    Correction limiting method based on the relative particle velocity.
+
+    This method limits the velocity correction to that of a rebound with a
+    coefficient of restitution \f$e\f$. The relative velocity of the particle
+    with respect to the mean value is used to calculate the direction and
+    magnitude of the limited velocity.
 
 SourceFiles
     relative.C
diff --git a/src/lagrangian/intermediate/submodels/MPPIC/PackingModels/Explicit/Explicit.H b/src/lagrangian/intermediate/submodels/MPPIC/PackingModels/Explicit/Explicit.H
index 054d53eaa5..15937ffb69 100644
--- a/src/lagrangian/intermediate/submodels/MPPIC/PackingModels/Explicit/Explicit.H
+++ b/src/lagrangian/intermediate/submodels/MPPIC/PackingModels/Explicit/Explicit.H
@@ -39,6 +39,7 @@ Description
         D Snider
         Journal of Computational Physics
         Volume 170, Issue 2, Pages 523-549, July 2001
+    \endverbatim
 
 SourceFiles
     Explicit.C
diff --git a/src/lagrangian/intermediate/submodels/MPPIC/ParticleStressModels/HarrisCrighton/HarrisCrighton.H b/src/lagrangian/intermediate/submodels/MPPIC/ParticleStressModels/HarrisCrighton/HarrisCrighton.H
index 3eb3702340..dd99caf882 100644
--- a/src/lagrangian/intermediate/submodels/MPPIC/ParticleStressModels/HarrisCrighton/HarrisCrighton.H
+++ b/src/lagrangian/intermediate/submodels/MPPIC/ParticleStressModels/HarrisCrighton/HarrisCrighton.H
@@ -29,12 +29,12 @@ Description
 
     The stress value takes the following form:
     \f[
-        \dfrac{P_s \alpha^\beta}{\max \left( \alpha_{pack} - \alpha , \epsilon
-        ( 1 - \alpha ) \right) }
+        \frac{P_s \alpha^\beta}{ \mathrm{max} \left( \alpha_{pack} - \alpha ,
+        \epsilon ( 1 - \alpha ) \right) }
     \f]
-    Here, $\alpha$ is the volume fraction of the dispersed phase, and the other
-    values are modelling constants. A small value $\epsilon$ is used to limit
-    the denominator to ensure numerical stability.
+    Here, \f$\alpha\f$ is the volume fraction of the dispersed phase, and the
+    other values are modelling constants. A small value \f$\epsilon\f$ is used
+    to limit the denominator to ensure numerical stability.
 
     Reference:
     \verbatim
diff --git a/src/lagrangian/intermediate/submodels/MPPIC/ParticleStressModels/Lun/Lun.H b/src/lagrangian/intermediate/submodels/MPPIC/ParticleStressModels/Lun/Lun.H
index c00787ce78..522acf3d25 100644
--- a/src/lagrangian/intermediate/submodels/MPPIC/ParticleStressModels/Lun/Lun.H
+++ b/src/lagrangian/intermediate/submodels/MPPIC/ParticleStressModels/Lun/Lun.H
@@ -29,13 +29,13 @@ Description
 
     The stress value takes the following form:
     \f[
-        \left[ \alpha \rho + \alpha^2 \rho (1 + e) \frac{3}{5}
-        \left( \frac{\alpha_{pack}}{\alpha_{pack} - \alpha}
-        \right)^\frac{1}{3} \right] \frac{1}{3} \sigma^2
+        \left( \alpha \rho + \alpha^2 \rho (1 + e) \frac{3}{5}
+        \left( 1 - \left( \frac{\alpha}{\alpha_{pack}} \right)^\frac{1}{3}
+        \right) \right) \frac{1}{3} \sigma^2
     \f]
-    Here, $\alpha$ is the volume fraction of the dispersed phase, $\rho$ is the
-    density of the dispersed phase, $e$ is a coefficient of restitution, and
-    $\sigma$ is the RMS velocityh fluctuation.
+    Here, \f$\alpha\f$ is the volume fraction of the dispersed phase,
+    \f$\rho\f$ is the density of the dispersed phase, \f$e\f$ is a coefficient
+    of restitution, and \f$\sigma\f$ is the RMS velocity fluctuation.
 
     Reference:
     \verbatim
diff --git a/src/lagrangian/intermediate/submodels/MPPIC/TimeScaleModels/nonEquilibrium/nonEquilibrium.H b/src/lagrangian/intermediate/submodels/MPPIC/TimeScaleModels/nonEquilibrium/nonEquilibrium.H
index 4eeb3ac482..1c48f9bed1 100644
--- a/src/lagrangian/intermediate/submodels/MPPIC/TimeScaleModels/nonEquilibrium/nonEquilibrium.H
+++ b/src/lagrangian/intermediate/submodels/MPPIC/TimeScaleModels/nonEquilibrium/nonEquilibrium.H
@@ -36,6 +36,7 @@ Description
         P O'Rourke and D Snider
         Chemical Engineering Science
         Volume 65, Issue 22, Pages 6014-6028, November 2010
+    \endverbatim
 
 SourceFiles
     nonEquilibrium.C