[GeoMechanicsApplication] Implement well constitutive behaviour for line thermal element

applications/GeoMechanicsApplication/custom_constitutive/thermal_filter_law.cpp

@@ -0,0 +1,45 @@
+// KRATOS___
+//     //   ) )
+//    //         ___      ___
+//   //  ____  //___) ) //   ) )
+//  //    / / //       //   / /
+// ((____/ / ((____   ((___/ /  MECHANICS
+//
+//  License:         geo_mechanics_application/license.txt
+//
+//  Main authors:    Mohamed Nabi
+//                   John van Esch
+//
+
+#include "custom_constitutive/thermal_filter_law.h"
+#include "geo_mechanics_application_variables.h"
+
+namespace Kratos
+{
+
+GeoThermalFilterLaw::GeoThermalFilterLaw() : mNumberOfDimensions{1} {}
+
+GeoThermalFilterLaw::GeoThermalFilterLaw(SizeType NumberOfDimensions)
+    : mNumberOfDimensions{NumberOfDimensions}
+{
+    KRATOS_ERROR_IF(mNumberOfDimensions != 1)
+        << "Got invalid number of dimensions. The dimension has to be 1, but got: " << mNumberOfDimensions << std::endl;
+}
+
+ConstitutiveLaw::Pointer GeoThermalFilterLaw::Clone() const
+{
+    return Kratos::make_shared<GeoThermalFilterLaw>(*this);
+}
+
+SizeType GeoThermalFilterLaw::WorkingSpaceDimension() { return mNumberOfDimensions; }
+
+Matrix GeoThermalFilterLaw::CalculateThermalFilterMatrix(const Properties& rProp, const ProcessInfo& rProcessInfo) const
+{
+    KRATOS_TRY
+
+    return ScalarMatrix(1, 1, rProp[THERMAL_CONDUCTIVITY_WATER]);
+
+    KRATOS_CATCH("")
+}
+
+} // Namespace Kratos

applications/GeoMechanicsApplication/custom_constitutive/thermal_filter_law.h

@@ -0,0 +1,59 @@
+// KRATOS___
+//     //   ) )
+//    //         ___      ___
+//   //  ____  //___) ) //   ) )
+//  //    / / //       //   / /
+// ((____/ / ((____   ((___/ /  MECHANICS
+//
+//  License:         geo_mechanics_application/license.txt
+//
+//  Main authors:    Mohamed Nabi
+//                   John van Esch
+//
+
+#pragma once
+
+#include "includes/constitutive_law.h"
+
+namespace Kratos
+{
+
+/**
+ * @class GeoThermalFilterLaw
+ * @ingroup GeoMechanicsApplication
+ * @brief This class defines the thermal filter for heat cases
+ */
+class KRATOS_API(GEO_MECHANICS_APPLICATION) GeoThermalFilterLaw : public ConstitutiveLaw
+{
+public:
+    /// Counted pointer of GeoThermalFilterLaw
+    KRATOS_CLASS_POINTER_DEFINITION(GeoThermalFilterLaw);
+
+    GeoThermalFilterLaw();
+
+    explicit GeoThermalFilterLaw(SizeType NumberOfDimensions);
+
+    ConstitutiveLaw::Pointer Clone() const override;
+
+    SizeType WorkingSpaceDimension() override;
+
+    Matrix CalculateThermalFilterMatrix(const Properties& rProp, const ProcessInfo& rProcessInfo) const;
+
+private:
+    SizeType mNumberOfDimensions;
+
+    friend class Serializer;
+
+    void save(Serializer& rSerializer) const override
+    {
+        KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, ConstitutiveLaw)
+        rSerializer.save("NumberOfDimensions", mNumberOfDimensions);
+    }
+
+    void load(Serializer& rSerializer) override
+    {
+        KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, ConstitutiveLaw)
+        rSerializer.load("NumberOfDimensions", mNumberOfDimensions);
+    }
+};
+} // namespace Kratos.

applications/GeoMechanicsApplication/geo_mechanics_application.cpp

@@ -354,6 +354,8 @@ void KratosGeoMechanicsApplication::Register() {
 
     KRATOS_REGISTER_CONSTITUTIVE_LAW("GeoThermalDispersion2DLaw",     mGeoThermalDispersion2DLaw)
     KRATOS_REGISTER_CONSTITUTIVE_LAW("GeoThermalDispersion3DLaw",     mGeoThermalDispersion3DLaw)
+    KRATOS_REGISTER_CONSTITUTIVE_LAW("GeoThermalFilter2DLaw",         mGeoThermalFilter2DLaw)
+    KRATOS_REGISTER_CONSTITUTIVE_LAW("GeoThermalFilter3DLaw",         mGeoThermalFilter3DLaw)
 
     //Register Variables
     KRATOS_REGISTER_VARIABLE( VELOCITY_COEFFICIENT )

applications/GeoMechanicsApplication/geo_mechanics_application.h

@@ -132,6 +132,7 @@
 #include "custom_constitutive/small_strain_umat_3D_interface_law.hpp"
 #include "custom_constitutive/small_strain_umat_3D_law.hpp"
 #include "custom_constitutive/thermal_dispersion_law.h"
+#include "custom_constitutive/thermal_filter_law.h"
 
 namespace Kratos {
 
@@ -601,6 +602,9 @@ class KRATOS_API(GEO_MECHANICS_APPLICATION) KratosGeoMechanicsApplication : publ
     const GeoThermalDispersionLaw mGeoThermalDispersion2DLaw{ConstitutiveLaw::SizeType(2)};
     const GeoThermalDispersionLaw mGeoThermalDispersion3DLaw{ConstitutiveLaw::SizeType(3)};
 
+    const GeoThermalFilterLaw mGeoThermalFilter2DLaw{std::size_t(1)};
+    const GeoThermalFilterLaw mGeoThermalFilter3DLaw{std::size_t(1)};
+
     ///@}
 
 }; // Class KratosGeoMechanicsApplication

applications/GeoMechanicsApplication/tests/cpp_tests/test_thermal_filter_law.cpp

@@ -0,0 +1,51 @@
+// KRATOS___
+//     //   ) )
+//    //         ___      ___
+//   //  ____  //___) ) //   ) )
+//  //    / / //       //   / /
+// ((____/ / ((____   ((___/ /  MECHANICS
+//
+//  License:         geo_mechanics_application/license.txt
+//
+//  Main authors:    Mohamed Nabi
+//
+
+#include "custom_constitutive/thermal_filter_law.h"
+#include "geo_mechanics_application.h"
+#include "includes/ublas_interface.h"
+#include "testing/testing.h"
+
+namespace Kratos::Testing {
+
+KRATOS_TEST_CASE_IN_SUITE(CalculateThermalFilterLawMatrix, KratosGeoMechanicsFastSuite)
+{
+    Model current_model;
+    auto& r_model_part = current_model.CreateModelPart("ModelPart");
+
+    auto cond_prop = r_model_part.CreateNewProperties(0);
+    cond_prop->SetValue(THERMAL_CONDUCTIVITY_WATER, 1000.0);
+
+    const SizeType dimension = 1;
+    GeoThermalFilterLaw geo_thermal_filter_law(dimension);
+    ProcessInfo info;
+
+    const Matrix thermal_filter_matrix = geo_thermal_filter_law.CalculateThermalFilterMatrix(*cond_prop, info);
+
+    Matrix expected_solution = ZeroMatrix(1, 1);
+    expected_solution(0, 0)  = 1000.0;
+
+    constexpr double tolerance{1.0e-6};
+
+     KRATOS_EXPECT_NEAR(thermal_filter_matrix(0, 0), expected_solution(0, 0), tolerance);
+
+}
+
+KRATOS_TEST_CASE_IN_SUITE(GetWorkingSpaceDimension_ReturnsFilterDimensionValue, KratosGeoMechanicsFastSuite)
+{
+    constexpr SizeType dimension = 1;
+    GeoThermalDispersionLaw geo_thermal_filter_law(dimension);
+
+    KRATOS_EXPECT_EQ(geo_thermal_filter_law.WorkingSpaceDimension(), dimension);
+}
+
+} // namespace Kratos::Testing