WettingPhaseVanGenuchten.h

Go to the documentation of this file.

 
#pragma once
 
#include "RelativePermeability.h"
 
namespace MaterialLib
{
namespace PorousMedium
{
class WettingPhaseVanGenuchten final : public RelativePermeability
{
public:
    WettingPhaseVanGenuchten(const double Sr, const double Smax, const double m,
                             const double krel_min)
        : RelativePermeability(Sr, Smax), _m(m), _krel_min(krel_min)
    {
    }

    std::string getName() const override
    {
        return "Wetting phase van Genuchten relative permeability model.";
    }

    double getValue(const double saturation) const override;

    double getdValue(const double saturation) const override;
 
private:
    const double _m;         
    const double _krel_min;  
};
 
}  // namespace PorousMedium
}  // namespace MaterialLib