Redistribution of the post-reaction internal energies in DSMC using quantum-kinetic model

The Direct Simulation Monte Carlo (DSMC) method has been largely adopted to analyze problems regarding hypersonic, non-equilibrium, and microscopic flows. In this study, we investigate the thermal-chemical effects on combustion at the microscopic scale using this particle collision-based method. It is realized that the existing Larsen-Borgnakke (L-B) model dealing with transfers of various internal energies cannot provide valid solutions for the reactions, and consequently the system fails to reach thermal equilibrium. To overcome this problem, we propose a modified quantum-kinetic (Q-K) model and corresponding redistribution algorithm to satisfy the required detailed balance based on the solver dsmcFoam+ in the open-source software OpenFOAM. This allows a more straightforward way to handle post-energy redistribution in chemical reactions in comparison to those of the other methods, thus reducing the computational cost and manipulation. To verify accuracy, spontaneous combustion of premixed H 2 / O 2 is simulated, which includes polyatomic reactions and non-equilibrium processes. Compared with the L-B redistribution method, substantial improvement and excellent solutions to the issues are demonstrated by using the new approach, paving the way for accurate and efficient studies of complex problems involving polyatomic chemical reactions and non-equilibrium processes. • A new method of DSMC established in OpenFOAM for simulation of flow fields involving thermal non-equilibrium and chemical reactions. • The approach is based on a modified quantum-kinetic model and corresponding redistribution algorithm. • An equilibrium between relative translational and internal energy modes at the final state is guaranteed while detailed balance is satisfied. • Accuracy verified via simulations of spontaneous combustion of premixed H 2 /O 2 , which includes polyatomic reactions and highly non-equilibrium processes. • More straightforward way of handling post-energy redistribution in chemical reactions as compared to other methods, thus reducing largely the computational cost and manipulation.