The Structure and Evolution of a Two-Dimensional H2/O2/Ar Cellular Detonation

This paper reports on two-dimensional numerical simulation of cellular detonation wave in a / / mixture with low initial pressure using a detailed chemical reaction model and high order WENO scheme. Before the final equilibrium structure is produced, a fairly regular but still non-equilibrium mode is observed during the early stage of structure formation process. The numerically tracked detonation cells show that the cell size always adapts to the channel height such that the cell ratio is fairly independent of the grid sizes and initial and boundary conditions. During the structural evolution in a detonation cell, even as the simulated detonation wave characteristics suggest the presence of an ordinary detonation, the evolving instantaneous detonation state indicates a mainly underdriven state. As a considerable region of the gas mixture in a cell is observed to be ignited by the incident wave and transverse wave, it is further suggested that these two said waves play an essential role in the detonation propagation.

Influence of the Processing Conditions on the Characteristics of the Clad Layers Produced with Laminar Plasma Technology

Laminar plasma technology was used to produce ceramic hardened layers of Al2O3-40% mass Ni composite powders on stainless steel substrates. In order to investigate the influences of processing conditions on the morphologies of the surface modified layers, two different powder-feeding methods were tested, one with carrier gas called the powder injection method, and the other without carrier gas called powder transfers method. The microscopic investigations demonstrate that the cross-section of the clad layers consists of two distinct microstructural regions, in which the Al2O3 phases exhibit different growth mechanisms. When the powder transfers method is adopted, the number density and volume fraction of the Al2O3 particles increase considerably and their distributions exhibit zonal periodical characteristics. When the powder-feeding rate increases, the microstructure of the Al2O3 phases changes from a small globular to a long needle shape. Finite element simulations show that the transient thermo-physical features of the pool substances, such as solidification rate and cooling rate, influence strongly the mechanisms of the nucleation and the directional growth of the Al2O3 phases in the thermal processing.

Determination of Proper Processing Conditions of Material Surface Heat Treatment Through Finite Element Method

A two-dimensional model has been developed based on the experimental results of stainless steel remelting with the laminar plasma technology to investigate the transient thermo-physical characteristics of the melt pool liquids. The influence of the temperature field, temperature gradient, solidification rate and cooling rate on the processing conditions has been investigated numerically. Not only have the appropriate processing conditions been determined according to the calculations, but also they have been predicted with a criterion established based on the concept of equivalent temperature area density (ETAD) that is actually a function of the processing parameters and material properties. The comparison between the resulting conditions shows that the ETAD method can better predict the optimum condition.

激光诱发活塞的热负荷

活塞的热疲劳性能对柴油发动机的全寿命至关重要。由于能量有限和可控性差等缺点，现有实验系统均不能满意地进行活塞热负荷模拟实验。为此，提出并建立了一套激光诱发活塞热负荷的实验系统。该系统通过对激光束的空间整形，使之投射到活塞表面后诱发的温度场分布满足特定要求。基于PROFIBUS-DP现场总线技术实现了系统集成和实验过程的全反馈控制。针对活塞的典型热负荷条件，即高周热疲劳和热冲击分别进行实验，以模拟正常工作循环和“启动一停车”等热负荷或转速突变工况。通过设置加热一冷却周期或上限下限温度，可以获得相应的热负荷模拟实验效果。研究结果表明，采用经光束整形的激光进行活塞热负荷模拟实验具有周期短、可控性好等优点。 

Concentration distribution in solution crystal growth: effect of moving interface conditions

The linear diffusion-reaction theory with finite interface kinetics is employed to describe the dissolution and the growth processes. The results show that it is imperative to consider the effect of the moving interfaces on the concentration distribution at the growth interface for some cases. For small aspect ratio and small gravity magnitude, the dissolution and the growth interfaces must be treated as the moving boundaries within an angle range of 0 degrees < gamma < 50 degrees in this work. For large aspect ratio or large gravity magnitude, the effect of the moving interfaces on the concentration distribution at the growth interface can be neglected except for gamma < - 50 degrees.

Wave dynamic processes induced by a supersonic projectile discharging from a shock tube

A numerical study on wave dynamic processes occurring in muzzle blast flows, which are created by a supersonic projectile released from the open-end of a shock tube into ambient air, is described in this paper. The Euler equations, assuming axisymmetric flows, are solved by using a dispersion-controlled scheme implemented with moving boundary conditions. Three test cases are simulated for examining friction effects on the muzzle flow. From numerical simulations, the wave dynamic processes, including two blast waves, two jet flows, the bow shock wave and their interactions in the muzzle blasts, are demonstrated and discussed in detail. The study shows that the major wave dynamic processes developing in the muzzle flow remain similar when the friction varies, but some wave processes, such as shock-shock interactions, shock-jet interactions and the contact surface instability, get more intensive, which result in more complex muzzle blast flows.

Experimental study on two-phase gas-liquid flow patterns at normal and reduced gravity conditions

Experimental studies have been performed for horizontal two-phase air-water flows at normal and reduced gravity conditions in a square cross-section channel. The experiments at reduced gravity are conducted on board the Russian IL-76 reduced gravity airplane. Four flow patterns, namely bubble, slug, slug-annular transition and annular flows, are observed depending on the liquid and gas superficial velocities at both conditions. Semi-theoretical Weber number model is developed to include the shape influence on the slug-annular transition. It is shown that its prediction is in reasonable agreement with the experimental slug-annular transition under both conditions. For the case of two-phase gas-liquid flow with large value of the Froude number, the drift-flux model can predict well the observed boundary between bubble and slug flows.

Statistical Simulation of Rarefied Gas Flows in Micro-Channels

Rarefied gas flows through micro-channels are simulated using particle approaches, named as the information preservation (IP) method and the direct simulation Monte Carlo (DSMC) method. In simulating the low speed flows in long micro-channels the DSMC method encounters the problem of large sample size demand and the difficulty of regulating boundary conditions at the inlet and outlet. Some important computational issues in the calculation of long micro-channel flows by using the IP method, such as the use the conservative form of the mass conservation equation to guarantee the adjustment of the inlet and outlet boundary conditions and the super-relaxation scheme to accelerate the convergence process, are addressed. Stream-wise pressure distributions and mass fluxes through micro-channels given by the IP method agree well with experimental data measured in long micro-channels by Pong et al. (with a height to length ratio of 1.2:3000), Shih et al. (l.2:4800), Arkilic et al. and Arkilic (l.3:7500), respectively. The famous Knudsen minimum of normalized mass flux is observed in IP and DSMC calculations of a short micro-channel over the entire flow regime from continuum to free molecular, whereas the slip Navier-Stokes solution fails to predict it.

Modelling of Laminar Plasma Jet Impinging on a Flatplate with Approximate Box Relaxation Method

Approximate Box Relaxation method was used t'o simulate a plasma jet flow impinging on a flatplate at atmospheric pressure, to achieve a better understanding of the characteristics of plasma jet in materials surface treating. The flow fields under different conditions were simulated and analyzed. The distributions of temperature, velocity and pressure were obtained by modelling. Computed results indicate that this numerical method is suitable for simulation of the flow characteristics of plasma jet: and is helpful for understanding of the mechanism of the plasma-material processing.

Damage analysis of tetragonal perovskite structure ceramics implicated by asymptotic field solutions and boundary conditions

Cracking of ceramics with tetragonal perovskite grain structure is known to appear at different sites and scale level. The multiscale character of damage depends on the combined effects of electromechanical coupling, prevailing physical parameters and boundary conditions. These detail features are exhibited by application of the energy density criterion with judicious use of the mode I asymptotic and full field solution in the range of r/a = 10(-4) to 10(-2) where r and a are, respectively, the distance to the crack tip and half crack length. Very close to the stationary crack tip, bifurcation is predicted resembling the dislocation emission behavior invoked in the molecular dynamics model. At the macroscopic scale, crack growth is predicted to occur straight ahead with two yield zones to the sides. A multiscale feature of crack tip damage is provided for the first time. Numerical values of the relative distances and bifurcation angles are reported for the PZT-4 ceramic subjected to different electric field to applied stress ratio and boundary conditions that consist of the specification of electric field/mechanical stress, electric displacement/mechanical strain, and mixed conditions. To be emphasized is that the multiscale character of damage in piezoceramics does not appear in general. It occurs only for specific combinations of the external and internal field parameters, elastic/piezoelectric/dielectric constants and specified boundary conditions. (C) 2002 Published by Elsevier Science Ltd.

Compact difference approximation with consistent boundary condition

For simulating multi-scale complex flow fields it should be noted that all the physical quantities we are interested in must be simulated well. With limitation of the computer resources it is preferred to use high order accurate difference schemes. Because of their high accuracy and small stencil of grid points computational fluid dynamics (CFD) workers pay more attention to compact schemes recently. For simulating the complex flow fields the treatment of boundary conditions at the far field boundary points and near far field boundary points is very important. According to authors' experience and published results some aspects of boundary condition treatment for far field boundary are presented, and the emphasis is on treatment of boundary conditions for the upwind compact schemes. The consistent treatment of boundary conditions at the near boundary points is also discussed. At the end of the paper are given some numerical examples. The computed results with presented method are satisfactory.

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.

Numerical study on rainfall infiltration in rock-soil slop

A mathematical model for the rain infiltration in the rock-soil slop has been established and solved by using the finite element method. The unsteady water infiltrating process has been simulated to get water content both in the homogeneous and heterogeneous media. The simulated results show that the rock blocks in the rock-soil slop can cause the wetting front moving fast. If the rain intensity is increased, the saturated region will be formed quickly while other conditions are the same. If the rain intensity keeps a constant, it is possible to accelerate the generation of the saturated region by properly increasing the vertical filtration rate of the rock-soil slop. However, if the vertical filtration rate is so far greater than the rain intensity, it will be difficult to form the saturated region in the rock-soil slop. The numerical method was verified by comparing the calculation results with the field test data.

Effects of entry conditions on cracked kerosene-fueled supersonic combustor performance

Supersonic combustion of thermally cracked kerosene was experimentally investigated in two model supersonic combustors with different entry cross-section areas. Effects of entry static pressure, entry Mach number, combustor entry geometry, and injection scheme on combustor performance were systematically investigated and discussed based on the measured static pressure distribution and specific thrust increment due to combustion. In addition, the methodology for characterizing flow rate and composition of cracked kerosene was detailed. Using a pulsed Schlieren system, the interaction of supercritical and cracked kerosene jet plumes with a Mach 2.5 crossflow was also visualized at different injection temperatures. The present experimental results suggest that the use of a higher combustor entry Mach number as well as a larger combustor duct height would suppress the boundary layer separation near the combustor entrance and avoid the problem of inlet un- start.