Heptane.h

Go to the documentation of this file.

//! @file Heptane.h
#ifndef TPX_HEPTANE_H
#define TPX_HEPTANE_H

#include "cantera/tpx/Sub.h"

namespace tpx
{
//! Pure species representation of heptane. Values and functions are
//! from "Thermodynamic Properties in SI" by W.C. Reynolds
class Heptane : public Substance
{
public:
    Heptane() {
        m_name = "Heptane";
        m_formula = "C7H16";
    }

    double MolWt();
    double Tcrit();
    double Pcrit();
    double Vcrit();
    double Tmin();
    double Tmax();

    //! Pressure. Equation P-2 in Reynolds.
    double Pp();

    /*!
     * internal energy.
     * See Reynolds eqn (15) section 2
     *  u = (the integral from T to To of co(T)dT) +
     *         sum from i to N ([C(i) - T*Cprime(i)] + uo
     */
    double up();

    //! Entropy. See Reynolds eqn (16) section 2
    double sp();

    //! Pressure at Saturation. Equation S-2 in Reynolds.
    double Psat();

private:
    //! liquid density. Equation D2 in Reynolds.
    double ldens();

    /*!
     * C returns a multiplier in each term of the sum
     * in P-2, used in conjunction with C in the function Pp
     * - j is used to represent which of the values in the summation to calculate
     * - j=0 is the second additive in the formula in reynolds
     * - j=1 is the third...
     */
    double C(int jm, double, double, double, double);

    //! derivative of C(i)
    double Cprime(int i,  double, double, double);

    /*!
     * I = integral from o-rho { 1/(rho^2) * H(i, rho) d rho }
     * ( see section 2 of Reynolds TPSI )
     */
    double I(int i,  double, double);

    /*!
     * H returns a multiplier in each term of the sum in P-2.
     * this is used in conjunction with C in the function Pp
     * this represents the product rho^n
     * - i=0 is the second additive in the formula in reynolds
     * - i=1 is the third ...
     */
    double H(int i, double egrho);
};

}

#endif // ! TPX_HEPTANE_H