DiffusionVelocity.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 "MassMoleFractions.h"
#include "PhaseTransitionData.h"
#include "Porosity.h"
#include "ProcessLib/Reflection/ReflectionData.h"
#include "Saturation.h"
 
namespace ProcessLib::TH2M
{
namespace ConstitutiveRelations
{
template <int DisplacementDim>
struct DiffusionVelocityData
{
    GlobalDimVector<DisplacementDim> d_CG;
    GlobalDimVector<DisplacementDim> d_WG;
    GlobalDimVector<DisplacementDim> d_CL;
    GlobalDimVector<DisplacementDim> d_WL;
 
    static auto reflect()
    {
        using Self = DiffusionVelocityData<DisplacementDim>;
        namespace R = ProcessLib::Reflection;
 
        return std::tuple{
            R::makeReflectionData("diffusion_velocity_gas_gas", &Self::d_CG),
            R::makeReflectionData("diffusion_velocity_vapour_gas", &Self::d_WG),
            R::makeReflectionData("diffusion_velocity_solute_liquid",
                                  &Self::d_CL),
            R::makeReflectionData("diffusion_velocity_liquid_liquid",
                                  &Self::d_WL)};
    }
};
 
template <int DisplacementDim>
struct DiffusionVelocityModel
{
    void eval(
        CapillaryPressureGradientData<DisplacementDim> const& grad_p_cap,
        GasPressureGradientData<DisplacementDim> const& grad_p_GR,
        MassMoleFractionsData const& mass_mole_fractions_data,
        PhaseTransitionData const& phase_transition_data,
        PorosityData const& porosity_data,
        SaturationData const& S_L_data,
        TemperatureGradientData<DisplacementDim> const& grad_T,
        DiffusionVelocityData<DisplacementDim>& diffusion_velocity_data) const;
};
 
extern template struct DiffusionVelocityModel<2>;
extern template struct DiffusionVelocityModel<3>;
}  // namespace ConstitutiveRelations
}  // namespace ProcessLib::TH2M