sootModel.h

Go to the documentation of this file.

#pragma once
 
#include "sootDefs.h"
#include "state.h"
 
#include "nucleationModels/nucleationModel.h"
#include "growthModels/growthModel.h"
#include "oxidationModels/oxidationModel.h"
#include "coagulationModels/coagulationModel.h"
 
#include <vector>
 
namespace soot {
 
 
class sootModel {
 
 
public:
 
    size_t            nsoot;    
 
    nucleationModel  *nucl;     
    growthModel      *grow;     
    oxidationModel   *oxid;     
    coagulationModel *coag;     
 
    bool    mechsNewedHere;     
 
    psdMech psdMechType;        
 
    //----------- for the SECT model interface
 
    std::vector<double> mBins;  
    sourceTerms sources;        
 
 
    virtual void setSourceTerms(state &state) = 0;
 
    void checkSpec();
 
    //----------- for the LOGN model interface
 
    virtual double pahSootCollisionRatePerDimer(const double mDimer) const {return 0;}
 
    //----------- for the SECT model interface
 
    virtual double get_M0_sectional(const state &state) {
        throw std::runtime_error("ERROR: sootModel::get_M0_sectional should not be called, but child sootModel_SECT::get_M0_sectional is fine");
        return 0;
    }
    virtual double get_M1_sectional(const state &state) {
        throw std::runtime_error("ERROR: sootModel::get_M1_sectional should not be called, but child sootModel_SECT::get_M1_sectional is fine");
        return 0;
    }
 
 
    sootModel(size_t            nsoot_,
              nucleationModel  *nucl_,
              growthModel      *grow_,
              oxidationModel   *oxid_,
              coagulationModel *coag_);
 
    sootModel(size_t          nsoot_,
              nucleationMech  Nmech,
              growthMech      Gmech,
              oxidationMech   Omech,
              coagulationMech Cmech);
 
    virtual ~sootModel() {
        if(mechsNewedHere){
            delete(nucl);
            delete(grow);
            delete(oxid);
            delete(coag);
        }
    }
 
};
} // namespace soot