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

Abstract: 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 C++ software for CFD. First, our dynamic load balancing model DLBFoam (Tekg\"ul et al., 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 linear algebra package LAPACK. We demonstrate the speed-up capabilities of this new library on various combustion problems. These test problems include a 2D turbulent reacting shear layer and 3D stratified combustion to highlight the favorable scaling aspects of the library on ignition/flame front initiation setups for dual-fuel combustion. Furthermore, two fundamental 3D demonstrations are provided on non-premixed and partially premixed flames, namely the ECN Spray A and the Sandia flame D experimental configurations. 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.