OpenFOAM is one of the most widely used open-source computational fluid dynamics tools and often employed for chemical engineering applications. However, there is no systematic assessment of OpenFOAM’s numerical accuracy and parallel performance for chemically reacting flows. For the first time, this work provides a direct comparison between OpenFOAM’s built-in flow solvers as well as its reacting flow extension EBIdnsFoam with four other, well established high-fidelity combustion codes. Quantification of OpenFOAM’s numerical accuracy is achieved with a benchmark suite that has recently been established by Abdelsamie et al. (Comput Fluids 223:104935, 2021. https://doi.org/10.1016/j.compfluid.2021.104935) for combustion codes. Fourth-order convergence can be achieved with OpenFOAM’s own cubic interpolation scheme and excellent agreement with other high-fidelity codes is presented for incompressible flows as well as more complex cases including heat conduction and molecular diffusion in multi-component mixtures. In terms of computational performance, the simulation of incompressible non-reacting flows with OpenFOAM is slower than the other codes, but similar performance is achieved for reacting flows with excellent parallel scalability. For the benchmark case of hydrogen flames interacting with a Taylor–Green vortex, differences between low-Mach and compressible solvers are identified which highlight the need for more investigations into reliable benchmarks for reacting flow solvers. The results from this work provide the first contribution of a fully implicit compressible combustion solver to the benchmark suite and are thus valuable to the combustion community. The OpenFOAM cases are publicly available and serve as guide for achieving the highest numerical accuracy as well as a basis for future developments.

Identifying combustion regimes in terms of premixed and non-premixed characteristics is an important task for understanding combustion phenomena and the structure of flames. A quasi-DNS database of the compositionally inhomogeneous partially premixed Sydney/Sandia flame in configuration FJ-5GP-Lr75-57 is used to directly compare different types of flame regime markers from literature. In the simulation of the flame, detailed chemistry and diffusion models are utilized and no turbulence and combustion models are used as the flame front and flow are fully resolved near the nozzle. This allows evaluating the regime markers as a post-processing step without modeling assumptions and directly comparing regime markers based on gradient alignment, drift term analysis and gradient free regime identification. The goal is not to find the correct regime marker, which might be impossible due to the different set of assumptions of every marker and the generally vague definition of the partially premixed regime itself, but to compare their behavior when applied to a resolved turbulent flame with partially premixed characteristics.