View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0001012 | OpenFOAM | Bug | public | 2013-09-13 14:03 | 2013-09-13 15:21 |
| Reporter | Assigned To | henry | |||
| Priority | normal | Severity | major | Reproducibility | always |
| Status | resolved | Resolution | fixed | ||
| Platform | Linux | OS | Ubuntu | OS Version | 9.10 |
| Summary | 0001012: powerSeriesReactionRate inconsistent with Chemkin definition | ||||
| Description | The 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 Reproduce | Change 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 Information | See equation 80 on page 36 of the Chemkin manual (attached). | ||||
| Tags | No tags attached. | ||||
|
2013-09-13 14:03
|
|
|
|
Thanks for the report. Resolved by commit 13b7ee65d421a08cb6f8ab7f1070a696a8fb8d0e |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2013-09-13 14:03 |
|
New Issue | |
| 2013-09-13 14:03 |
|
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 |
