Density100MPa.cpp

Go to the documentation of this file.

 
#include "Density100MPa.h"
 
namespace
{
// NaX_HighP_polyfrac_CC.pickle
// date extracted 2015-06-23 15:38:35 file mtime 2015-06-23 15:19:57
const double Density100MPa_c[] = {
    0.3490302932983226,     /* a0 */
    -0.0014061345691831226, /* a1 */
    -0.0007399303393402753, /* a2 */
    5.129318840267485e-09,  /* a3 */
    5.243619689772646e-07,  /* a4 */
    6.347011955956523e-10,  /* a5 */
    -9.919599580166727e-11  /* a6 */
};
 
}  // namespace
 
namespace Adsorption
{
double Density100MPa::getAdsorbateDensity(const double T_Ads) const
{
    return -0.0013 * T_Ads * T_Ads + 0.3529 * T_Ads + 1049.2;
}
 
// Thermal expansivity model for water found in the works of Hauer
double Density100MPa::getAlphaT(const double T_Ads) const
{
    const double rho = -0.0013 * T_Ads * T_Ads + 0.3529 * T_Ads + 1049.2;
    const double drhodT = -0.0026 * T_Ads + 0.3529;
 
    return -drhodT / rho;
}
 
// Characteristic curve. Return W (A)
double Density100MPa::characteristicCurve(const double A) const
{
    double W = curvePolyfrac(Density100MPa_c, A);  // cm^3/g
 
    if (W < 0.0)
    {
        W = 0.0;  // TODO [CL] debug output
    }
 
    return W / 1.e3;  // m^3/kg
}
 
double Density100MPa::dCharacteristicCurve(const double A) const
{
    return dCurvePolyfrac(Density100MPa_c, A);
}
 
}  // namespace Adsorption