TH2M/ConstitutiveRelations/TransportPorosity.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 "Base.h"
#include "Biot.h"
#include "Bishops.h"
#include "Porosity.h"
#include "ProcessLib/ConstitutiveRelations/MechanicalStrainData.h"
#include "ProcessLib/Reflection/ReflectionData.h"
#include "Saturation.h"
#include "SolidCompressibility.h"
#include "SolidThermalExpansion.h"
 
namespace ProcessLib::TH2M
{
namespace ConstitutiveRelations
{
struct TransportPorosityData
{
    double phi = nan;
 
    static auto reflect()
    {
        using Self = TransportPorosityData;
        namespace R = ProcessLib::Reflection;
 
        return std::tuple{
            R::makeReflectionData("transport_porosity", &Self::phi)};
    }
};
 
template <int DisplacementDim>
struct TransportPorosityModel
{
    void 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;
};
 
extern template struct TransportPorosityModel<2>;
extern template struct TransportPorosityModel<3>;
 
}  // namespace ConstitutiveRelations
}  // namespace ProcessLib::TH2M