ReactionMechanismGenerator / ReactionMechanismGenerator/RMG-database
Some MultiPDepArrhenius rates in libraries don't validate in Cantera
Nobody has claimed this yet.
- Dominant language
- Python
- Stars
- 127
- Forks
- 160
- Avg merge
- 6d 2h
- Merged PRs (30d)
- 2
Description
The NOx2018 library index = 258, label = "C2H3 + CH2O <=> C2H4 + HCO", has a MultiPDepArrhenius rate expression that comes from (and generates) DUPLICATE sets of PLOG rates in Chemkin format. When translated to Cantera, it chokes with
***********************************************************************
CanteraError thrown by Plog::validate:
Invalid rate coefficient for reaction 'C2H3 + CH2O <=> C2H4 + HCO'
at P = 101.325, T = 200
***********************************************************************
The automatic validation checks an array of T and P, and reveals that at 200K the net rate (the sum of the two expressions) is negative.
The source of the NOx2018 library is
Peter Glarborg, James A. Miller, Branko Ruscic, Stephen J. Klippenstein, Modeling nitrogen chemistry in combustion, Progress in Energy and Combustion Science, 67, 31-68, (2018)
https://doi.org/10.1016/j.pecs.2018.01.002.
The chemkin file contains a comment pointing to a @cfgoldsmith paper on the second (negative) term. This term is exactly the same expression as the following reaction which has the same reactants but 3 products C2H4+H+CO. I guess the trick is to subtract the "hot" route direct to H+CO from the HCO channel and add in the stand-alone reaction. Makes sense. But at low temperatures, this expression is faster than the original, and the net rate ends up negative.
The original they used appears in several other models, including AramcoMech 2, though I'm not sure who used it first.
C2H3+CH2O=C2H4+HCO 1.4E10 0.209 3934 ! 1 atm
PLOG/ 0.001 1.1E07 1.090 1807/
PLOG/ 0.01 2.5E07 0.993 1995/
PLOG/ 0.1 2.5E08 0.704 2596/
PLOG/ 1 1.4E10 0.209 3934/
PLOG/ 10 3.5E13 -0.726 6944/
PLOG/ 100 3.3E14 -0.866 10966/
PLOG/ 1000 1.7E01 3.170 9400/
DUPLICATE
C2H3+CH2O=C2H4+HCO -8.5E19 -2.264 23862 ! 1 atm
PLOG/ 1.000E-03 -2.3E16 -1.269 20617 /
PLOG/ 1.000E-02 -5.2E16 -1.366 20805 /
PLOG/ 1.000E-01 -1.5E18 -1.769 22524 /
PLOG/ 1.000E00 -8.5E19 -2.264 23862 /
PLOG/ 1.000E01 -4.4E23 -3.278 27795 /
PLOG/ 1.000E02 -4.2E24 -3.418 31817 /
PLOG/ 1.000E03 -2.1E11 0.618 30251 /
DUPLICATE
! C.F. Goldsmith et al., J. Phys. Chem. A 116 (2012) 3325�3346.
!
C2H3+CH2O=C2H4+H+CO 8.5E19 -2.264 23862 ! 1 atm
PLOG/ 1.000E-03 2.3E16 -1.269 20617 /
PLOG/ 1.000E-02 5.2E16 -1.366 20805 /
PLOG/ 1.000E-01 1.5E18 -1.769 22524 /
PLOG/ 1.000E00 8.5E19 -2.264 23862 /
PLOG/ 1.000E01 4.4E23 -3.278 27795 /
PLOG/ 1.000E02 4.2E24 -3.418 31817 /
PLOG/ 1.000E03 2.1E11 0.618 30251 /
Any suggestions, @cfgoldsmith? I can think of...
- just revert to the AramcoMech rate and not use the hot route.
- change Cantera to not do these sanity checks? (or not do them at 200K)
- ...?
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start by reproducing the NOx2018 library index 258 translation and Cantera validation at 200 K and 101.325 kPa. Inspect the duplicate PLOG rate sets and the related C2H3 + CH2O = C2H4 + H + CO reaction described in the issue. Done means an agreed correction or validation policy is implemented and the translated rates no longer fail validation.
Written by the indexing model from the issue text.
Assessment
- Domain
- databases
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Mostly clear
- Newbie friendliness
- 35/100