diff --git a/Documentation/bibliography/other.bib b/Documentation/bibliography/other.bib
index 1d855d9ae13408f589224cdcd2fcd33115d80b10..93600dc78789790ddc49b2b193927532de21f027 100644
--- a/Documentation/bibliography/other.bib
+++ b/Documentation/bibliography/other.bib
@@ -209,3 +209,14 @@ URL = {https://doi.org/10.1680/geot.2008.58.3.157}
   year={1985},
   publisher={Elsevier}
 }
+
+@article{lenhard1987model,
+  title={A model for hysteretic constitutive relations governing multiphase flow: 2. Permeability-saturation relations},
+  author={Lenhard, RJ and Parker, JC},
+  journal={Water Resources Research},
+  volume={23},
+  number={12},
+  pages={2197--2206},
+  year={1987},
+  publisher={Wiley Online Library}
+}
diff --git a/MaterialLib/MPL/Properties/RelativePermeability/RelPermNonWettingPhaseVanGenuchtenMualem.h b/MaterialLib/MPL/Properties/RelativePermeability/RelPermNonWettingPhaseVanGenuchtenMualem.h
index 2e1014b1e3ecaa03221b33a51bd0e5bb3210a54a..8c35d7fa4b3ad4b305a055a75cd7fcca01f5fb8b 100644
--- a/MaterialLib/MPL/Properties/RelativePermeability/RelPermNonWettingPhaseVanGenuchtenMualem.h
+++ b/MaterialLib/MPL/Properties/RelativePermeability/RelPermNonWettingPhaseVanGenuchtenMualem.h
@@ -20,6 +20,20 @@ namespace MaterialPropertyLib
 {
 class Medium;
 
+/**
+ *  van Genuchten-Mualem relative permeability function for non-wetting phase
+ *  in terms of effective wetting-phase saturation \cite lenhard1987model :
+ *
+ *   \f[k_{rel}^n= (1 - S_e)^{1/2} (1 - S_e^{1/m})^{2m}\f]
+ *   with
+ *   \f[S_e=\frac{S^L-S^L_r}{S^L_{\mbox{max}}-S^L_r}\f]
+ *   where
+ *    \f{eqnarray*}{
+ *       &S^L_r&            \mbox{residual saturation of wetting phase,}\\
+ *       &S^L_{\mbox{max}}& \mbox{maximum saturation of wetting phase,}\\
+ *       &m\, \in (0, 1) &    \mbox{ exponent.}\\
+ *    \f}
+ */
 class RelPermNonWettingPhaseVanGenuchtenMualem final : public Property
 {
 public: