Search Results

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.

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.

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.

The Homogeneous Charge Compression Ignition (HCCI) is the third alternative for combustion in the reciprocating engine. Here, a homogeneous charge is used as in a spark ignited engine, but the charge is compressed to auto-ignition as in a diesel. The main difference compared with the Spark Ignition (SI) engine is the lack of flame propagation and hence the independence from turbulence. Compared with the diesel engine, HCCI has a homogeneous charge and hence no problems associated with soot and NOX formation. Earlier research on HCCI showed high efficiency and very low amounts of NOX, but HC and CO were higher than in SI mode. It was not possible to achieve high IMEP values with HCCI, the limit being 5 bar. Supercharging is one way to dramatically increase IMEP. The influence of supercharging on HCCI was therefore experimentally investigated. Three different fuels were used during the experiments: iso-octane, ethanol and natural gas.

Cycle-resolved end-gas temperatures were measured using dual-broadband rotational CARS in a single-cylinder spark-ignition engine. Simultaneous cylinder pressure measurements were used as an indicator for knock and as input data to numerical calculations. The chemical processes in the end-gas have been analysed with a detailed kinetic mechanism for mixtures of iso-octane and n-heptane at different Research Octane Numbers (RON'S). The end-gas is modelled as a homogeneous reactor that is compressed or expanded by the piston movement and the flame propagation in the cylinder. The calculated temperatures are in agreement with the temperatures evaluated from CARS measurements. It is found that calculations with different RON'S of the fuel lead to different levels of radical concentrations in the end-gas. The apperance of the first stage of the autoignition process is marginally influenced by the RON, while the ignition delay of the second stage is increased with increasing RON.

A detailed analysis of the end-gas temperature and pressure in gasoline engines has been performed. This analysis leads to a simplified zero-dimensional model, that considers both, the compression and the expansion of the end-gas by the piston movement, and the compression by the flame front. If autoignition occurs in the end-gas the sudden rise of the pressure and the heat release is calculated. The rate form of the first law of thermodynamics for a control volume combined with the mass conservation equation for an unsteady and a uniform-flow process are applied. The heat of formation in the end-gas due to the chemical activity has been taken into account. In addition, a chemical kinetic model has been applied in order to study the occurrence of autoignition and prediction of knock.

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.

The formation of exothermic centers was modeled with a Stochastic Reactor Model (SRM) to investigate their impact on HCCI combustion. By varying the exhaust valve temperature, and thus assigning more realistic wall temperatures, the formation of exothermic centers and the ignition timing was shifted in time. To be able to study the exothermic centers, their formation and their distribution, Scatter plots, standard deviation plots and Probability Density Function (PDF) plots were constructed on the basis of the data the SRM calculations provided. The standard deviation for the particle temperatures was found to be an useful indicator of the degree of homogeneity within the combustion chamber, and thus of how efficient the combustion process was. It was observed that when the standard deviation of the temperature was higher, the emissions of CO and of hydrocarbons present at the end of the closed cycle were higher.

The ability to predict cyclic variations is certainly useful in studying engine operating regimes, especially under unstable operating conditions where one single cycle may differ from another substantially and a single simulation may give rather misleading results. PDF based models such as Stochastic Reactor Models (SRM) are able to model cyclic variations, but these may be overpredicted if discretization is too coarse. The range of cyclic variations and the dependence of the ability to correctly assess their mean values on the number of cycles simulated were investigated. In most cases, the average values were assessed correctly on the basis of as few as 10 cycles, but assessing the complete range of cyclic variations could require a greater number of cycles. In studying average values, variations due too coarse discretization being employed are smaller than variations originating from changes in physical parameters, such as heat transfer and mixing parameters.

Operating the HCCI engine with dual fuels with a large difference in auto-ignition characteristics (octane number) is one way to control the HCCI operation. The effect of octane number on combustion, emissions and engine performance in a 6 cylinder SCANIA truck engine, fuelled with n-heptane and isooctane, and running in HCCI mode, are investigated numerically and compared with measurements taken from Olsson et al. [SAE 2000-01-2867]. To correctly simulate the HCCI engine operation, we implement a probability density function (PDF) based stochastic reactor model (including detailed chemical kinetics and accounting for inhomogeneities in composition and temperature) coupled with GT-POWER, a 1-D fluid dynamics based engine cycle simulator. Such a coupling proves to be ideal for the understanding of the combustion phenomenon as well as the gas dynamics processes intrinsic to the engine cycle.

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.

This paper focuses on the performance and efficiency of an HCCI (Homogenous Charge Compression Ignition) engine system running on natural gas or landfill gas for stationary applications. Zero dimensional modeling and simulation of the engine, turbo, inlet and exhaust manifolds and inlet air conditioner (intercooler/heater) are used to study the effect of compression ratio and exhaust turbine size on maximum mean effective pressure and efficiency. The extended Zeldovich mechanism is used to estimate NO-formation in order to determine operation limits. Detailed chemical kinetics is used to predict ignition timing. Simulation of the in-cylinder process gives a minimum λ-value of 2.4 for natural gas, regardless of compression ratio. This is restricted by the NO formation for richer mixtures. Lower compression ratios allow higher inlet pressure and hence higher load, but it also reduces indicated efficiency.

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.

Large two-stroke marine Diesel engines have special injector geometries, which differ substantially from the configurations used in most other Diesel engine applications. One of the major differences is that injector orifices are distributed in a highly non-symmetric fashion affecting the spray characteristics. Earlier investigations demonstrated the dependency of the spray morphology on the location of the spray orifice and therefore on the resulting flow conditions at the nozzle tip. Thus, spray structure is directly influenced by the flow formation within the orifice. Following recent Large Eddy Simulation resolved spray primary breakup studies, the present paper focuses on spray secondary breakup modelling of asymmetric spray structures in Euler-Lagrangian framework based on previously obtained droplet distributions of primary breakup.

The three-way catalytic converter (TWC) is the most common catalyst for gasoline engine exhaust gas after treatment. The reduction of carbon monoxide (CO), nitrogen oxides (NOx) and unburned hydrocarbons (HC) is achieved via oxidation of carbon monoxide and hydrocarbons, and reduction of nitrogen oxides. These conversion effects were simulated in previous works using single-channel approaches and detailed kinetic models. In addition to the single-channel model multiple representative catalyst channels are used in this work to take heat transfer between the channels into account. Furthermore, inlet temperature distribution is considered. Each channel is split into a user given number of cells and each cell is treated like a perfectly stirred reactor (PSR). The simulation is validated against an experimental four-stroke engine setup with emission outputs fed into a TWC.

We numerically simulate a Homogeneous Charge Compression Ignition (HCCI) engine fuelled with a blend of ethanol and diethyl ether by means of a stochastic reactor model (SRM). A 1D CFD code is employed to calculate gas flow through the engine, whilst the SRM accounts for combustion and convective heat transfer. The results of our simulations are compared to experimental measurements obtained using a Caterpillar CAT3401 single-cylinder Diesel engine modified for HCCI operation. We consider emissions of CO, CO2 and unburnt hydrocarbons as functions of the crank angle at 50% heat release. In addition, we establish the dependence of ignition timing, combustion duration, and emissions on the mixture ratio of the two fuel components. Good qualitative agreement is found between our computations and the available experimental data.

The influence of variable air-fuel ratio inside a spark ignition engine is examined by the use of an ionization sensor. The measured ion currents are used for predicting the local air-fuel ratio in the vicinity of the spark plug. In order to support the results, a theoretical analysis has been made. An instationary chemical kinetic model burning a mixture of iso-octane and n-heptane is used for the calculations. The results are used to reconstruct the crank angle resolved ion current that has been measured in an engine. This technique has been developed in order to offer a supplementary low-cost facility of controlling the air-fuel ratio within the combustion chamber of an engine.

In this work, we present an unsteady flamelet progress variable approach for diesel engine CFD combustion modeling. The progress variable is based on sensible enthalpy integrated over the flamelet and describes the transient flamelet ignition process. By using an unsteady flamelet library for the progress variable, the impact of local effects, for example variations in the turbulence field, effects of wall heat transfer etc. on the autoignition chemistry can be considered on a cell level. The coupling between the unsteady flamelet library and the transport equation for total enthalpy follows the ideas of the representative interactive flamelet approach. Since the progress variable gives a direct description of the state in the flamelet, the method can be compared to having a flamelet in each computational cell in the CFD grid.

We present a computational tool to develop an exhaust gas recirculation (EGR) - air-fuel ratio (AFR) operating range for homogeneous charge compression ignition (HCCI) engines. A single cylinder Ricardo E-6 engine running in HCCI mode, with external EGR is simulated using an improved probability density function (PDF) based engine cycle model. For a base case, the in-cylinder temperature and unburned hydrocarbon emissions predicted by the model show a satisfactory agreement with measurements [Oakley et al., SAE Paper 2001-01-3606]. Furthermore, the model is applied to develop the operating range for various combustion parameters, emissions and engine parameters with respect to the air-fuel ratio and the amount of EGR used. The model predictions agree reasonably well with the experimental results for various parameters over the entire EGR-AFR operating range thus proving the robustness of the PDF based model.

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.