NonWettingPhaseBrooksCoreyOilGas.h

Go to the documentation of this file.

 
#pragma once
 
#include "RelativePermeability.h"
 
namespace MaterialLib
{
namespace PorousMedium
{
class NonWettingPhaseBrooksCoreyOilGas final : public RelativePermeability
{
public:
    NonWettingPhaseBrooksCoreyOilGas(const double Snr, const double Snmax,
                                     const double m, const double krel_min)
        : RelativePermeability(1. - Snmax, 1. - Snr), _m(m), _krel_min(krel_min)
    {
    }
 
    std::string getName() const override
    {
        return "Non-wetting phase Brook-Corey relative permeability model.";
    }
 
    double getValue(const double saturation_w) const override;
 
    double getdValue(const double saturation_w) const override;
 
private:
    const double _m;         
    const double _krel_min;  
};
 
}  // namespace PorousMedium
}  // namespace MaterialLib