[GeoMechanicsApplication] Implement well constitutive behaviour for line pressure element

applications/GeoMechanicsApplication/custom_elements/README.md

@@ -175,5 +175,54 @@ where
 The superscripts $^l$ and $^r$ for Robin boundary condition indicate the left hands side (matrix) and right hand side (vector), respectively. The superscripts $^e$ and $^{ep}$ for $\Omega$ and $\Gamma$ indicate values in the element volume and perpendicular to element boundaries, respectively.  
 
 
+# Pressure Filter Line Element
+For laminar flow along the axis of the well and uniformly distributed storage over the length of the well bore (Diersch, 2014)
+
+```math
+\pi R^2 \left( \frac{1}{l_w} + \rho_0 g \beta \right) \frac{\partial h}{\partial t} - \pi R^2 K_w \frac{\partial}{\partial y} \left[ f_{\mu} \left( \frac{\partial h}{\partial y} + \chi e \right)\right] = -Q_w \delta \left( y - y_w \right)
+```
+
+where
+
+```math
+K_w = \frac{R^2 \rho_0 g}{8 \mu_0} \quad \quad \quad \quad f_u = \frac{\mu_0}{\mu} \quad \quad \quad \quad \chi = \frac{\rho - \rho_0}{\rho_0} \quad \quad \quad \quad h = \frac{p}{\rho_0 g} + y
+```
+
+- $Q_w$			= pumping rate sink $\mathrm{\left[ m^3/s \right]}$
+- $t$				= time $\mathrm{\left[ s \right]}$
+- $y$				= vertical coordinate $\mathrm{\left[ m \right]}$
+- $y_w$				= location of the discharge point $\mathrm{\left[ m \right]}$
+- $h$				= hydraulic head in the well $\mathrm{\left[ m \right]}$
+- $l_w$				= total length of liquid filled well bore $\mathrm{\left[ m \right]}$
+- $R$				= radius of the well casing $\mathrm{\left[ m \right]}$
+- $K_w$				= Hagen-Poiseuille permeability $\mathrm{\left[ m \right]}$
+- $\delta$			= Dirichlet delta function $\mathrm{\left[ - \right]}$
+- $\beta$			= compressibility of the liquid $\mathrm{\left[ m^2/N \right]}$
+- $f_{\mu}$			= viscosity relation function of liquid $\mathrm{\left[ - \right]}$
+- $\chi$			= buoyancy coefficient $\mathrm{\left[ - \right]}$
+- $e$				= gravitational unit vector $\mathrm{\left[ - \right]}$
+- $g$				= gravitational acceleration $\mathrm{\left[ m/s^2 \right]}$
+- $\rho_0$			= reference density of the fluid $\mathrm{\left[ kg/m^3 \right]}$
+- $\mu_0$			= reference viscosity of the fluid $\mathrm{\left[ Pas \right]}$
+- $p$				= liquid pressure $\mathrm{\left[ Pa \right]}$
+
+reformulated in pressure and preserving mass
+
+```math
+\rho^w \left( \frac{1}{\rho^w g l_w} + \beta \right) \frac{\partial p}{\partial t} - \rho^w \frac{\partial}{\partial y} \left[ \frac{R^2}{8 \mu} \left( \frac{\partial p}{\partial y} - \rho^w g \right)\right] = -\frac{\rho^w Q_w}{\pi R^2}
+```
+
+The one dimensionalpressure equation with $\alpha = 1$ and $n = 1$ captures the well flow equation.
+
+```math
+\rho^w \beta^* \frac{\partial p}{\partial t} - \frac{\partial}{\partial y} \left[ \frac{\rho^w K^*}{\mu} \left( \frac{\partial p}{\partial y} - \rho^w g \right) \right] = - \rho^w Q^*_w
+```
+
+where the intrinsic permeability and boundary condition are given by
+
+```math
+\beta^* = \frac{1}{\rho^w g l_w} + \beta \quad \quad \quad \quad K^* = \frac{R^2}{8} \quad \quad \quad \quad Q^*_w = \frac{Q_w}{\pi R^2}
+```
+
 ## Bibliography
 Diersch, H.-J. G., 2014. FEFLOW; Finite Element Modeling of Flow, Mass and Heat Transport in Porous and Fractured Media. Springer.

applications/GeoMechanicsApplication/custom_elements/filter_compressibility_calculator.cpp

@@ -68,10 +68,8 @@ Matrix FilterCompressibilityCalculator::CalculateCompressibilityMatrix() const
 double FilterCompressibilityCalculator::CalculateElasticCapacity(double ProjectedGravity) const
 {
     const auto&  r_properties = mInputProvider.GetElementProperties();
-    const double result =
-        1.0 / (r_properties[DENSITY_WATER] * ProjectedGravity * r_properties[FILTER_LENGTH]) +
-        1.0 / r_properties[BULK_MODULUS_FLUID];
-    return result;
+    return 1.0 / (r_properties[DENSITY_WATER] * ProjectedGravity * r_properties[FILTER_LENGTH]) +
+           1.0 / r_properties[BULK_MODULUS_FLUID];
 }
 
 const Properties& FilterCompressibilityCalculator::InputProvider::GetElementProperties() const

applications/GeoMechanicsApplication/custom_elements/transient_Pw_line_element.h

@@ -281,9 +281,8 @@ class KRATOS_API(GEO_MECHANICS_APPLICATION) TransientPwLineElement : public Elem
     {
         const auto number_integration_points = GetGeometry().IntegrationPointsNumber(GetIntegrationMethod());
         GeometryType::JacobiansType J_container{number_integration_points};
-        for (std::size_t i = 0; i < number_integration_points; ++i) {
-            J_container[i].resize(GetGeometry().WorkingSpaceDimension(),
-                                  GetGeometry().LocalSpaceDimension(), false);
+        for (auto& container : J_container) {
+            container.resize(GetGeometry().WorkingSpaceDimension(), GetGeometry().LocalSpaceDimension(), false);
         }
         GetGeometry().Jacobian(J_container, this->GetIntegrationMethod());
 

applications/GeoMechanicsApplication/tests/test_pressure_line_element/test_oblique_filter_element2D3N/MaterialParameters.json

@@ -7,14 +7,14 @@
         "name"                             : "GeoLinearElasticPlaneStrain2DLaw"
       },
       "Variables" : {
-        "IGNORE_UNDRAINED"                 : false,
         "DENSITY_SOLID"                    : 2.650000e+03,
         "DENSITY_WATER"                    : 1.000000e+03,
         "POROSITY"                         : 1.000000e-01,
         "BULK_MODULUS_SOLID"               : 9.000000e+19,
         "BULK_MODULUS_FLUID"               : 1.000000e+20,
         "PERMEABILITY_XX"                  : 9.084000e-06,
         "DYNAMIC_VISCOSITY"                : 1.000000e-03,
+        "BIOT_COEFFICIENT"                 : 1.000000e+00,
         "RETENTION_LAW"                    : "PressureFilterLaw",
         "CROSS_AREA"                       : 1.0,
         "FILTER_LENGTH"                    : 3.0

applications/GeoMechanicsApplication/tests/test_pressure_line_element/test_oblique_filter_element3D3N/MaterialParameters.json

@@ -7,14 +7,14 @@
         "name"                             : "GeoIncrementalLinearElastic3DLaw"
       },
       "Variables" : {
-        "IGNORE_UNDRAINED"                 : false,
         "DENSITY_SOLID"                    : 2.650000e+03,
         "DENSITY_WATER"                    : 1.000000e+03,
         "POROSITY"                         : 1.000000e-01,
         "BULK_MODULUS_SOLID"               : 9.000000e+19,
         "BULK_MODULUS_FLUID"               : 1.000000e+20,
         "PERMEABILITY_XX"                  : 9.084000e-06,
         "DYNAMIC_VISCOSITY"                : 1.000000e-03,
+        "BIOT_COEFFICIENT"                 : 1.000000e+00,
         "RETENTION_LAW"                    : "PressureFilterLaw",
         "CROSS_AREA"                       : 1.0,
         "FILTER_LENGTH"              : 3.0