The cellular structure of a two-dimensional H2/O2/Ar detonation wave

In this paper, the cellular structure of a two-dimensional detonation wave in a low pressure H2/O2/Ar mixture calculated with a detailed chemical reaction model, high order scheme and high resolution grids is investigated. The regular cellular structure is produced about 1 ms after introducing perturbations in the reaction zone of a steady one-dimensional detonation wave. It is found from the present resolution study that the discrepancies concerning the structure type arising from the coarser grid employed can be resolved using a sufficiently fine grid size of 0.05 mm and below and shows a double-Mach-like strong-type configuration. During the structure evolution process, the structure configuration does not change much in the periods before and after the triple point collision. Through the triple point collision, three regular collision processes are observed and are followed by a quick change to the double-Mach-like configuration. The simulated structure tracks show that there are three different tracks associated with different triple points or the kink on the transverse wave. Comparisons with previous work and experiments indicate the presence of a strong structure for an ordinary detonation.

Effects of Translational Nonequilibrium on Performance of Flowing Chemical Oxygen-Iodine Laser

The effect of the translational nonequilibrium on performance modeling of flowing chemical oxygen-iodine lasers (COIL) is emphasized in this paper. The spectral line broadening (SLB) model is a basic factor for predicting the performances of flowing COIL. The Voigt profile function is a well-known SLB model and is usually utilized. In the case of gas pressure in laser cavity less than 5 torr, a low pressure limit expression of the Voigt profile function is used. These two SLB models imply that ail lasing particles can interact with monochromatic laser radiation. Basically, the inhomogeneous broadening effects are not considered in these two SLB models and they cannot predict the spectral content. The latter requires consideration of finite translational relaxation rate. Unfortunately, it is rather difficult to solve simultaneously the Navier-Stokes (NS) equations and the conservation equations of the number of lasing particles per unit volume and per unit frequency interval. In the operating condition of flowing COIL, it is possible to obtain a perturbational solution of the conservational equations for lasing particles and deduce a new relation between the gain and the optical intensity, i.e., a new gain-saturation relation. By coupling the gain-saturation relation with other governing equations (such as the NS equations, chemical reaction equations and the optical model of gain-equal-loss), We have numerically calculated the performances of flowing COIL. The present results are compared with those obtained by the common rate-equation (RE) model, in which the Voigt profile function and its low pressure limit expression are used. The difference of different model's results is great. For instance, in the case of lasing frequency coinciding with the central frequency of line profile and very low gas pressure, the gain-saturation relation of the present model is quite different with that of the RE model.

Size effects on the melting of nickel nanowires: a molecular dynamics study

The melting process of nickel nanowires are simulated by using molecular dynamics with the quantum Sutten-Chen many-body force field. The wires studied were approximately cylindrical in cross-section and periodic boundary conditions were applied along their length; the atoms were arranged initially in a face-centred cubic structure with the [0 0 1] direction parallel to the long axis of the wire. The size effects of the nanowires on the melting temperatures are investigated. We find that for the nanoscale regime, the melting temperatures of Ni nanowires are much lower than that of the bulk and are linear with the reciprocal of the diameter of the nanowire. When a nanowire is heated up above the melting temperature, the neck of the nanowire begins to arise and the diameter of neck decreases rapidly with the equilibrated running time. Finally, the breaking of nanowire arises, which leads to the formation of the spherical clusters. (C) 2004 Elsevier B.V. All rights reserved.

Gain saturation effects of a flowing chemical oxygen-iodine laser

A new oxygen-iodine medium gain model is developed to include pumping and deactivation of the upper laser levels of the iodine atoms, hyperfine and translation relaxation, as well as the flowing effect. The rate equations for gain of a supersonic flowing cw oxygen-iodine laser (COIL) are described when the medium is stimulated by a single-mode field. The general solution of the self-consistency integral equation is obtained. The result shows that the saturation behaviour in low pressure of the COIL differs from both the inhomogeneous and homogeneous broadening, and exhibits an 'anomalous' saturation phenomenon.

Theoretical and experimental studies of the bell-jar-top inductively coupled plasma

The present paper describes a systematic study of argon plasmas in a bell-jar inductively coupled plasma (ICP) source over the range of pressure 5-20 mtorr and power input 0.2-0.5 kW, Experimental measurements as well as results of numerical simulations are presented. The models used in the study include the well-known global balance model (or the global model) as well as a detailed two-dimensional (2-D) fluid model of the system, The global model is able to provide reasonably accurate values for the global electron temperature and plasma density, The 2-D model provides spatial distributions of various plasma parameters that make it possible to compare with data measured in the experiments, The experimental measurements were obtained using a tuned Langmuir double-probe technique to reduce the RF interference and obtain the light versus current (I-V) characteristics of the probe. Time-averaged electron temperature and plasma density were measured for various combinations of pressure and applied RF power, The predictions of the 2-D model were found to be in good qualitative agreement with measured data, It was found that the electron temperature distribution T-e was more or less uniform in the chamber, It was also seen that the electron temperature depends primarily on pressure, but is almost independent of the power input, except in the very low-pressure regime. The plasma density goes up almost linearly with the power input.

Two-dimensional simulation of inductive plasma sources with self-consistent power deposition

A time averaged two-dimensional fluid model including an electromagnetic module with self-consistent power deposition was developed to simulate the transport of a low pressure radio frequency inductively coupled plasma source. Comparsions with experiment and previous simulation results show, that the fluid model is feasible in a certain range of gas pressure. In addition, the effects of gas pressure and power input have been discussed.

Parameter-estimation of a spherical top inductively-coupled plasma source with global-model

Using spatially averaged global model, we succeed in obtaining some plasma parameters for a low pressure inductively coupled plasma source of our laboratory. As far as the global balance is concerned, the models can give reasonable results of the parameters, such as the global electron temperature and the ion impacting energy, etc. It is found that the ion flow is hardly affected by the neutral gas pressure. Finally, the magnetic effects are calculated by means of the method. The magnetic field can play an important role to increase plasma density and ion current.

Behavior of ar plasma formed in a high-density plasma source - an ecr reactor

In order to develop the ultra-large scale integration(ULSI), low pressure and high density plasma apparatus are required for etching and deposit of thin films. To understand critical parameters such as the pressure, temperature, electrostatic potential and energy distribution of ions impacting on the wafer, it is necessary to understand how these parameters are influenced by the power input and neutral gas pressure. In the present work, a 2-D hybrid electron fluid-particle ion model has been developed to simulate one of the high density plasma sources-an Electron Cyclotron Resonance (ECR) plasma system with various pressures and power inputs in a non-uniform magnetic field. By means of numerical simulation, the energy distributions of argon ion impacting on the wafer are obtained and the plasma density, electron temperature and plasma electrostatic potential are plotted in 3-D. It is concluded that the plasma density depends mainly on both the power input and neutral gas pressure. However, the plasma potential and electron temperature can hardly be affected by the power input, they seem to be primarily dependent on the neutral gas pressure. The comparison shows that the simulation results are qualitatively in good agreement with the experiment measurements.

The positive floating potential of high-temperature particles in plasma

It has been predicted that the floating potential of particles in plasma may become positive when the particle surface temperature is high enough, but, to our knowledge, no positive floating potential has been obtained yet. In the present paper the floating potential theory of high-temperature particles in plasma is developed to cover the positive potential range for the first time, and a general approximate analytical formula for the positive floating potential with a thin plasma sheath and subsonic plasma flow is derived from the new model recently proposed by the authors. The results show that when the floating potential is positive, the net flux of charge incident on the particle approaches a constant similar to the 'electron saturation' phenomena in the case of the electric probes.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

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.

Numerical and Experimental Investigation On The Prevention of Co Deflagration

The exhaust gases from industrial furnaces contain a huge amount of heat and chemical enthalpy. However, it is hard to recover this energy since exhaust gases invariably contain combustible components such as carbon monoxide (CC). If the CO is unexpectedly ignited during the heat recovery process, deflagration or even detonation could occur, with serious consequences such as complete destruction of the equipment. In order to safely utilize the heat energy contained in exhaust gas, danger of its explosion must be fully avoided. The mechanism of gas deflagration and its prevention must therefore be studied. In this paper, we describe a numerical and experimental investigation of the deflagration process in a semi-opened tube. The results show that, upon ignition, a low-pressure wave initially spreads within the tube and then deflagration begins. For the purpose of preventing deflagration, an appropriate amount of nitrogen was injected into the tube at a fixed position. Both simulation and experimental results have shown that the injection of inert gas can successfully interrupt the deflagration process. The peak value of the deflagration pressure can thereby be reduced by around 50%. (C) 2008 Elsevier Ltd. All rights reserved.

Characterization of m-plane ZnO thin film on gamma-LiAlO2 (100) substrate by metal-organic chemical vapor deposition

Non-polar (1 (1) over bar 00)m-plane ZnO thin film has been prepared on gamma-LiAlO2 (100)substrate via the low pressure metal organic chemical vapor deposition. Obvious intensity variation of the E-2 mode in the polarized Raman spectra and the absorption edge shift in the polarized optical transmission spectra indicate that the m-plane film exhibits optical anisotropy, which have applications in certain optical devices, such as the UV modulator and polarization-dependent beam switch. From the atomic force microscopy images, highly-oriented uniform-sized grains of rectangular shape were observed. (c) 2008 Elsevier B.V. All rights reserved.

真空中有机污染物对TiO2/SiO2介质膜损伤阈值的影响

本文主要研究了真空中有机污染物对薄膜阈值的影响。采用电子枪蒸镀的方式制备了由TiO2/SiO2两种材料组合而成的1064nm高反膜，分别在大气和真空中对其损伤阈值进行了测量。发现在真空中污染物易黏附于光学元件表面，导致损伤；真空中的氧分子会加速损伤过程，根据破斑形貌的分析认为是由于在激光作用下氧分子与有机污染物发生反应造成的；真空中的有机污染物被去除后，真空和大气中的损伤阈值差别不大，破斑形貌相似。结果表明有机污染物是真空中损伤阈值降低的主要原因。