Search Results

A simulation method is presented for the analysis of combustion in spark ignition (SI) engines operated at elevated exhaust gas recirculation (EGR) level and employing multiple spark plug technology. The modeling is based on a zero-dimensional (0D) stochastic reactor model for SI engines (SI-SRM). The model is built on a probability density function (PDF) approach for turbulent reactive flows that enables for detailed chemistry consideration. Calculations were carried out for one, two, and three spark plugs. Capability of the SI-SRM to simulate engines with multiple spark plug (multiple ignitions) systems has been verified by comparison to the results from a three-dimensional (3D) computational fluid dynamics (CFD) model. Numerical simulations were carried for part load operating points with 12.5%, 20%, and 25% of EGR. At high load, the engine was operated at knock limit with 0%, and 20% of EGR and different inlet valve closure timing.

The use of complex reaction schemes is accompanied by high computational cost in 3D CFD simulations but is particularly important to predict pollutant emissions in internal combustion engine simulations. One solution to tackle this problem is to solve the chemistry prior the CFD run and store the chemistry information in look-up tables. The approach presented combines pre-tabulated progress variable-based source terms for auto-ignition as well as soot and NOx source terms for emission predictions. The method is coupled to the 3D CFD code CONVERGE v2.4 via user-coding and tested over various speed and load passenger-car Diesel engine conditions. This work includes the comparison between the combustion progress variable (CPV) model and the online chemistry solver in CONVERGE 2.4. Both models are compared by means of combustion and emission parameters. A detailed n-decane/α-methyl-naphthalene mechanism, comprising 189 species, is used for both online and tabulated chemistry simulations.

This paper reports on a fast predictive combustion tool employing detailed chemistry. The model is a stochastic reactor based, discretised probability density function model, without spatial resolution. Employing detailed chemistry has the potential of predicting emissions, but generally results in very high CPU costs. Here it is shown that CPU times of a couple of minutes per cycle can be reached when applying detailed chemistry, and CPU times below 10 seconds per cycle can be reached when using reduced chemistry while still catching in-cylinder in-homogeneities. This makes the tool usable for efficient engine performance mapping and optimisation. To meet CPU time requirements, automatically load balancing parallelisation was included in the model. This allowed for an almost linear CPU speed-up with number of cores available.

This paper reports on a turbulent flame propagation model combined with a zero-dimensional two-zone stochastic reactor model (SRM) for efficient predictive SI in-cylinder combustion calculations. The SRM is a probability density function based model utilizing detailed chemistry, which allows for accurate knock prediction. The new model makes it possible to - in addition - study the effects of fuel chemistry on flame propagation, yielding a predictive tool for efficient SI in-cylinder calculations with all benefits of detailed kinetics. The turbulent flame propagation model is based on a recent analytically derived formula by Kolla et al. It was simplified to better suit SI engine modelling, while retaining the features allowing for general application. Parameters which could be assumed constant for a large spectrum of situations were replaced with a small number of user parameters, for which assumed default values were found to provide a good fit to a range of cases.

In one-dimensional engine simulation programs the simulation of engine performance is mostly done by parameter fitting in order to match simulations with experimental data. The extensive fitting procedure is especially needed for emissions formation - CO, HC, NO, soot - simulations. An alternative to this approach is, to calculate the emissions based on detailed kinetic models. This however demands that the in-cylinder combustion-flow interaction can be modeled accurately, and that the CPU time needed for the model is still acceptable. PDF based stochastic reactor models offer one possible solution. They usually introduce only one (time dependent) parameter - the mixing time - to model the influence of flow on the chemistry. They offer the prediction of the heat release, together with all emission formation, if the optimum mixing time is given.

An experimental study of the Homogeneous Charge Compression Ignition (HCCI) combustion process has been conducted by using chemiluminescence imaging. The major intent was to characterize the flame structure and its transient behavior. To achieve this, time resolved images of the naturally emitted light were taken. Emitted light was studied by recording its spectral content and applying different filters to isolate species like OH and CH. Imaging was enabled by a truck-sized engine modified for optical access. An intensified digital camera was used for the imaging. Some imaging was done using a streak-camera, capable of taking eight arbitrarily spaced pictures during a single cycle, thus visualizing the progress of the combustion process. All imaging was done with similar operating conditions and a mixture of n-heptane and iso-octane was used as fuel. Some 20 crank angles before Top Dead Center (TDC), cool flames were found to exist.

Interactions between in-cylinder combustion and emission aftertreatment need to be understood for optimizing the overall powertrain system. Numerical investigations can aid this process. For this purpose, simple and numerically fast, but still accurate models are needed for in-cylinder combustion and exhaust aftertreatment. The chemical processes must be represented in sufficient detail to predict engine power, fuel consumption, and tailpipe emission levels of NOx, soot, CO and unburned hydrocarbons. This paper reports on a new transient one-dimensional catalyst model. This model makes use of a detailed kinetic mechanism to describe the catalytic reactions. A single-channel or a set of representative channels are used in the presented approach. Each channel is discretized into a number of cells. Each cell is treated as a perfectly stirred reactor (PSR) with a thin film layer for washcoat treatment. Heat and mass transport coefficients are calculated from Nusselt and Sherwood laws.

A zero-dimensional, three-zone model is developed in order to study the gas thermodynamic characteristics and its relation to knock in SI engines. The first zone is the zone behind the flame front, i.e. the burned gas products. The second zone is the unburned gas ahead of the flame front. The end gas adjacent to the wall, in the boundary layer, is not included in the second zone but it is treated as a separate zone, i.e. the third zone. A detailed analysis of the gas thermodynamic state, including heat transfer analysis between the zones and the walls and mass transfer analysis between the zones combined with a detailed chemical kinetic mechanism in each zone have been performed. The effects of piston movement, flame propagation and transient behavior of the thermal boundary layer are modeled. A sudden rise of pressure and temperature and associated heat release in the end gas are calculated if autoignition occurs.

The solution mapping method Adaptive Polynomial Tabulation (APT) for complex chemistry is presented. The method has the potential of reducing the computational time required for stochastic reactor model simulations of the HCCI combustion process. In this method the solution of the initial value chemical rate equation system is approximated in real-time with zero, first and second order polynomial expressions. These polynomials are algebraic functions of a progress variable, pressure and total enthalpy. The chemical composition space is divided a priori into block-shaped regions (hypercubes) of the same size. Each hypercube may be divided in real-time into adaptive hypercubes of different sizes. During computations, initial conditions are stored in the adaptive hypercubes. Two concentric Ellipsoids of Accuracy (EOA) are drawn around each stored initial condition.

Today numerical models are a major part of the diesel engine development. They are applied during several stages of the development process to perform extensive parameter studies and to investigate flow and combustion phenomena in detail. The models are divided by complexity and computational costs since one has to decide what the best choice for the task is. 0D models are suitable for problems with large parameter spaces and multiple operating points, e.g. engine map simulation and parameter sweeps. Therefore, it is necessary to incorporate physical models to improve the predictive capability of these models. This work focuses on turbulence and mixing modeling within a 0D direct injection stochastic reactor model. The model is based on a probability density function approach and incorporates submodels for direct fuel injection, vaporization, heat transfer, turbulent mixing and detailed chemistry.

Mixing in wall-jets was investigated in an optical heavy-duty diesel engine with several injector configurations and injection pressures. Laser-induced fluorescence (LIF) was employed in non-reacting conditions in order to quantitatively measure local equivalence ratios in colliding wall-jets. A novel laser diagnostic technique, Structured Laser Illumination Planar Imaging (SLIPI), was successfully implemented in an optical engine and permits to differentiate LIF signal from multiply scattered light. It was used to quantitatively measure local equivalence ratio in colliding wall-jets under non-reacting conditions. Mixing phenomena in wall-jets were analyzed by comparing the equivalence ratio in the free part of the jet with that in the recirculation zone where two wall-jets collide. These results were then compared to φ predictions for free-jets. It was found that under the conditions tested, increased injection pressure did not increase mixing in the wall-jets.

The novel wave-shaped bowl piston geometry design with protrusions has been proved in previous studies to enhance late-cycle mixing and therefore significantly reduce soot emissions and increase engine thermodynamic efficiency. The wave-shaped piston is characterized by the introduction of evenly spaced protrusions around the inner wall of the bowl, with a matching number with the number of injection holes, i.e., flames. The interactions between adjacent flames strongly affect the in-cylinder flow and the wave shape is designed to guide the near-wall flow. The flow re-circulation produces a radial mixing zone (RMZ) that extends towards the center of the piston bowl, where unused air is available for oxidation promotion. The waves enhance the flow re-circulation and thus increase the mixing intensity of the RMZ.

The soot distribution as function of ambient O₂ mole fraction in a heavy-duty diesel engine was investigated at low load (6 bar IMEP) with laser-induced incandescence (LII) and natural luminosity. A Multi-YAG laser system was utilized to create time-resolved LII using 8 laser pulses with a spacing of one CAD with detection on an 8-chip framing camera. It is well known that the engine-out smoke level increases with decreasing oxygen fraction up to a certain level where it starts to decrease again. For the studied case the peak occurred at an O₂ fraction of 11.4%. When the oxygen fraction was decreased successively from 21% to 9%, the initial soot formation moved downstream in the jet. At the lower oxygen fractions, below 12%, no soot was formed until after the wall interaction. At oxygen fractions below 11% the first evidence of soot is in the recirculation zone between two adjacent jets.

A heavy duty diesel engine operating case producing no engine-out smoke was studied using combined simultaneous optical diagnostics. The case was close to a typical low load modern diesel operating point without EGR. Parallels were drawn to the conceptual model by Dec and results from high-pressure combustion vessels. Optical results revealed that no soot was present in the upstream part of the jet cross-section. Soot was only observed in the recirculation zones close to the bowl perimeter. This indicated very slow soot formation and was explained by a significantly higher air entrainment rate than in Dec's study. The local fuel-air equivalence ratio, Φ, at the lift-off length was estimated to be 40% of the value in Dec's study. The lower Φ in the jet produced a different Φ -T-history, explaining the soot results. The increased air entrainment rate was mainly due to smaller nozzle holes and increased TDC density.

The development of future gasoline engines is dominated by the study of new technologies aimed at reducing the engine negative environmental impact and increase its thermal efficiency. One common trend is to develop smaller engines able to operate in stoichiometric conditions across the whole engine map for better efficiency, lower fuel consumption, and optimal conversion rate of the three-way catalyst (TWC). Water injection is one promising technique, as it significantly reduces the engine knock tendency and avoids fuel enrichment for exhaust temperature mitigation at high power operation. With the focus on reducing the carbon footprint of the automotive sector, another vital topic of research is the investigation of new alternative CO2-neutral fuels or so-called eFuels. Several studies have already shown how these new synthetic fuels can be produced by exploiting renewable energy sources and can significantly reduce engine emissions.

The local equivalence ratio, Φ, was measured in fuel jets using laser-induced spontaneous Raman scattering in an optical heavy duty diesel engine. The measurements were performed at 1200 rpm and quarter load (6 bar IMEP). The objective was to study factors influencing soot formation, such as gas entrainment and lift-off position, and to find correlations with engine-out particulate matter (PM) levels. The effects of nozzle hole size, injection pressure, inlet oxygen concentration, and ambient density at TDC were studied. The position of the lift–off region was determined from OH chemiluminescence images of the flame. The liquid penetration length was measured with Mie scattering to ensure that the Raman measurement was performed in the gaseous part of the spray. The local Φ value was successfully measured inside a fuel jet. A surprisingly low correlation coefficient between engine-out PM and the local Φ in the reaction zone were observed.

Water injection can be applied to spark ignited gasoline engines to increase the Knock Limit Spark Advance and improve the thermal efficiency. The Knock Limit Spark Advance potential of 6 °CA to 11 °CA is shown by many research groups for EN228 gasoline fuel using experimental and simulation methods. The influence of water is multi-layered since it reduces the in-cylinder temperature by vaporization and higher heat capacity of the fresh gas, it changes the chemical equilibrium in the end gas and increases the ignition delay and decreases the laminar flame speed. The aim of this work is to extend the analysis of water addition to different octane ratings. The simulation method used for the analysis consists of a detailed reaction scheme for gasoline fuels, the Quasi-Dimensional Stochastic Reactor Model and the Detonation Diagram. The detailed reaction scheme is used to create the dual fuel laminar flame speed and combustion chemistry look-up tables.

A method for automatic reduction of detailed reaction mechanisms using simultaneous sensitivity, reaction flow and lifetime analysis has been developed and applied to a two-zone model of an SI engine fuelled with Primary Reference Fuel (PRF). Species which are less relevant for the occurrence of autoignition in the end gas are declared redundant. They are identified and eliminated for different pre-set minimum levels of reaction flow and sensitivity. The resulting skeletal mechanism is valid in the ranges of initial and boundary values for which the analyses have been performed. A measure of species lifetime is calculated from the chemical source terms, and the species with the lifetime shorter than and mass-fraction less than specified limits are selected for removal.

Today’s car manufactures inevitably have to focus on the reduction of fuel consumption while maintaining high performance standards. In this respect, the downsized turbocharged DISI (Direct Injection Spark Ignition) engine represents an appealing solution. However, downsizing is limited because of knocking phenomena occurring at high- and full-load conditions due to autoignition of the unburned mixture ahead the flame front. A common way of reducing knock tendencies is provided by Exhaust Gas Recirculation (EGR). However, EGR modifies the chemical composition of the cylinder charge and recirculated species like nitric oxide (NO) or unburned Hydrocarbons (HC) particularly increase the reactivity of the unburned mixture. In other words, the EGR influences the Octane Number (ON) of the in-cylinder gases.

Laser-Induced Incandescence (LII) has traditionally been considered a straightforward and reliable optical diagnostic technique for in-cylinder soot measurements. As a result, it is nowadays even possible to buy turn-key LII measurement systems. During recent years, however, attention has been drawn to a number of unresolved challenges with LII. Many of these are relevant mostly for particle sizing using time-resolved LII, but also two-dimensional soot volume fraction measurements are affected, especially in regions with high soot concentrations typically found in combustion engines. In this work the focus is on the specific challenges involved in performing high-repetition rate measurements with LII in diesel engines. All the mentioned issues might not be possible to overcome but they should nevertheless be known and their potential impact should be considered.