Search Results

Recent research on thermal reciprocating engines has focused on the influence of lubricant oil on the combustion process, which can lead to highly undesired super-knock events. Low-Speed Pre-Ignition (LSPI) events severely limit the further development of Direct Injection Spark Ignition Engines (DISI), preventing high efficiencies from being achieved. However, there is still a lack of knowledge about the fundamental mechanisms leading to LSPI, due to the complex phenomena involved and the interaction between lubricant oil and fuel. Understanding how the presence of lubricant oil traces affects gasoline chemical reactivity is an essential step for performing successful numerical simulations aimed at predicting the onset of LSPI phenomena. Reaction mechanisms able to predict oil-fuel interaction have been proposed, but they are computationally demanding.

Computational fluid dynamics of gas-fueled large-bore spark ignition engines with pre-chamber ignition can speed up the design process of these engines provided that 1) the reliability of the results is not affected by poor meshing and 2) the time cost of the meshing process does not negatively compensate for the advantages of running a computer simulation. In this work a flame propagation model that runs with arbitrary hybrid meshes was developed and coupled with the KIVA4-MHI CFD solver, in order to address these aims. The solver follows the G-Equation level-set method for turbulent flame propagation by Tan and Reitz, and employs improved numerics to handle meshes featuring different cell types such as hexahedra, tetrahedra, square pyramids and triangular prisms. Detailed reaction kinetics from the SpeedCHEM solver are used to compute the non-equilibrium composition evolution downstream and upstream of the flame surface, where chemical equilibrium is instead assumed.

Improvements to combustion models for modeling spark ignition engines using the G-equation flame propagation model and detailed chemical kinetics have been performed. The improvements include revision of a PRF chemistry mechanism, precise calculation of “primary heat release” based on the sub-grid scale unburned/burnt volumes of flame-containing cells, modeling flame front quenching in highly stratified mixtures, introduction of a Damkohler model for assessing the combustion regime of flame-containing cells, and a better method of modeling the effects of the local residual value on the burning velocity. The validation of the revised PRF mechanism shows that the calculated ignition delay matches shock tube data very well. The improvements to the “primary heat release” model based on the cell unburned/burnt volumes more precisely consider the chemical kinetics heat release in unburned regions, and thus are thought to be physically reasonable.

Fuel spray atomization and breakup processes within a direct-injection spark-ignition (DISI) engine and outside the engine were modeled using a modified KIVA-3V code with improved spray models. The structures of the predicted sprays were qualitatively compared with planar images. The considered sprays were created by a prototype pressure-swirl injector and the planar images were obtained by laser sheet imaging in an optical DISI engine. In the out-of-engine case, the spray was injected into atmospheric air, and was modeled in a two dimensional bomb. In the engine case, the injection started from 270° ATDC, and full 3-D computations in the same engine were performed. In both cases, two liquid injection pressure conditions were applied, that is, 3.40 MPa and 6.12 MPa. The model gives good prediction of the tip penetration, and external spray shape, but the internal structure prediction has relatively lower accuracy, especially near the spray axis.

In this paper, knock in a Ford single cylinder direct-injection spark-ignition (DISI) engine was modeled and investigated using the KIVA-3V code with a G-equation combustion model coupled with detailed chemical kinetics. The deflagrative turbulent flame propagation was described by the G-equation combustion model. A 22-species, 42-reaction iso-octane (iC8H18) mechanism was adopted to model the auto-ignition process of the gasoline/air/residual-gas mixture ahead of the flame front. The iso-octane mechanism was originally validated by ignition delay tests in a rapid compression machine. In this study, the mechanism was tested by comparing the simulated ignition delay time in a constant volume mesh with the values measured in a shock tube under different initial temperature, pressure and equivalence ratio conditions, and acceptable agreements were obtained.

Homogeneous charge compression ignition (HCCI) combustion is considered to be an attractive alternative to traditional internal combustion engine operation because of its extremely low levels of pollutant emissions. However, there are several difficulties that must be overcome for HCCI practical use, such as difficult ignition timing controllability. Indeed, too early or too late ignition can occur with obvious drawbacks. In addition, the increase in cyclic variation caused by the ignition timing uncertainty can lead to uneven engine operation. As a way to solve the combustion phasing control problem, dual-fuel combustion has been proposed. It consists of a diesel pilot injection used to ignite a pre-mixture of gasoline (or other high octane fuel) and air. Although dual-fuel combustion is an attractive way to achieve controllable HCCI operation, few studies are available to help the understanding of its in-cylinder combustion behavior.

A continuous multi-component fuel evaporation model has been integrated with an improved G-equation combustion and detailed chemical kinetics model. The integrated code has been successfully used to simulate a gasoline direct injection engine. In the multi-component fuel model, the theory of continuous thermodynamics is used to model the properties and composition of multi-component fuels such as gasoline. In the improved G-equation combustion model a flamelet approach based on the G-equation is used that considers multi-component fuel effects. To precisely calculate the local and instantaneous residual which has a great effect on the laminar flame speed, a “transport equation residual” model is used. A Damkohler number criterion is used to determine the combustion mode in flame containing cells.

An experimental study was conducted on an air cooled high-speed, direct-injection diesel generator that investigated the use of gasoline in a dual fuel PCCI strategy. The single-speed generator used in the study has an effective compression ratio of 17 and runs at 3600 rev/min. Varying amounts of gasoline were introduced into the combustion chamber through a port injection system. The generator uses an all-mechanical diesel fuel injection system that has a fixed injection timing. The experiments explored variable intake temperatures and fuel split quantities to investigate different combustion phasing regimes. Results from the study showed low combustion efficiency at low load. Low load operation was also characterized by high levels of HC and CO (in excess of 20 g/kwh and 50 g/kwh respectively). Operation at 75% load was more efficient, and displayed three different combustion regimes that are possible with PIG (port injected gasoline) dual fuel PCCI.

A numerical study was performed to compare the formation of nitric oxide (NO) and nitrogen dioxide (NO₂), collectively termed NOx, resulting from biodiesel and diesel combustion in an internal combustion engine. It has been shown that biodiesel tends to increase NOx compared to diesel, and to-date, there is no widely accepted explanation. Many factors can lead to increased NOx formation and it was of interest to determine if fuel chemistry plays a significant role. Therefore, in order to isolate the fuel chemistry from mixing processes typical in a compression ignition engine, sprays were not considered in the present investigation. The current study compares the NOx formation of surrogates for biodiesel (as represented by methyl butanoate and n-heptane) and diesel (n-heptane) under completely homogeneous conditions. Combustion of each fuel was simulated using the Senkin code for both an adiabatic, constant volume reactor, and an adiabatic, single-zone HCCI engine model.

Experiments show that intake flow details have a significant influence on High-Speed Direct-Injection (HSDI) diesel engine soot emissions. Four different intake modes were simulated using the combination of the CFD codes, STAR-CD and KIVA-3V, to investigate spray-intake flow-emission interaction characteristics. The simulation results were compared to steady-state flow bench data and engine experimental data. It was found that it is difficult to accurately predict the timing of the small pilot and main combustion events, simultaneously, with current simplified ignition models. NOx emissions were predicted well, however, an insensitivity of the soot emissions to the details of the intake process was found, mainly due to the deficiencies in predicting the ignition delay. The results show that a strong swirling flow causes the formed soot to remain within the bowl, leading to high soot emissions.

A spray model for pressure-swirl atomizers that is based on a linearized instability analysis of liquid sheets has been combined with an ignition and combustion model for stratified charge spark ignition engines. The ignition model has been advanced, such that the presence of dual spark plugs can now be accounted for. Independent validation of the spray model is achieved by investigating a pressure-swirl injector inside a pressure bomb containing air at ambient temperature. In a second step, the complete model is used to estimate the performance of a small marine DISI Two-Stroke engine operating in stratified charge mode. Simulation results and experimental data are compared for several different injection timings and the agreement is generally good such that there is confidence in the predictive quality of the combustion model. Finally the model is applied in a conceptual study to investigate possible benefits of split injections.

Detailed chemical kinetics was implemented in the KIVA-3V multidimensional CFD code to study the combustion process in Homogeneous Charge Compression Ignition (HCCI) engines. The CHEMKIN code was implemented such that the chemistry and flow solutions were coupled. Detailed reaction mechanisms were used to simulate the fuel chemistry of ignition and combustion. Effects of turbulent mixing on the reaction rates were also considered. The model was validated using the experimental data from two modified heavy-duty diesel engines, including a Volvo engine and a Caterpillar engine operated at the HCCI mode. The results show that good levels of agreement were obtained using the present KIVA/CHEMKIN model for a wide range of engine conditions, including various fuels, injection systems, engine speeds, and EGR levels. Ignition timings were predicted well without the need to adjust any kinetic constants.

The present study investigated HCCI combustion in a heavy-duty diesel engine both experimentally and numerically. The engine was equipped with a hollow-cone pressure-swirl injector using gasoline direct injection. Characteristics of HCCI combustion were obtained by very early injection with a heated intake charge. Experimental results showed an increase in NOx emission and a decrease in UHC as the injection timing was retarded. It was also found that optimization can be achieved by controlling the intake temperature together with the start-of-injection timing. The experiments were modeled by using an engine CFD code with detailed chemistry. The CHEMKIN code was implemented into KIVA-3V such that the chemistry and flow solutions were coupled. The model predicted ignition timing, cylinder pressure, and heat release rates reasonably well. The NOx emissions were found to increase as the injection timing was retarded, in agreement with experimental results.

An improved discrete particle ignition kernel (DPIK) model and the G-equation combustion model have been developed and implemented in KIVA-3V. In the ignition model, the spark ignition kernel growth is tracked by Lagrangian markers and the spark discharge energy and flow turbulence effects on the ignition kernel growth are considered. The predicted ignition kernel size was compared with the available measurements and good agreement was obtained. Once the ignition kernel grows to a size where the turbulent flame is fully developed, the level set method (G-equation) is used to track the mean turbulent flame propagation. It is shown that, by ignoring the detailed turbulent flame brush structure, fine numerical resolution is not needed, thus making the models suitable for use in multidimensional modeling of SI engine combustion.

The present study uses a numerical model to investigate the effects of flow turbulence on premixed iso-octane HCCI engine combustion. Different levels of in-cylinder turbulence are generated by using different piston geometries, namely a disc-shape versus a square-shape bowl. The numerical model is based on the KIVA code which is modified to use CHEMKIN as the chemistry solver. A detailed reaction mechanism is used to simulate the fuel chemistry. It is found that turbulence has significant effects on HCCI combustion. In the current engine setup, the main effect of turbulence is to affect the wall heat transfer, and hence to change the mixture temperature which, in turn, influences the ignition timing and combustion duration. The model also predicts that the combustion duration in the square bowl case is longer than that in the disc piston case which agrees with the measurements.

A level set G-equation model has been developed to model the combustion process in spark ignition engines. The spark ignition process was modeled using an improved version of the Discrete Particle Ignition Kernel (DPIK) model. The two models were implemented into the KIVA-3V code to simulate SI engine combustion under both premixed and direct injection conditions. In the ignition model, the ignition kernel growth is tracked by Lagrangian markers, and spark discharge energy and flow turbulence effects on the kernel growth are considered. Once the ignition kernel grows to a size where the turbulent flame is fully developed, the G-equation model is used to track the mean turbulent flame evolution. When combined with a characteristic time scale combustion model, the models were also used to simulate stratified combustion in DISI engines, where the triple flame structure must be considered.

The combustion mechanism in a Compression Ignition Homogeneous Charge (CIHC) engine was studied. Previous experiments done on a four-stroke CIHC engine were modeled using the KIVA-II code with modifications to the combustion, heat transfer, and crevice flow submodels. A laminar and turbulence characteristic time combustion model that has been used for spark-ignited engine studies was extended to allow predictions of ignition. The rate of conversion from one chemical species to another is modeled using a characteristic time which is the sum of a laminar (high temperature) chemistry time, an ignition (low temperature) chemistry time, and a turbulence mixing time. The ignition characteristic time was modeled using data from elementary initiation reactions and has the Arrhenius form. It was found to be possible to match all engine test cases reasonably well with one set of combustion model constants.

Combustion and flow were calculated in a spark-ignited two-stroke crankcase-scavenged engine using a laminar and turbulent characteristic-time combustion submodel in the three-dimensional KIVA code. Both premixed-charge and fuel-injected cases were examined. A multi-cylinder engine simulation program was used to specify initial and boundary conditions for the computation of the scavenging process. A sensitivity study was conducted using the premixed-charge engine data. The influence of different port boundary conditions on the scavenging process was examined. At high delivery ratios, the results were insensitive to variations in the scavenging flow or residual fraction details. In this case, good agreement was obtained with the experimental data using an existing combustion submodel, previously validated in a four-stroke engine study.

In automotive industry it has been a challenge to retain diesel-like thermal efficiency while maintaining low emissions. Numerous studies have shown significant progress in achieving low emissions through the introduction of common-rail injection systems, multiple injections and exhaust gas recirculation and by using a high octane number fuel, like gasoline, to achieve adequate premixing. On the other hand, low temperature combustion strategies, like HCCI and PCCI, have also shown promising results in terms of reducing both NOx and soot emissions simultaneously. With the increasing capacity of computers, multi-dimensional CFD engine modeling enables a reasonably good prediction of combustion characteristics and pollutant emissions, which is the motivation behind the present research. The current research effort presents an optimization study of light-duty compression ignition engine performance, while meeting the emission regulation targets.

Due to concerns regarding pollutant and CO2 emissions, advanced combustion modes that can simultaneously reduce exhaust emissions and improve thermal efficiency have been widely investigated. The main characteristic of the new combustion strategies, such as HCCI and LTC, is that the formation of a homogenous mixture or a controllable stratified mixture is required prior to ignition. The major issue with these approaches is the lack of a direct method for the control of ignition timing and combustion rate, which can be only indirectly controlled using high EGR rates and/or lean mixtures. Homogeneous Charge Progressive Combustion (HCPC) is based on the split-cycle principle. Intake and compression phases are performed in a reciprocating external compressor, which drives the air into the combustor cylinder during the combustion process, through a transfer duct. A transfer valve is positioned between the compressor cylinder and the transfer duct.