MMCollisionInt.h

Go to the documentation of this file.

/**
 * @file MMCollisionInt.h
 *  Monchick and Mason collision integrals
 */
 
// This file is part of Cantera. See License.txt in the top-level directory or
// at https://cantera.org/license.txt for license and copyright information.
 
#ifndef CT_MMCOLLISIONINT_H
#define CT_MMCOLLISIONINT_H
 
#include "cantera/base/ct_defs.h"
 
namespace Cantera
{
 
//! Calculation of Collision integrals
/*!
 * This class provides functions that interpolate the tabulated collision
 * integrals in Monchick and Mason @cite monchick1961.
 *
 * @ingroup tranprops
 */
class MMCollisionInt
{
public:
    MMCollisionInt() {}
    virtual ~MMCollisionInt() {}
 
    //! Initialize the object for calculation
    /*!
     *  @param tsmin       Minimum value of Tstar to carry out the fitting
     *  @param tsmax       Maximum value of Tstar to carry out the fitting
     *  @param loglevel    Set the loglevel for the object. The default
     *                     loglevel is zero, indicating no output.
     */
    void init(double tsmin, double tsmax, int loglevel = 0);
 
    double omega22(double ts, double deltastar);
    double astar(double ts, double deltastar);
    double bstar(double ts, double deltastar);
    double cstar(double ts, double deltastar);
    void fit(int degree, double deltastar, double* astar, double* bstar, double* cstar);
    void fit_omega22(int degree, double deltastar, double* om22);
    double omega11(double ts, double deltastar) {
        return omega22(ts, deltastar)/astar(ts, deltastar);
    }
 
private:
    double fitDelta(int table, int ntstar, int degree, double* c);
 
    vector<vector<double>> m_o22poly;
    vector<vector<double>> m_apoly;
    vector<vector<double>> m_bpoly;
    vector<vector<double>> m_cpoly;
 
    static double delta[8];
    static double tstar22[37];
 
    //! Table of omega22 values from MM
    static double omega22_table[37*8];
 
    //! table of tstar values
    static double tstar[39];
 
    //! astar table from MM
    static double astar_table[39*8];
 
    //! bstar table from MM
    static double bstar_table[39*8];
 
    //! cstar table from MM
    static double cstar_table[39*8];
 
    //! Log temp
    vector<double> m_logTemp;
 
    int m_nmin;
    int m_nmax;
 
    //! loglevel
    int m_loglevel;
};
}
#endif