View Issue Details

IDProjectCategoryView StatusLast Update
0001012OpenFOAMBugpublic2013-09-13 15:21
Reporteruser503Assigned Tohenry  
PrioritynormalSeveritymajorReproducibilityalways
Status resolvedResolutionfixed 
PlatformLinuxOSUbuntuOS Version9.10
Summary0001012: powerSeriesReactionRate inconsistent with Chemkin definition
DescriptionThe calculation of reaction rate in powerSeriesReactionRate is inconsistent with the Chemkin power series reaction rate definition, creating errors when reading from a Chemkin reaction input file.

The Chemkin power series reaction exponential argument is b0/T + b1/T^2 + b2/T^3 + b3/T^4 while the OpenFOAM implementation indexes the power from 0 giving b0 + b1/T + b2/T^2 + b3/T^3.

Changing line 91 of powerSeriesReactionRateI.H from 'expArg += coeffs_[n]/pow(T,n);' to 'expArg += coeffs_[n]/pow(T,n+1);' should fix it. This is the same in 2.1.x and 2.2.x.
Steps To ReproduceChange the Arrhenius reactions in one of the tutorials to powerSeries setting b0 to -Ta and b1-3 to 0 (which should be the same).
Additional InformationSee equation 80 on page 36 of the Chemkin manual (attached).
TagsNo tags attached.

Activities

user503

2013-09-13 14:03

 

chemkinCK.pdf (594,154 bytes)

henry

2013-09-13 15:21

manager   ~0002491

Thanks for the report.

Resolved by commit 13b7ee65d421a08cb6f8ab7f1070a696a8fb8d0e

Issue History

Date Modified Username Field Change
2013-09-13 14:03 user503 New Issue
2013-09-13 14:03 user503 File Added: chemkinCK.pdf
2013-09-13 15:21 henry Note Added: 0002491
2013-09-13 15:21 henry Status new => resolved
2013-09-13 15:21 henry Resolution open => fixed
2013-09-13 15:21 henry Assigned To => henry