From 6638102829cd866f0342bf5d503422fa9a43f2ac Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Tue, 10 Nov 2020 14:36:15 +0100
Subject: [PATCH] [MPL/StrainDependentPermeability] Added a description of the
 model

---
 Documentation/bibliography/ogs.bib            | 18 ++++++++++
 .../Properties/StrainDependentPermeability.h  | 36 +++++++++++++++++++
 2 files changed, 54 insertions(+)

diff --git a/Documentation/bibliography/ogs.bib b/Documentation/bibliography/ogs.bib
index 24711e036db..033ebfbc598 100644
--- a/Documentation/bibliography/ogs.bib
+++ b/Documentation/bibliography/ogs.bib
@@ -194,3 +194,21 @@
   publisher={Springer},
   category={Processes},
 }
+
+@inproceedings{xu2011simulation,
+  title={Simulation of dilatancy-controlled gas migration process in saturated argillaceous rock},
+  author={Xu, WJ and Shao, H and Hesser, J and Wang, WQ and Kolditz, O and Popp, T},
+  booktitle={Computational geomechanics, COMGEO II. Proceedings of the 2nd international symposium on computational geomechanics},
+  pages={693--703},
+  year={2011}
+}
+
+@article{xu2013coupled,
+  title={Coupled multiphase flow and elasto-plastic modelling of in-situ gas injection experiments in saturated claystone (Mont Terri Rock Laboratory)},
+  author={Xu, WJ and Shao, H and Hesser, J and Wang, W and Schuster, K and Kolditz, O},
+  journal={Engineering Geology},
+  volume={157},
+  pages={55--68},
+  year={2013},
+  publisher={Elsevier}
+}
diff --git a/MaterialLib/MPL/Properties/StrainDependentPermeability.h b/MaterialLib/MPL/Properties/StrainDependentPermeability.h
index dc2549c8721..e2ef82ee83b 100644
--- a/MaterialLib/MPL/Properties/StrainDependentPermeability.h
+++ b/MaterialLib/MPL/Properties/StrainDependentPermeability.h
@@ -23,6 +23,42 @@ struct Parameter;
 
 namespace MaterialPropertyLib
 {
+/**
+ * \brief A strain dependent intrinsic permeability model.
+ *
+ *  The model was proposed
+ *  in \cite xu2011simulation and it was further investigated
+ *  in \cite xu2013coupled .
+ *
+ *   The model takes the form of
+ *  \f[ \mathbf{k} =f(\epsilon_v) e^{b_1  {\bar\epsilon}^p}\mathbf{k}_0
+ *  \f]
+ *  with
+ *  \f[ f(\epsilon_v)=
+ *      \begin{cases}
+ *       10^{b_2  \epsilon_v}, & \epsilon_v <=0\\
+ *      10^{b_3  \epsilon_v}, & \epsilon_v >0
+ *      \end{cases}
+ *  \f]
+ *   where
+ *   <table>
+ *   <tr><td>\f$ \epsilon_v \f$  <td> the volumetric strain,
+ *   <tr><td> \f$ {\bar\epsilon}^p\f$ <td> the equivalent plastic
+ *    strain,
+ *   <tr><td>\f$\mathbf{k}_0\f$  <td> the initial intrinsic permeability,
+ *   <tr><td>\f$b_1,\,b_2,\,b_3\f$  <td> the three parameters.
+ * </table>
+ *
+ *  * Note: In \cite xu2011simulation  and \cite xu2013coupled ,
+ * from the point of view of experiment of permeability change, the symbols of
+ * \f$ \Delta \epsilon_v \f$ and \f$\Delta {\bar\epsilon}^p\f$  are used for \f$
+ * \epsilon_v \f$ and \f$
+ * {\bar\epsilon}^p\f$, respectively. That means that the symbols of
+ * \f$ \Delta \epsilon_v \f$ and \f$\Delta {\bar\epsilon}^p\f$
+ *  refer to the increment of the volumetric strain and the increment of the
+ * equivalent plastic strain, respectively, from the beginning of the
+ * experiment.
+ */
 template <int DisplacementDim>
 class StrainDependentPermeability final : public Property
 {
-- 
GitLab