NonWettingPhaseVanGenuchten.h

Go to the documentation of this file.

#pragma once
 
#include "RelativePermeability.h"
 
namespace MaterialLib
{
namespace PorousMedium
{
class NonWettingPhaseVanGenuchten final : public RelativePermeability
{
public:
    NonWettingPhaseVanGenuchten(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 van Genuchten 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