CapillaryPressureSaturation.h

Go to the documentation of this file.

 
#pragma once
 
#include <limits>
#include <string>
 
namespace MaterialLib
{
namespace PorousMedium
{
class CapillaryPressureSaturation
{
public:
    CapillaryPressureSaturation(const double Sr, const double Sg_r,
                                const double Smax, const double Pc_max)
        : _saturation_r(Sr),
          _saturation_nonwet_r(Sg_r),
          _saturation_max(Smax),
          _pc_max(Pc_max)
    {
    }
    virtual ~CapillaryPressureSaturation() = default;
 
    virtual std::string getName() const = 0;
 
    virtual double getCapillaryPressure(const double saturation) const = 0;
    virtual double getSaturation(const double capillary_pressure) const = 0;
 
    virtual double getdPcdS(const double saturation) const = 0;
 
    virtual double getd2PcdS2(const double saturation) const = 0;
 
protected:
    const double _saturation_r;         
    const double _saturation_nonwet_r;  
    const double _saturation_max;       
    const double _pc_max;               
 
    const double _minor_offset = std::numeric_limits<double>::epsilon();
};
 
}  // namespace PorousMedium
}  // namespace MaterialLib