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

• Published in: Physics of fluids (1994) 2022-02, Vol.34 (2)
• Main Authors: Morev, Ilya, Tekgül, Bulut, Gadalla, Mahmoud, Shahanaghi, Ali, Kannan, Jeevananthan, Karimkashi, Shervin, Kaario, Ossi, Vuorinen, Ville
• Format: Article
• Language: English
• Subjects: Chemistry; Computational fluid dynamics; Differential equations; Dynamic loads; Flame propagation; Flow simulation; Fluid dynamics; Fluid flow; Fuel combustion; Libraries; Linear algebra; Load balancing; Mathematical models; Multiprocessing; Ordinary differential equations; Physics; Premixed flames; Reacting flow; Reaction kinetics; Shear layers; Simulation
• ISSN: 1070-6631; 1089-7666
• DOI: 10.1063/5.0077437

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.