SaturationLuMcCartney.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#ifdef __GNUC__
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wmaybe-uninitialized"
#endif
#include <boost/math/differentiation/autodiff.hpp>
#ifdef __GNUC__
#pragma GCC diagnostic pop
#endif
 
#include "MaterialLib/MPL/Properties/GibbsFreeEnergy/DimensionLessGibbsFreeEnergyRegion1.h"
#include "MaterialLib/MPL/Property.h"
 
namespace MaterialPropertyLib
{
class Medium;
class Phase;
class Component;
class SaturationLuMcCartney final : public Property
{
public:
    SaturationLuMcCartney(std::string name, std::string const& material);
 
    void checkScale() const override
    {
        if (!std::holds_alternative<Medium*>(scale_))
        {
            OGS_FATAL(
                "The property 'SaturationLuMcCartney' 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;
    PropertyDataType d2Value(VariableArray const& variable_array,
                             Variable const variable1, Variable const variable2,
                             ParameterLib::SpatialPosition const& /*pos*/,
                             double const /*t*/,
                             double const /*dt*/) const override;
 
private:
    template <typename PressureType, typename TemperatureType>
    boost::math::differentiation::promote<PressureType, TemperatureType>
    water_content(PressureType p_cap, TemperatureType T,
                  VariableArray const& variable_array) const;
    template <typename PressureType, typename TemperatureType>
    boost::math::differentiation::promote<PressureType, TemperatureType>
    adsorptive_water_content(PressureType p_cap, TemperatureType T,
                             VariableArray const& variable_array) const;
    template <typename PressureType, typename TemperatureType>
    boost::math::differentiation::promote<PressureType, TemperatureType>
    capillary_water_content(PressureType p_cap, TemperatureType T,
                            VariableArray const& variable_array) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> theta_a_max(
        TemperatureType T, VariableArray const& variable_array) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> psi_max(
        TemperatureType T) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> c(
        TemperatureType T) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> CEC(
        TemperatureType T) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> psi_c(
        TemperatureType T, VariableArray const& variable_array) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> A_H(
        TemperatureType T) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> chi(
        TemperatureType T) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> delta_h(
        TemperatureType T) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> alpha(
        TemperatureType T) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> epsilon_w(
        TemperatureType T) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> zeta_s(
        TemperatureType T) const;
    template <typename TemperatureType>
    boost::math::differentiation::promote<TemperatureType> rho_w(
        TemperatureType T) const;
 
    std::string const material_;
 
    const MaterialLib::Fluid::DimensionLessGibbsFreeEnergyRegion1
        gibbs_free_energy_;
    // constants
    static constexpr double ref_T_ = 1386;  
    static constexpr double ref_p_ = 1.653e7;

    double const nu_w = 1.8e-5;  
    double const A = 0.92;  
    double const B = 0.9;   
    double const T1 = 373.15;  
    double const T2 =
        1273.15;  
    double const nu_e =
        2.45e9;                
    double const Tr = 293.15;  
    double const C1 = -0.00151;  
    double const n_w = 1.33;     

    double E1_minus_EL;  
    double nu_mr;        
    double CEC_max;      
    double bw;
    double M;  

    double SSA;        
    double epsilon_s;  
    double n_s;        
    double eta_alpha;  
    double
        alpha_0_inv;  
    double chi_r;
    double N;  
};
}  // namespace MaterialPropertyLib