OGS
WaterVapourEnthalpyIAPWSIF97Region4.cpp
Go to the documentation of this file.
1
12
#include "
WaterVapourEnthalpyIAPWSIF97Region4.h
"
13
14
#include <cmath>
15
16
#include "
BaseLib/Error.h
"
17
#include "
MaterialLib/Fluid/GibbsFreeEnergy/DimensionlessGibbsFreeEnergyRegion2.h
"
18
#include "
MaterialLib/MPL/Properties/WaterSaturationCurveIAPWSIF97Region4.h
"
19
#include "
MaterialLib/PhysicalConstant.h
"
20
21
namespace
MaterialPropertyLib
22
{
23
PropertyDataType
WaterVapourEnthalpyIAPWSIF97Region4::value
(
24
VariableArray
const
& variable_array,
25
ParameterLib::SpatialPosition
const
&
/*pos*/
,
double
const
/*t*/
,
26
double
const
/*dt*/
)
const
27
{
28
double
const
p = variable_array.
liquid_phase_pressure
;
29
36
if
((p < 611.213) || (p > 22.064e6))
37
{
38
WARN
(
39
"Pressure is out of the range for the water vapour saturation "
40
"enthalpy."
);
41
}
42
43
static
constexpr
double
ref_T_
= 540;
44
static
constexpr
double
ref_p_
= 1.e6;
45
46
double
const
T_s =
IAPWSIF97Region4::waterSaturationTemperature
(p);
47
double
const
tau =
ref_T_
/ T_s;
48
double
const
pi = p /
ref_p_
;
49
50
double
dgamma_dtau =
51
MaterialLib::Fluid::DimensionlessGibbsFreeEnergyRegion2::getdGammadTau
(
52
tau, pi);
53
54
return
tau *
55
MaterialLib::PhysicalConstant::SpecificGasConstant::WaterVapour
*
56
T_s * dgamma_dtau;
57
}
58
59
PropertyDataType
WaterVapourEnthalpyIAPWSIF97Region4::dValue
(
60
VariableArray
const
&
/*variable_array*/
,
Variable
const
/*variable*/
,
61
ParameterLib::SpatialPosition
const
&
/*pos*/
,
double
const
/*t*/
,
62
double
const
/*dt*/
)
const
63
{
64
OGS_FATAL
(
65
"WaterVapourEnthalpyIAPWSIF97Region4::dValue is not implemented."
);
66
}
67
68
}
// namespace MaterialPropertyLib
DimensionlessGibbsFreeEnergyRegion2.h
Error.h
OGS_FATAL
#define OGS_FATAL(...)
Definition
Error.h:26
WARN
void WARN(fmt::format_string< Args... > fmt, Args &&... args)
Definition
Logging.h:40
PhysicalConstant.h
WaterSaturationCurveIAPWSIF97Region4.h
WaterVapourEnthalpyIAPWSIF97Region4.h
MaterialPropertyLib::Property::value
virtual PropertyDataType value() const
Definition
Property.cpp:76
MaterialPropertyLib::VariableArray
Definition
VariableType.h:99
MaterialPropertyLib::VariableArray::liquid_phase_pressure
double liquid_phase_pressure
Definition
VariableType.h:181
MaterialPropertyLib::WaterVapourEnthalpyIAPWSIF97Region4::dValue
PropertyDataType dValue(VariableArray const &variable_array, Variable const variable, ParameterLib::SpatialPosition const &pos, double const t, double const dt) const override
Definition
WaterVapourEnthalpyIAPWSIF97Region4.cpp:59
ParameterLib::SpatialPosition
Definition
SpatialPosition.h:27
MaterialLib::Fluid::DimensionlessGibbsFreeEnergyRegion2::getdGammadTau
double getdGammadTau(const double tau, const double pi)
Definition
DimensionlessGibbsFreeEnergyRegion2.cpp:85
MaterialLib::PhysicalConstant::SpecificGasConstant::WaterVapour
constexpr double WaterVapour
Specific gas constant for water vapour.
Definition
PhysicalConstant.h:95
MaterialPropertyLib::IAPWSIF97Region4::waterSaturationTemperature
double waterSaturationTemperature(double const pressure)
Definition
WaterSaturationCurveIAPWSIF97Region4.h:23
MaterialPropertyLib
Definition
ChemicalSolverInterface.h:21
MaterialPropertyLib::ref_p_
static constexpr double ref_p_
reference pressure in Pa.
Definition
WaterEnthalpyIAPWSIF97Region1.cpp:27
MaterialPropertyLib::Variable
Variable
Definition
VariableType.h:30
MaterialPropertyLib::ref_T_
static constexpr double ref_T_
reference temperature in K.
Definition
WaterEnthalpyIAPWSIF97Region1.cpp:26
MaterialPropertyLib::PropertyDataType
std::variant< double, Eigen::Matrix< double, 2, 1 >, Eigen::Matrix< double, 3, 1 >, Eigen::Matrix< double, 2, 2 >, Eigen::Matrix< double, 3, 3 >, Eigen::Matrix< double, 4, 1 >, Eigen::Matrix< double, 6, 1 >, Eigen::MatrixXd > PropertyDataType
Definition
Property.h:31
MaterialLib
MPL
Properties
Enthalpy
WaterVapourEnthalpyIAPWSIF97Region4.cpp
Generated by
1.12.0