50 if (temperature > H2O::criticalTemperature())
51 temperature = H2O::criticalTemperature();
53 Scalar Tr = temperature/H2O::criticalTemperature();
56 static const Scalar c[6] = {
57 1.99274064, 1.09965342, -0.510839303,
58 -1.75493479,-45.5170352, -6.7469445e5
60 static const Scalar d[6] = {
62 16/3.0, 43/3.0, 110/3.0
64 static const Scalar q = -0.023767;
68 for (
int i = 0; i < 6; ++i) {
69 f += c[i]*pow(tau, d[i]);
78 exp((H2O::tripleTemperature() - temperature)/100);
81 return exp(exponent)*H2O::vaporPressure(temperature);
Scalar henryIAPWS(Scalar E, Scalar F, Scalar G, Scalar H, Scalar temperature)
The Henry constants in liquid water using the IAPWS 2004 formulation.
Definition henryiapws.hh:40