On the conservative property of particle-based Fokker-Planck method for rarefied gas flows

The Fokker-Planck-type approximation of the full Boltzmann equation has aroused intense research interest due to its potential for the stochastic particle simulation of rarefied gas flows. The ellipsoidal statistical Fokker-Planck (ES-FP) model treats the evolution of molecular velocity as a continuous stochastic process, and it satisfies the basic requirements for a proper gas-kinetic model including the H-theorem and an adjustable Prandtl number. The ES-FP model can be numerically implemented with computational particles in a Monte Carlo manner. Two different particle ES-FP schemes are presented. The first scheme utilizes the exact stochastic integral solution of the Langevin equations corresponding to the ES-FP equation and couples free-molecular moves and intermolecular collisions. The second scheme is designed to enforce the conservation of momentum and energy during the numerical simulation based on the decoupled algorithm and the analysis of the specific conditions for the conservative property. Numerical tests are conducted to demonstrate the performances of different schemes. In the simulation of a homogeneous gas system, the ES-FP scheme without enforcement of conservation leads to unphysical variation in the momentum and loss in energy, whereas the conservative ES-FP scheme strictly maintains the momentum and energy of the system. For the Mach 6 flows over the leading edge of a flat plate and over a round-nosed blunt body, the non-conservative ES-FP scheme underestimates the shock angle and the shock standoff distance, makes inaccurate predictions of aerodynamic force and heating, and produces low-temperature anomalies in front of the shock waves. In comparison with the results given by the direct simulation Monte Carlo method, the results of the conservative ES-FP simulations show satisfactory accuracy for the flow fields as well as the distributions of pressure, friction, and heat flux on the wall surfaces.