A detailed chemistry multi-cycle simulation of a gasoline fueled HCCI engine operated with NVO

A previously developed Stochastic Reactor Model (SRM) is used to simulate combustion in a four cylinder in-line four-stroke naturally aspirated direct injection Spark Ignition (SI) engine modified to run in Homogeneous Charge Compression Ignition (HCCI) mode with a Negative Valve Overlap (NVO). A portion of the fuel is injected during NVO to increase the cylinder temperature and enable HCCI combustion at a compression ratio of 12:1. The model is coupled with GT-Power, a one-dimensional engine simulation tool used for the open valve portion of the engine cycle. The SRM is used to model in-cylinder mixing, heat transfer and chemistry during the NVO and main combustion. Direct injection is simulated during NVO in order to predict heat release and internal Exhaust Gas Recycle (EGR) composition and mass. The NOx emissions and simulated pressure profiles match experimental data well, including the cyclic fluctuations. The model predicts combustion characteristics at different fuel split ratios and injection timings. The effect of fuel reforming on ignition timing is investigated along with the causes of cycle to cycle variations and unstable operation. A detailed flux analysis during NVO unearths interesting results regarding the effect of NOx on ignition timing compared with its effect during the main combustion. © 2009 SAE International.

Convective heat dissipation with lattice-frame materials

This paper presents experimental results on heat transfer and pressure drop for a compact heat sink made of fully triangulated, lightweight (porosity∼0.938), aluminum lattice-frame materials (LFMs). Due to the inherent structural anisotropy of the LFMs, two mutually perpendicular orientations were selected for the measurements. Constant heat flux was applied to the heat sink under steady state conditions, and dissipated by forced air convection. The experimental data were compared with those predicted from an analytical model based on fin analogy. The experimental results revealed that pressure drop is strongly dependent upon the orientation of the structure, due mainly to the flow blockage effect. For heat transfer measurements, typical local temperature distributions on the substrate under constant heat flux conditions were captured with infrared camera. The thermal behavior of LFMs was found to follow closely that of cylinder banks, with early transition Reynolds number (based on strut diameter) equal to about 300. The Nusselt number prediction from the fin-analogy correlates well with experimental measurements, except at low Reynolds numbers where a slightly underestimation is observed. Comparisons with empty channels and commonly used heat exchanger media show that the present LFM heat sink can remove heat approximately seven times more efficient than an empty channel and as efficient as a bank of cylinders at the same porosity level. The aluminum LFMs are extremely stiff and strong, making them ideal candidates for multifunctional structures requiring both heat dissipation and mechanical load carrying capabilities. © 2003 Elsevier Ltd. All rights reserved.

The ensemble statistics of the energy of a random system subjected to harmonic excitation

This paper is concerned with the ensemble statistics of the response to harmonic excitation of a single dynamic system such as a plate or an acoustic volume. Random point process theory is employed, and various statistical assumptions regarding the system natural frequencies are compared, namely: (i) Poisson natural frequency spacings, (ii) statistically independent Rayleigh natural frequency spacings, and (iii) natural frequency spacings conforming to the Gaussian orthogonal ensemble (GOE). The GOE is found to be the most realistic assumption, and simple formulae are derived for the variance of the energy of the system under either point loading or rain-on-the-roof excitation. The theoretical results are compared favourably with numerical simulations and experimental data for the case of a mass loaded plate. © 2003 Elsevier Ltd. All rights reserved.

Mechanical and magnetic properties of metal fibre networks, with and without a polymeric matrix

Bonded networks of metal fibres are highly porous, permeable materials, which often exhibit relatively high strength. Material of this type has been produced, using melt-extracted ferritic stainless steel fibres, and characterised in terms of fibre volume fraction, fibre segment (joint-to-joint) length and fibre orientation distribution. Young's moduli and yield stresses have been measured. The behaviour when subjected to a magnetic field has also been investigated. This causes macroscopic straining, as the individual fibres become magnetised and tend to align with the applied field. The modeling approach of Markaki and Clyne, recently developed for prediction of the mechanical and magneto-mechanical properties of such materials, is briefly summarised and comparisons are made with experimental data. The effects of filling the inter-fibre void with compliant (polymeric) matrices have also been explored. In general the modeling approach gives reliable predictions, particularly when the network architecture has been characterised using X-ray tomography. © 2005 Published by Elsevier Ltd.

Modelling soot formation in a DISI engine

In this work, the formation of soot in a Direct Injection Spark Ignition (DISI) engine is simulated using the Stochastic Reactor Model (SRM) engine code. Volume change, convective heat transfer, turbulent mixing, direct injection and flame propagation are accounted for. In order to simulate flame propagation, the cylinder is divided into an unburned, entrained and burned zone, with the rate of entrainment being governed by empirical equations but combustion modelled with chemical kinetics. The model contains a detailed chemical mechanism as well as a highly detailed soot formation model, however computation times are relatively short. The soot model provides information on the morphology and chemical composition of soot aggregates along with bulk quantities, including soot mass, number density, volume fraction and surface area. The model is first calibrated by simulating experimental data from a Gasoline Direct Injection (GDI) Spark Ignition (SI) engine. The model is then used to simulate experimental data from the literature, where the numbers, sizes and derived mass particulate emissions from a 1.83 L, 4-cylinder, 4 valve production DISI engine were examined. Experimental results from different injection and spark timings are compared with the model and the qualitative trends in aggregate size distribution and emissions match the exhaust gas measurements well. © 2010 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Study of lateral mode SOI-MEMS resonators for reduced anchor loss

MEMS resonators fabricated in silicon-on-insulator (SOI) technology must be clamped to the substrate via anchoring stems connected either from within the resonator or through the sides, with the side-clamped solution often employed due to manufacturing constraints. This paper examines the effect of two types of commonly used side-clamped, anchoring-stem geometries on the quality factor of three different laterally-driven resonator topologies. This study employs an analytical framework which considers the relative distribution of strain energies between the resonating body and clamping stems. The ratios of the strain energies are computed using ANSYS FEA and used to provide an indicator of the expected anchor-limited quality factors. Three MEMS resonator topologies have been fabricated and characterized in moderate vacuum. The associated measured quality factors are compared against the computed strain energy ratios, and the trends are shown to agree well with the experimental data. © 2011 IOP Publishing Ltd.

Dynamics of collapsing fluidized beds and its application in the simulation of pulsed fluidized beds

A new method to study collapsing process of fluidized bed was proposed. The method is based on the analysis of the pressure variation during collapse. A model is proposed to describe the pressure variation on any location of the bed during collapse. Three kind of particles were fluidized by air and the pressure variation after shutting down the gas supply was measured by pressure transducers and simulated by the proposed model. The simulated results were in good agreement with the experimental data. The parameters of the fluidized bed, such as bubble fraction alpha(b), particle fraction in dense phase alpha(p), bubble velocity u(b) and slip velocity of particle u(i), can be estimated from the measured pressure data by using the proposed model. The model was also used to calculate the pressure variation in pulsed fluidized beds with low frequency.

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.

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.

How important are shock waves to single-bubble sonoluminescence?

By solving numerically the full set of hydrodynamic equations governing the pulsation of a bubble,we show that shock waves are often absent in a stable sonoluminescing bubble. Nevertheless, for a wide range of physical parameters, a continuous compressional wave emerges and heats up the bubble, and the resulting black-body radiations have pulse heights and widths that agree with experimental data. Shock waves, being much less robust, are not essential for stable single-bubble sonoluminescence.

Study Of The Damage-Induced Anisotropy Of Quasi-Brittle Materials Using The Component Assembling Model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Characteristics Of Energy Dissipation In Hyperconcentrated Flows

An equilibrium equation for the turbulence energy in sediment-laden flows was derived on the basis of solid-liquid two-phase flow theory. The equation was simplified for two-dimensional, uniform, steady and fully developed turbulent hyperconcentrated flows. An energy efficiency coefficient of suspended-load motion was obtained from the turbulence energy equation, which is defined as the ratio of the sediment suspension energy to the turbulence energy of the sediment-laden flows. Laboratory experiments were conducted to investigate the characteristics of energy dissipation in hyperconcentrated flows. A total of 115 experimental runs were carried out, comprising 70 runs with natural sediments and 45 runs with cinder powder. Effects of sediment concentration on sediment suspension energy and flow resistance were analyzed and the relation between the energy efficiency coefficient of suspended-load motion and sediment concentration was established on the basis of experimental data. Furthermore, the characteristics of energy dissipation in hyperconcentrated flows were identified and described. It was found that the high sediment concentration does not increase the energy dissipation; on the contrary, it decreases flow resistance.

High Aspect Ratio (L/D) Riser Viv Prediction Using Wake Oscillator Model

A two-dimensional (2-D) vortex-induced vibration (VIV) prediction model for high aspect ratio (LID) riser subjected to uniform and sheared flow is studied in this paper. The nonlinear structure equations are considered. The near wake dynamics describing the fluctuating nature of vortex shedding is modeled using classical van der Pol equation. A new approach was applied to calibrate the empirical parameters in the wake oscillator model. Compared the predicted results with the experimental data and computational fluid dynamic (CFD) results. Good agreements are observed. It can be concluded that the present model can be used as simple computational tool in predicting some aspects of VIV of long flexible structures. (C) 2008 Elsevier Ltd. All rights reserved.

The Head-On Colliding Process Of Binary Liquid Droplets At Low Velocity: High-Speed Photography Experiments And Modeling

The experimental and theoretical studies are reported in this paper for the head-on collisions of a liquid droplet with another of the same fluid resting on a solid substrate. The droplet on the hydrophobic polydimethylsiloxane (PDMS) substrate remains in a shape of an approximately spherical segment and is isometric to an incoming droplet. The colliding process of the binary droplets was recorded with high-speed photography. Head-on collisions saw four different types of response in our experiments: complete rebound, coalescence, partial rebound With conglutination, and coalescence accompanied by conglutination. For a complete rebound, both droplets exhibited remarkable elasticity and the contact time of the two colliding droplets was found to be in the range of 10-20 ms. With both droplets approximately considered as elastic bodies, Hertz contact theory was introduced to estimate the contact time for the complete rebound case. The estimated result Was found to be on the same order of magnitude as the experimental data, which indicates that the present model is reasonable. (C) 2008 Elsevier Inc. All rights reserved.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.