BishopsPowerLaw.cpp

Go to the documentation of this file.

 
#include "BishopsPowerLaw.h"
 
namespace MaterialPropertyLib
{
BishopsPowerLaw::BishopsPowerLaw(std::string name, double const exponent)
    : m_(exponent)
{
    name_ = std::move(name);
}
 
void BishopsPowerLaw::checkScale() const
{
    if (!std::holds_alternative<Medium*>(scale_))
    {
        OGS_FATAL(
            "The property 'BishopsPowerLaw' is implemented on the 'media' "
            "scale only.");
    }
}
 
PropertyDataType BishopsPowerLaw::value(
    VariableArray const& variable_array,
    ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/,
    double const /*dt*/) const
{
    auto const S_L = variable_array.liquid_saturation;
 
    return std::pow(S_L, m_);
}
 
PropertyDataType BishopsPowerLaw::dValue(
    VariableArray const& variable_array, Variable const variable,
    ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/,
    double const /*dt*/) const
{
    if (variable != Variable::liquid_saturation)
    {
        OGS_FATAL(
            "BishopsPowerLaw::dValue is implemented for derivatives with "
            "respect to liquid saturation only.");
    }
 
    auto const S_L = variable_array.liquid_saturation;
 
    return m_ * std::pow(S_L, m_ - 1.);
}
}  // namespace MaterialPropertyLib