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.

A backward-simulation based Rao-Blackwellized particle smoother for conditionally linear Gaussian models

In this article, we develop a new Rao-Blackwellized Monte Carlo smoothing algorithm for conditionally linear Gaussian models. The algorithm is based on the forward-filtering backward-simulation Monte Carlo smoother concept and performs the backward simulation directly in the marginal space of the non-Gaussian state component while treating the linear part analytically. Unlike the previously proposed backward-simulation based Rao-Blackwellized smoothing approaches, it does not require sampling of the Gaussian state component and is also able to overcome certain normalization problems of two-filter smoother based approaches. The performance of the algorithm is illustrated in a simulated application. © 2012 IFAC.

深水机器人低成本导航系统的位置估计方法研究

依据6000米自治水下机器人及其长基线声学定位系统现有的导航设备,将测距声信标和机器人载体携带的低成本导航传感器:涡轮式计程仪,压力传感器以及TCM2电子罗盘测量的导航数据相融合,分别提出两种基于EKF的导航数据融合算法,对机器人的位置以及水流参数进行估计,解决复杂环境下的深水机器人位置估计问题.蒙特卡洛仿真实验和湖上试验数据后处理表明,设计的位置估计算法收敛快,精度高,计算时间小,能够满足深水机器人的导航需要.

IN-SITU TEM PLASMA CHIP NANOFABRICATION AND CHARACTERIZATION

A silicon-based microcell was fabricated with the potential for use in in-situ transmission electron microscopy (TEM) of materials under plasma processing. The microcell consisted of 50 nm-thick film of silicon nitride observation window with 60μm distance between two electrodes. E-beam scattering Mont Carlo simulation showed that the silicon nitride thin film would have very low scattering effect on TEM primary electron beam accelerated at 200 keV. Only 4.7% of primary electrons were scattered by silicon nitride thin film and the Ar gas (60 μm thick at 1 atm pressure) filling the space between silicon nitride films. Theoretical calculation also showed low absorption of high-energy e-beam electrons. Because the plasma cell needs to survive the high vacuum TEM chamber while holding 1 atm internal pressure, a finite element analysis was performed to find the maximum stress the low-stress silicon nitride thin film experienced under pressure. Considering the maximum burst stress of low-stress silicon nitride thin film, the simulation results showed that the 50 nm silicon nitride thin film can be used in TEM under 1 atm pressure as the observation window. Ex-situ plasma generation experiment demonstrated that air plasma can be ignited at DC voltage of 570. A Scanning electron microscopy (SEM) analysis showed that etching and deposition occurred during the plasma process and larger dendrites formed on the positive electrode.

The generalized Onsager model for the secondary flow in a high-speed rotating cylinder

The generalizations of the Onsager model for the radial boundary layer and the Carrier-Maslen model for the end-cap axial boundary layer in a high-speed rotating cylinder are formulated for studying the secondary gas flow due to wall heating and due to insertion of mass, momentum and energy into the cylinder. The generalizations have wider applicability than the original Onsager and Carrier-Maslen models, because they are not restricted to the limit A >> 1, though they are restricted to the limit R e >> 1 and a high-aspect-ratio cylinder whose length/diameter ratio is large. Here, the stratification parameter A = root m Omega(2)R(2)/2k(B)T). This parameter A is the ratio of the peripheral speed, Omega R, to the most probable molecular speed, root 2k(B)T/m, the Reynolds number Re = rho w Omega R(2)/mu, where m is the molecular mass, Omega and R are the rotational speed and radius of the cylinder, k(B) is the Boltzmann constant, T is the gas temperature, rho(w) is the gas density at wall, and mu is the gas viscosity. In the case of wall forcing, analytical solutions are obtained for the sixth-order generalized Onsager equations for the master potential, and for the fourth-order generalized Carrier-Maslen equation for the velocity potential. For the case of mass/momentum/energy insertion into the flow, the separation-of-variables procedure is used, and the appropriate homogeneous boundary conditions are specified so that the linear operators in the axial and radial directions are self-adjoint. The discrete eigenvalues and eigenfunctions of the linear operators (sixth-order and second-order in the radial and axial directions for the Onsager equation, and fourth-order and second-order in the axial and radial directions for the Carrier-Maslen equation) are determined. These solutions are compared with direct simulation Monte Carlo (DSMC) simulations. The comparison reveals that the boundary conditions in the simulations and analysis have to be matched with care. The commonly used `diffuse reflection' boundary conditions at solid walls in DSMC simulations result in a non-zero slip velocity as well as a `temperature slip' (gas temperature at the wall is different from wall temperature). These have to be incorporated in the analysis in order to make quantitative predictions. In the case of mass/momentum/energy sources within the flow, it is necessary to ensure that the homogeneous boundary conditions are accurately satisfied in the simulations. When these precautions are taken, there is excellent agreement between analysis and simulations, to within 10 %, even when the stratification parameter is as low as 0.707, the Reynolds number is as low as 100 and the aspect ratio (length/diameter) of the cylinder is as low as 2, and the secondary flow velocity is as high as 0.2 times the maximum base flow velocity. The predictions of the generalized models are also significantly better than those of the original Onsager and Carrier-Maslen models, which are restricted to thin boundary layers in the limit of high stratification parameter.

The generalized Onsager model for a binary gas mixture

The Onsager model for the secondary flow field in a high-speed rotating cylinder is extended to incorporate the difference in mass of the two species in a binary gas mixture. The base flow is an isothermal solid-body rotation in which there is a balance between the radial pressure gradient and the centrifugal force density for each species. Explicit expressions for the radial variation of the pressure, mass/mole fractions, and from these the radial variation of the viscosity, thermal conductivity and diffusion coefficient, are derived, and these are used in the computation of the secondary flow. For the secondary flow, the mass, momentum and energy equations in axisymmetric coordinates are expanded in an asymptotic series in a parameter epsilon = (Delta m/m(av)), where Delta m is the difference in the molecular masses of the two species, and the average molecular mass m(av) is defined as m(av) = (rho(w1)m(1) + rho(w2)m(2))/rho(w), where rho(w1) and rho(w2) are the mass densities of the two species at the wall, and rho(w) = rho(w1) + rho(w2). The equation for the master potential and the boundary conditions are derived correct to O(epsilon(2)). The leading-order equation for the master potential contains a self-adjoint sixth-order operator in the radial direction, which is different from the generalized Onsager model (Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159), since the species mass difference is included in the computation of the density, viscosity and thermal conductivity in the base state. This is solved, subject to boundary conditions, to obtain the leading approximation for the secondary flow, followed by a solution of the diffusion equation for the leading correction to the species mole fractions. The O(epsilon) and O(epsilon(2)) equations contain inhomogeneous terms that depend on the lower-order solutions, and these are solved in a hierarchical manner to obtain the O(epsilon) and O(epsilon(2)) corrections to the master potential. A similar hierarchical procedure is used for the Carrier-Maslen model for the end-cap secondary flow. The results of the Onsager hierarchy, up to O(epsilon(2)), are compared with the results of direct simulation Monte Carlo simulations for a binary hard-sphere gas mixture for secondary flow due to a wall temperature gradient, inflow/outflow of gas along the axis, as well as mass and momentum sources in the flow. There is excellent agreement between the solutions for the secondary flow correct to O(epsilon(2)) and the simulations, to within 15 %, even at a Reynolds number as low as 100, and length/diameter ratio as low as 2, for a low stratification parameter A of 0.707, and when the secondary flow velocity is as high as 0.2 times the maximum base flow velocity, and the ratio 2 Delta m/(m(1) + m(2)) is as high as 0.5. Here, the Reynolds number Re = rho(w)Omega R-2/mu, the stratification parameter A = root m Omega R-2(2)/(2k(B)T), R and Omega are the cylinder radius and angular velocity, m is the molecular mass, rho(w) is the wall density, mu is the viscosity and T is the temperature. The leading-order solutions do capture the qualitative trends, but are not in quantitative agreement.

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.

Vaporization of heated materials into discharge plasmas

The vaporization of condensed materials in contact with high-current discharge plasmas is considered. A kinetic numerical method named direct simulation Monte Carlo (DSMC) and analytical kinetic approaches based on the bimodal distribution function approximation are employed. The solution of the kinetic layer problem depends upon the velocity at the outer boundary of the kinetic layer which varies from very small, corresponding to the high-density plasma near the evaporated surface, up to the sound speed, corresponding to evaporation into vacuum. The heavy particles density and temperature at the kinetic and hydrodynamic layer interface were obtained by the analytical method while DSMC calculation makes it possible to obtain the evolution of the particle distribution function within the kinetic layer and the layer thickness.

Thickness And Component Distributions Of Yttrium-Titanium Alloy Films In Electron-Beam Physical Vapor Deposition

Thickness and component distributions of large-area thin films are an issue of international concern in the field of material processing. The present work employs experiments and direct simulation Monte Carlo (DSMC) method to investigate three-dimensional low-density, non-equilibrium jets of yttrium and titanium vapor atoms in an electron-beams physical vapor deposition (EBPVD) system furnished with two or three electron-beams, and obtains their deposition thickness and component distributions onto 4-inch and 6-inch mono-crystal silicon wafers. The DSMC results are found in excellent agreement with our measurements, such as evaporation rates of yttrium and titanium measured in-situ by quartz crystal resonators, deposited film thickness distribution measured by Rutherford backscattering spectrometer (RBS) and surface profilometer and deposited film molar ratio distribution measured by RBS and inductively coupled plasma atomic emission spectrometer (ICP-AES). This can be taken as an indication that a combination of DSMC method with elaborate measurements may be satisfactory for predicting and designing accurately the transport process of EBPVD at the atomic level.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Rarefied gas dynamics: fundamentals, simulations and micro flows

Kinetic Study Of The Rayleigh-Benard Flows

The present paper employs the direct simulation Monte Carlo (DSMC) method to study the Rayleigh-Benard flows, where the temperature ratio of the upper to lower plate is fixed to 0.1. For a Knudsen number (Kn) of 0.01, as the Rayleigh number (Ra) increases, the flow changes from the thermal conductive state to the convective state at about Ra=1700, and the calculated relation of heat flux through the lower plate versus Ra is in good agreement with classical experimental and theoretical results. For Kn=0.05, the thermal conductive state remains stable, and the increase of Ra cannot trigger thermal instability.