IP Simulation of Micro Gas Flows under 3-D Head Sliders

Gas film lubrication of a three-dimensional flat read-write head slider is calculated using the information preservation (IP) method and the direct simulation Monte Carlo (DSMC) method, respectively. The pressure distributions on the head slider surface at different velocities and flying heights obtained by the two methods are in excellent agreement. IP method is also employed to deal with head slider with three-dimensional complex configuration. The pressure distribution on the head slider surface and the net lifting force obtained by the IP method also agree well with those of DSMC method. Much less (of the order about 10(2) less) computational time (the sum of the time used to reach a steady stage and the time used in sampling process) is needed by the IP method than the DSMC method and such an advantage is more remarkable as the gas velocity decreases.

Simulation of gas flow over a micro cylinder

Gas flow over a micro cylinder is simulated using both a compressible Navier-Stokes solver and a hybrid continuum /particle approach. The micro cylinder flow has low Reynolds number because of the small length scale and the low speed, which also indicates that the rarefied gas effect exists in the flow. A cylinder having a diameter of 20 microns is simulated under several flow conditions where the Reynolds number ranges from 2 to 50 and the Mach number varies from 0.1 to 0.8. It is found that the low Reynolds number flow can be compressible even when the Mach number is less than 0.3, and the drag coefficient of the cylinder increases when the Reynolds number decreases. The compressible effect will increase the pressure drag coefficient although the friction coefficient remains nearly unchanged. The rarefied gas effect will reduce both the friction and pressure drag coefficients, and the vortex in the flow may be shrunk or even disappear.

Development of an efficient particle approach for micro-scale gas flow simulations

The micro-scale gas flows are usually low-speed flows and exhibit rarefied gas effects. It is challenging to simulate these flows because traditional CFD method is unable to capture the rarefied gas effects and the direct simulation Monte Carlo (DSMC) method is very inefficient for low-speed flows. In this study we combine two techniques to improve the efficiency of the DSMC method. The information preservation technique is used to reduce the statistical noise and the cell-size relaxed technique is employed to increase the effective cell size. The new cell-size relaxed IP method is found capable of simulating micro-scale gas flows as shown by the 2D lid-driven cavity flows.

Statistical analysis of microscale gas flows and stock price changes

First, recent studies on the information preservation (IP) method, a particle approach for low-speed micro-scale gas flows, are reviewed. The IP method was validated for benchmark issues such as Couette, Poiseuille and Rayleigh flows, compared well with measured data for typical internal flows through micro-channels and external flows past micro flat plates, and combined with the Navier-Stokes equations to be a hybrid scheme for subsonic, rarefied gas flows. Second, the focus is moved to the microscopic characteristic of China stock market, particularly the price correlation between stock deals. A very interesting phenomenon was found that showed a reverse transition behaviour between two neighbouring price changes. This behaviour significantly differs from the transition rules for atomic and molecular energy levels, and it is very helpful to understand the essential difference between stock markets and nature.

Kinetic Studies of Gas Flows

Our recent studies on kinetic behaviors of gas flows are reviewed in this paper. These flows have a wide range of background, but share a common feature that the flow Knudsen number is larger than 0.01. Thus kinetic approaches such as the direct simulation Monte Carlo method are required for their description. In the past few years, we studied several micro/nano-scale flows by developing novel particle simulation approach, and investigated the flows in low-pressure chambers and at high altitude. In addition, the microscopic behaviors of a couple of classical flow problems were analyzed, which shows the potential for kinetic approaches to reveal the microscopic mechanism of gas flows.

3D Direct Simulation Monte Carlo Modelling of the Inner Gas Coma of Comet 67P/Churyumov-Gerasimenko: A Parameter Study

Direct Simulation Monte Carlo (DSMC) is a powerful numerical method to study rarefied gas flows such as cometary comae and has been used by several authors over the past decade to study cometary outflow. However, the investigation of the parameter space in simulations can be time consuming since 3D DSMC is computationally highly intensive. For the target of ESA's Rosetta mission, comet 67P/Churyumov-Gerasimenko, we have identified to what extent modification of several parameters influence the 3D flow and gas temperature fields and have attempted to establish the reliability of inferences about the initial conditions from in situ and remote sensing measurements. A large number of DSMC runs have been completed with varying input parameters. In this work, we present the simulation results and conclude on the sensitivity of solutions to certain inputs. It is found that among cases of water outgassing, the surface production rate distribution is the most influential variable to the flow field.

Couette Gas Flow Between Cylinders with Different Temperature

Петър Господинов, Добри Данков, Владимир Русинов, Стефан Стефанов - Изследвано е стационарно течение на Кует на разреден газ в случая на въртене на вътрешния цилиндър и неподвижен външен цилиндър чрез използване на DSMC метод и числено решение на уравненията на Навие–Стокс за относително малка (дозвукова) скорост на въртене. Изследвани са различни случаи при промяна на температурата на въртящият се цилиндър и числото на Кнудсен. Целта на изследването е да се установи влиянието на малки скорости на въртене върху макрохарактеристиките – плътността, скоростта и температурата на газа. Установено е добро съвпадение на резултатите получени по двата метода за Kn = 0.02. Получените резултати са важни при решаването на неравнинни, задачи от микрофлуидиката с отчитане на ефектите на кривината. Ключови думи: механика на флуидите, кинетична теория, разреден газ, DSMC.

Non-linear couette flow with heat transfer using the B-G-K model

In recent years a large number of investigators have devoted their efforts to the study of flow and heat transfer in rarefied gases, using the BGK [1] model or the Boltzmann kinetic equation. The velocity moment method which is based on an expansion of the distribution function as a series of orthogonal polynomials in velocity space, has been applied to the linearized problem of shear flow and heat transfer by Mott-Smith [2] and Wang Chang and Uhlenbeck [3]. Gross, Jackson and Ziering [4] have improved greatly upon this technique by expressing the distribution function in terms of half-range functions and it is this feature which leads to the rapid convergence of the method. The full-range moments method [4] has been modified by Bhatnagar [5] and then applied to plane Couette flow using the B-G-K model. Bhatnagar and Srivastava [6] have also studied the heat transfer in plane Couette flow using the linearized B-G-K equation. On the other hand, the half-range moments method has been applied by Gross and Ziering [7] to heat transfer between parallel plates using Boltzmann equation for hard sphere molecules and by Ziering [83 to shear and heat flow using Maxwell molecular model. Along different lines, a moment method has been applied by Lees and Liu [9] to heat transfer in Couette flow using Maxwell's transfer equation rather than the Boltzmann equation for distribution function. An iteration method has been developed by Willis [10] to apply it to non-linear heat transfer problems using the B-G-K model, with the zeroth iteration being taken as the solution of the collisionless kinetic equation. Krook [11] has also used the moment method to formulate the equivalent continuum equations and has pointed out that if the effects of molecular collisions are described by the B-G-K model, exact numerical solutions of many rarefied gas-dynamic problems can be obtained. Recently, these numerical solutions have been obtained by Anderson [12] for the non-linear heat transfer in Couette flow,

Direct numerical test of the B-G-K model equation by the DSMC method

In the present paper the rarefied gas how caused by the sudden change of the wall temperature and the Rayleigh problem are simulated by the DSMC method which has been validated by experiments both in global flour field and velocity distribution function level. The comparison of the simulated results with the accurate numerical solutions of the B-G-K model equation shows that near equilibrium the BG-K equation with corrected collision frequency can give accurate result but as farther away from equilibrium the B-G-K equation is not accurate. This is for the first time that the error caused by the B-G-K model equation has been revealed.

Examination of the LBM in simulation of microchannel flow in transitional regime

The gas flows in micro-electro-mechanical systems possess relatively large Knudsen number and usually belong to the slip flow and transitional flow regimes. Recently the lattice Boltzmann method (LBM) was proposed by Nie et al. in Journal of Statistical Physics, vol. 107, pp. 279-289, in 2002 to simulate the microchannel and microcavity flows in the transitional flow regime. The present article intends to test the feasibility of doing so. The results of using the lattice Boltzmann method and the direct simulation Monte Carlo method show good agreement between them for small Kn (Kn = 0.0194), poor agreement for Kn = 0.194, and large deviation for Kn = 0.388 in simulating microchannel flows. This suggests that the present version of the lattice Boltzmann method is not feasible to simulate the transitional channel flow.

矩形微槽道气体流动的速度分布

矩形微槽道的各个流向截面可以局部近似为平面Poiseuille流动，应用信息保存（IP）方法和直接模拟Monte Carlo（DSMC）方法计算了从连续介质区到自由分子流区的平面Poiseuille流动，利用其结果对Beskok-Karniadakis公式和质量流率动理论因子进行修正和重新拟合，给出在整个稀薄气体流动领域都适用的微槽道气体流动速度分布.

认识稀薄气体动力学

以通俗易懂的方式介绍了空气动力学当气体间断分子效应显著时发展起来的特殊分支－－稀薄气体动力学。讨论了非平衡现象与稀薄气体动力学的。通过与8速度气体模型的间断Boltzmann方程的对比，解释了Boltzmann方程碰撞项的物理意义和数学困难，简要综述了其一般解法。讨论了分在物体表面的反射和问题的边界条件，着重介绍了直接模拟Monte Carlo (DSMC)方法和为克服低速稀薄流动（如MEMS中流动）中模拟困难的信息保存（IP）方法。

第24届国际稀薄气体动力学会议简介

第24届国际稀薄气体动力学会议于2004年7月11日至16日在意大利巴里市召开.会议主席为巴里大学化学系主任马里奥?卡彼特利教授.

Statistical Simulation of Thin-Film Bearings

The information preservation (IP) method and the direct simulation Monte Carlo (DSMC) method are used to simulate the gas flows between the write/read head and the platter of the disk drive (the slider bearing problem). The results of both methods are in good agreement with numerical solution of the Reynolds equation in the cases studied. However, the DSMC method owing to the problem of large sample size demand and the difficulty in regulating boundary conditions at the inlet and outlet was able to simulate only short bearings, while IP simulates the bearing of authentic length ~1000 m ? and can provide more detailed flow information.

Information Preservation (IP) Method in Simulation

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.