henryiapws.hh

Go to the documentation of this file.

// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
// vi: set et ts=4 sw=4 sts=4:
//
// SPDX-FileCopyrightInfo: Copyright © DuMux Project contributors, see AUTHORS.md in root folder
// SPDX-License-Identifier: GPL-3.0-or-later
//
#ifndef DUMUX_HENRY_IAPWS_HH
#define DUMUX_HENRY_IAPWS_HH
 
#include <dumux/material/components/h2o.hh>
 
namespace Dumux {
template <class Scalar>
inline Scalar henryIAPWS(Scalar E,
                         Scalar F,
                         Scalar G,
                         Scalar H,
                         Scalar temperature)
{
    using H2O = Dumux::Components::H2O<Scalar>;
 
    // regularizing temperature helps for stability.
    // Results are unphysical!
    if (temperature > H2O::criticalTemperature())
        temperature = H2O::criticalTemperature();
 
    Scalar Tr = temperature/H2O::criticalTemperature();
    Scalar tau = 1 - Tr;
 
    static const Scalar c[6] = {
        1.99274064, 1.09965342, -0.510839303,
        -1.75493479,-45.5170352, -6.7469445e5
    };
    static const Scalar d[6] = {
        1/3.0, 2/3.0, 5/3.0,
        16/3.0, 43/3.0, 110/3.0
    };
    static const Scalar q = -0.023767;
 
    Scalar f = 0;
    using std::pow;
    for (int i = 0; i < 6; ++i) {
        f += c[i]*pow(tau, d[i]);
    }
 
    Scalar exponent =
        q*F +
        E/temperature*f +
        (F +
         G*pow(tau, 2.0/3) +
         H*tau)*
        exp((H2O::tripleTemperature() - temperature)/100);
 
    using std::exp;
    return exp(exponent)*H2O::vaporPressure(temperature);
}
}
 
#endif