Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

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 m ? 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 ( 1000 L m ? = ) 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 ( 5 L m ? = ) 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.

Rarefied Gas Flows in Micro-Channels

利用基于分子模型的统计模拟方法--信息保存方法(IP)统计模拟了实验条件下微槽道气体流动,仔细讨论了用IP方法模拟长槽道稀薄气流时遇到的问题,并给出了解决的方法,即采取守恒形式的控制方程避免质量流量计算误差积累,并利用超松弛方法使收敛过程加速。将IP计算结果与压力分布和质量流量实验数据进行了比较。

Application of Electrical Resistance Tomography System to Monitor Gas/Liquid Two-Phase Flow in a Horizontal Pipe

The Electrical Resistance Tomography (ERT) technique possesses great potential in monitoring widely exiting industrial two/multi-phase flow. For vertical pipe flow and inclined pipe flow, some application studies with exciting results have been reported, but there is rarely a paper regarding the application of ERT to horizontal gas/liquid pipe flow. This paper addresses this issue and proposes a smart method, Liquid Level Detection method, to conventional ERT system. The enhanced ERT system using the new method can monitor horizontal pipe flow effectively and its application is no longer restricted by the flow conditions. Some experimental results from monitoring an air/water slug pipe flow are presented.

Characteristics of Micronozzle Gas Flows

The fluid characteristics of gas flows in the micronozzle whose throat height is 20 μm were investigated by the direct simulation Monte Carlo (DSMC) method. In a series of cases, the dependence of mass flux on the pressure difference was gained, and the DSMC's results show good agreement with the experimental data. The comparison of mass flux and the Mach number contours between the DSMC and Navier-Stokes equations adding slip boundary also reveals quantitatively that the continuum model will be invalid gradually even when the average Knudsen number is smaller than 0.01. As one focus of the present paper, the phenomenon of the multiple expansion-compression waves that comes from the nozzle's divergent part was analyzed in detailed.

Study of Gas Combustion in a Vented Cylindrical Vessel

Complicated interaction of a flame front with a turbulent flow induced by venting is studied during combustion of the stoichiometric propane/air mixture in a relatively large vented cylindrical vessel. Flame position, its shape, and combustion pressure were measured as a function of time and vent parameters. The experimental data were used to verify numerical simulation of the combustion process. The proposed numerical model satisfactorily simulates the main features of combustion in a closed and vented vessel such as flame configuration, flow and temperature fields, and pressure variation pattern. Simulated velocity and temperature distribution are very useful pieces of information because they are not available from experiments.

Experimental study and simulation principles of an oil-gas multiphase transportation system

Presented is an experimental study on the performance of an oil-gas multiphase transportation system, especially on the multiphase flow patterns, multiphase pumping and multiphase metering of the system. A dynamic simulation analysis is conducted to deduce simulation parameters of the system and similarity criteria under simplified conditions are obtained. The reliability and feasibility of two-phase flow experiment with oil and natural gas simulated by water and air are discussed by using the similarity criteria.

A Generalized Soft-Sphere for Monte Carlo Simulation

A new collision model, called the generalized soft-sphere (GSS) model, is introduced. It has the same total cross section as the generalized hard-sphere model [Phys. Fluids A 5, 738 (1993)], whereas the deflection angle is calculated by the soft-sphere scattering model [Phys. Fluids A 3, 2459 (1991)]. In virtue of a two-term formula given to fit the numerical solutions of the collision integrals for the Lennard-Jones (6-12) potential and for the Stockmayer potential, the parameters involved in the GSS model are determined explicitly that may fully reproduce the transport coefficients under these potentials. Coefficients of viscosity, self-diffusion and diffusion for both polar and nonpolar molecules given by the GSS model and experiment are in excellent agreement over a wide range of temperature from low to high.

Numerical Simulation of Nox Formation in Coal Combustion with Inlet Natural Gas Burning

A full two-fluid model of reacting gas-particle flows and coal combustion is used to simulate coal combustion with and without inlet natural gas added in the inlet. The simulation results for the case without natural gas burning is in fair agreement with the experimental results reported in references. The simulation results of different natural gas adding positions indicate that the natural gas burning can form lean oxygen combustion enviroment at the combustor inlet region and the NOz concentration is reduced. The same result can be obtained from chemical equilibrium analysis.

Gas exchange by bubbles in waves

A simple and feasible model feet the calculation of the gas transfer by bubble clouds is proposed in this article. N-2, O-2, and CO2 transferred by bubble clouds are obtained. At wind speed of 10 m/s, the calculated supersaturation of dissolved oxygen is 1.93-3.89% in agreement with the measurement.

Flow properties of a dusty-gas point source in a supersonic free stream

By using Lagrangian method, the flow properties of a dusty-gas point source in a supersonic free stream were studied and the particle parameters in the near-symmetry-axis region were obtained. It is demonstrated that fairly inertial particles travel along oscillating and intersecting trajectories between the bow and termination shock waves. In this region,formation of "multi-layer structure" in particle distribution with alternating low- and highdensity layers is revealed. Moreover, sharp accumulation of particles occurs near the envelopes of particle trajectories.

Saltation in wind blown sand

The Boltzmann equation of the sand particle velocity distribution function in wind-blown sand two-phase flow is established based on the motion equation of single particle in air. And then, the generalized balance law of particle property in single phase granular flow is extended to gas-particle two-phase flow. The velocity distribution function of particle phase is expanded into an infinite series by means of Grad's method and the Gauss distribution is used to replace Maxwell distribution. In the case of truncation at the third-order terms, a closed third-order moment dynamical equation system is constructed. The theory is further simplified according to the measurement results obtained by stroboscopic photography in wind tunnel tests.

Conjugate Model for Heat and Mass Transfer of Porous Wall in the High Temperature Gas Flow

Heat and mass transfer of a porous permeable wall in a high temperature gas dynamical flow is considered. Numerical simulation is conducted on the ground of the conjugate mathematical model which includes filtration and heat transfer equations in a porous body and boundary layer equations on its surface. Such an approach enables one to take into account complex interaction between heat and mass transfer in the gasdynamical flow and in the structure subjected to this flow. The main attention is given to the impact of the intraporous heat transfer intensity on the transpiration cooling efficiency.

Shock wave diffraction by a square cavity filled with dusty gas

A simple two-dimensional square cavity model is used to study shock attenuating effects of dust suspension in air. The GRP scheme for compressible flows was extended to simulate the fluid dynamics of dilute dust suspensions, employing the conventional two-phase approximation. A planar shock of constant intensity propagated in pure air over Aat ground and diffracted into a square cavity filled with a dusty quiescent suspension. Shock intensities were M-s = 1.30 and M-s = 2.032, dust loading ratios were alpha = 1 and alpha = 5, and particle diameters were d = 1, 10 and 50 mum. It was found that the diffraction patterns in the cavity were decisively attenuated by the dust suspension, particularly for the higher loading ratio. The particle size has a pronounced effect on the flow and wave pattern developed inside the cavity. Wall pressure historics were recorded for each of the three cavity walls, showing a clear attenuating effect of the dust suspension.

Predictions for Partial-Dislocation-Mediated Processes in Nanocrystalline Ni by Generalized Planar Fault Energy Curves: An Experimental Evaluation

Generalized planar fault energy (GPFE) curves have been used to predict partial-dislocation-mediated processes in nanocrystalline materials, but their validity has not been evaluated experimentally. We report experimental observations of a large quantity of both stacking faults and twins in nc Ni deformed at relatively low stresses in a tensile test. The experimental findings indicate that the GPFE curves can reasonably explain the formation of stacking faults, but they alone were not able to adequately predict the propensity of deformation twinning.

A generalized self-consistent method accounting for fiber section shape

A three-phase confocal elliptical cylinder model is proposed for fiber-reinforced composites, in terms of which a generalized self-consistent method is developed for fiber-reinforced composites accounting for variations in fiber section shapes and randomness in fiber section orientation. The reasonableness of the fiber distribution function in the present model is shown. The dilute, self-consistent, differential and Mori-Tanaka methods are also extended to consider randomness in fiber section orientation in a statistical sense. A full comparison is made between various micromechanics methods and with the Hashin and Shtrikman's bounds. The present method provides convergent and reasonable results for a full range of variations in fiber section shapes (from circular fibers to ribbons), for a complete spectrum of the fiber volume fraction (from 0 to 1, and the latter limit shows the correct asymptotic behavior in the fully packed case) and for extreme types of the inclusion phases (from voids to rigid inclusions). A very different dependence of the five effective moduli on fiber section shapes is theoretically predicted, and it provides a reasonable explanation on the poor correlation between previous theory and experiment in the case of longitudinal shear modulus.