Search Results

It has been known from multi-zone simulations that HCCI combustion can be significantly affected by temperature stratification of the in-cylinder gas. With the same combustion timing (i.e. crank angles at 50% heat release, denoted as CA50), large temperature stratification tends to prolong the combustion duration and lower down the in-cylinder pressure-rise-rate. With low pressure-rise-rate HCCI engines can be operated at high load, therefore it is of practical importance to look into more details about how temperature stratification affects the auto-ignition process. It has been realized that multi-zone simulations can not account for the effects of spatial structures of the stratified temperature field, i.e. how the size of the hot and cold spots in the temperature field could affect the auto-ignition process. This question is investigated in the present work by large eddy simulation (LES) method which is capable of resolving the in-cylinder turbulence field in space and time.

The paper describes a novel method for calculating the residual gas fraction and the trapping efficiency in a 2 stroke engine. Assuming one dimensional compressible flow through the inlet and exhaust ports, the method estimates the instantaneous mass flowing in and out from the combustion chamber; later the residual gas fraction and trapping efficiency are estimated combining together the perfect displacement and perfect mixing scavenging models. It is assumed that when the intake port opens, the fresh mixture is pushing out the burned charge without any mixing and after a multiple of the time needed for the largest eddy to perform one rotation, the two gasses are instantly mixed up together and expelled. The result is a very simple algorithm that does not require much computational time and is able to estimate with high level of precision the trapping efficiency and the residual gas fraction in 2 stroke engines.

In-cylinder flow measurements, conventional swirl measurements and CFD-simulations have been performed and then compared. The engine studied is a single cylinder version of a Scania D12 that represents a modern heavy-duty truck size engine. Bowditch type optical access and flat piston is used. The cylinder head was also measured in a steady-flow impulse torque swirl meter. From the two-dimensional flow-field, which was measured in the interval from -200° ATDC to 65° ATDC at two different positions from the cylinder head, calculations of the vorticity, turbulence and swirl were made. A maximum in swirl occurs at about 50° before TDC while the maximum vorticity and turbulence occurs somewhat later during the compression stroke. The swirl centre is also seen moving around and it does not coincide with the geometrical centre of the cylinder. The simulated flow-field shows similar behaviour as that seen in the measurements.

In this work, constant volume combustion is studied using a zero-dimensional FORTRAN code, which is a wide-ranging chemical kinetic simulation that allows a closed system of gases to be described on the basis of a set of initial conditions. The model provides an engine- or reactor-like environment in which the engine simulations allow for a variable system volume and heat transfer both to and from the system. The combustion chamber is divided into two zones as burned and unburned ones, which are separated by an assumed thin flame front in the combustion model used for this work. Equilibrium assumptions have been adopted for the modeling of the thermal ionization, where Saha's equation was derived for singly ionized molecules. The investigation is focused on the thermal ionization of NO as well as for other species. The outputs generated by the model are temperature profiles, species concentration profiles, ionization degree and an electron density for each zone.

This paper presents a new method to derive the tire forces for simultaneous braking and cornering, by combining empirical models for pure braking and cornering using brush-model tire mechanics. The method is aimed at simulation of vehicle handling, and is of intermediate complexity such that it may be implemented and calibrated by the end user. The brush model states that the contact patch between the tire and the road is divided into an adhesion region where the rubber is gripping the road and a sliding region where the rubber slides on the road surface. The total force generated by the tire is then composed of components from these two regions. In the proposed method the adhesion and the sliding forces are extracted from an empirical pure-slip tire model and then scaled to account for the combined-slip condition. The combined-slip self-aligning torque is described likewise.

An experimental and computational study of the near-wall combustion in a Homogeneous Charge Compression Ignition (HCCI) engine has been conducted by applying laser based diagnostic techniques in combination with numerical modeling. Our major intent was to characterize the combustion in the velocity- and thermal boundary layers. The progress of the combustion was studied by using fuel tracer LIF, the result of which was compared with LDA measurements of the velocity boundary layer along with numerical simulations of the reacting boundary layer. Time resolved images of the PLIF signal were taken and ensemble averaged images were calculated. In the fuel tracer LIF experiments, acetone was seeded into the fuel as a tracer. It is clear from the experiments that a proper set of backgrounds and laser profiles are necessary to resolve the near-wall concentration profiles, even at a qualitative level.

SI Engine knock is caused by autoignition in the unburnt part of the mixture (end-gas) ahead of the propagating flame. Autoignition of the end-gas occurs when the temperature and pressure exceeds a critical limit when comparatively slow reactions-releasing moderate amounts of heat-transform into ignition and rapid heat release. In this paper the difference in the heat released in the end-gas-by low temperature chemistry-between lean, rich, stochiometric, and stoichiometric mixtures diluted with cooled EGR was examined by measuring the temperature in the end-gas with Dual Broadband Rotational CARS. The measured temperature history was compared with an isentropic temperature calculated from the cylinder pressure trace. The experimentally obtained values for knock onset were compared with results from a two-zone thermodynamic model including detailed chemistry modeling of the end-gas reactions.

A comparison of scanning mobility particle sizer (SMPS) and laser-induced incandescence (LII) measurements of diesel particulate matter (PM) was performed. The results reveal the significance of the aggregate nature of diesel PM on interpretation of size and volume fraction measurements obtained with an SMPS, and the accuracy of primary particle size measurements by LII. Volume fraction calculations based on the mobility diameter measured by the SMPS substantially over-predict the space-filling volume fraction of the PM. Correction algorithms for the SMPS measurements, to account for the fractal nature of the aggregate morphology, result in a substantial reduction in the reported volume. The behavior of the particulate volume fraction, mean and standard deviation of the mobility diameter, and primary particle size are studied as a function of the EGR for a range of steady-state engine speeds and loads for a turbocharged direct-injection diesel engine.

This paper presents a large eddy simulation study of the liquid spray mixing with hot ambient gas in a constant volume vessel under engine-like conditions with the injection pressure of 1500 bar, ambient density 22.8 kg/m₃, ambient temperature of 900 K and an injector nozzle of 0.09 mm. The simulation results are compared with the experiments carried out by Pickett et al., under similar conditions. Under modern direct injection diesel engine conditions, it has been argued that the liquid core region is small and the droplets after atomization are fine so that the process of spray evaporation and mixing with the air is controlled by the heat and mass transfer between the ambient hot gas and central fuel flow. To examine this hypothesis a simple spray breakup model is tested in the present LES simulation. The simulations are performed using an open source compressible flow solver, in OpenFOAM.

Gasoline direct injection engines can greatly improve the fuel economy, but the idea mixture distribution cannot be easily controlled. In this paper, the linearized instability sheet atomization (LISA) and large eddy simulation (LES) implemented into KIVA-3V code were used to study the gasoline hollow cone spray process for gasoline direct injection (GDI) in a constant volume vessel. The three-dimensional results show that the LISA model can effectively simulate the gasoline hollow cone spray and obtain the string structure compared to the experiment data. And the velocity interpolation method can reduce the grid dependency of spray simulation. Using dense grid (about 8 million cells) in LES and RANS all can obtain the good spray tip penetration and width. Unlike diesel spray, for gasoline spray there are not big difference between the results using LES and RANS. In additional the ambient pressure significantly influence the gasoline spray shape.

In this work, large eddy simulation (LES) with a K-equation subgrid turbulent kinetic energy model is implemented into the CFD code KIVA3V to study the features of liquid fuel spray and combustion using gradually varying grid in a constant volume chamber. The characteristic time-scale combustion model (CTC) incorporating a turbulent timescale is adopted to predict the combustion process and the SHELL auto-ignition model is used to predict auto-ignition. Combustion is also simulated using Parallel Detailed Chemistry with Lu's n-heptane reduced mechanism (58 species), which has been added into the KIVA3V-LES code. The computational results are compared with Sandia experimental data for non-reacting and reacting cases. As a result, LES can capture the complex structure of the spray and temperature distribution as well as the trend of ignition delay and flame lift-off length variations. Better results are obtained using the Parallel Detailed Chemistry than the CTC model.

In this paper a fast NOx model is presented which can be used for engine optimization, aftertreatment control or virtual mapping. A cylinder pressure trace is required as input data. High calculation speed is obtained by using table interpolation to calculate equilibrium temperatures and species concentrations. Test data from a single-cylinder engine and from a complete six-cylinder engine have been used for calibration and validation of the model. The model produces results of good agreement with emission measurements using approximately 50 combustion product zones and a calculation time of one second per engine cycle. Different compression ratios, EGR rates, injection timing, inlet pressures etc. were used in the validation tests.

Early injection timing is an effective measure of pre-mixture formation for diesel low-temperature combustion. Three algebraic subgrid models (Smagorinsky model, dynamic Smagorinsky model and WALE model) and one-equation kinetic energy turbulent model using modified TAB breakup model (MTAB model) have been implemented into KIVA3V code to make a detailed large eddy simulation of the atomization and evaporation processes of early injection timing in a constant volume chamber and a Ford high-speed direct-injection diesel engine. The results show that the predictive vapor mass fraction and liquid penetration using LES is in good agreement with the experiment results. In combustion chamber, the sub-grid turbulent kinetic energy and viscosity using LES are less than with the RANS models, and following the increasing time, the sub-grid turbulent kinetic energy and viscosity also increase and are concentrated on the spray area.

A comprehensive computer simulation has been developed to predict the transient performance of a rotary engine. An adiabatic throttling process has been assumed to account for crevice leakage. The simulation has been tuned to fit experimental data obtained from tests upon a NORTON NR 601 engine. The potential of the model as an aid to engine design and development is illustrated by means of predictions of the effects of changes in some of the design parameters.

A comprehensive computer simulation has been developed to predict the transient performance of a spark ignition engine. The model has been tuned to fit experimental data obtained from tests upon a FORD SIERRA engine. The potential of the model as an aid to engine design and development is illustrated by means of predictions of the effects of changes in moment of inertia and ignition advance.

In this paper the combustion chamber wall temperature was measured by the use of thermographic phosphor. The temperature was monitored over a large time window covering a load transient. Wall temperature measurement provide helpful information in all engines. This temperature is for example needed when calculating heat losses to the walls. Most important is however the effect of the wall temperature on combustion. The walls can not heat up instantaneously and the slowly increasing wall temperature following a load transient will affect the combustion events sucseeding the transient. The HCCI combustion process is, due to its dependence on chemical kinetics more sensitive to wall temperature than Otto or Diesel engines. In depth knowledge about transient wall temperature could increase the understanding of transient HCCI control. A “black box” state space model was derived which is useful when predicting transient wall temperature.

An air hybrid is a vehicle with an ICE modified to also work as an air compressor and air motor. The engine is connected to two air reservoirs, normally the atmosphere and a high pressure tank. The main benefit of such a system is the possibility to make use of the kinetic energy of the vehicle otherwise lost when braking. The main difference between the air hybrid developed in this paper and earlier air hybrid concepts is the introduction of a pressure tank that substitutes the atmosphere as supplier of low air pressure. By this modification, a very high torque can be achieved in compressor mode as well as in air motor mode. A model of an air hybrid with two air tanks was created using the engine simulation code GT-Power. The results from the simulations were combined with a driving cycle to estimate the reduction in fuel consumption.

As part of on-going research into ground effect aerodynamics at the University of Southampton, attempts have been made to shed light on variables that may influence the characteristic shape of a typical multi element airfoil downforce curve while varying ride height. To achieve the stated goal, a commercial CFD software package was used to perform a comparative aerodynamic analysis study. The height of a double element airfoil above the ground was varied, while the values of lift and drag obtained were recorded to provide baseline information. The angle of attack of the flap and the main element were then changed in order to discern any effects on the lift curve. Also investigated was the effect that the relative sizes of the main element and flap had on the lift and drag curves, since modern racing car wings vary in this manner across their span.

This paper reports an one–dimensional modeling procedure of the hot spot autoignition with a detailed chemistry and multi–species transport in the end gas in an SI engine. The governing equations for continuity of mass, momentum, energy and species for an one–dimensional, unsteady, compressible, laminar, reacting flow and thermal fields are discretized and solved by a fully implicit method. A chemical kinetic mechanism is used for the primary reference fuels n–heptane and iso–octane. This mechanism contains 510 chemical reactions and 75 species. The change of the cylinder pressure is calculated from both flame propagation and piston movement. The turbulent velocity of the propagating flame is modeled by the Wiebe function. Adiabatic conditions, calculated by minimizing Gibb's free energy at each time step, are assumed behind the flame front in the burned gas.