From afc8ace232ed9396524bee58ddb5832284c9191b Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Dennis=20Gl=C3=A4ser?= Date: Wed, 13 Apr 2022 14:26:08 +0000 Subject: [PATCH 1/3] [doc][elastic] cleanup --- dumux/geomechanics/elastic/model.hh | 25 ++++++++++++++----------- 1 file changed, 14 insertions(+), 11 deletions(-) diff --git a/dumux/geomechanics/elastic/model.hh b/dumux/geomechanics/elastic/model.hh index 1b7ca4fea0..7b83074241 100644 --- a/dumux/geomechanics/elastic/model.hh +++ b/dumux/geomechanics/elastic/model.hh @@ -21,27 +21,30 @@ * \ingroup Elastic * \brief A geomechanical model. * - * This model describe the deformation of the solid body, using theory of linear elasticity: + * This model describes the deformation of a solid body using the theory of linear elasticity. + * The momentum balance equation of the solid body can be expressed by \f[ - \boldsymbol{\nabla}\cdot\boldsymbol{\sigma_s} + \mathbf{F} + \rho_s \mathbf{g} = \rho\ddot{\mathbf{u}} + \nabla\cdot\boldsymbol{\sigma} + \rho \mathbf{g} = \rho\ddot{\mathbf{u}}, \f] - * We assume the quasi-static conditions, so the acceleration term \f$\rho\ddot{\mathbf{u}} \approx 0\f$. + * where \f$\boldsymbol{\sigma} \f$ is the stress tensor, \f$\rho \f$ is the density of + * the solid, and \f$\mathbf{u} = \mathbf{x} - \mathbf{x}_{\mathrm{initial}} \f$ is the displacement, + * defined as the difference in material points \f$\mathbf{x} \f$ and \f$\mathbf{x}_{\mathrm{initial}} \f$ + * in the deformed and undeformed (initial) state, respectively. The model assumes quasi-static conditions, + * that is, the above momentum balance equation is solved under the assumption that the acceleration term + * \f$\rho\ddot{\mathbf{u}} \approx 0\f$. * - * For isotropic materials, the stress tensor \f$\boldsymbol{\sigma_s} \f$ can be calculated after Hookes' Law + * Per default, Hookes' Law is used for expressing the stress tensor \f$\boldsymbol{\sigma} \f$ as a function of the + * displacement: \f[ - {\boldsymbol {\sigma_s}}= \lambda \mathrm{tr}(\varepsilon) \boldsymbol{\mathrm I}+ 2G \varepsilon, + \boldsymbol{\sigma} = \lambda\mathrm{tr}(\boldsymbol{\varepsilon}) \mathbf{I} + 2G \boldsymbol{\varepsilon}, \f] * with \f[ - {\boldsymbol {\varepsilon }}={\frac {1}{2}}\left[{\boldsymbol {\nabla }}\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )^{\mathrm {T} }\right]. + \boldsymbol{\varepsilon} = \frac{1}{2} \left[ \nabla\mathbf{u} + (\nabla\mathbf{u})^{\mathrm{T}} \right]. \f] * - * Gravity can be enabled or disabled via the property system. - * - * The equations are discretized using a vertex-centered finite volume (box) scheme as spatial discretization. The time discretization is not needed due to the quasi-static conditions. + * The equations are discretized using a vertex-centered finite volume (box) scheme as spatial discretization. * PrimaryVariables are the displacements in each direction \f$\mathbf{u} \f$. - * - * Pay attention to the sign convention, where tension has a positive! */ #ifndef DUMUX_GEOMECHANICS_ELASTIC_MODEL_HH #define DUMUX_GEOMECHANICS_ELASTIC_MODEL_HH -- GitLab From 458f96e1ccd7c80f69ef3b53c9c7f00d1b33f2fc Mon Sep 17 00:00:00 2001 From: Timo Koch Date: Fri, 22 Apr 2022 14:21:32 +0000 Subject: [PATCH 2/3] Fixup doc --- dumux/geomechanics/elastic/model.hh | 1 - 1 file changed, 1 deletion(-) diff --git a/dumux/geomechanics/elastic/model.hh b/dumux/geomechanics/elastic/model.hh index 7b83074241..cc1db6ccf0 100644 --- a/dumux/geomechanics/elastic/model.hh +++ b/dumux/geomechanics/elastic/model.hh @@ -43,7 +43,6 @@ \boldsymbol{\varepsilon} = \frac{1}{2} \left[ \nabla\mathbf{u} + (\nabla\mathbf{u})^{\mathrm{T}} \right]. \f] * - * The equations are discretized using a vertex-centered finite volume (box) scheme as spatial discretization. * PrimaryVariables are the displacements in each direction \f$\mathbf{u} \f$. */ #ifndef DUMUX_GEOMECHANICS_ELASTIC_MODEL_HH -- GitLab From 9072d64603d929a20d34b69a3ff381c5f46602e9 Mon Sep 17 00:00:00 2001 From: Yue Wang Date: Mon, 25 Apr 2022 19:55:50 +0200 Subject: [PATCH 3/3] [geomech] change the model descriptions --- dumux/geomechanics/elastic/model.hh | 9 ++++-- dumux/geomechanics/poroelastic/model.hh | 37 +++++++++++++++---------- 2 files changed, 28 insertions(+), 18 deletions(-) diff --git a/dumux/geomechanics/elastic/model.hh b/dumux/geomechanics/elastic/model.hh index cc1db6ccf0..bce39074e9 100644 --- a/dumux/geomechanics/elastic/model.hh +++ b/dumux/geomechanics/elastic/model.hh @@ -24,10 +24,11 @@ * This model describes the deformation of a solid body using the theory of linear elasticity. * The momentum balance equation of the solid body can be expressed by \f[ - \nabla\cdot\boldsymbol{\sigma} + \rho \mathbf{g} = \rho\ddot{\mathbf{u}}, + \nabla\cdot\boldsymbol{\sigma} + \rho \mathbf{g} + \mathbf{f} = \rho\ddot{\mathbf{u}}, \f] * where \f$\boldsymbol{\sigma} \f$ is the stress tensor, \f$\rho \f$ is the density of - * the solid, and \f$\mathbf{u} = \mathbf{x} - \mathbf{x}_{\mathrm{initial}} \f$ is the displacement, + * the solid, \f$\mathbf{f} \f$ in \f$\mathrm{N/m^3} \f$ is the external force acting on the body per unit volume (e.g. magnetism), + * and \f$\mathbf{u} = \mathbf{x} - \mathbf{x}_{\mathrm{initial}} \f$ is the displacement, * defined as the difference in material points \f$\mathbf{x} \f$ and \f$\mathbf{x}_{\mathrm{initial}} \f$ * in the deformed and undeformed (initial) state, respectively. The model assumes quasi-static conditions, * that is, the above momentum balance equation is solved under the assumption that the acceleration term @@ -43,7 +44,9 @@ \boldsymbol{\varepsilon} = \frac{1}{2} \left[ \nabla\mathbf{u} + (\nabla\mathbf{u})^{\mathrm{T}} \right]. \f] * - * PrimaryVariables are the displacements in each direction \f$\mathbf{u} \f$. + * Primary variables are the displacements in each direction \f$\mathbf{u} \f$. + * Gravity can be enabled or disabled via a runtime parameter. + * */ #ifndef DUMUX_GEOMECHANICS_ELASTIC_MODEL_HH #define DUMUX_GEOMECHANICS_ELASTIC_MODEL_HH diff --git a/dumux/geomechanics/poroelastic/model.hh b/dumux/geomechanics/poroelastic/model.hh index 852ccd6bd5..98d590e4f4 100644 --- a/dumux/geomechanics/poroelastic/model.hh +++ b/dumux/geomechanics/poroelastic/model.hh @@ -21,32 +21,39 @@ * \ingroup PoroElastic * \brief A poroelastic geomechanical model * - * The adapted geomechanical model describes the deformation of the solid matrix filled with fluids, using linear poro-elasticity: + * This model describes the deformation of a porous medium using the theory of linear poroelasticity. + * The momentum balance equation of a porous medium can be expressed by \f[ - \boldsymbol{\nabla}\cdot\boldsymbol{\sigma} + \mathbf{F} + \rho \mathbf{g} = \rho\ddot{\mathbf{u}} + \nabla\cdot\boldsymbol{\sigma_{\mathrm{eff}}} + \rho \mathbf{g} + \mathbf{f} = \rho\ddot{\mathbf{u}}, \f] - * where the density of the matrix \f$\rho = \phi \sum_{\beta} S_{\beta} \rho_{\beta} + (1-\phi) \rho_s\f$. + * where \f$\boldsymbol{\sigma_{\mathrm{eff}}} \f$ is the effective stress tensor, + * \f$\rho = (1 - \phi) \rho_s + \phi \rho_f \f$ is the average density of solids and fluids within the porous medium, + * \f$\mathbf{f} \f$ in \f$\mathrm{N/m^3} \f$ is the external force acting on the body per unit volume (e.g. magnetism), + * and \f$\mathbf{u} = \mathbf{x} - \mathbf{x}_{\mathrm{initial}} \f$ is the displacement, + * defined as the difference in material points \f$\mathbf{x} \f$ and \f$\mathbf{x}_{\mathrm{initial}} \f$ + * in the deformed and undeformed (initial) state, respectively. The model assumes quasi-static conditions, + * that is, the above momentum balance equation is solved under the assumption that the acceleration term + * \f$\rho\ddot{\mathbf{u}} \approx 0\f$. * - * We assume the quasi-static conditions, so the acceleration term \f$\rho\ddot{\mathbf{u}} \approx 0\f$. - * - * For isotropic materials, the stress tensor \f$\boldsymbol{\sigma} \f$ can be calculated after Hookes' Law, including the effective fluid pressure. + * Using the concept of the effective stress, the effective stress tensor \f$\boldsymbol{\sigma_{\mathrm{eff}}} \f$ is + * determined by the stress tensor \f$\boldsymbol{\sigma} \f$ , the effective pore pressure \f$p_{\mathrm{eff}} \f$ and the Biot's coefficient \f$\alpha \f$ : \f[ - {\boldsymbol {\sigma}}= \lambda \mathrm{tr}(\varepsilon) \boldsymbol{\mathrm I}+ 2G \varepsilon - \alpha p_{\mathrm{eff}} \boldsymbol{\mathrm I}, + \boldsymbol{\sigma_{\mathrm{eff}}} = \boldsymbol{\sigma} - \alpha p_{\mathrm{eff}} \mathbf{I} \f] - * with + * + * Per default, Hookes' Law is used for expressing the stress tensor \f$\boldsymbol{\sigma} \f$ as a function of the + * displacement: \f[ - {\boldsymbol {\varepsilon }}={\frac {1}{2}}\left[{\boldsymbol {\nabla }}\mathbf {u} +({\boldsymbol {\nabla }}\mathbf {u} )^{\mathrm {T} }\right], + \boldsymbol{\sigma} = \lambda\mathrm{tr}(\boldsymbol{\varepsilon}) \mathbf{I} + 2G \boldsymbol{\varepsilon}, \f] - * and + * with \f[ - p_{\mathrm{eff}} = \sum_{\beta} S_{\beta} p_{\beta}. + \boldsymbol{\varepsilon} = \frac{1}{2} \left[ \nabla\mathbf{u} + (\nabla\mathbf{u})^{\mathrm{T}} \right]. \f] - * Gravity can be enabled or disabled via a runtime parameter. * - * The equations are discretized using a vertex-centered finite volume (box) scheme as spatial discretization. The time discretization is not needed due to the quasi-static conditions. - * PrimaryVariables are the displacements in each direction \f$\mathbf{u} \f$. + * Primary variables are the displacements in each direction \f$\mathbf{u} \f$. + * Gravity can be enabled or disabled via a runtime parameter. * - * Pay attention to the sign convention, where tension has a positive! */ #ifndef DUMUX_GEOMECHANICS_POROELASTIC_MODEL_HH #define DUMUX_GEOMECHANICS_POROELASTIC_MODEL_HH -- GitLab