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.

Combustion information from three combustion chamber geometries was analyzed: Pancake and horseshoe geometry on a single-cylinder research engine, and pentroof geometry in a turbocharged four-cylinder production engine. Four different fuels were used. In the horseshoe configuration, the cylinder pressure traces from the burnt gas and from the end-gas pocket were evaluated. It is shown that the characteristics of knock are to a large degree a function of the combustion chamber geometry and that they are influenced strongly by the transducer position. It is shown for pentroof geometry that the number of cycles required to properly describe the knock population is a function of the knock intensity. A large error potential is shown for samples smaller than about 100 - 200 consecutive cycles. Good agreement between knock description based on accelerometer data and based on pressure data was found.

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.

A single cylinder, naturally aspirated, four-stroke and camless (Otto) engine was operated in homogeneous charge compression ignition (HCCI) mode with commercial gasoline. The valve timing could be adjusted during engine operation, which made it possible to optimize the HCCI engine operation for different speed and load points in the part-load regime of a 5-cylinder 2.4 liter engine. Several tests were made with differing combinations of speed and load conditions, while varying the valve timing and the inlet manifold air pressure. Starting with conventional SI combustion, the negative valve overlap was increased until HCCI combustion was obtained. Then the influences of the equivalence ratio and the exhaust valve opening were investigated. With the engine operating on HCCI combustion, unthrottled and without preheating, the exhaust valve opening, the exhaust valve closing and the intake valve closing were optimized next.

Future requirements for emission reduction from combustion engines in ground vehicles might be met by using the HCCI combustion concept. In this study, negative valve overlap (NVO) and low lift, short duration, camshaft profiles, were used to initiate HCCI combustion by increasing the internal exhaust gas recirculation (EGR) and thus retaining sufficient thermal energy for chemical reactions to occur when a pilot injection was introduced prior to TDC, during the NVO. One of the crucial parameters to control in HCCI combustion is the combustion phasing and one way of doing this is to vary the relative ratio of fuel injected in pilot and main injections. The combustion phasing is also influenced by the total amount of fuel supplied to the engine, the combustion phasing is thus affected when the load is changed. This study focuses on the reactions that occur in the highly diluted environment during the NVO when load and pilot to main ratio are changed.

The compression ratio is an important engine design parameter. It determines to a large extend engine properties like the achievable efficiency, the heat losses from the combustion chamber and the exhaust losses. The same properties are affected by insulation of the combustion chamber. It is therefore especially important to know the compression ratio when doing experiments with thermal barrier coatings (TBC). In case of porous TBCs, the standard methods to measure the compression ratio can give wrong results. When measuring the compression ratio by volume, using a liquid, it is uncertain if the liquid fills the total porous volume of the coating. And for a thermodynamic compression ratio estimation, a model for the heat losses is needed, which is not available when doing experiments with insulation. The subject of this paper is the evaluation of an alternative method to assess the compression ratio.

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.

Heavy-duty vehicles are primarily powered by diesel fuel, emitting CO2 emissions regardless of the exhaust after-treatment system. Contrastingly, a hydrogen engine has the potential to decarbonize the transportation sector as hydrogen is a carbon free, renewable fuel. In this study, a multi-physics 1D simulation tool (GT-Power) is used to model the gas exchange process and performance prediction of a two-stroke hydrogen engine. The aim is to establish a maximum torque-level for a four-stroke hydrogen engine and then utilize different methods for two-stroke modeling to achieve similar torque by optimizing the gas exchange process. A camless engine is used as base, enabling the flexibility to utilize approximately square valve lift profiles. The preliminary step is the GT-Power model validation, which has been done using diesel and hydrogen engines (single-cylinder heavy-duty) experiments at different operating points (871 rpm, 1200 rpm, 1259 rpm, and 1508 rpm).

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.

A Spray-Guided (SG) Direct-Injection (DI) Spark-Ignition (SI) engine is widely recognized to be a promising technology capable for substantially reducing fuel consumption and carbon dioxide emissions. Accordingly, there is a strong need for developing models of some effects specific to stratified turbulent burning under conditions of elevated and rapidly varying pressure. Two such effects were addressed in the present work by performing unsteady three-dimensional URANS simulations of stratified turbulent combustion in a SG DISI engine. First, a simple method of evaluation equilibrium combustion temperature, implemented into the CFD code OpenFOAM®, was improved in order to take into account the dissociation of the combustion products. Second, stratified turbulent combustion is affected by fluctuations in mixture composition. A widely used approach to modeling this effect consists of invoking a presumed Probability Density Function (PDF) for mixture fraction f.

Projected stringent emissions legislation will make tough demands on engine development. For diesel engines, in which combustion and emissions formation are governed by the spray formation and mixing processes, fuel injection plays a major role in the future development of cleaner engines. It is therefore important to study the fundamental features of the fuel injection process. In an engine the fuel is injected at high pressure into a pressurized and hot environment of air, which causes droplet formation and fuel evaporation. The injected fuel then forms a gaseous phase surrounding the liquid phase. The amount of evaporated fuel in relation to the total amount of injected fuel is of importance for engine performance, i.e. ignition delay and mixing rate. In this paper, the fraction of evaporated fuel was determined for sprays, using different orifice diameters ranging from 0.100 mm up to 0.227 mm, with the aid of a high-pressure spray chamber.

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.

Future pressures to reduce the fuel consumption of passenger cars may require the exploitation of alternative combustion strategies for gasoline engines to replace, or use in combination with the conventional stoichiometric spark ignition (SSI) strategy. Possible options include homogeneous lean charge spark ignition (HLCSI), stratified charge spark ignition (SCSI) and homogeneous charge compression ignition (HCCI), all of which are intended to reduce pumping and thermal losses. In the work presented here four different combustion strategies were evaluated using the same engine: SSI, HLCSI, SCSI and HCCI. HLCSI was achieved by early injection and operating the engine lean, close to its stability limits. SCSI was achieved using the spray-guided technique with a centrally placed multi-hole injector and spark-plug. HCCI was achieved using a negative valve overlap to trap hot residuals and thus generate auto-ignition temperatures at the end of the compression stroke.

Experiments were performed using a Light-Duty, single-cylinder, research engine in which the emissions, fuel consumption and combustion characteristics of two Fischer-Tropsch (F-T) Diesel fuels derived from natural gas and two conventional Diesel fuels (Swedish low sulfur Diesel and European EN 590 Diesel) were compared. Due to their low aromatic contents combustion with the F-T Diesel fuels resulted in lower soot emissions than combustion with the conventional Diesel fuels. The hydrocarbon emissions were also significantly lower with F-T fuel combustion. Moreover the F-T fuels tended to yield lower CO emissions than the conventional Diesel fuels. The low emissions from the F-T Diesel fuels, and the potential for producing such fuels from biomass, are powerful reason for future interest and research in this field.

Global warming driven by “greenhouse gas” emissions is an increasingly serious concern of both the public and legislators. A potentially potent way to reduce these emissions and conserve fossil fuel resources is to use n-butanol, iso-butanol or octanol (2-ethylhexanol) from renewable sources as alternative fuels in diesel engines. The effects of adding these substances to diesel fuel were therefore tested in a single-cylinder heavy duty diesel engine operated using factory settings. These alcohols have better calorific values, flash points, lubricity, cetane numbers and solubility in diesel than shorter-chain alcohols. However, they have lower cetane numbers than diesel, so either hydrotreated vegetable oil (HVO) or Di-tertiary-butyl peroxide (DTBP) was added to the diesel-alcohol mixtures to generate blends with the same Cetane Number (CN) as diesel.

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.