An accurate moving wall boundary algorithm for Direct Simulation of Monte Carlo in unsteady rarefied flow

TL;DR

This paper introduces a precise moving wall boundary algorithm for DSMC simulations of unsteady rarefied flows, improving accuracy in particle-wall collision predictions by accounting for wall motion during collisions.

Contribution

The proposed algorithm accurately models particle collisions with moving walls in DSMC, handling complex geometries and unsteady flows, surpassing previous approximate methods.

Findings

The algorithm aligns with Maxwellian solutions and conserves particles.

It accurately predicts collision dynamics in complex, moving geometries.

Unsteady simulations capture hysteresis effects missed by steady-state models.

Abstract

An accurate algorithm is proposed to improve the prediction of a particle in collision with a moving wall within the direct simulation Monte Carlo (DSMC) framework for the simulation of unsteady rarefied flows. This algorithm is able to predict the particle-wall collision in a coupled manner by removing the assumption employed by the approximate algorithm, in that the wall is frozen during the collision. The trajectory equation of the particle is theoretically constructed in a moving object coordinate system. It can accurately describe the geometries of the collision between a particle and an arbitrarily shaped object of which the motion incorporates both translation and rotation, thus allowing to deal with complex problems. In contrast, the approximate algorithm ignores the effect of the moving wall on the particle movement during the collision, and therefore induces error that is an increasing function of the wall velocity. Four rarefied flow problems are applied to validate the accurate algorithm. It is shown that the algorithm can produce results perfectly consistent with the Maxwellian theoretical solutions and ensure particle conservation to avoid gas leakage. It is also shown in a three-dimensional case of a re-entry module that the steady simulation fails to reproduce the hysteresis effect while the unsteady simulation using the accurate algorithm can do that, indicating that the unsteady simulation with an appropriate algorithm as proposed in the present work is essentially required in such applications.