BrooksCoreyCapillaryPressureSaturation.h

Go to the documentation of this file.

 
#pragma once
 
#include <limits>  // std::numeric_limits
 
#include "CapillaryPressureSaturation.h"
 
#include "BaseLib/Error.h"
 
namespace MaterialLib
{
namespace PorousMedium
{
class BrooksCoreyCapillaryPressureSaturation final
    : public CapillaryPressureSaturation
{
public:
    BrooksCoreyCapillaryPressureSaturation(const double pb, const double Sr,
                                           const double Sg_r, const double Smax,
                                           const double m, const double Pc_max)
        : CapillaryPressureSaturation(Sr, Sg_r, Smax, Pc_max), _pb(pb), _m(m)
    {
    }
 
    std::string getName() const override
    {
        return "Brooks-Corey capillary pressure saturation model.";
    }
 
    double getCapillaryPressure(const double saturation) const override;
 
    double getSaturation(const double capillary_pressure) const override;
 
    double getdPcdS(const double saturation) const override;
 
    double getd2PcdS2(const double /*saturation*/) const override
    {
        OGS_FATAL(
            "Second derivative of BrooksCorey d^2p_cap / dS^2 is not "
            "implemented.");
    }
 
private:
    const double _pb;  
    const double _m;   
};
 
}  // namespace PorousMedium
}  // namespace MaterialLib