A variable time step self-consistent mean field DSMC model for three-dimensional environments
Author(s) -
Otto Schullian,
Hanne Antila,
Brianna R. Heazlewood
Publication year - 2022
Publication title -
the journal of chemical physics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.071
H-Index - 357
eISSN - 1089-7690
pISSN - 0021-9606
DOI - 10.1063/5.0083033
Subject(s) - hagen–poiseuille equation , direct simulation monte carlo , mean free path , statistical physics , flow (mathematics) , limit (mathematics) , monte carlo method , range (aeronautics) , mechanics , collision frequency , free molecular flow , work (physics) , physics , computational physics , thermodynamics , materials science , mathematics , plasma , mathematical analysis , dynamic monte carlo method , statistics , optics , quantum mechanics , scattering , composite material
A self-consistent mean field direct simulation Monte Carlo (SCMFD) algorithm was recently proposed for simulating collision environments for a range of one-dimensional model systems. This work extends the one-dimensional SCMFD approach to three dimensions and introduces a variable time step (3D-vt-SCMFD), enabling the modeling of a considerably wider range of different collision environments. We demonstrate the performance of the augmented method by modeling a varied set of test systems: ideal gas mixtures, Poiseuille flow of argon, and expansion of gas into high vacuum. For the gas mixtures, the 3D-vt-SCMFD method reproduces the properties (mean free path, mean free time, collision frequency, and temperature) in excellent agreement with theoretical predictions. From the Poiseuille flow simulations, we extract flow profiles that agree with the solution to the Navier–Stokes equations in the high-density limit and resemble free molecular flow at low densities, as expected. The measured viscosity from 3D-vt-SCMF is ∼15% lower than the theoretical prediction from Chapman–Enskog theory. The expansion of gas into vacuum is examined in the effusive regime and at the hydrodynamic limit. In both cases, 3D-vt-SCMDF simulations produce gas beam density, velocity, and temperature profiles in excellent agreement with analytical models. In summary, our tests show that 3D-vt-SCMFD is robust and computationally efficient, while also illustrating the diversity of systems the SCMFD model can be successfully applied to.
Accelerating Research
Robert Robinson Avenue,
Oxford Science Park, Oxford
OX4 4GP, United Kingdom
Address
John Eccles HouseRobert Robinson Avenue,
Oxford Science Park, Oxford
OX4 4GP, United Kingdom