VanGenuchtenCapillaryPressureSaturation.h

Go to the documentation of this file.

 
#pragma once
 
#include <cmath>
#include <limits>
#include "CapillaryPressureSaturation.h"
 
namespace MaterialLib
{
namespace PorousMedium
{
class VanGenuchtenCapillaryPressureSaturation final
    : public CapillaryPressureSaturation
{
public:
    VanGenuchtenCapillaryPressureSaturation(const double pb, const double Sr,
                                            const double Sg_r,
                                            const double Smax, const double m,
                                            const double Pc_max,
                                            const bool has_regularized)
        : CapillaryPressureSaturation(Sr, Sg_r, Smax, Pc_max),
          _pb(pb),
          _m(m),
          _has_regularized(has_regularized)
    {
    }
 
    std::string getName() const override
    {
        return "van Genuchten water retention 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;
 
private:
    const double _pb;             
    const double _m;              
    const bool _has_regularized;  
    const double _xi = 1e-5;  
 
private:
    double getPcBarvGSg(double Sg) const;
    double getSBar(double Sg) const;
    double getPcvGSg(double Sg) const;
    double getdPcdSvGBar(double Sg) const;
    double getdPcdSvG(const double Sg) const;
};
 
}  // namespace PorousMedium
}  // namespace MaterialLib