[PL] Constitutive data provide reflection meta data

ProcessLib/ThermoRichardsMechanics/Constitutive/Base.h

@@ -13,6 +13,7 @@
 #include "MaterialLib/MPL/Medium.h"
 #include "MathLib/KelvinVector.h"
 #include "ParameterLib/SpatialPosition.h"
+#include "ProcessLib/Reflection/ReflectionData.h"
 
 namespace ProcessLib::ThermoRichardsMechanics
 {
@@ -116,6 +117,13 @@ struct StrainData
 {
     // TODO Move initialization to the local assembler.
     KelvinVector<DisplacementDim> eps = KVzero<DisplacementDim>();
+
+    static auto reflect()
+    {
+        using Self = StrainData<DisplacementDim>;
+
+        return ProcessLib::Reflection::reflectWithName("epsilon", &Self::eps);
+    }
 };
 
 }  // namespace ProcessLib::ThermoRichardsMechanics

ProcessLib/ThermoRichardsMechanics/Constitutive/DarcyLaw.h

@@ -11,7 +11,6 @@
 #pragma once
 
 #include "LiquidDensity.h"
-#include "LiquidViscosity.h"
 #include "Permeability.h"
 #include "ThermoOsmosis.h"
 
@@ -21,6 +20,14 @@ template <int DisplacementDim>
 struct DarcyLawData
 {
     Eigen::Vector<double, DisplacementDim> v_darcy;
+
+    static auto reflect()
+    {
+        using Self = DarcyLawData<DisplacementDim>;
+
+        return ProcessLib::Reflection::reflectWithName("velocity",
+                                                       &Self::v_darcy);
+    }
 };
 
 template <int DisplacementDim>

ProcessLib/ThermoRichardsMechanics/Constitutive/LiquidDensity.h

@@ -19,6 +19,12 @@ struct LiquidDensityData
     double rho_LR;
     double drho_LR_dp;
     double drho_LR_dT;
+
+    static auto reflect()
+    {
+        return ProcessLib::Reflection::reflectWithName(
+            "liquid_density", &LiquidDensityData::rho_LR);
+    }
 };
 
 template <int DisplacementDim>

ProcessLib/ThermoRichardsMechanics/Constitutive/LiquidViscosity.h

@@ -10,13 +10,19 @@
 
 #pragma once
 
-#include "Porosity.h"
+#include "ProcessLib/ThermoRichardsMechanics/Constitutive/Base.h"
 
 namespace ProcessLib::ThermoRichardsMechanics
 {
 struct LiquidViscosityData
 {
     double viscosity;
+
+    static auto reflect()
+    {
+        return ProcessLib::Reflection::reflectWithName(
+            "viscosity", &LiquidViscosityData::viscosity);
+    }
 };
 
 template <int DisplacementDim>

ProcessLib/ThermoRichardsMechanics/Constitutive/Porosity.h

@@ -11,7 +11,6 @@
 #pragma once
 
 #include "Bishops.h"
-#include "MathLib/KelvinVector.h"
 #include "Saturation.h"
 #include "SolidCompressibility.h"
 
@@ -20,6 +19,12 @@ namespace ProcessLib::ThermoRichardsMechanics
 struct PorosityData
 {
     double phi;
+
+    static auto reflect()
+    {
+        return ProcessLib::Reflection::reflectWithName("porosity",
+                                                       &PorosityData::phi);
+    }
 };
 
 struct TransportPorosityData
@@ -28,8 +33,8 @@ struct TransportPorosityData
 
     static auto reflect()
     {
-        return std::make_tuple(
-            ReflectionData{"transport_porosity", &TransportPorosityData::phi});
+        return ProcessLib::Reflection::reflectWithName(
+            "transport_porosity", &TransportPorosityData::phi);
     }
 };
 

ProcessLib/ThermoRichardsMechanics/Constitutive/Saturation.h

@@ -17,6 +17,12 @@ namespace ProcessLib::ThermoRichardsMechanics
 struct SaturationData
 {
     double S_L;
+
+    static auto reflect()
+    {
+        return ProcessLib::Reflection::reflectWithName("saturation",
+                                                       &SaturationData::S_L);
+    }
 };
 
 struct SaturationDataDeriv

ProcessLib/ThermoRichardsMechanics/Constitutive/SolidDensity.h

@@ -19,6 +19,12 @@ struct SolidDensityData
 {
     double rho_SR;
     double dry_density_solid;
+
+    static auto reflect()
+    {
+        return ProcessLib::Reflection::reflectWithName(
+            "dry_density_solid", &SolidDensityData::dry_density_solid);
+    }
 };
 
 template <int DisplacementDim>

ProcessLib/ThermoRichardsMechanics/Constitutive/SolidMechanics.h

@@ -25,6 +25,15 @@ struct SolidMechanicsDataStateful
     // TODO it seems fragile that some data have to be initialized that way.
     KelvinVector<DisplacementDim> sigma_eff = KVzero<DisplacementDim>();
     KelvinVector<DisplacementDim> eps_m = KVzero<DisplacementDim>();
+
+    static auto reflect()
+    {
+        using Self = SolidMechanicsDataStateful<DisplacementDim>;
+
+        // TODO add eps_m?
+        return ProcessLib::Reflection::reflectWithName("sigma",
+                                                       &Self::sigma_eff);
+    }
 };
 
 template <int DisplacementDim>

ProcessLib/ThermoRichardsMechanics/Constitutive/Swelling.h

@@ -20,6 +20,14 @@ template <int DisplacementDim>
 struct SwellingDataStateful
 {
     KelvinVector<DisplacementDim> sigma_sw = KVzero<DisplacementDim>();
+
+    static auto reflect()
+    {
+        using Self = SwellingDataStateful<DisplacementDim>;
+
+        return ProcessLib::Reflection::reflectWithName("swelling_stress",
+                                                       &Self::sigma_sw);
+    }
 };
 
 template <int DisplacementDim>

ProcessLib/ThermoRichardsMechanics/ConstitutiveSetting.h

@@ -30,6 +30,7 @@
 
 namespace ProcessLib::ThermoRichardsMechanics
 {
+
 /// Data whose state must be tracked by the TRM process.
 template <int DisplacementDim>
 struct StatefulData
@@ -40,6 +41,19 @@ struct StatefulData
     StrainData<DisplacementDim> eps_data;
     SwellingDataStateful<DisplacementDim> swelling_data;
     SolidMechanicsDataStateful<DisplacementDim> s_mech_data;
+
+    static auto reflect()
+    {
+        using Self = StatefulData<DisplacementDim>;
+
+        return ProcessLib::Reflection::reflectWithoutName(
+            &Self::S_L_data,
+            &Self::poro_data,
+            &Self::transport_poro_data,
+            &Self::eps_data,
+            &Self::swelling_data,
+            &Self::s_mech_data);
+    }
 };
 
 /// Data that is needed for output purposes, but not directly for the assembly.
@@ -50,6 +64,16 @@ struct OutputData
     LiquidDensityData rho_L_data;
     LiquidViscosityData mu_L_data;
     SolidDensityData rho_S_data;
+
+    static auto reflect()
+    {
+        using Self = OutputData<DisplacementDim>;
+
+        return ProcessLib::Reflection::reflectWithoutName(&Self::darcy_data,
+                                                          &Self::rho_L_data,
+                                                          &Self::mu_L_data,
+                                                          &Self::rho_S_data);
+    }
 };
 
 /// Data that is needed for the equation system assembly.