ThermoOsmoticCoefficient.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 "BaseLib/Error.h"
#include "MaterialLib/MPL/Medium.h"
#include "MaterialLib/MPL/Utils/FormEigenTensor.h"
#include "ParameterLib/SpatialPosition.h"
 
namespace ProcessLib
{
template <int GlobalDim>
Eigen::Matrix<double, GlobalDim, GlobalDim> getThermoOsmoticCoefficient(
    MaterialPropertyLib::Medium const& medium,
    MaterialPropertyLib::VariableArray const& variable_array,
    ParameterLib::SpatialPosition const& pos,
    double const t,
    double const dt,
    Eigen::Matrix<double, GlobalDim, GlobalDim> const& intrinsic_permeability,
    double const liquid_dynamic_viscosity)
{
    if (medium.hasProperty(
            MaterialPropertyLib::PropertyType::thermal_osmosis_permeability))
    {
        if (liquid_dynamic_viscosity <= 0.)
        {
            OGS_FATAL(
                "Liquid dynamic viscosity must be > 0 when using "
                "thermal_osmosis_permeability, but is {:g} in {:s}.",
                liquid_dynamic_viscosity, medium.description());
        }
 
        auto const epsilon_T = medium
                                   .property(MaterialPropertyLib::PropertyType::
                                                 thermal_osmosis_permeability)
                                   .value<double>(variable_array, pos, t, dt);
        return epsilon_T * intrinsic_permeability / liquid_dynamic_viscosity;
    }
 
    if (medium.hasProperty(
            MaterialPropertyLib::PropertyType::thermal_osmosis_coefficient))
    {
        return MaterialPropertyLib::formEigenTensor<GlobalDim>(
            medium
                .property(MaterialPropertyLib::PropertyType::
                              thermal_osmosis_coefficient)
                .value(variable_array, pos, t, dt));
    }
 
    // Neither thermal_osmosis_permeability nor thermal_osmosis_coefficient is
    // defined.
    return Eigen::Matrix<double, GlobalDim, GlobalDim>::Zero();
}
}  // namespace ProcessLib