chemFoam solver gets wrong results when using second order time scheme
Summary
When using the chemFoam solver to to run the corresponding tutorial case, we found that using "Euler" as the ddt scheme leads to correct results, which is the same as the CHEMKIN result (the final temperature is 2660K). While using "backward" or the other second order schemes leads to a much lower final temperature. We are not sure if this issue still exists in reactingFoam, which has been widely used in turbulent combustion simulations by OpenFOAM. This question is very important for us.
Steps to reproduce
-
Run the tutorial case with default settings and plot the results. The case dictionary is: /tutorials/combustion/chemFoam/gri
-
Change the ddt scheme from "Euler" into "backward" and run the case again.
-
Plot the results and compare them with the CHEMKIN result provided in the tutorial case: /tutorials/combustion/chemFoam/gri/chemkin/senk.out
(The OpenFOAM version we use is v2112)
Example case
What is the current bug behaviour?
The final temperature does not agree with the CHEMKIN result for second order ddt schemes.
What is the expected correct behavior?
The final temperature agree with the CHEMKIN result for different ddt schemes.
Relevant logs and/or images
Environment information
- OpenFOAM version :v2112
- Operating system :Linux
- Hardware info :
- Compiler :gcc