Non-equilibrium parameter for a hybrid Fokker-Planck/DSMC scheme

Abstract

In this work we investigate the kinetic Fokker-Planck (FP) model [1, 2] concerning its applicability to simulate hyper-sonic, rarefied gas flow. As in the direct simulation Monte Carlo (DSMC) algorithm, the FP model simulates the stochastic motion of a set of particles through the computational domain which results in a simple coupling of both methods. No collisions must be calculated in the FP approach and so larger cell and time step sizes than by DSMC calculations can be used. Because of this, the FP model holds the potential to be more efficient than DSMC where the Knudsen number is small. Further we present and investigate a hybridisation scheme for the FP and DSMC algorithms. We derive a non-equilibrium parameter from the Fokker-Planck- and Boltzmann-operator which indicates the validity of the FP model. In hybrid simulations this parameter is used to partition the computational domain in DSMC and FP regions. The particle motion is processed similarly in both regions, only the procedure for assigning new particle velocities differs. To achieve an equal number of particles per cell, weighting factors depended on the cell size are used. The hybridisation scheme is tested for one-dimensional shock and two-dimensional cylinder flows and gives results in agreement to pure DSMC simulations. An efficiency study shows, that the hybrid scheme is roughly five times faster than a pure DSMC simulation for a hypersonic cylinder flow test case.