A dynamic load balancing model coupled with DAC and ISAT for a stochastic turbulent combustion model

Due to the composition-dependent stiffness of chemistry, simulations of reactive turbulent flows may present computational load imbalance among parallel processes when spatial decomposition is used for parallelisation, causing high CPU idle time and waste of computational resources. To increase computational efficiency, a dynamic load balancing (DLB) model is proposed to redistribute computational load among computing cores. The DLB model exploits a decomposition in the mixture fraction space with two dynamic adjusting decomposition strategies to realise load redistribution. The DLB model is suitable for the integration of chemistry on stochastic particles in hybrid Eulerian/Lagrangian turbulent combustion models in which the Eulerian field is conventionally decomposed statically in physical space in a way that balances the computational load for the solution of the Navier-Stokes equation but which does not generally lead to balanced load for the computation of the composition fields. Here it is tested using an OpenFOAM-based platform, mmcFoam, which is a comprehensive object-orientated C++ library for stochastic turbulent combustion modelling. Apart from direct integration (DI) for chemistry, the DLB model is also coupled with dynamic adaptive chemistry (DAC) and in situ adaptive tabulation (ISAT), which allows for extra speedup. The performance of the coupled models is validated and assessed for two laboratory flame conditions that exhibit different levels of computational load imbalance. Overall, the DLB model effectively balances the computational load distribution and increases the effective usage of computing power, shortening the simulation wall time required. Moreover, a strong scaling test is carried out using up to 512 cores. Although all approaches have sub-ideal scalability, the scalability of each with DLB is significantly better than without DLB. While DLB-ISAT has relatively poor scalability compared to the DI- and DAC-based methods, DLB-ISAT still ranks the fastest among the algorithms in all scaling trials