[GeoMechanicsApplication] Implement well constitutive behaviour for line thermal element

applications/GeoMechanicsApplication/custom_constitutive/thermal_filter_law.cpp

@@ -31,18 +31,14 @@ ConstitutiveLaw::Pointer GeoThermalFilterLaw::Clone() const
     return Kratos::make_shared<GeoThermalFilterLaw>(*this);
 }
 
-SizeType GeoThermalFilterLaw::WorkingSpaceDimension() { return mNumberOfDimensions; }
+SizeType GeoThermalFilterLaw::WorkingSpaceDimension() { return 1; }
 
-Matrix GeoThermalFilterLaw::CalculateThermalFilterMatrix(const Properties&  rProp,
-                                                                 const ProcessInfo& rProcessInfo) const
+Matrix GeoThermalFilterLaw::CalculateThermalFilterMatrix(const Properties& rProp, const ProcessInfo& rProcessInfo) const
 {
     KRATOS_TRY
 
-    Matrix result = ZeroMatrix(mNumberOfDimensions, mNumberOfDimensions);
-
-    const auto x = static_cast<int>(indexThermalFlux::X);
-    result(x, x) = rProp[THERMAL_CONDUCTIVITY_WATER];
-
+    Matrix result = ZeroMatrix(1, 1);
+    result(0, 0) = rProp[THERMAL_CONDUCTIVITY_WATER];
     return result;
 
     KRATOS_CATCH("")

applications/GeoMechanicsApplication/tests/test_thermal_element/test_thermal_filter_element/README.md

@@ -5,8 +5,7 @@
 **Source files:** [Thermal filter element with fixed temperature](https://github.com/KratosMultiphysics/Kratos/tree/master/applications/GeoMechanicsApplication/tests/test_thermal_element/test_thermal_filter_element)
 
 ## Case Specification
-In this thermal test case, a $2 \mathrm{[m]} \times 1 \mathrm{[m]}$ soil block is considered, with everywhere set to 0 $\mathrm{[^\circ C]}$. A well/filter passes vertically through middle of the domain. The temperature along the filter is kept constant at 0 $\mathrm{[^\circ C]}$. A sudden jump at the right and left boundary temperatures are given for the soil. These are both 25 $\mathrm{[^\circ C]}$. The simulation spans 100 days to allow for a transition from the initial to a linear temperature profile between the two sides. This test is conducted for various configurations, including 2D3N and 2D3N soil elements, and corresponding 2D2N and 2D3N filter elements. The temperature distribution is then evaluated with its own result.
-The boundary conditions are shown below:
+In this thermal test case, a $2 \mathrm{[m]} \times 1 \mathrm{[m]}$ soil block is considered, with everywhere set to 0 $\mathrm{[^\circ C]}$. A well/filter passes vertically through middle of the domain. The temperature along the filter is kept constant at 0 $\mathrm{[^\circ C]}$. A sudden jump at the right and left boundary temperatures are given for the soil. These are both 25 $\mathrm{[^\circ C]}$. The simulation spans 100 days to allow for a transition from the initial to a linear temperature profile between the two sides. This test is conducted for various configurations, including 2D3N and 2D6N soil elements, and corresponding 2D2N and 2D3N filter elements. The temperature distribution is then evaluated with its analytical result. The boundary conditions are shown below:
 
 <img src="../documentation_data/test_thermal_filter_element.svg" alt="Visualization of the Boundary conditions" title="Visualization of the Boundary conditions" width="600">