A parallel cell‐based DSMC method on unstructured adaptive meshes

A parallel DSMC method based on a cell‐based data structure is developed for the efficient simulation of rarefied gas flows on PC‐clusters. Parallel computation is made by decomposing the computational domain into several subdomains. Dynamic load balancing between processors is achieved based on the number of simulation particles and the number of cells allocated in each subdomain. Adjustment of cell size is also made through mesh adaptation for the improvement of solution accuracy and the efficient usage of meshes. Applications were made for a two‐dimensional supersonic leading‐edge flow, the axi‐symmetric Rothe's nozzle, and the open hollow cylinder flare flow for validation. It was found that the present method is an efficient tool for the simulation of rarefied gas flows on PC‐based parallel machines. Copyright © 2004 John Wiley & Sons, Ltd.     

低速二维流动的粒子仿真研究

本文运用信息保存法对低速二维的流动现象进行模拟，考察了低速条件下的有限平板绕流以及微槽道气体流动问题。研究表明：在对低速流动的模拟过程中，运用IP法在能够获得较好的结果的同时，具有比DSMC方法更高的计算效率。     

双组元姿控发动机反流特性的DSMC模拟

为准确获得双组元姿控发动机的羽流及其反流特性,采用直接模拟蒙特卡罗(DSMC)方法结合可变软球(VSS)模型数值模拟以N2O4/MMH为推进剂的真空羽流。喷口参数由喷管流动的DSMC法计算得到。与文献比较表明本文计算程序可信。模拟结果表明,分子质量越小的气体越容易进入反流区;本文的计算条件下,在偏离X轴正向116°~140°范围内,组分数密度和质量流率以及温度均出现了最大值。     

Beyond direct simulation Monte Carlo (DSMC) modelling of collision environments

ABSTRACT

Direct simulation Monte Carlo (DSMC) models have been successfully adopted and adapted to describe gas flows in a wide range of environments since the method was first introduced by Bird in the 1960s. We propose a new approach to modelling collisions between gas-phase particles in this work – operating in a similar way to the DSMC model, but with one key difference. Particles move in a mean field, generated by all previously propagated particles, which removes the requirement that all particles be propagated simultaneously. This yields a significant reduction in computation effort and lends itself to applications for which DSMC becomes intractable, such as when a species of interest is only a minor component of a large gas mixture.     

DSMC方法下给定热流边界条件的实施

本文详细讨论了DSMC方法中流体温度、能量及边界热流的统计方法,发展了一种从边界热流求得与壁面碰撞分子反射速度的方法。该方法被称为逆温度抽样算法(ITS,Inverse Temperature Sampling)方法。在此基础上,本文发展了DSMC方法中壁面处给定热流边界条件的实施方法。计算结果表明: ITS方法能准确抽样反射分子的特征温度,进而求得分子反射速度。基于该方法的给定热流边界条件可以准确求得壁面处温度分布,以及流场内的压力、速度。     

A hybrid DSMC/Navier–Stokes frame to solve mixed rarefied/nonrarefied hypersonic flows over nano‐plate and micro‐cylinder

We extend a hybrid DSMC/Navier–Stokes (NS) approach to unify the DSMC and the NS simulators in one framework capable of solving the mixed non‐equilibrium and near‐equilibrium flow regions efficiently. Furthermore, we use a one‐way state‐based coupling (Dirichlet–Dirichlet boundary‐condition coupling) to transfer the required information from the continuum region to the rarefied one. The current hybrid DSMC–NS frame is applied to the hypersonic flows over nanoflat plate and microcylinder cases. The achieved solutions are compared with the pure DSMC and NS solutions. The results show that the current hybrid approach predicts the surface heat transfer rate and shear stress magnitudes very accurately. Some important conclusions can be drawn from this study. For example, although the shock wave region would be a non‐equilibrium region, it is not necessary to use a pure DSMC simulator to solve it entirely. This is important when the researchers wish to predict the surface properties such as velocity slip, temperature jump, wall heat flux rate, and friction drag magnitudes accurately. Our investigation showed that our hybrid solution time would be at least 40% (for the flat plate) and 35% (for the cylinder) of the time that must be spent by a pure DSMC solver to attain the same accuracy.Copyright © 2013 John Wiley & Sons, Ltd.     

DSMC approach to nonstationary Mach reflection of strong incoming shock waves using a smoothing technique

The direct simulation Monte Carlo (DSMC) method is applied to simulation of nonstationary Mach reflection of strong shock waves. Normally the DSMC method is very time consuming in solving unsteady flow field problems especially for high Mach numbers, because of the necessity of iterative calculations to eliminate the inherent statistical fluctuation caused by a finite sample size. A central weighted smoothing technique is introduced to process the DSMC results, so that the iteration time can be significantly reduced. In spite of some relaxations of the shock wave structure, the smoothing technique is verified to be useful to estima te the flow fields qualitatively and even quantitatively by using a relatively small sample size. The comparison between the present approach and the kineticmodel approach (Xu et al. 1991a, 1991b) on the application to unsteady rarefied flow fields was also carried out.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

高速流场气动光学RANS/DSMC混合数值模拟算法研究

提出了适用于高速流场气动光学数值模拟的RANS/DSMC混合算法.通过RANS对全局时均流场进行数值模拟,再对其中局部流场的脉动量采用DSMC进行数值模拟,以Maxwellian速度分布实现RANS宏观量信息向DSMC微观量信息的传递.采用超声速环境下尖劈模型对该混合算法进行校验,对比试验结果论证了算法的正确性.在计算...     

Searching for a near neighbor particle in DSMC cells using pseudo-subcells

LeBeau et al. (2003) [4] introduced the ‘virtual-subcell’ (VSC) method of finding a collision partner for a given DSMC particle in a cell; all potential collision partners in the cell are examined to find the nearest neighbor, which becomes the collision partner. Here I propose a modification of the VSC method, the ‘pseudo-subcell’ (PSC) method, whereby the search for a collision partner stops whenever a ‘near-enough’ particle is found, i.e. whenever another particle is found within the ‘pseudo-subcell’ of radius δ centered on the first particle. The radius of the pseudo-subcell is given by δ = Fd_n, where d_n is the expected distance to the nearest neighbor and F is a constant which can be adjusted to give a desired trade-off between CPU time and accuracy as measured by a small mean collision separation (MCS). For 3D orthogonal cells, of various aspect ratios, d_n/L ≈ 0.746/N^0.383 where N is the number of particles in the cell and L is the cube root of the cell volume. There is a good chance that a particle will be found in the pseudo-subcell and there is a good chance that such a particle is in fact the nearest neighbor. If no particle is found within the pseudo-subcell the closest particle becomes the collision partner.     

微尺度气体滑动轴承的Monte Carlo模拟与性能分析

改进了直接模拟Monte Carlo (DSMC)方法并模拟研究了轴承构型、滑动速度、壁温及环境压力对微气体滑动轴承内部压力分布及承载能力的影响.结果表明:轴承的几何构型和滑动速度对轴承的性能影响很大,对于相同长度的轴承,气体的峰值压力与轴承的承载能力随轴承出口尺寸的减小和滑动速度的增大而增大;当轴承的几何形状和滑动速度固定时,通道壁温和环境压力是影响轴承性能的重要因素,壁温越高,轴承的承载能力越强;环境压力不同,轴承性能亦有所不同.