RelativePermeabilityCurve.h

Go to the documentation of this file.

 
#pragma once
 
#include <algorithm>
#include <memory>
 
#include "MathLib/InterpolationAlgorithms/PiecewiseLinearInterpolation.h"
#include "RelativePermeability.h"
 
namespace MaterialLib
{
namespace PorousMedium
{
class RelativePermeabilityCurve final : public RelativePermeability
{
public:
    explicit RelativePermeabilityCurve(
        std::unique_ptr<MathLib::PiecewiseLinearInterpolation>&& curve_data)
        : RelativePermeability(curve_data->getSupportMin(),
                               curve_data->getSupportMax()),
          _curve_data(std::move(curve_data))
    {
    }
 
    std::string getName() const override
    {
        return "Relative permeability curve.";
    }
 
    double getValue(const double saturation) const override
    {
        const double S = std::clamp(saturation, _saturation_r + _minor_offset,
                                    _saturation_max - _minor_offset);
 
        return _curve_data->getValue(S);
    }
 
    double getdValue(const double saturation) const override
    {
        const double S = std::clamp(saturation, _saturation_r + _minor_offset,
                                    _saturation_max - _minor_offset);
 
        return _curve_data->getDerivative(S);
    }
 
private:
    std::unique_ptr<MathLib::PiecewiseLinearInterpolation> _curve_data;
};
}  // namespace PorousMedium
}  // namespace MaterialLib