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.

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.

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 dynamics: fundamentals, simulations and micro flows

IP方法模拟硬盘读写头气膜润滑问题

本论文中采用IP方法模拟硬盘读写头的气膜润滑问题，并选择耗时但精准的DSMC方法作为该问题中IP方法的检验标准，IP方法计算得到的读写头表面压力分布及读写头所受净升力为读写头的设计提供可靠的参考数据。 首先，将读写头复杂表面简化为平面并假定尺寸为微米（远小于真实的毫米），微米尺寸时可以得到DSMC的模拟结果用于比较。IP法对该问题的模拟结果与DSMC的模拟结果完全吻合。 其次，选定一个真实的读写头表面几何，但仍然假定尺寸为微米。IP法在该问题的应用过程中遇到四点困难，我们提出或采取了相应的解决办法。在复杂表面几何的读写头气膜润滑问题中，IP法的模拟结果与DSMC的模拟结果也相符一致。 最后完成真实几何、真实尺寸（毫米）的读写头气膜润滑问题的模拟，这也是本论文的最终研究目标。IP法模拟该问题时需要采用大网格，为此我们构造二维检验模型验证了大网格的合理性。 真实几何读写头气膜润滑问题的模拟结果表明，在流场的局部区域压力沿着高度方向变化显著，这意味着读写头问题的传统计算方法（即概括性雷诺方程）在该局部区域不成立，故用它计算得到的结果值得认真检验，但文献中尚未给出过这类检验。这正是该论文的研究背景和采用新方法（IP方法）的原因。 本论文研究的问题属于过渡领域、低速流动问题，故本论文还介绍了该领域另外一种粒子模拟方法（LVDSMC方法）以及我们对该方法提出的一些改进。 关键词：硬盘读写头，气膜润滑，概括性雷诺方程，IP方法，DSMC方法，过渡领域，LVDSMC方法

EXPERIMENTS ON AQUATIC ECOSYSTEM REGIME SHIFT IN ENCLOSURES NEAR LAKE DIANCHI, CHINA

To discover how a lake converts from a turbid state to clean state, and what drives this process, we constructed controlled enclosure ecosystems and used the ecological remediation method to force ecosystems to convert from the turbid state to the clean state. Our results show that the driving forces include temperature., macrophyte, silver carp and mussel, which form a combined force to drive the controlled ecosystem to switch. There is a threshold existing in treated enclosure ecosystem during the conversion from turbid to clean state. When TP <0.09 mg.L-1, Chl-a <0.036 mg.L-1, transparency >62 cm, TN <2.15 mg.L-1, CODMn <13.7 mg.L-1, tubidity <10, and the number of algal cells <10(6) cells.L-1, the treated ecosystem changes sharply from turbid to clean state. The conversion process can be divided into three phases: turbid state, clean-turbid transitional state as well as clean state, and described with the power function Y = a*X-b (where Y is water parameter, X is time, a and b are constants), which indicates that the shift in the enclosure ecosystem from turbid to clean state is discontinuous.

Arc Voltage Fluctuation in DC Laminar and Turbulent Plasma Jets Generat Ion

Arc voltage fluctuations in a direct current (DC) non-transferred arc plasma generator are experimentally studied, in generating a jet in the laminar, transitional and turbulent regimes. The study is with a view toward elucidating the mechanism of the fluctuations and their relationship with the generating parameters, arc root movement and flow regimes. Results indicate that the existence of a 300 Hz alternating current (AC) component in the power supply ripples does not cause the transition of the laminar plasma jet into a turbulent state. There exists a high frequency fluctuation at 4 kHz in the turbulent jet regime. It may be related to the rapid movement of the anode attachment point of the arc.

Experimental Study on the Thermal Argon Plasma Generation and Jet Length Change Characteristics at Atmospheric Pressure

The generation, jet length and flow-regime change characteristics of argon plasma issuing into ambient air have been experimentally examined. Different torch structures have been used in the tests. Laminar plasma jets can be generated within a rather wide range of working-gas flow rates, and an unsteady transitional flow state exists between the laminar and turbulent flow regimes. The high-temperature region length of the laminar plasma jet can be over an order longer than that of the turbulent plasma jet and increases with increasing argon flow rate or arc current, while the jet length of the turbulent plasma is less influenced by the generating parameters. The flow field of the plasma jet has very high radial gradients of plasma parameters, and a Reynolds number alone calculated in the ordinary manner may not adequately serve as a criterion for transition. The laminar plasma jet can have a higher velocity than that of an unsteady or turbulent jet. The long laminar plasma jet has good stiffness to withstand the impact of laterally injected cold gas and particulate matter. It could be used as a rather ideal object for fundamental studies and be applied to novel materials processing due to its attractive stable and adjustable properties.

Configuration optimization of hypersonic vehicles under transitional flow conditions

The rarefied gas effects on several configurations are investigated under hypersonic flow conditions using the direct simulation Mont Carlo method. It is found that the Knudsen number, the Mach number, and the angle of attack all play a mixed role in the aerodynamics of a flat plate. The ratio of lift to drag decreases as the Knudsen number increases. Studies on 3D delta wings show that the ratio of lift to drag could be increased by decreasing the wing thickness and/or by increasing the wing span. It is also found that the waveriders could produce larger ratio of lift to drag as compared with the delta wing having the same length, wing span, and cross section area.

A scaling law of photoionization in few-cycle regime

In this paper, a scaling law of photoionization of atoms irradiated by intense, few- cycle laser pulses is established. The scaling law sets a relation to the phase- dependent ionization with the kinetic energy of photoelectrons, the duration and peak intensity of short pulses, and the ionization potential of the target atoms. We find that it will be advantageous to manifest the phase- dependent photoionization by choosing the target atoms with larger ionization potential, using laser with smaller carrier- frequency, and increasing the pulse intensity. (c) 2007 Optical Society of America.

Electron injection by a nanowire in the bubble regime

The triggering of wave-breaking in a three-dimensional laser plasma wake (bubble) is investigated. The Coulomb potential from a nanowire is used to disturb the wake field to initialize the wave-breaking. The electron acceleration becomes more stable and the laser power needed for self-trapping is lowered. Three-dimensional particle-in-cell simulations were performed. Electrons with a charge of about 100 pC can be accelerated stably to energy about 170 MeV with a laser energy of 460 mJ. The first step towards tailoring the electron beam properties such as the energy, energy spread, and charge is discussed. (C) 2007 American Institute of Physics.

Quantum path control in few-optical-cycle regime

We theoretically show that selection of a single quantum path in high-order harmonics generation can be realized in a few-optical-cycle regime with two-color schemes. We also demonstrate, in theory as well, the generation of spectrally smooth and ultrabroad extreme ultraviolet supercontinuum in argon gas which can produce single similar to 79 as pulses with currently available ultrafast laser sources. Our finding can be beneficial for generating isolated sub-100 as extreme ultraviolet pulses.

Bubble regime for ion acceleration in a laser-driven plasma

Proton trapping and acceleration by an electron bubble-channel structure in laser interaction with high-density plasma is investigated by using three-dimensional particle-in-cell simulations. It is shown that protons can be trapped, bunched, and efficiently accelerated for appropriate laser and plasma parameters, and the proton acceleration is enhanced if the plasma consists mainly of heavier ions such as tritium. The observed results are analyzed and discussed in terms of a one-dimensional analytical three-component-plasma wake model.