Search Results

The possibility to operate current diesel engines in dual-fuel mode with the addition of an alternative fuel is fundamental to accelerate the energy transition to achieve carbon neutrality. The simulation of the dual- fuel combustion process with 0D/1D combustion models is fundamental for the performance prediction, but still particularly challenging, due to chemical interactions of the mixture. The authors defined a novel data-driven workflow for the development of combustion reaction mechanisms and used it to generate a dual-fuel mechanism for Ammonia and Diesel Primary Reference Fuels (DPRF) suitable for efficient combustion simulations in heavy duty engines, with variable cetane number Diesel fuels. A baseline reaction mechanism was created by merging the detailed ammonia mechanism by Glarborg et al. with reaction pathways for n- hexadecane and 2,2,4,4,6,8,8-heptamethylnonane from a well-established multi-component fuel mechanism.

We propose a novel dual-fuel combustion model for simulating heavy-duty engines with the G-Equation. Dual-Fuel combustion strategies in such engines features direct injection of a high-reactivity fuel into a lean, premixed chamber which has a high resistance to autoignition. Distinct combustion modes are present: the DI fuel auto-ignites following chemical ignition delay after spray vaporization and mixing; a reactive front is formed on its surroundings; it develops into a well-structured turbulent flame, which propagates within the premixed charge. Either direct chemistry or the flame-propagation approach (G- Equation), taken alone, do not produce accurate results. The proposed Dual-Fuel model decides what regions of the combustion chamber should be simulated with either approach, according to the local flame state; and acts as a “kernel” model for the G- Equation model. Direct chemistry is run in the regions where a premixed front is not present.

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.

Reactivity-controlled compression ignition (RCCI) has been shown to be capable of providing improved engine efficiencies coupled with the benefit of low emissions via in-cylinder fuel blending. Much of the previous body of work has studied the benefits of RCCI operation using high injection pressures (e.g., 500 bar or greater) with common rail injection (CRI) hardware. However, low-pressure fueling technology is capable of providing significant cost savings. Due to the broad market adoption of gasoline direct injection (GDI) fueling systems, a market-type prototype GDI injector was selected for this study. Single-cylinder light-duty engine experiments were undertaken to examine the performance and emissions characteristics of the RCCI combustion strategy with low-pressure GDI technology and compared against high injection pressure RCCI operation. Gasoline and diesel were used as the low-reactivity and high-reactivity fuels, respectively.

In a recent study, quantitative measurements were presented of in-cylinder spatial distributions of mixture equivalence ratio in a single-cylinder light-duty optical diesel engine, operated with a non-reactive mixture at conditions similar to an early injection low-temperature combustion mode. In the experiments a planar laser-induced fluorescence (PLIF) methodology was used to obtain local mixture equivalence ratio values based on a diesel fuel surrogate (75% n-heptane, 25% iso-octane), with a small fraction of toluene as fluorescing tracer (0.5% by mass). Significant changes in the mixture's structure and composition at the walls were observed due to increased charge motion at high swirl and injection pressure levels. This suggested a non-negligible impact on wall heat transfer and, ultimately, on efficiency and engine-out emissions.

The paper presents the development of a novel approach to the solution of detailed chemistry in internal combustion engine simulations, which relies on the analytical computation of the ordinary differential equations (ODE) system Jacobian matrix in sparse form. Arbitrary reaction behaviors in either Arrhenius, third-body or fall-off formulations can be considered, and thermodynamic gas-phase mixture properties are evaluated according to the well-established 7-coefficient JANAF polynomial form. The current work presents a full validation of the new chemistry solver when coupled to the KIVA-4 code, through modeling of a single cylinder Caterpillar 3401 heavy-duty engine, running in two-stage combustion mode.

An investigation of high speed direct injection (DI) compression ignition (CI) engine combustion fueled with gasoline injected using a triple-pulse strategy in the low temperature combustion (LTC) regime is presented. This work aims to extend the operation ranges for a light-duty diesel engine, operating on gasoline, that have been identified in previous work via extended controllability of the injection process. The single-cylinder engine (SCE) was operated at full load (16 bar IMEP, 2500 rev/min) and computational simulations of the in-cylinder processes were performed using a multi-dimensional CFD code, KIVA-ERC-Chemkin, that features improved sub-models and the Chemkin library. The oxidation chemistry of the fuel was calculated using a reduced mechanism for primary reference fuel combustion chosen to match ignition characteristics of the gasoline fuel used for the SCE experiments.

Many recent studies have shown that the Reactivity Controlled Compression Ignition (RCCI) combustion strategy can achieve high efficiency with low emissions. However, it has also been revealed that RCCI combustion is difficult at high loads due to its premixed nature. To operate at moderate to high loads with gasoline/diesel dual fuel, high amounts of EGR or an ultra low compression ratio have shown to be required. Considering that both of these approaches inherently lower thermodynamic efficiency, in this study natural gas was utilized as a replacement for gasoline as the low-reactivity fuel. Due to the lower reactivity (i.e., higher octane number) of natural gas compared to gasoline, it was hypothesized to be a better fuel for RCCI combustion, in which a large reactivity gradient between the two fuels is beneficial in controlling the maximum pressure rise rate.

A generalized re-normalization group (RNG) turbulence model based on the local "dimensionality" of the flow field is proposed. In this modeling approach the model coefficients C₁, C₂, and C₃ are all constructed as functions of flow strain rate. In order to further validate the proposed turbulence model, the generalized RNG closure model was applied to model the backward facing step flow (a classic test case for turbulence models). The results indicated that the modeling of C₂ in the generalized RNG closure model is reasonable, and furthermore, the predictions of the generalized RNG model were in better agreement with experimental data than the standard RNG turbulence model. As a second step, the performance of the generalized RNG closure was investigated for a complex engine flow.

The simulation of combustion chemistry in internal combustion engines is challenging due to the need to include detailed reaction mechanisms to describe the engine physics. Computational times needed for coupling full chemistry to CFD simulations are still too computationally demanding, even when distributed computer systems are exploited. For these reasons the present paper proposes a time scale separation approach for the integration of the chemistry differential equations and applies it in an engine CFD code. The time scale separation is achieved through the estimation of a characteristic time for each of the species and the introduction of a sampling timestep, wherein the chemistry is subcycled during the overall integration. This allows explicit integration of the system to be carried out, and the step size is governed by tolerance requirements.

The role of the fluid motion in a diesel engine on mixing and combustion was investigated using the CFD code Kiva-3v. The study considered pre-mixed charge compression ignition (PCCI) combustion that is a hybrid combustion system characterized by early injection timings and high amounts of EGR dilution to delay the start and lower the temperature of combustion. The fuel oxidizer mixture is not homogeneous at the start of combustion and therefore requires further mixing for complete combustion. PCCI combustion systems are characterized by relatively high CO and UHC emissions. This work investigates attenuating CO emissions by enhancing mixing processes through non-uniform flowfield motions. The fluid motion was characterized by the amount of average angular rotation about the cylindrical axis (swirl ratio) and the amount of non-uniform motion imparted by the relative amounts of mass inducted through tangential and helical intake ports in a 0.5L HSDI diesel engine.

The present work proposes a methodology for diesel engine development using multi-dimensional CFD and computer optimization. A multi-objective genetic algorithm coupled with the KIVA3V Release 2 code was used to optimize a high speed direct injection (HSDI) diesel engine for passenger car applications. The simulations were conducted using high-throughput computing with the CONDOR system. An automated grid generator was used for efficient mesh generation with 11 variable piston bowl geometry parameters. The first step in the procedure was to search for an optimal nozzle and piston bowl design. In this case, spray targeting, swirl ratio, and piston bowl shape were optimized separately for two full-load cases using simpler efficient combustion models (the characteristic time scale model and the shell ignition model). The optimal designs from the two optimizations were then validated using a combustion model with detailed chemistry (KIVA-CHEMKIN model and ERC n-heptane mechanism).

Two different types of mesh used for diesel combustion with the KIVA-4 code are compared. One is a well established conventional KIVA-3 type polar mesh. The other is a non-polar mesh with uniform size throughout the piston bowl so as to reduce the number of cells and to improve the quality of the cell shapes around the cylinder axis which can contain many fuel droplets that affect prediction accuracy and the computational time. This mesh is specialized for the KIVA-4 code which employs an unstructured mesh. To prevent dramatic changes in spray penetration caused by the difference in cell size between the two types of mesh, a recently developed spray model which reduces mesh dependency of the droplet behavior has been implemented. For the ignition and combustion models, the Shell model and characteristic time combustion (CTC) model are employed.

A study of partially premixed combustion (PPC) with non-oxygenated 91 pump octane number1 (PON) commercially available gasoline was performed using a heavy-duty (HD) compression-ignition (CI) 2.44 l Caterpillar 3401E single-cylinder oil test engine (SCOTE). The experimental conditions selected were a net indicated mean effective pressure (IMEP) of 11.5 bar, an engine speed of 1300 rev/min, an intake temperature of 40°C with intake and exhaust pressures of 200 and 207 kPa, respectively. The baseline case for all studies presented had 0% exhaust gas recirculation (EGR), used a dual injection strategy a -137 deg ATDC pilot SOI and a -6 deg ATDC main start-of-injection (SOI) timing with a 30/70% pilot/main fuel split for a total of 5.3 kg/h fueling (equating to approximately 50% load). Combustion and emissions characteristics were explored relative to the baseline case by sweeping main and pilot SOI timings, injection split fuel percentage, intake pressure, load and EGR levels.

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.

A transport equation residual model incorporating refined G-equation and detailed chemical kinetics combustion models has been developed and implemented in the ERC KIVA-3V release2 code for Gasoline Direct Injection (GDI) engine simulations for better predictions of flame propagation. In the transport equation residual model a fictitious species concept is introduced to account for the residual gases in the cylinder, which have a great effect on the laminar flame speed. The residual gases include CO2, H2O and N2 remaining from the previous engine cycle or introduced using EGR. This pseudo species is described by a transport equation. The transport equation residual model differentiates between CO2 and H2O from the previous engine cycle or EGR and that which is from the combustion products of the current engine cycle.

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.

The lift-off length plays a significant role in spray combustion as it influences the air entrainment upstream of the lift-off location and hence the soot formation. Accurate prediction of lift-off length thus becomes a prerequisite for accurate soot prediction in lifted flames. In the present study, KIVA-3v coupled with CHEMKIN, as developed at the Engine Research Center (ERC), is used as the CFD model. Experimental data from the Sandia National Labs. is used for validating the model predictions of n-heptane lift-off lengths and soot formation details in a constant volume combustion chamber. It is seen that the model predictions, in terms of lift-off length and soot mass, agree well with the experimental results for low ambient density (14.8 kg/m3) cases with different EGR rates (21% O2 - 8% O2). However, for high density cases (30 kg/m3) with different EGR rates (15% O2 - 8% O2) disagreements were found.

A computational study was performed to evaluate the effects of bowl geometry, fuel spray targeting and swirl ratio under highly diluted, low-temperature combustion conditions in a heavy-duty diesel engine. This study is used to examine aspects of low-temperature combustion that are affected by mixing processes and offers insight into the effect these processes have on emissions formation and oxidation. The foundation for this exploratory study stems from a large data set which was generated using a genetic algorithm optimization methodology. The main results suggest that an optimal combination of spray targeting, swirl ratio and bowl geometry exist to simultaneously minimize emissions formation and improve soot and CO oxidation rates. Spray targeting was found to have a significant impact on the emissions and fuel consumption performance, and was furthermore found to be the most influential design parameter explored in this study.

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.