Revision e1aa6fc73e9dc533aa874ad45a7ea8ae3ceef455 authored by Lars Bilke on 21 May 2021, 08:55:35 UTC, committed by Lars Bilke on 21 May 2021, 08:55:35 UTC
Small fixes See merge request ogs/ogs!3631
Density100MPa.cpp
/**
* \file
* \copyright
* Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*/
#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
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...