sootModel.cc

Go to the documentation of this file.

#include "sootModel.h"
#include <iostream>
 
#include "nucleationModels/nucleationModel_NONE.h"
#include "nucleationModels/nucleationModel_LL.h"
#include "nucleationModels/nucleationModel_LIN.h"
#include "nucleationModels/nucleationModel_LINA1.h"
#include "nucleationModels/nucleationModel_PAH.h"
 
#include "growthModels/growthModel_NONE.h"
#include "growthModels/growthModel_LL.h"
#include "growthModels/growthModel_LIN.h"
#include "growthModels/growthModel_HACA.h"
 
#include "oxidationModels/oxidationModel_NONE.h"
#include "oxidationModels/oxidationModel_LL.h"
#include "oxidationModels/oxidationModel_HACA.h"
#include "oxidationModels/oxidationModel_LEE_NEOH.h"
#include "oxidationModels/oxidationModel_NSC_NEOH.h"
#include "oxidationModels/oxidationModel_OPTJ.h"
#include "oxidationModels/oxidationModel_OPTG.h"
 
#include "coagulationModels/coagulationModel_NONE.h"
#include "coagulationModels/coagulationModel_FM.h"
#include "coagulationModels/coagulationModel_CONTINUUM.h"
#include "coagulationModels/coagulationModel_HM.h"
#include "coagulationModels/coagulationModel_FUCHS.h"
 
using namespace std;
using namespace soot;
 
 
sootModel::sootModel(size_t            nsoot_,
                     nucleationModel  *nucl_,
                     growthModel      *grow_,
                     oxidationModel   *oxid_,
                     coagulationModel *coag_) :
        nsoot(nsoot_), nucl(nucl_), grow(grow_), oxid(oxid_), coag(coag_), 
        mechsNewedHere(false), sources(nsoot_) {
        checkSpec();
        nucl->SM = this;
        grow->SM = this;
        oxid->SM = this;
        coag->SM = this;
}
 
 
 
sootModel::sootModel(size_t          nsoot_,
                     nucleationMech  Nmech,
                     growthMech      Gmech,
                     oxidationMech   Omech,
                     coagulationMech Cmech) : nsoot(nsoot_), 
                                              mechsNewedHere(true),
                                              sources(nsoot_) {
 
    //---------- set nucleation model
 
    switch (Nmech) {
        case nucleationMech::NONE  : nucl = new nucleationModel_NONE();   break;
        case nucleationMech::LL    : nucl = new nucleationModel_LL();     break;
        case nucleationMech::LIN   : nucl = new nucleationModel_LIN();    break;
        case nucleationMech::LINA1 : nucl = new nucleationModel_LINA1();  break;
        case nucleationMech::PAH   : nucl = new nucleationModel_PAH();    break;
        default: throw domain_error("Invalid nucleation model requested");
    }
 
    //---------- set growth model
 
    switch (Gmech) {
        case growthMech::NONE : grow = new growthModel_NONE(); break;
        case growthMech::LL   : grow = new growthModel_LL();   break;
        case growthMech::LIN  : grow = new growthModel_LIN();  break;
        case growthMech::HACA : grow = new growthModel_HACA(); break;
        default: throw domain_error("Invalid growth model requested");
    }
 
    //---------- set oxidation model
 
    switch (Omech) {
        case oxidationMech::NONE     : oxid = new oxidationModel_NONE();     break;
        case oxidationMech::LL       : oxid = new oxidationModel_LL();       break;
        case oxidationMech::HACA     : oxid = new oxidationModel_HACA();     break;
        case oxidationMech::LEE_NEOH : oxid = new oxidationModel_LEE_NEOH(); break;
        case oxidationMech::NSC_NEOH : oxid = new oxidationModel_NSC_NEOH(); break;
        case oxidationMech::OPTJ     : oxid = new oxidationModel_OPTJ();     break;
        case oxidationMech::OPTG     : oxid = new oxidationModel_OPTG();     break;
        default: throw domain_error("Invalid oxidation model requested");
    }
 
    //---------- set coagulation model
 
    switch (Cmech) {
        case coagulationMech::NONE      : coag = new coagulationModel_NONE();      break;
        case coagulationMech::FM        : coag = new coagulationModel_FM();        break;
        case coagulationMech::CONTINUUM : coag = new coagulationModel_CONTINUUM(); break;
        case coagulationMech::HM        : coag = new coagulationModel_HM();        break;
        case coagulationMech::FUCHS     : coag = new coagulationModel_FUCHS();     break;
        default: throw domain_error("Invalid coagulation model requested");
    }
 
    //----------- 
 
    nucl->SM = this;
    grow->SM = this;
    oxid->SM = this;
    coag->SM = this;
 
    //----------- 
 
    checkSpec();
}
 
 
void sootModel::checkSpec() {
 
    if((grow->mechType==growthMech::HACA    && oxid->mechType!=oxidationMech::HACA) ||
       (oxid->mechType==oxidationMech::HACA && grow->mechType!=growthMech::HACA))
        cerr << endl
             << "**********************************************************************" << endl
             << "Warning: HACA surface growth should pair with HACA oxid and vice-versa" << endl
             << "**********************************************************************" << endl << endl;
        //throw domain_error("Invalid model combination: HACA surface growth must pair with HACA oxid and vice-versa");
}