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.

Numerical Simulation of the Flow Field and Concentration Distribution in the Bulk Growth of Silicon Carbide Crystals

The physical vapor transport (PVT) method is being widely used to grow large-size single SiC crystals. The growth process is associated with heat and mass transport in the growth chamber, chemical reactions among multiple species as well as phase change at the crystal/gas interface. The current paper aims at studying and verifying the transport mechanism and growth kinetics model by demonstrating the flow field and species concentration distribution in the growth system. We have developed a coupled model, which takes into account the mass transport and growth kinetics. Numerical simulation is carried out by employing an in-house developed software based on finite volume method. The results calculated are in good agreement with the experimental observation.

Processing Simulation of Virtual Manufacture in Laser Surface Treatmenthang

Processing simulation is at the bottom of the coral technology of VM and is also difficult due to the complexity of mechanism and diversity of parameters. Previously much research has been mainly carried out on the geometrical simulation or physical simulation respectively. The aim of this paper is to study the processing simulation in laser surface treatment based on the mechanism, put forward the architecture of the whole processing simulation and give the models of the processing. As a result the data structure layers in the whole simulation is presented.

Development of a Cell Size Relaxed Scheme for the Direct Simulation Monte Carlo Method

The conventional direct simulation Monte Carlo (DSMC) method has a strong restriction on the cell size because simulated particles are selected randomly within the cell for collisions. Cells with size larger than the molecular mean free path are generally not allowed in correct DSMC simulations. However, the cell-size induced numerical error can be controlled if the gradients of flow properties are properly involved during collisions. In this study, a large cell DSMC scheme is proposed to relax the cell size restriction. The scheme is applied to simulate several test problems and promising results are obtained even when the cell size is greater than 10 mean free paths of gas molecules. However, it is still necessary, of course, that the cell size be small with respect to the flow field structures that must be resolved.

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.

Computation of turbulence-generated noise by large-eddy simulation

The hybrid method of large eddy simulation (LES) and the Lighthill analogy is being developed to compute the sound radiated from turbulent flows. The results obtained from the hybrid method are often contaminated by the absence of small scales in LES, since the energy level of sound is much smaller than that of turbulent flows. Previous researches investigate the effects of subgrid sacle (SGS) eddies on the frequency spectra of sound radiated by isotropic turbulence and suggest a SGS noise model to represent the SGS contributions to the frequency spectra. Their investigations are conducted in physical space and are unavoidably influenced by boundary conditions. In this paper, we propose to perform such calculations in Fourier space so that the effects of boundary conditions can be correctly treated. Posteriori tests are carried out to investigate the SGS contribution to the sound. The results obtained recover the -7/2 law within certain wave-number ranges, but under-estimate the amplitudes of the frequency spectra. The reason for the underestimation is also discussed.

Prediction of particle distribution in isotropic turbulence by large-eddy simulation

The recent application of large-eddy simulation (LES) to particle-laden turbulence requires that the LES with a subgrid scale (SGS) model could accurately predict particle distributions. Usually, a SGS particle model is used to recover the small-scale structures of velocity fields. In this study, we propose a rescaling technique to recover the effects of small-scale motions on the preferential concentration of inertial particles. The technique is used to simulate particle distribution in isotropic turbulence by LES and produce consistent results with direct numerical simulation (DNS). Key words: particle distribution, particle-laden turbulence, large-eddy simulation, subgrid scale model.

An algorithm for coarse particle sedimentation simulation by Stokesian dynamics

Coarse Particle sedimentation is studied by using an algorithm with no adjustable parameters based on stokesian dynamics. Only inter-particle interactions of hydrodynamic force and gravity are considered. The sedimentation of a simple cubic array of spheres is used to verify the computational results. The scaling and parallelism with OpenMP of the method are presented. Random suspension sedimentation is investigated with Mont Carlo simulation. The computational results are shown in good agreement with experimental fitting at the lower computational cost of O(N In N).

Size-dependent elastic properties of Ni nanofilms by molecular dynamics simulation

Size-dependent elastic properties of Ni nanofilms are investigated by molecular dynamics ( MD) simulations with embedded atom method (EAM). The surface effects are considered by calculating the surface relaxation, surface energy, and surface stress. The Young's modulus and yield stress are obtained as functions of thickness and crystallographic orientation. It is shown that the surface relaxation has important effects on the the elastic properties at nanoscale. When the surface relaxation is outward, the Young's modulus decreases with the film thickness decreasing, and vice versa. The results also show that the yield stresses of the films increase with the films becoming thinner. With the thickness of the nanofilms decreasing, the surface effects on the elastic properties become dominant.

Molecular dynamics simulation of adhesion forces in a dipalmitoylphosphatidylcholine membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.