WaterViscosityIAPWS.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 Phase;

class WaterViscosityIAPWS final : public Property
{
public:
    explicit WaterViscosityIAPWS(std::string name) { name_ = std::move(name); }
    void checkScale() const override
    {
        if (!std::holds_alternative<Phase*>(scale_))
        {
            OGS_FATAL(
                "The property 'WaterViscosityIAPWS' is "
                "implemented on the 'Phase' 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:
    static constexpr double ref_T_ = 647.096;  
    static constexpr double ref_rho_ =
        322.0;  
    static constexpr double ref_mu_ = 1.0e-6;  
 
    // Coefficients Hi and Hij are given in two static arrays in the cpp file.
};
}  // namespace MaterialPropertyLib