coagulationModel.h

Go to the documentation of this file.

#pragma once
 
#include "state.h"
 
namespace soot {
 
class sootModel;
 
 
class coagulationModel {
 
 
public:
 
    sootModel *SM;                
    coagulationMech mechType;     
 
protected:
 
    double FM_multiplier;         
 
 
public:
 
    virtual double getCoagulationSootRate(const state& state, double m1, double m2) const = 0;
 
    virtual double getKfm(const state& state) const {
        throw std::runtime_error("ERROR: coagulationModel::getKfm should not be called, but child coagulationModel_FM::getKfm is fine");
        return -1.0;}
 
    virtual double getKc(const state& state) const {
        throw std::runtime_error("ERROR: coagulationModel::getKc should not be called, but child coagulationModel_CONTINUUM::getKc is fine");
        return -1.0;}
 
    virtual double getKcp(const state& state) const {
        throw std::runtime_error("ERROR: coagulationModel::getKcp should not be called, but child coagulationModel_CONTINUUM::getKcp is fine");
        return -1.0;}
 
    void set_FM_multiplier(double FM_multiplier_) {
        FM_multiplier = FM_multiplier_;
    }
 
 
    coagulationModel() = default;
    virtual ~coagulationModel() = default;      // Virtual deconstructor since interface will be used in the form of pointers
 
};
} // namespace soot