DensityLegacy.cpp

Go to the documentation of this file.

 
#include "DensityLegacy.h"
 
namespace
{
// parameters from least squares fit (experimental data)
const double DensityLegacy_c[] = {
    0.34102920966608297,    -0.0013106032830951296, -0.00060754147575378876,
    3.7843404172683339e-07, 4.0107503869519016e-07, 3.1274595098338057e-10,
    -7.610441241719489e-11};
const double DensityLegacy_c[] = {…};
 
}  // namespace
namespace {…}
 
namespace Adsorption
{
double DensityLegacy::getAdsorbateDensity(const double T_Ads) const
{
    // set reference state for adsorbate EOS in Hauer
    const double T0 = 293.15;
    const double rho0 = 998.084;
    const double alpha0 = 2.06508e-4;  // K; kg/m^3; 1/K
 
    return (rho0 * (1. - alpha0 * (T_Ads - T0)));  // in kg/m^3
}
double DensityLegacy::getAdsorbateDensity(const double T_Ads) const {…}
 
// Thermal expansivity model for water found in the works of Hauer
double DensityLegacy::getAlphaT(const double T_Ads) const
{
    // set reference state for adsorbate EOS in Hauer
    const double T0 = 293.15;
    const double alpha0 = 2.06508e-4;  // K; 1/K
 
    return (alpha0 / (1. - alpha0 * (T_Ads - T0)));  // in 1/K
}
double DensityLegacy::getAlphaT(const double T_Ads) const {…}
 
// Characteristic curve. Return W (A)
double DensityLegacy::characteristicCurve(const double A) const
{
    double W = curvePolyfrac(DensityLegacy_c, A);  // cm^3/g
 
    if (W < 0.0)
    {
        W = 0.0;  // TODO [CL] debug output
    }
 
    return W / 1.e3;  // m^3/kg
}
double DensityLegacy::characteristicCurve(const double A) const {…}
 
double DensityLegacy::dCharacteristicCurve(const double A) const
{
    return dCurvePolyfrac(DensityLegacy_c, A);
}
double DensityLegacy::dCharacteristicCurve(const double A) const {…}
 
}  // namespace Adsorption