A second-order particle Fokker-Planck model for rarefied gas flows

The direct simulation Monte Carlo (DSMC) method has become a powerful tool for studying rarefied gas flows. However, for the DSMC method to be effective, the cell size must be smaller than the mean free path, and the time step smaller than the mean collision time. These constraints make it difficult to use the DSMC method in multiscale rarefied gas flows. Over the past decade, the particle Fokker-Planck (FP) method has been studied to address computational cost issues in the near-continuum regime. To capture the main features of the Boltzmann equation, various FP models have been proposed, such as the quadratic entropic FP (Quad-EFP) and the ellipsoidal statistical FP (ESFP). Nevertheless, few studies have clearly demonstrated that the FP method offers a computational advantage over the DSMC method without sacrificing accuracy. This is because conventional particle FP methods have employed first-order accuracy schemes. The present study proposes a unified stochastic particle ESFP (USP-ESFP) model. This model improves the accuracy of shear stress and heat flux predictions. Additionally, a spatial interpolation scheme is introduced to the particle FP method. The numerical test cases include relaxation problem, Couette flows, Poiseuille flows, velocity perturbation, and hypersonic flows around a cylinder. The results show that the USP-ESFP model agrees well with both analytical and DSMC results. Furthermore, the USP-ESFP model is found to be less sensitive to cell size and time step than the DSMC method, resulting in a factor of four speed-up for the considered hypersonic flow around a cylinder.