state.h

Go to the documentation of this file.

#pragma once
 
#include "sootDefs.h"
 
#include <vector>
 
namespace soot {
 
 
class state {
 
 
    public:
 
        double              T = 0;         
        double              P = 0;         
        double              rhoGas = 0;    
        double              MWGas = 0;     
        double              muGas = 0;     
        std::vector<double> yGas;          
        std::vector<double> yPah;          
 
        int                 nsoot;         
        std::vector<double> sootVar;       
        std::vector<double> absc;          
        std::vector<double> wts;           
        std::vector<double> sootScales;    
        double              cMin = 100;    
 
 
 
        void setState(double T_, double P_, double rhoGas_, double muGas_,
                      std::vector<double> yGas_, std::vector<double> yPAH_, 
                      std::vector<double> sootVar_, int nsoot_, double cMin_ = 100);
 
        double getGasSpC(gasSp sp)  const { return rhoGas * yGas[(int)sp] / gasSpMW[(int)sp]; }
 
        double getGasSpP(gasSp sp)  const { return yGas[(int)sp] * MWGas / gasSpMW[(int)sp] * P; }
 
        double getGasMeanFreePath() const { return muGas/rhoGas*sqrt(M_PI*MWGas/(2.0*Rg*T)); }
 
        double get_pahSpC(pahSp sp) const { return rhoGas * yPah[(int)sp] / pahSpMW[(int)sp]; }
 
        double get_pahSpP(pahSp sp) const { return yPah[(int)sp] * MWGas  / pahSpMW[(int)sp] * P; }
 
        void setSootScales(std::vector<double> &sootScales_) { sootScales = sootScales_; }
 
 
    state(size_t nsoot_=0) :
        nsoot(nsoot_), 
        sootVar(std::vector<double>(   nsoot_,           0.0)),
        sootScales(std::vector<double>(nsoot_,           1.0)),
        absc(std::vector<double>(      nsoot_/2,         0.0)),
        wts(std::vector<double>(       nsoot_/2,         0.0)),
        yGas(std::vector<double>(      (int)gasSp::size, 0.0)),
        yPah(std::vector<double>(      (int)pahSp::size, 0.0)) { }
    ~state() = default;
 
};
}     // namespace soot