chemFoam solver gets wrong results when using second order time scheme
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)
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
- OpenFOAM version :v2112
- Operating system :Linux
- Hardware info :
- Compiler :gcc