New Treatment of Pressure Boundary Conditions for DSMC Method in Micro-Channel Flow Simulations

Using DSMC to simulate micro flows in micro-channels, the numerical treatment of boundary conditions is very important. In this paper, several previous numerical treatments of boundary conditions are discussed with their merits and demerits, and a new treatment method based on the assumption of certain pressure distribution in the cells for boundary conditions is proposed. As comparable validity tests, it is applied in the DSMC simulations for the Poiseuille micro flows in micro-channels with four types of classical pressure boundary conditions. The dimensionless velocity profiles are shown and compared with analytical solutions derived from the Navier-Stokes equations with slip boundary conditions. The pressure distributions along the centerline of the micro-channel with the different boundary conditions are presented, and the simulation solutions agree well with the slip analytical solutions. As the Knudsen number increased, a strong linearity of the pressure distribution can be evidently predicted by the new method. Compared with the inlet and outlet velocity distribution, it is shown that the new method has better efficiency than the previous methods in the convergence.