TH2M/ConstitutiveRelations/PermeabilityModel.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "PermeabilityModel.h"
 
#include "MaterialLib/MPL/Utils/FormEigenTensor.h"
 
namespace ProcessLib::TH2M
{
namespace ConstitutiveRelations
{
template <int DisplacementDim>
void 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
{
    namespace MPL = MaterialPropertyLib;
 
    static constexpr int kelvin_vector_size =
        MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
    using Invariants = MathLib::KelvinVector::Invariants<kelvin_vector_size>;
 
    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 = media_data.relative_permeability_nonwetting_phase_prop
                      .template value<double>(variables, x_t.x, x_t.t, x_t.dt);
 
    out.dk_rel_G_dS_L = media_data.relative_permeability_nonwetting_phase_prop
                            .template dValue<double>(
                                variables, MPL::Variable::liquid_saturation,
                                x_t.x, x_t.t, x_t.dt);
 
    out.k_rel_L = media_data.relative_permeability_prop.template value<double>(
        variables, x_t.x, x_t.t, x_t.dt);
 
    out.dk_rel_L_dS_L =
        media_data.relative_permeability_prop.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>(
        media_data.permeability_prop.value(variables, x_t.x, x_t.t, x_t.dt));
}
 
template struct PermeabilityModel<2>;
template struct PermeabilityModel<3>;
}  // namespace ConstitutiveRelations
}  // namespace ProcessLib::TH2M