CapillaryPressureSaturationCurve.h

Go to the documentation of this file.

 
#pragma once
 
#include <algorithm>
#include <memory>
 
#include "CapillaryPressureSaturation.h"
#include "MathLib/Curve/PiecewiseLinearMonotonicCurve.h"
 
namespace MaterialLib
{
namespace PorousMedium
{
class CapillaryPressureSaturationCurve final
    : public CapillaryPressureSaturation
{
public:
    explicit CapillaryPressureSaturationCurve(
        std::unique_ptr<MathLib::PiecewiseLinearMonotonicCurve>&& curve_data)
        : CapillaryPressureSaturation(
              curve_data->getSupportMin(),
              1.0 - curve_data->getSupportMax(),
              curve_data->getSupportMax(),
              curve_data->getValue(curve_data->getSupportMin())),
          _curve_data(std::move(curve_data))
    {
    }
 
    std::string getName() const override
    {
        return "Capillary pressure saturation curve.";
    }
 
    double getCapillaryPressure(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 getSaturation(const double capillary_pressure) const override
    {
        const double pc =
            std::clamp(capillary_pressure, _minor_offset, _pc_max);
        return _curve_data->getInverseVariable(pc);
    }
 
    double getdPcdS(const double saturation) const override
    {
        const double S = std::clamp(saturation, _saturation_r + _minor_offset,
                                    _saturation_max - _minor_offset);
 
        return _curve_data->getDerivative(S);
    }
 
    double getd2PcdS2(const double /*saturation*/) const override { return 0; }
 
private:
    std::unique_ptr<MathLib::PiecewiseLinearMonotonicCurve> _curve_data;
};
 
}  // namespace PorousMedium
}  // namespace MaterialLib