Fast reactive flow simulations using analytical Jacobian and dynamic load balancing in OpenFOAM

Ilya Morev*, Bulut Tekgül, Mahmoud Gadalla, Ali Shahanaghi, Jeevananthan Kannan, Shervin Karimkashi, Ossi Kaario, Ville Vuorinen

*Corresponding author for this work

Research output: Contribution to journalArticleScientificpeer-review

11 Citations (Scopus)
95 Downloads (Pure)


Detailed chemistry-based computational fluid dynamics (CFD) simulations are computationally expensive due to the solution of the underlying chemical kinetics system of ordinary differential equations (ODEs). Here, we introduce a novel open-source library aiming at speeding up such reactive flow simulations using OpenFOAM, an open-source software for CFD. First, our dynamic load balancing model by Tekgül et al. [“DLBFoam: An open-source dynamic load balancing model for fast reacting flow simulations in OpenFOAM,” Comput. Phys. Commun. 267, 108073 (2021)] is utilized to mitigate the computational imbalance due to chemistry solution in multiprocessor reactive flow simulations. Then, the individual (cell-based) chemistry solutions are optimized by implementing an analytical Jacobian formulation using the open-source library pyJac, and by increasing the efficiency of the ODE solvers by utilizing the standard linear algebra package. We demonstrate the speed-up capabilities of this new library on various combustion problems. These test problems include a two-dimensional (2D) turbulent reacting shear layer and three-dimensional (3D) stratified combustion to highlight the favorable scaling aspects of the library on ignition and flame front initiation setups for dual-fuel combustion. Furthermore, two fundamental 3D demonstrations are provided on non-premixed and partially premixed flames, viz., the Engine Combustion Network Spray A and the Sandia flame D experimental configurations, which were previously considered unfeasible using OpenFOAM. The novel model offers up to two orders of magnitude speed-up for most of the investigated cases. The openly shared code along with the test case setups represent a radically new enabler for reactive flow simulations in the OpenFOAM framework.
Original languageEnglish
Article number021801
Number of pages14
Issue number2
Publication statusPublished - 1 Feb 2022
MoE publication typeA1 Journal article-refereed


Dive into the research topics of 'Fast reactive flow simulations using analytical Jacobian and dynamic load balancing in OpenFOAM'. Together they form a unique fingerprint.

Cite this