TH2M/ConstitutiveRelations/TransportPorosity.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "TransportPorosity.h"
 
namespace ProcessLib::TH2M
{
namespace ConstitutiveRelations
{
 
template <int DisplacementDim>
void TransportPorosityModel<DisplacementDim>::eval(
    SpaceTimeData const& x_t, MediaData const& media_data,
    SaturationData const& S_L_data,
    PrevState<SaturationData> const& S_L_prev_data,
    CapillaryPressureData const& p_cap, GasPressureData const& p_GR,
    BishopsData const& chi_S_L, PrevState<BishopsData> const& chi_S_L_prev,
    SolidCompressibilityData const& solid_compressibility,
    MechanicalStrainData<DisplacementDim> const& eps_m_data,
    PrevState<MechanicalStrainData<DisplacementDim>> const& eps_m_prev_data,
    PrevState<TransportPorosityData> const& transport_porosity_prev_data,
    PorosityData const& poro_data,
    TransportPorosityData& transport_porosity_data) const
{
    static int const KelvinVectorSize =
        MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
    using Invariants = MathLib::KelvinVector::Invariants<KelvinVectorSize>;
 
    MaterialPropertyLib::VariableArray variables;
    MaterialPropertyLib::VariableArray variables_prev;
 
    variables.grain_compressibility = solid_compressibility();
 
    variables.liquid_saturation = S_L_data.S_L;
    variables_prev.liquid_saturation = S_L_prev_data->S_L;
 
    variables.effective_pore_pressure =
        (1 - chi_S_L.chi_S_L) * p_GR.pG +
        chi_S_L.chi_S_L * (p_GR.pG - p_cap.pCap);
 
    // Used in MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp
    // and MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp
    variables_prev.effective_pore_pressure =
        (1 - chi_S_L_prev->chi_S_L) * p_GR.pG_prev +
        chi_S_L_prev->chi_S_L * (p_GR.pG_prev - p_cap.pCap_prev);
 
    variables.volumetric_mechanical_strain =
        Invariants::trace(eps_m_data.eps_m);
    variables_prev.volumetric_mechanical_strain =
        Invariants::trace(eps_m_prev_data->eps_m);
 
    variables_prev.transport_porosity = transport_porosity_prev_data->phi;
    variables.porosity = poro_data.phi;
    assert(media_data.transport_porosity_prop != nullptr);
    transport_porosity_data.phi =
        media_data.transport_porosity_prop->template value<double>(
            variables, variables_prev, x_t.x, x_t.t, x_t.dt);
}
 
template struct TransportPorosityModel<2>;
template struct TransportPorosityModel<3>;
}  // namespace ConstitutiveRelations
}  // namespace ProcessLib::TH2M