LinearSaturationSwellingStress.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "LinearSaturationSwellingStress.h"
 
#include "MaterialLib/MPL/Medium.h"
#include "MaterialLib/MPL/VariableType.h"
 
namespace MaterialPropertyLib
{
LinearSaturationSwellingStress::LinearSaturationSwellingStress(
    std::string name, double const coefficient,
    double const reference_saturation)
    : coefficient_(coefficient), reference_saturation_(reference_saturation)
{
    name_ = std::move(name);
}
 
PropertyDataType LinearSaturationSwellingStress::value(
    const VariableArray& /*variable_array*/,
    const ParameterLib::SpatialPosition& /*pos*/, const double /*t*/,
    const double /*dt*/) const
{
    OGS_FATAL(
        "LinearSaturationSwellingStress value call requires previous time step "
        "values.");
}
 
PropertyDataType LinearSaturationSwellingStress::value(
    const VariableArray& variable_array,
    const VariableArray& variable_array_prev,
    const ParameterLib::SpatialPosition& /*pos*/, const double /*t*/,
    const double /*dt*/) const
{
    // Sl <= S_max is guaranteed by the saturation property or
    // the saturation calculation.
    const double Sl = variable_array.liquid_saturation;
 
    if (Sl < reference_saturation_)
    {
        return 0.0;
    }
 
    const double Sl_prev = variable_array_prev.liquid_saturation;
 
    return coefficient_ * (Sl - Sl_prev);
}
 
PropertyDataType LinearSaturationSwellingStress::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(
            "LinearSaturationSwellingStress::dValue is implemented for "
            "derivatives with respect to liquid saturation only.");
    }
 
    // Sl <= S_max is guaranteed by the saturation property or
    // the saturation calculation.
    const double Sl = variable_array.liquid_saturation;
 
    return Sl < reference_saturation_ ? 0.0 : coefficient_;
}
 
}  // namespace MaterialPropertyLib