Experimental investigation of a control method for SI-HCCI-SI transition in a multi-cylinder gasoline engine

In HCCI engines, the Air/Fuel Ratio (AFR) and Residual Gas Fraction (RGF) are difficult to control during the SI-HCCI-SI transition, and this may result in incomplete combustion and/or high pressure raise rates. As a result, there may be undesirably high engine load fluctuations. The objectives of this work are to further understand this process and develop control methods to minimize these load fluctuations. This paper presents data on instantaneous AFR and RGF measurements, both taken by novel experimental techniques. The data provides an insight into the cyclic AFR and RGF fluctuations during the switch. These results suggest that the relatively slow change in the intake Manifold Air Pressure (MAP) and actuation time of the Variable Valve Timing (VVT) are the main causes of undesired AFR and RGF fluctuations, and hence an unacceptable Net IMEP (NIMEP) fluctuation. We also found large cylinder-to-cylinder AFR variations during the transition. Therefore, besides throttle opening control and VVT shifting, cyclic and individual cylinder fuel injection control is necessary to achieve a smooth transition. The control method was developed and implemented in a test engine, and the result was a considerably reduced NIMEP fluctuation during the mode switch. The instantaneous AFR and RGF measurements could furthermore be adopted to develop more sophisticated control methods for SI-HCCI-SI transitions. © 2010 SAE International.

Application of Three-Dimensional Discrete Element Face-to-Face Contact Model with Fissure Water Pressure to Stability Analysis of Landslide in Panluo Iron Mine

Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical phi numerically calculated is less than the one calculated by use of the limit equilibrium method for the same C. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.

Pressure and thermal effects on elastic properties of single crystal alpha-Al2O3 from Monte-Carlo simulation

A rectangular structural unit cell of a-Al2O3 is generated from its hexagonal one. For the rectangular structural crystal with a simple interatomic potential [Matsui, Mineral Mag. 58A, 571 (1994)], the relations of lattice constants to homogeneous pressure and temperature are calculated by using Monte-Carlo method at temperature 298K and 0 GPa, respectively. Both numerical results agree with experimental ones fairly well. By comparing pair distribution function, the crystal structure of a-Al2O3 has no phase transition in the range of systematic parameters. Based on the potential model, pressure dependence of isothermal bulk moduli is predicted. Under variation of general strains, which include of external and internal strains, elastic constants of a-Al2O3 in the different homogeneous load are determined. Along with increase of pressure, axial elastic constants increase appreciably, but nonaxial elastic constants are slowly changed.

An Approximate Method for Evaluating the Shear Band Thickness in Saturated Soils

formula for the thickness of a shear band formed in saturated soils under a simple shear or a combined stress state has been proposed. It is shown that the shear band thickness is dependent on the pore pressure properties of the material and the dilatancy rate, but is independent of the details of the combined stress state. This is in accordance with some separate experimental observations.

Elastic Behaviour and Microstructural Characteristics of Nd60Al10Fe20Co10 Bulk Metallic Glass Investigated by Ultrasonic Measurement under Pressure

An ultrasonic pulse-echo method was used to measure the transit time of longitudinal and transverse (10 MHz) elastic waves in a Nd60Al10Fe20Co10 bulk metallic glass (BMG). The measurements were carried out under hydrostatic pressure up to 0.5 GPa at room temperature. On the basis of experimental data for the sound velocities and density, the elastic moduli and Debye temperature of the BMG were derived as a function of pressure. Murnaghan's equation of state is obtained. The normal behaviour of the positive pressure dependence of the ultrasonic velocities was observed for this glass. Moreover, the compression curve, the elastic constants, and the Debye temperature of the BMG are calculated on the basis of the similarity between their physical properties in the glassy state and those in corresponding crystalline state. These results confirm qualitatively the theoretical predictions concerning the features of the microstructure and interatomic bonding in the Nd60Al10Fe20Co10 BMG.

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.

The Pressure Transient Analysis Of Deformation Of Fractal Medium

The assumption of constant rock properties in pressure-transient analysis of stress-sensitive reservoirs can cause significant errors in the estimation of temporal and spatial variation of pressure. In this article, the pressure transient response of the fractal medium in stress-sensitive reservoirs was studied by using the self-similarity solution method and the regular perturbation method. The dependence of permeability on pore pressure makes the flow equation strongly nonlinear. The nonlinearities associated with the governing equation become weaker by using the logarithm transformation. The perturbation solutions for a constant pressure production and a constant rate production of a linear-source well were obtained by using the self-similarity solution method and the regular perturbation method in an infinitely large system, and inquire into the changing rule of pressure when the fractal and deformation parameters change. The plots of typical pressure curves were given in a few cases, and the results can be applied to well test analysis.

A thermal insulation method for a piezoelectric transducer

This study deals with the sources of signal distortion of a piezoelectric transducer heated by measured gas flow. These signal distortions originate from both unloading of preload on a piezocrystal because of expansion of a diaphragm in the test apparatus and the pyroelectric effect of a heated piezoelectric crystal. A plastic film on the diaphragm of the transducer can effectively insulate the diaphragm and the piezocrystal within transducer from heating by gas flow, eliminating the sources of distortion. A method for evaluating the thickness of the film is proposed.

Compact finite difference - Fourier spectral method for three-dimensional incompressible Navier-Stokes equations

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.

Self-consistent kinetic modeling of low-pressure inductively coupled radio-frequency discharges

An efficient method for solving the spatially inhomogeneous Boltzmann equation in a two-term approximation for low-pressure inductively coupled plasmas has been developed. The electron distribution function (EDF), a function of total electron energy and two spatial coordinates, is found self-consistently with the static space-charge potential which is computed from a 2D fluid model, and the rf electric field profile which is calculated from the Maxwell equations. The EDF and the spatial distributions of the electron density, potential, temperature, ionization rate, and the inductive electric field are calculated and discussed. (C) 1996 American Institute of Physics.

Non-die explosive forming of spherical pressure-vessels

This paper presents a newly developed method of manufacturing spherical pressure vessels based on the technology of non-die explosive forming. Compared with the traditional method, this technology does not need any dies and pressing equipment, so that the cost of the production process can be greatly reduced, especially for vessels of less than 100 m3 capacity.

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

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.

Growth of Silicon Carbide Bulk Crystals by Physical Vapor Transport Method and Modeling Efforts in the Process Optimization

Silicon carbide bulk crystals were grown in an induction-heating furnace using the physical vapor transport method. Crystal growth modeling was performed to obtain the required inert gas pressure and temperatures for sufficiently large growth rates. The SiC crystals were expanded by designing a growth chamber having a positive temperature gradient along the growth interface. The obtained 6H-SiC crystals were cut into wafers and characterized by Raman scattering spectroscopy and X-ray diffraction, and the results showed that most parts of the crystals had good crystallographic structures.