ClausiusClapeyron.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 "BaseLib/ConfigTree-fwd.h"
#include "MaterialLib/MPL/Property.h"
 
namespace MaterialPropertyLib
{
class Medium;
class Phase;
class Component;
class ClausiusClapeyron final : public Property
{
public:
    explicit ClausiusClapeyron(std::string name,
                               const double triple_temperature,
                               const double triple_pressure,
                               const double critical_temperature,
                               const double critical_pressure,
                               const double ref_temperature,
                               const double ref_pressure);
 
    void checkScale() const override;
 
    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:
    double T_triple_;
    double p_triple_;
    double T_critical_;
    double p_critical_;
    double T_ref_;
    double p_ref_;
 
    double molarMass(std::variant<Medium*, Phase*, Component*> const scale,
                     VariableArray const& variable_array,
                     ParameterLib::SpatialPosition const& pos, double const t,
                     double const dt) const;
};
 
}  // namespace MaterialPropertyLib