SaturationDependentSwelling.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 "MaterialLib/MPL/Property.h"
 
namespace ParameterLib
{
struct CoordinateSystem;
}
 
namespace MaterialPropertyLib
{
class Medium;
class Phase;
class Component;

class SaturationDependentSwelling final : public Property
{
private:
    std::array<double, 3> const p_;
    std::array<double, 3> const lambda_;
    double const S_min_;  //< Lower saturation limit.
    double const S_max_;  //< Upper saturation limit.
    ParameterLib::CoordinateSystem const* const local_coordinate_system_;
 
public:
    SaturationDependentSwelling(
        std::string name,
        std::array<double, 3>
            swelling_pressures,
        std::array<double, 3>
            exponents,
        double const lower_saturation_limit,
        double const upper_saturation_limit,
        ParameterLib::CoordinateSystem const* const local_coordinate_system);
 
    void checkScale() const override;
 
    PropertyDataType value(VariableArray const& variable_array,
                           VariableArray const& variable_array_prev,
                           ParameterLib::SpatialPosition const& pos,
                           double const t, double const dt) const override;
    PropertyDataType dValue(VariableArray const& variable_array,
                            VariableArray const& variable_array_prev,
                            Variable const variable,
                            ParameterLib::SpatialPosition const& pos,
                            double const t, double const dt) const override;
};
 
}  // namespace MaterialPropertyLib