CapillaryPressureRegularizedVanGenuchten.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "MaterialLib/MPL/Property.h"
 
namespace MaterialPropertyLib
{
class Medium;

class CapillaryPressureRegularizedVanGenuchten final : public Property
{
public:
    CapillaryPressureRegularizedVanGenuchten(
        double const residual_liquid_saturation,
        double const maximum_liquid_saturation,
        double const exponent,
        double const p_b);
 
    void checkScale() const override
    {
        if (!std::holds_alternative<Medium*>(scale_))
        {
            OGS_FATAL(
                "The property 'CapillaryPressureRegularizedVanGenuchten' is "
                "implemented on the 'media' scale only.");
        }
    }

    PropertyDataType value(VariableArray const& variable_array,
                           ParameterLib::SpatialPosition const& pos,
                           double const t,
                           double const dt) const override;

    PropertyDataType dValue(VariableArray const& variable_array,
                            Variable const variable,
                            ParameterLib::SpatialPosition const& pos,
                            double const t,
                            double const dt) const override;
 
private:
    double const Sg_r_;    
    double const Sg_max_;  
    double const m_;       
    double const p_b_;
    static constexpr double xi_ = 1e-5;
 
    double const PcBarvGSg_Sg_max_;
    double const dPcdSvGBarSg_max_;

    double getPcBarvGSg(double const Sg) const;
    double getSBar(double const Sg) const;
    double getPcvGSg(double const Sg) const;
    double getdPcdSvGBar(double const Sg) const;
    double getdPcdSvG(double const Sg) const;
};
 
}  // namespace MaterialPropertyLib