sootModel_MOMIC.h

Go to the documentation of this file.

#pragma once
 
#include "sootModel.h"
#include "sootDefs.h"
#include "state.h"
 
namespace soot {
 
 
class sootModel_MOMIC : public sootModel {
 
 
private:
    size_t Nmom;              
 
    std::vector<double> Mp6;  
    std::vector<double> Mq6;  
    std::vector<size_t> np;   
    std::vector<size_t> nq;   
 
    std::vector<std::vector<double> > diffTable;     
 
    double mDn46,                                    
           mDn26, mD26;
    double mDn36,                                    
           mDn16, mD16,  mD36,  mD56,  mD76,
           mD96,  mD116, mD136, mD156, mD176, mD196;
 
 
public:
 
    virtual void setSourceTerms(state &state);
 
    virtual double pahSootCollisionRatePerDimer(const state &state, const double mDimer);
 
private:
 
    void                downselectIfNeeded(std::vector<double> &M);
    double              f_grid(int x, int y);
    double              g_grid(int y);
    std::vector<double> MOMICCoagulationRates(const state& state, std::vector<double>& M);
    double              Mr(const double r);                              // fractional moments using polynomial interp among log M
    void                set_diffTable(const std::vector<double> &M);  // needed by Mr
    void                set_fractional_moments_Mp6_Mq6();
    void                set_mDimerPowers();
 
 
public:
 
    sootModel_MOMIC(size_t            nsoot_,
                    nucleationModel  *nucl_,
                    growthModel      *grow_,
                    oxidationModel   *oxid_,
                    coagulationModel *coag_);
 
    sootModel_MOMIC(size_t           nsoot_,
                    nucleationMech   Nmech,
                    growthMech       Gmech,
                    oxidationMech    Omech,
                    coagulationMech  Cmech);
 
    virtual ~sootModel_MOMIC() {};
 
};
} // namspace soot