coagulationModel_FM.cc

Go to the documentation of this file.

#include "coagulationModels/coagulationModel_FM.h"
#include "sootDefs.h"
 
using namespace std;
using namespace soot;
 
 
double coagulationModel_FM::getCoagulationSootRate(const state& state,
                                                   double m1,
                                                   double m2) const {
 
    if (m1 <= mmin || m2 <= mmin) return 0.0;
 
    double m13_m13 = pow(m1, onethird) + pow(m2, onethird);
 
    return getKfm(state) * sqrt(1/m1 + 1/m2)*m13_m13*m13_m13;
}