TH2M/ConstitutiveRelations/ConstitutiveData.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "Advection.h"
#include "Biot.h"
#include "Bishops.h"
#include "CEquation.h"
#include "ConstitutiveDensity.h"
#include "DarcyVelocity.h"
#include "DiffusionVelocity.h"
#include "ElasticTangentStiffnessData.h"
#include "Enthalpy.h"
#include "EquivalentPlasticStrainData.h"
#include "FluidDensity.h"
#include "Gravity.h"
#include "InternalEnergy.h"
#include "MassMoleFractions.h"
#include "MechanicalStrainModel.h"
#include "PermeabilityData.h"
#include "PhaseTransitionData.h"
#include "Porosity.h"
#include "ProcessLib/ConstitutiveRelations/EffectiveStressData.h"
#include "ProcessLib/ConstitutiveRelations/StrainData.h"
#include "ProcessLib/ConstitutiveRelations/StressData.h"
#include "ProcessLib/Reflection/ReflectionData.h"
#include "PureLiquidDensity.h"
#include "Saturation.h"
#include "SolidCompressibility.h"
#include "SolidDensity.h"
#include "SolidHeatCapacity.h"
#include "SolidMechanics.h"
#include "SolidThermalExpansion.h"
#include "Swelling.h"
#include "TEquation.h"
#include "ThermalConductivity.h"
#include "TotalStress.h"
#include "TransportPorosity.h"
#include "UEquation.h"
#include "VapourPartialPressure.h"
#include "Viscosity.h"
#include "WEquation.h"
 
namespace ProcessLib::TH2M
{
namespace ConstitutiveRelations
{
template <int DisplacementDim>
struct StatefulData
{
    SaturationData S_L_data;
    SwellingDataStateful<DisplacementDim> swelling_data;
    BishopsData chi_S_L;
    ProcessLib::ConstitutiveRelations::EffectiveStressData<DisplacementDim>
        eff_stress_data;
    MechanicalStrainData<DisplacementDim> mechanical_strain_data;
    PureLiquidDensityData rho_W_LR;
    ConstituentDensityData constituent_density_data;
    InternalEnergyData internal_energy_data;
    PorosityData porosity_data;
    TransportPorosityData transport_porosity_data;
 
    static auto reflect()
    {
        using Self = StatefulData<DisplacementDim>;
 
        return Reflection::reflectWithoutName(&Self::S_L_data,
                                              &Self::swelling_data,
                                              &Self::chi_S_L,
                                              &Self::eff_stress_data,
                                              &Self::porosity_data,
                                              &Self::transport_porosity_data);
    }
};
 
template <int DisplacementDim>
struct StatefulDataPrev
{
    PrevState<SaturationData> S_L_data;
    PrevState<SwellingDataStateful<DisplacementDim>> swelling_data;
    PrevState<BishopsData> chi_S_L;
    PrevState<
        ProcessLib::ConstitutiveRelations::EffectiveStressData<DisplacementDim>>
        eff_stress_data;
    PrevState<MechanicalStrainData<DisplacementDim>> mechanical_strain_data;
    PrevState<PureLiquidDensityData> rho_W_LR;
    PrevState<ConstituentDensityData> constituent_density_data;
    PrevState<InternalEnergyData> internal_energy_data;
    PrevState<PorosityData> porosity_data;
    PrevState<TransportPorosityData> transport_porosity_data;
 
    StatefulDataPrev<DisplacementDim>& operator=(
        StatefulData<DisplacementDim> const& state)
    {
        S_L_data = state.S_L_data;
        swelling_data = state.swelling_data;
        chi_S_L = state.chi_S_L;
        eff_stress_data = state.eff_stress_data;
        mechanical_strain_data = state.mechanical_strain_data;
        rho_W_LR = state.rho_W_LR;
        constituent_density_data = state.constituent_density_data;
        internal_energy_data = state.internal_energy_data;
        porosity_data = state.porosity_data;
        transport_porosity_data = state.transport_porosity_data;
 
        return *this;
    }
};

template <int DisplacementDim>
struct OutputData
{
    PermeabilityData<DisplacementDim> permeability_data;
    FluidEnthalpyData fluid_enthalpy_data;
    SolidEnthalpyData solid_enthalpy_data;
    MassMoleFractionsData mass_mole_fractions_data;
    FluidDensityData fluid_density_data;
    VapourPartialPressureData vapour_pressure_data;
    SolidDensityData solid_density_data;
    DiffusionVelocityData<DisplacementDim> diffusion_velocity_data;
    DarcyVelocityData<DisplacementDim> darcy_velocity_data;
    ProcessLib::ConstitutiveRelations::StrainData<DisplacementDim> eps_data;
    TotalStressData<DisplacementDim> total_stress_data;
 
    static auto reflect()
    {
        using Self = OutputData<DisplacementDim>;
 
        return Reflection::reflectWithoutName(&Self::permeability_data,
                                              &Self::fluid_enthalpy_data,
                                              &Self::solid_enthalpy_data,
                                              &Self::mass_mole_fractions_data,
                                              &Self::fluid_density_data,
                                              &Self::vapour_pressure_data,
                                              &Self::solid_density_data,
                                              &Self::diffusion_velocity_data,
                                              &Self::darcy_velocity_data,
                                              &Self::eps_data,
                                              &Self::total_stress_data);
    }
};

template <int DisplacementDim>
struct ConstitutiveData
{
    SolidMechanicsDataStateless<DisplacementDim> s_mech_data;
};

template <int DisplacementDim>
struct ConstitutiveTempData
{
    SwellingDataStateless<DisplacementDim> swelling_data;
    ElasticTangentStiffnessData<DisplacementDim> C_el_data;
    BiotData biot_data;
    SolidCompressibilityData beta_p_SR;
    SolidThermalExpansionData<DisplacementDim> s_therm_exp_data;
    EquivalentPlasticStrainData equivalent_plastic_strain_data;
    ViscosityData viscosity_data;
    PhaseTransitionData phase_transition_data;
    SolidHeatCapacityData solid_heat_capacity_data;
    ThermalConductivityData<DisplacementDim> thermal_conductivity_data;
    EffectiveVolumetricEnthalpy effective_volumetric_enthalpy_data;
    AdvectionData<DisplacementDim> advection_data;
    VolumetricBodyForce<DisplacementDim> volumetric_body_force;
    FC1Data<DisplacementDim> fC_1;
    FC2aData fC_2a;
    FC3aData fC_3a;
    FC4LCpGData<DisplacementDim> fC_4_LCpG;
    FC4LCpCData<DisplacementDim> fC_4_LCpC;
    FC4LCTData<DisplacementDim> fC_4_LCT;
    FC4MCpGData fC_4_MCpG;
    FC4MCpCData fC_4_MCpC;
    FC4MCTData fC_4_MCT;
    FC4MCuData fC_4_MCu;
 
    FW1Data<DisplacementDim> fW_1;
    FW2Data fW_2;
    FW3aData fW_3a;
    FW4LWpGData<DisplacementDim> fW_4_LWpG;
    FW4LWpCData<DisplacementDim> fW_4_LWpC;
    FW4LWTData<DisplacementDim> fW_4_LWT;
    FW4MWpGData fW_4_MWpG;
    FW4MWpCData fW_4_MWpC;
    FW4MWTData fW_4_MWT;
    FW4MWuData fW_4_MWu;
 
    FT1Data fT_1;
    FT2Data<DisplacementDim> fT_2;
    FT3Data<DisplacementDim> fT_3;
 
    FU2KUpCData fu_2_KupC;
};

template <int DisplacementDim>
struct DerivativesData
{
    SaturationDataDeriv dS_L_dp_cap;
    AdvectionDerivativeData<DisplacementDim> advection_d_data;
    PorosityDerivativeData porosity_d_data;
    ThermalConductivityDerivativeData<DisplacementDim>
        thermal_conductivity_d_data;
    SolidDensityDerivativeData solid_density_d_data;
    EffectiveVolumetricInternalEnergyDerivatives
        effective_volumetric_internal_energy_d_data;
    EffectiveVolumetricEnthalpyDerivatives effective_volumetric_enthalpy_d_data;
    FC2aDerivativeData dfC_2a;
    FC3aDerivativeData dfC_3a;
    FC4LCpGDerivativeData<DisplacementDim> dfC_4_LCpG;
    FC4LCpCDerivativeData<DisplacementDim> dfC_4_LCpC;
    FC4MCpGDerivativeData dfC_4_MCpG;
    FC4MCTDerivativeData dfC_4_MCT;
    FC4MCuDerivativeData dfC_4_MCu;
    FW2DerivativeData dfW_2;
    FW3aDerivativeData dfW_3a;
    FW4LWpGDerivativeData<DisplacementDim> dfW_4_LWpG;
    FW4LWpCDerivativeData<DisplacementDim> dfW_4_LWpC;
    FT1DerivativeData dfT_1;
    FT2DerivativeData<DisplacementDim> dfT_2;
    FU1KUTDerivativeData<DisplacementDim> dfu_1_KuT;
    FU2KUpCDerivativeData dfu_2_KupC;
};
 
}  // namespace ConstitutiveRelations
}  // namespace ProcessLib::TH2M