ProcessLib::TH2M::ConstitutiveRelations::PermeabilityModel< DisplacementDim > Struct Template Reference

Detailed Description

template<int DisplacementDim>
struct ProcessLib::TH2M::ConstitutiveRelations::PermeabilityModel< DisplacementDim >

Definition at line 25 of file PermeabilityModel.h.

#include <PermeabilityModel.h>

Public Member Functions

void

eval (SpaceTimeData const &x_t, MediaData const &media_data, SaturationData const &S_L_data, GasPressureData const &p_g, CapillaryPressureData const &p_cap, TemperatureData const &T_data, TransportPorosityData const &transport_poro_data, TotalStressData< DisplacementDim > const &total_stress_data, MechanicalStrainData< DisplacementDim > const &mechanical_strain_data, StrainData< DisplacementDim > const &eps_data, EquivalentPlasticStrainData const &equivalent_plastic_strain, PermeabilityData< DisplacementDim > &out) const

Member Function Documentation

eval()

template<int DisplacementDim>

void ProcessLib::TH2M::ConstitutiveRelations::PermeabilityModel< DisplacementDim >::eval	(	SpaceTimeData const &	x_t,
		MediaData const &	media_data,
		SaturationData const &	S_L_data,
		GasPressureData const &	p_g,
		CapillaryPressureData const &	p_cap,
		TemperatureData const &	T_data,
		TransportPorosityData const &	transport_poro_data,
		TotalStressData< DisplacementDim > const &	total_stress_data,
		MechanicalStrainData< DisplacementDim > const &	mechanical_strain_data,
		StrainData< DisplacementDim > const &	eps_data,
		EquivalentPlasticStrainData const &	equivalent_plastic_strain,
		PermeabilityData< DisplacementDim > &	out ) const

Definition at line 19 of file PermeabilityModel.cpp.

{
    namespace MPL = MaterialPropertyLib;
 
    static constexpr int kelvin_vector_size =
        MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
    using Invariants = MathLib::KelvinVector::Invariants<kelvin_vector_size>;
 
    auto const& medium = media_data.medium;
 
    MPL::VariableArray variables;
    variables.liquid_saturation = S_L_data.S_L;
    variables.temperature = T_data.T;
    variables.capillary_pressure = p_cap.pCap;
    variables.gas_phase_pressure = p_g.pG;
    variables.transport_porosity = transport_poro_data.phi;
 
    out.k_rel_G =
        medium
            .property(MPL::PropertyType::relative_permeability_nonwetting_phase)
            .template value<double>(variables, x_t.x, x_t.t, x_t.dt);
 
    out.dk_rel_G_dS_L =
        medium[MPL::PropertyType::relative_permeability_nonwetting_phase]
            .template dValue<double>(variables,
                                     MPL::Variable::liquid_saturation, x_t.x,
                                     x_t.t, x_t.dt);
 
    out.k_rel_L = medium.property(MPL::PropertyType::relative_permeability)
                      .template value<double>(variables, x_t.x, x_t.t, x_t.dt);
 
    out.dk_rel_L_dS_L = medium[MPL::PropertyType::relative_permeability]
                            .template dValue<double>(
                                variables, MPL::Variable::liquid_saturation,
                                x_t.x, x_t.t, x_t.dt);
 
    // For stress dependent permeability.
    variables.total_stress =
        MathLib::KelvinVector::kelvinVectorToSymmetricTensor(
            total_stress_data.sigma_total);
 
    variables.mechanical_strain.emplace<KelvinVector<DisplacementDim>>(
        mechanical_strain_data.eps_m);
 
    variables.equivalent_plastic_strain = equivalent_plastic_strain();
 
    variables.volumetric_strain = Invariants::trace(eps_data.eps);
 
    out.Ki = MPL::formEigenTensor<DisplacementDim>(
        medium.property(MPL::PropertyType::permeability)
            .value(variables, x_t.x, x_t.t, x_t.dt));
}

References MaterialPropertyLib::VariableArray::capillary_pressure, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::dk_rel_G_dS_L, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::dk_rel_L_dS_L, ProcessLib::ConstitutiveRelations::SpaceTimeData::dt, ProcessLib::ConstitutiveRelations::StrainData< DisplacementDim >::eps, ProcessLib::ConstitutiveRelations::MechanicalStrainData< DisplacementDim >::eps_m, MaterialPropertyLib::VariableArray::equivalent_plastic_strain, MaterialPropertyLib::formEigenTensor(), MaterialPropertyLib::VariableArray::gas_phase_pressure, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::k_rel_G, ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::k_rel_L, MathLib::KelvinVector::kelvin_vector_dimensions(), MathLib::KelvinVector::kelvinVectorToSymmetricTensor(), ProcessLib::TH2M::ConstitutiveRelations::PermeabilityData< DisplacementDim >::Ki, MaterialPropertyLib::liquid_saturation, MaterialPropertyLib::VariableArray::liquid_saturation, MaterialPropertyLib::VariableArray::mechanical_strain, ProcessLib::TH2M::ConstitutiveRelations::MediaData::medium, ProcessLib::TH2M::ConstitutiveRelations::CapillaryPressureData::pCap, MaterialPropertyLib::permeability, ProcessLib::TH2M::ConstitutiveRelations::GasPressureData::pG, ProcessLib::TH2M::ConstitutiveRelations::TransportPorosityData::phi, MaterialPropertyLib::relative_permeability, MaterialPropertyLib::relative_permeability_nonwetting_phase, ProcessLib::TH2M::ConstitutiveRelations::SaturationData::S_L, ProcessLib::TH2M::ConstitutiveRelations::TotalStressData< DisplacementDim >::sigma_total, ProcessLib::TH2M::ConstitutiveRelations::TemperatureData::T, ProcessLib::ConstitutiveRelations::SpaceTimeData::t, MaterialPropertyLib::VariableArray::temperature, MaterialPropertyLib::VariableArray::total_stress, MaterialPropertyLib::VariableArray::transport_porosity, MaterialPropertyLib::VariableArray::volumetric_strain, and ProcessLib::ConstitutiveRelations::SpaceTimeData::x.

---

The documentation for this struct was generated from the following files:

• ProcessLib/TH2M/ConstitutiveRelations/PermeabilityModel.h
• ProcessLib/TH2M/ConstitutiveRelations/PermeabilityModel.cpp