“…There are different types of continuum breakdown parameter used in the literatures. 7,24,[45][46][47][48][49] Most of the parameters proposed and used in these works are based on the calculation of flow gradients. In a statistical method, calculation of gradients is always subject to noise, which can lead to significant computational errors.…”
Summary
A new 2D parallel multispecies polyatomic particle–based hybrid flow solver is developed by coupling the Direct Simulation Monte Carlo (DSMC) method with a novel Dynamic Collision Limiter (DCL) approach to solve multiscale transitional flows. The hybrid DSMC‐DCL solver can solve nonequilibrium multiscale flows with length scales ranging from continuum to rarefied. The DCL method, developed in this work, dynamically assigns different number of collisions in cells, which is based on the local value of K‐S parameter such that the number of collisions per time step is limited in near‐equilibrium flow regions. Present hybrid solver uses the Kolmogorov‐Smirnov statistical test as the continuum breakdown parameter, based on which, the solution domain is decomposed into near‐equilibrium and nonequilibrium flow regions. Direct Simulation Monte Carlo is used where nonequilibrium flow regions are encountered, while the DCL method is used where flow regions are found to be in near‐equilibrium state. In this work, we have studied hypersonic flow of nitrogen over a blunt body with an aerospike and supersonic flow of argon through a micronozzle. The results obtained by the hybrid DSMC‐DCL solver are compared and shown to agree well with the experimental data and with those obtained from DSMC, with significant savings in the computational cost.
“…There are different types of continuum breakdown parameter used in the literatures. 7,24,[45][46][47][48][49] Most of the parameters proposed and used in these works are based on the calculation of flow gradients. In a statistical method, calculation of gradients is always subject to noise, which can lead to significant computational errors.…”
Summary
A new 2D parallel multispecies polyatomic particle–based hybrid flow solver is developed by coupling the Direct Simulation Monte Carlo (DSMC) method with a novel Dynamic Collision Limiter (DCL) approach to solve multiscale transitional flows. The hybrid DSMC‐DCL solver can solve nonequilibrium multiscale flows with length scales ranging from continuum to rarefied. The DCL method, developed in this work, dynamically assigns different number of collisions in cells, which is based on the local value of K‐S parameter such that the number of collisions per time step is limited in near‐equilibrium flow regions. Present hybrid solver uses the Kolmogorov‐Smirnov statistical test as the continuum breakdown parameter, based on which, the solution domain is decomposed into near‐equilibrium and nonequilibrium flow regions. Direct Simulation Monte Carlo is used where nonequilibrium flow regions are encountered, while the DCL method is used where flow regions are found to be in near‐equilibrium state. In this work, we have studied hypersonic flow of nitrogen over a blunt body with an aerospike and supersonic flow of argon through a micronozzle. The results obtained by the hybrid DSMC‐DCL solver are compared and shown to agree well with the experimental data and with those obtained from DSMC, with significant savings in the computational cost.
“…To reduce them, the initial stages for some regimes (with a dense initial cloud, e.g. Θ > 10 for b = 10) were calculated using the 'equilibrium' modification of the DSMC method [46].…”
A theoretical study of the time-of-flight (TOF) distributions of neutral species produced by low-fluence pulsed laser evaporation in a vacuum has been performed. A database of TOF distributions has been calculated by the direct simulation Monte Carlo method. The calculated TOF signals have been fitted by a modification of the shifted Maxwell-Boltzmann distribution function with the only parameter being the shift velocity. The dependence of the shift velocity on the number of evaporated monolayers has been described by an analytical expression, derived based on the gas-dynamic analysis of pulsed gas expansion into a vacuum. This expression allows the derivation of a new formula for TOF distributions for neutral particles with the only parameter being the surface temperature. This new formula has been used for the analysis of experimental data on pulsed laser ablation of niobium, copper, graphite and gold. The evaluated surface temperature agrees well with the results of the thermal model calculations and available experimental data with an error of 10% instead of 50-150% for commonly used formulas.
“…Titov and Levin [6] found, in collision-limited DSMC, that two collisions per time step per particle are sufficient for the computed nonequilibrium distribution to relax to one differing negligibly from the corresponding Maxwell-Boltzmann equilibrium distribution. This distribution is valid for a gas in local thermal equilibrium, and is derived assuming an infinite time has elapsed for relaxation, and that no spatial gradients in density, bulk velocity or temperature exist to disturb the equilibrium.…”
Section: Validity Of the Assumption Of Infinite Collision Ratementioning
The Quiet Direct Simulation (QDS) scheme is a numerical method for modelling gas flows, based on kinetic theory, with some similarities to the Lattice Boltzmann Method (LBM). It differs from LBM notably in that the discrete molecular velocities are not constant but are reset each timestep according to local values of bulk velocity and temperature. For this reason it performs well in highly compressible flows. Two features of the scheme limit its accuracy in low Mach number flows. QDS assumes a Maxwell distribution of molecular velocities. The validity of this assumption may be tested by calculating the gradient Knudsen number and average number of collisions per timestep. The separation of collision and streaming leads to excessive diffusion of momentum, leading to a very high effective viscosity of the modelled gas when the grid spacing is larger than the mean free path. This numerical dissipation is different in character from the dissipation due to the finite order of the spatial reconstruction, common to all finite volume methods, which is also present. The effective viscosity is quantified for simple shear flows and tested in models of a 2D channel flow. A crude model of intermolecular collision during streaming is implemented and shown to reduce the effective viscosity.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.