OGS
Adsorption::DensityNunez Class Reference

Detailed Description

Definition at line 17 of file DensityNunez.h.

#include <DensityNunez.h>

Inheritance diagram for Adsorption::DensityNunez:
[legend]
Collaboration diagram for Adsorption::DensityNunez:
[legend]

Public Member Functions

double getAdsorbateDensity (const double T_Ads) const override
 
double getAlphaT (const double T_Ads) const override
 
double characteristicCurve (const double A) const override
 
double dCharacteristicCurve (const double A) const override
 
- Public Member Functions inherited from Adsorption::AdsorptionReaction
double getEquilibriumLoading (const double p_Ads, const double T_Ads, const double M_Ads) const override
 
double getEnthalpy (const double p_Ads, const double T_Ads, const double M_Ads) const override
 
double getReactionRate (const double p_Ads, const double T_Ads, const double M_Ads, const double loading) const override
 
- Public Member Functions inherited from Adsorption::Reaction
virtual ~Reaction ()=default
 

Additional Inherited Members

- Static Public Member Functions inherited from Adsorption::AdsorptionReaction
static double getEvaporationEnthalpy (const double T_Ads)
 
static double getEquilibriumVapourPressure (const double T_Ads)
 
static double getMolarFraction (double xm, double M_this, double M_other)
 
static double dMolarFraction (double xm, double M_this, double M_other)
 
static double getLoading (const double rho_curr, const double rho_dry)
 
- Static Public Member Functions inherited from Adsorption::Reaction
static std::unique_ptr< ReactionnewInstance (BaseLib::ConfigTree const &conf)
 

Member Function Documentation

◆ characteristicCurve()

double Adsorption::DensityNunez::characteristicCurve ( const double  A) const
overridevirtual

Implements Adsorption::AdsorptionReaction.

Definition at line 72 of file DensityNunez.cpp.

73 {
74  double W = curvePolyfrac(DensityNunez_c, A); // cm^3/g
75 
76  if (W < 0.0)
77  {
78  W = 0.0; // TODO [CL] debug output
79  }
80 
81  return W / 1.e3; // m^3/kg
82 }
double curvePolyfrac(const double *coeffs, const double x)
Definition: Adsorption.h:52

References Adsorption::curvePolyfrac(), and anonymous_namespace{DensityNunez.cpp}::DensityNunez_c.

◆ dCharacteristicCurve()

double Adsorption::DensityNunez::dCharacteristicCurve ( const double  A) const
overridevirtual

Implements Adsorption::AdsorptionReaction.

Definition at line 84 of file DensityNunez.cpp.

85 {
86  return dCurvePolyfrac(DensityNunez_c, A);
87 }
double dCurvePolyfrac(const double *coeffs, const double x)
Definition: Adsorption.h:61

References Adsorption::dCurvePolyfrac(), and anonymous_namespace{DensityNunez.cpp}::DensityNunez_c.

◆ getAdsorbateDensity()

double Adsorption::DensityNunez::getAdsorbateDensity ( const double  T_Ads) const
overridevirtual

Implements Adsorption::AdsorptionReaction.

Definition at line 30 of file DensityNunez.cpp.

31 {
32  // TODO admissible T range: 273.16 K <= T_Ads <= 633.15 K
33  const double a[] = {1.0644e3, -8.01905, 1.445348e-2, -4.19589e-6,
34  -4.5294e-9};
35  const double b[] = {-8.039e-3, 1.8698e-5, -2.3015e-8, 2.3809e-11,
36  -1.388e-14};
37  const double u =
38  a[0] + T_Ads * (a[1] + T_Ads * (a[2] + T_Ads * (a[3] + T_Ads * a[4])));
39  const double v =
40  1.0 +
41  T_Ads *
42  (b[0] +
43  T_Ads * (b[1] + T_Ads * (b[2] + T_Ads * (b[3] + T_Ads * b[4]))));
44  return u / v;
45 }

◆ getAlphaT()

double Adsorption::DensityNunez::getAlphaT ( const double  T_Ads) const
overridevirtual

Implements Adsorption::AdsorptionReaction.

Definition at line 48 of file DensityNunez.cpp.

49 {
50  // TODO admissible T range: 273.16 K <= T_Ads <= 633.15 K
51  const double a[] = {1.0644e3, -8.01905, 1.445348e-2, -4.19589e-6,
52  -4.5294e-9};
53  const double b[] = {-8.039e-3, 1.8698e-5, -2.3015e-8, 2.3809e-11,
54  -1.388e-14};
55  const double u =
56  a[0] + T_Ads * (a[1] + T_Ads * (a[2] + T_Ads * (a[3] + T_Ads * a[4])));
57  const double v =
58  1.0 +
59  T_Ads *
60  (b[0] +
61  T_Ads * (b[1] + T_Ads * (b[2] + T_Ads * (b[3] + T_Ads * b[4]))));
62  const double du =
63  a[1] + T_Ads * (2.0 * a[2] + T_Ads * (3.0 * a[3] + T_Ads * 4.0 * a[4]));
64  const double dv =
65  b[0] + T_Ads * (2.0 * b[1] +
66  T_Ads * (3.0 * b[2] +
67  T_Ads * (4.0 * b[3] + T_Ads * 5.0 * b[4])));
68  return dv / v - du / u;
69 }

The documentation for this class was generated from the following files: