From 8e31407c61e3deb4709ab6077c143d6773368c9f Mon Sep 17 00:00:00 2001 From: Yeongdo Park Date: Wed, 7 Nov 2018 07:25:15 -0500 Subject: [PATCH] dynamic cast to derived reaction classes and print Arrhenius parameters and third body efficiencies --- ct2foam.C | 27 +++++++++++++++++++++++++++ 1 file changed, 27 insertions(+) diff --git a/ct2foam.C b/ct2foam.C index b798379..dc5e158 100644 --- a/ct2foam.C +++ b/ct2foam.C @@ -137,14 +137,41 @@ int main(int argc, char *argv[]) { case Cantera::ELEMENTARY_RXN: Info << "ELEMENTARY_RXN" << endl; + + Info << std::dynamic_pointer_cast(r)->rate.preExponentialFactor() << tab + << std::dynamic_pointer_cast(r)->rate.temperatureExponent() << tab + << std::dynamic_pointer_cast(r)->rate.activationEnergy_R() << endl; break; case Cantera::THREE_BODY_RXN: Info << "THREE_BODY_RXN" << endl; + + Info << std::dynamic_pointer_cast(r)->rate.preExponentialFactor() << tab + << std::dynamic_pointer_cast(r)->rate.temperatureExponent() << tab + << std::dynamic_pointer_cast(r)->rate.activationEnergy_R() << endl; + + for (label l = 0; l < nsp; l++) + { + Info << gas_.speciesName(l) << tab + << std::dynamic_pointer_cast(r)->third_body.efficiency(gas_.speciesName(l)) << endl; + } + break; case Cantera::FALLOFF_RXN: Info << "FALLOFF_RXN" << endl; + + Info << std::dynamic_pointer_cast(r)->low_rate.preExponentialFactor() << tab; + Info << std::dynamic_pointer_cast(r)->low_rate.temperatureExponent() << tab; + Info << std::dynamic_pointer_cast(r)->low_rate.activationEnergy_R() << endl; + + Info << std::dynamic_pointer_cast(r)->high_rate.preExponentialFactor() << tab; + Info << std::dynamic_pointer_cast(r)->high_rate.temperatureExponent() << tab; + Info << std::dynamic_pointer_cast(r)->high_rate.activationEnergy_R() << endl; + + Info << std::dynamic_pointer_cast(r)->falloff->getType() << tab + << std::dynamic_pointer_cast(r)->falloff->nParameters() << endl; + break; case Cantera::PLOG_RXN: