WaterDensityIAPWSIF97Region1.h

Go to the documentation of this file.

 
#pragma once
 
#include "MaterialLib/Fluid/GibbsFreeEnergy/DimensionLessGibbsFreeEnergyRegion1.h"
#include "MaterialLib/MPL/Property.h"
 
namespace MaterialPropertyLib
{
class Phase;
 
class WaterDensityIAPWSIF97Region1 final : public Property
{
public:
    explicit WaterDensityIAPWSIF97Region1(std::string name)
        : gibbs_free_energy_(
              MaterialLib::Fluid::DimensionLessGibbsFreeEnergyRegion1())
    {
        name_ = std::move(name);
    }
    void checkScale() const override
    {
        if (!std::holds_alternative<Phase*>(scale_))
        {
            OGS_FATAL(
                "The property 'WaterDensityIAPWSIF97Region1' 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:
    const MaterialLib::Fluid::DimensionLessGibbsFreeEnergyRegion1
        gibbs_free_energy_;
 
    static constexpr double ref_T_ = 1386;     
    static constexpr double ref_p_ = 1.653e7;  
};
}  // namespace MaterialPropertyLib