Deflections and curvatures of a film-substrate structure with the presence of gradient stress in MEMS applications

The close form solutions of deflections and curvatures for a film-substrate composite structure with the presence of gradient stress are derived. With the definition of more precise kinematic assumption, the effect of axial loading due to residual gradient stress is incorporated in the governing equation. The curvature of film-substrate with the presence of gradient stress is shown to be nonuniform when the axial loading is nonzero. When the axial loading is zero, the curvature expressions of some structures derived in this paper recover the previous ones which assume the uniform curvature. Because residual gradient stress results in both moment and axial loading inside the film-substrate composite structure, measuring both the deflection and curvature is proposed as a safe way to uniquely determine the residual stress state inside a film-substrate composite structure with the presence of gradient stress.

Drag on an ellipsoid particle in free-molecular plasma flow

Based on the analysis of molecular gas dynamics, the drag and moment acting on an ellipsoid particle of revolution X-2/a(2) + Y-2/a(2) + Z(2)/c(2) = 1, as an example of nonspherical particles, are studied under the condition of free-molecular plasma flow with thin plasma sheaths. A nonzero moment which causes nonspherical particle self-oscillation and self-rotation around its own axis in the plasma flow-similar to the pitching moment in aerodynamics-is discovered for the first time. When the ratio of axis length c/a is unity, the moment is zero and the drag formula are reduced to the well-known results of spherical particles. The effects of the particle-plasma relative velocity, the plasma temperature, and the particle materials on the drag and moment are also investigated.

Distribution, stability, bifurcations and catastrophe of steady rotation of a symmetrical heavy gyroscope with viscous-liquid-filled cavity

The number, the angles of orientation and the stability in Rumyantsev Movchan's sense of oblique steady rotations of a symmetric heavy gyroscope with a cavity completely filled with a uniform viscous liquid, possessing a fixed point 0 on its symmetric axis. are given for various values of the parameters. By taking the square of the upright component of the angular momentum M2 as a control parameter, three types of bifurcation diagrams of the steady rotations, two types of jumps and two kinds of local catastrophes, one being the symmetric reduced cusp type and the other being of the symmetric reduced butterfly type, are obtained. By taking account of the M2-damping owing to the moment of unavoidable faint friction, two different modes for the gyroscope, initially in a stable quasi-steady upright rotation with a nutation angle theta(s) equal to zero, to topple over are found.

On dilatational plastic constitutive equation of ductile materials and plastic loading paths at bifurcation

The dilatational plastic constitutive equation presented in this paper is proved to be in a form of generality. Based on this equation, the constitutive behaviour of materials at the moment of bifurcation is demonstrated to follow a loading path with the response as "soft" as possible.

Statistical-theory of multiphoton excitation of polyatomic-molecules based on the wigner function

This paper is aimed at establishing a statistical theory of rotational and vibrational excitation of polyatomic molecules by an intense IR laser. Starting from the Wigner function of quantum statistical mechanics, we treat the rotational motion in the classical approximation; the vibrational modes are classified into active ones which are coupled directly with the laser and the background modes which are not coupled with the laser. The reduced Wigner function, i.e., the Wigner function integrated over all background coordinates should satisfy an integro-differential equation. We introduce the idea of ``viscous damping'' to handle the interaction between the active modes and the background. The damping coefficient can be calculated with the aid of the well-known Schwartz–Slawsky–Herzfeld theory. The resulting equation is solved by the method of moment equations. There is only one adjustable parameter in our scheme; it is introduced due to the lack of precise knowledge about the molecular potential. The theory developed in this paper explains satisfactorily the recent absorption experiments of SF6 irradiated by a short pulse CO2 laser, which are in sharp contradiction with the prevailing quasi-continuum theory. We also refined the density of energy levels which is responsible for the muliphoton excitation of polyatomic molecules.

A closure theory of intermittency of turbulence

The variational approach to the closure problem of turbulence theory, proposed in an earlier article [Phys. Fluids 26, 2098 (1983); 27, 2229 (1984)], is extended to evaluate the flatness factor, which indicates the degree of intermittency of turbulence. Since the flatness factor is related to the fourth moment of a turbulent velocity field, the corresponding higher-order terms in the perturbation solution of the Liouville equation have to be considered. Most closure methods discard these higher-order terms and fail to explain the intermittency phenomenon. The computed flatness factor of the idealized model of infinite isotropic turbulence ranges from 9 to 15 and has the same order of magnitude as the experimental data of real turbulent flows. The intermittency phenomenon does not necessarily negate the Kolmogorov k−5/3 inertial range spectrum. The Kolmogorov k−5/3 law and the high degree of intermittency can coexist as two consistent consequences of the closure theory of turbulence. The Kolmogorov 1941 theory [J. Fluid Mech. 62, 305 (1974)] cannot be disqualified merely because the energy dissipation rate fluctuates.

Chemical equilibrium analysis and CFD simulation of coal combustion and NOx formation

A full two-fluid model of reacting gas-particle flows with an algebraic unified second-order moment (AUSM) turbulence-chemistry model is used to simulate Beijing coal combustion and NOx formation. The sub-models are the k-epsilon-kp two-phase turbulence model, the EBU-Arrhenius volatile and CO combustion model, the six-flux radiation model, coal devolatilization model and char combustion model. The blocking effect on NOx formation is discussed. In addition, the chemical equilibrium analysis is used to predict NOx concentration at different temperature. Results of CID simulation and chemical equilibrium analysis show that, optimizing air dynamic parameters can delay the NOx formation and decrease NOx emission, but it is effective only in a restricted range. In order to decrease NOx emission near to zero, the re-burning or other chemical methods must be used.

Research progress on high-enthalpy and hypersonic flows

The research progress on high-enthalpy and hypersorlic flows having been achieved in the Institute of Mechanics, Chinese Academy of Sciences, is reported in this paper. The paper consists of three main parts: The first part is on the techniques to develop advanced hypersonic test facilities, in which the detonation-driven shock-reflected tunnel and the detonation-driven shock-expanded tube are introduced. The shock tunnel can be used for generating hypersonic flows of a Mach number ranging from 10 to 20, and the expansion tube is applicable to simulate the flows with a speed of 7 similar to 10km/s. The second part is dedicated to the shock tunnel nozzle flow diagnosis to examine properties of the hypersonic flows thus created. The third part is on experiments and numerical simulations. The experiments include measuring the aerodynamic pitching moment and heat transfer in hypersonic flows, and the numerical work reports nozzle flow simulations and flow non-equilibrium effects on the possible experiments that may be carried out on the above-mentioned hypersonic test facilities.

郭永怀文集

基于多矩的有限体积浸入边界法研究

浸入边界法（Immersed Boundary Method）是计算流体力学中求解具有复杂、移动边界流动问题的一类有效途径，该方法在笛卡尔坐标系上离散求解流体控制方程，并通过在控制方程中添加相应源相来代表浸入边界。尽管浸入边界法借助其简单、高效的显著特点在计算流体力学应用中显示出极强的生命力，特别是针对复杂的实际流动及动边界流动问题有着无可比拟的优势，但仍有许多问题需要进一步的研究。 本论文基于浸入边界方法及多矩VSIAM3（Volume/Surface Integrated Average Multi-Moment Method）格式提出了一种不可压缩流体求解数值格式。不可压N-S方程使用VSIAM3格式进行法进行离散，引入浸入边界法处理复杂、移动流动边界条件，使用虚拟网格方法计算动量方程修正项，同时还考虑了对连续方程的修正。VSIAM3格式是一种基于多矩的有限体积法，在方程的离散中总是使用两种或两种以上的矩，如：VIA（Volume Integrated Average）和SIA（Surface Integrated Average）。而不同的矩在求解过程中依据不同形式的控制方程使用不同的离散方法进行更新。VSIAM3格式更多的局地自由度及同时使用交错网格和同位网格的特点使浸入边界法的实施更加便利、高效。研究中，浸入边界法不仅应用于处理动力边界条件，同样可以处理热动力边界条件。 研究中对大量经典算例进行了数值实验，包括一维线性初始问题、方腔流问题、二维绕静止及振荡圆柱流动、三维绕球流动及热对流问题等。数值结果同实验值及其它计算结果保持一致，该算法可准确、高效处理具有复杂、移动边界及存在热对流的不可压流动问题，为实际应用打下了基础。

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

AUTOROTATION OF TWO TANDEM TRIANGULAR CYLINDERS

The autorotation of two tandem triangular cylinders at different gap distances is investigated by numerical simulations. At the Reynolds number of 200， three distinct regimes are observed with the increase of gap distance: namely, angular oscillation, quasi-periodic autorotation and ‘chaotic’ autorotation. For various gap distances, the characteristic of vortex shedding and vortex interaction are discussed. The phase graphs (angular acceleration vs. angular velocity) and the power spectra of moment are analyzed to characterize the motion of the cylinder. The Lyapunov exponent is also calculated to identify the existence of chaos.

Motion simulation of annular polishing

Based on the Coulomb friction model, the frictional motion model of workpiece relating to the polishing pad was presented in annular polishing. By the dynamic analysis software, the model was simulated and analysed. The conclusions from the results were that the workpiece did not rotate steadily. When the angular velocity of ring and the direction were the same as that of the polishing pad, the angular velocity of workpiece hoicked at the beginning and at the later stage were the same as that of the polishing pad before contacting with the ring. The angular velocity of workpiece vibrated at the moment of contacting with the ring. After that the angular velocity of workpiece increased gradually and fluctuated at a given value, while the angular velocity of ring decreased gradually and also fluctuated at a given value. Since the contact between the workpiece and the ring was linear, their linear velocities and directions should be the same. But the angular velocity of workpiece was larger than that of the polishing pad on the condition that the radius of the workpiece was less than that of the ring. This did not agree with the pure translation principle and the workpiece surface could not be flat, either. Consequently, it needed to be controlled with the angular velocity of ring and the radii of the ring and the workpiece, besides friction to make the angular velocity of workpiece equal to that of the polishing pad for obtaining fine surface flatness of the workpiece. Copyright © 2007 Inderscience Enterprises Ltd.}

Synthesis and magnetic properties of Ni-doped zinc oxide powders

Polycrystalline Zn1-xNixO diluted magnetic semiconductors have been successfully synthesized by an auto-combustion method. X-ray diffraction measurements indicated that the 5 at% Ni-cloped ZnO had the pure wurtzite structure. Refinements of cell parameters from powder diffraction data revealed that the cell parameters of Zn0.95Ni0.05O were a little bit larger than ZnO. Transmission electron microscopy observation showed that the as-synthesized powders were of the size similar to 60 nm. Magnetic investigations showed that the nanocystalline Zn0.95Ni0.05O possessed room temperature ferromagnetisin with the saturation magnetic moment of 0.1 emu/g (0.29 mu(B)/Ni2+). (c) 2005 Elsevier B.V. All rights reserved.

Synthesis and characterization of room-temperature ferromagnetism in Fe- and Ni-co-doped In2O3

Observation of room-temperature ferromagnetisin in Fe- and Ni-co-doped In2O3 samples (In0.9Fe0.1-xNix)(2)O-3 (0 <= x <= 0.1) prepared by citric acid sol-gel auto-igniting method is reported. All of the samples with intermediate x values are ferromagnetic at roomtemperature. The highest saturation magnetization (0.453 mu B/Fe + Ni ions) moment is reached in the sample with x = 0.04. The highest solubility of Fe and Ni ions in the In2O3 lattice is around 10 and 4 at%, respectively. The 10 at% Fe-doped sample is found to be weakly ferromagnetic, while the 10at% Ni-doped sample is paramagnetic. Extensive structure including Extended X-ray absorption fine structure (EXAFS), magnetic and magneto-transport including Hall effects studies on the samples indicate the observed ferromagnetism is intrinsic rather than from the secondary impurity phases. (c) 2007 Elsevier Ltd. All rights reserved.