Search Results

The 4th Workshop of the Engine Combustion Network (ECN) was held September 5-6, 2015 in Kyoto, Japan. This manuscript presents a summary of the progress in experiments and modeling among ECN contributors leading to a better understanding of soot formation under the ECN “Spray A” configuration and some parametric variants. Relevant published and unpublished work from prior ECN workshops is reviewed. Experiments measuring soot particle size and morphology, soot volume fraction (fv), and transient soot mass have been conducted at various international institutions providing target data for improvements to computational models. Multiple modeling contributions using both the Reynolds Averaged Navier-Stokes (RANS) Equations approach and the Large-Eddy Simulation (LES) approach have been submitted. Among these, various chemical mechanisms, soot models, and turbulence-chemistry interaction (TCI) methodologies have been considered.

In-cylinder flow and combustion processes simulated with the standard k-ε turbulence model and with an alternative model-employing a non-linear, quadratic equation for the turbulent stresses-are contrasted for both motored and fired engine operation at two loads. For motored operation, the differences observed in the predictions of mean flow development are small and do not emerge until expansion. Larger differences are found in the spatial distribution and magnitude of turbulent kinetic energy. The non-linear model generally predicts lower energy levels and larger turbulent time scales. With fuel injection and combustion, significant differences in flow structure and in the spatial distribution of soot are predicted by the two models. The models also predict considerably different combustion efficiencies and NOx emissions.

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.

In a recent experimental study the in-cylinder spatial distribution of mixture equivalence ratio was quantified under non-combusting conditions by planar laser-induced fluorescence (PLIF) of a fuel tracer (toluene). The measurements were made in a single-cylinder, direct-injection, light-duty diesel engine at conditions matched to an early-injection low-temperature combustion mode. A fuel amount corresponding to a low load (3.0 bar indicated mean effective pressure) operating condition was introduced with a single injection at -23.6° ATDC. The data were acquired during the mixture preparation period from near the start of injection (-22.5° ATDC) until the crank angle where the start of high-temperature heat release normally occurs (-5° ATDC). In the present study the measured in-cylinder images are compared with a fully resolved three-dimensional CFD model, namely KIVA3V-RANS simulations.

Premixed charge compression ignition (PCI) strategies offer the potential for simultaneously low NOx and soot emissions with diesel-like efficiency. However, these strategies are generally confined to low loads due to inadequate control of combustion phasing and heat-release rate. One PCI strategy, dual-fuel reactivity-controlled compression ignition (RCCI), has been developed to control combustion phasing and rate of heat release. The RCCI concept uses in-cylinder blending of two fuels with different auto-ignition characteristics to achieve controlled high-efficiency clean combustion. This study explores fuel reactivity stratification as a method to control the rate of heat release for PCI combustion. To introduce fuel reactivity stratification, the research engine is equipped with two fuel systems. A low-pressure (100 bar) gasoline direct injector (GDI) delivers iso-octane, and a higher-pressure (600 bar) common-rail diesel direct-injector delivers n-heptane.

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.

The influence of the constant C3, which multiplies the mean flow divergence term in the model equation for the turbulent kinetic energy dissipation, is examined in a motored diesel engine for three different swirl ratios and three different spatial locations. Predicted temporal histories of turbulence energy and its dissipation are compared with experimentally-derived estimates. A “best-fit” value of C3 = 1.75, with an approximate uncertainty of ±0.3 is found to minimize the error between the model predictions and the experiments. Using this best-fit value, model length scale behavior corresponds well with that of measured velocity-correlation integral scales during compression. During expansion, the model scale grows too rapidly. Restriction of the model assessment to the expansion stroke suggests that C3 = 0.9 is more appropriate during this period.

Two-Color imaging optics were developed and used to observe soot emission processes in a modern heavy-duty diesel engine. The engine was equipped with a common rail, electronically-controlled, high-pressure fuel injection system that is capable of up to four injection pulses per engine cycle. The engine was instrumented with an endoscope system for optical access for the combustion visualization. Multidimensional combustion and soot modeling results were used for comparisons to enhance the understanding and interpretation of the experimental data. Good agreement between computed and measured cylinder pressures, heat release and soot and NOx emissions was achieved. In addition, good qualitative agreement was found between in-cylinder soot concentration (KL) and temperature fields obtained from the endoscope images and those obtained from the multidimensional modeling.

Consideration of the heat transfer effects in low-heat-rejection engines has prompted further study into engine heat transfer phenomena. In a previous study, an approximate solution of the one-dimensional energy equation was acquired for transient, compressible, low-Mach number, turbulent boundary layers typical of those found in engines. The current study shows that an approximate solution of the one-dimensional energy equation with arbitrarily-distributed heat release can also be obtained. Using this model, the effects of high temperature walls, combustion, and autoignition on heat transfer can be studied. In the case of high temperature walls, the model predicts the expected behavior unless the quench distance gets very small. For combustion, the reaction must occur close to the wall for a direct effect on the heat transfer to be observed. With autoignition, instantaneous values of heat flux reach levels as high as 6 MW/m2, and oscillate in phase with the pressure wave.

High-speed video imaging in a swirl-supported (Rs = 1.7), direct-injection heavy-duty diesel engine operated with moderate-to-high EGR rates reveals a distinct correlation between the spatial distribution of luminous soot and mean flow vorticity in the horizontal plane. The temporal behavior of the experimental images, as well as the results of multi-dimensional numerical simulations, show that this soot-vorticity correlation is caused by the presence of a greater amount of soot on the windward side of the jet. The simulations indicate that while flow swirl can influence pre-ignition mixing processes as well as post-combustion soot oxidation processes, interactions between the swirl and the heat release can also influence mixing processes. Without swirl, combustion-generated gas flows influence mixing on both sides of the jet equally. In the presence of swirl, the heat release occurs on the leeward side of the fuel sprays.

Engine design is a time consuming process in which many costly experimental tests are usually conducted. With increasing prediction ability of engine simulation tools, engine design aided by CFD software is being given more attention by both industry and academia. It is also of much interest to be able to use design information gained from an existing engine design of one size in the design of engines of other sizes to reduce design time and costs. Therefore it is important to study size-scaling relationships for engines over wide range of operating conditions. This paper presents CFD studies conducted for two production diesel engines - a light-duty GM-Fiat engine (0.5L displacement) and a heavy-duty Caterpillar engine (2.5L displacement). Previously developed scaling arguments, including an equal spray penetration scaling model and an extended, equal flame lift-off length scaling model were employed to explore the parametric scaling connections between the two engines.

Various single and split injection schemes are studied to provide a better understanding of fuel distribution during cold starting in DI diesel engines. Improved spray-wall interaction, fuel film and multicomponent vaporization models are used to analyze the combustion processes. Better combustion characteristics are obtained for the split injection schemes than with a single injection. An analysis of the fuel impingement processes identifies the mechanisms involved in producing the differences in vaporization and combustion of the fuel. A greater amount of splashing occurred for the split injections compared to a single injection. This behavior is attributed to the decreased film thickness (less dissipation of impingement energy), the decreased impingement area (obtained by increasing the impingement Weber number), and most importantly, the reduced frequency of drop impingement.

The effect of injection rate shape on spray evolution and emission characteristics is investigated and a methodology for active control of fuel injection is proposed. Extensive validation of advanced vaporization and primary jet breakup models was performed with experimental data before studying the effects of systematic changes of injection rate shape. Excellent agreement with the experiments was obtained for liquid and vapor penetration lengths, over a broad range of gas densities and temperatures. Also the predicted flame lift-off lengths of reacting diesel fuel sprays were in good agreement with the experiments. After the validation of the models, well-defined rate shapes were used to study the effect of injection rate shape on liquid and vapor penetration, flame lift-off lengths and emission characteristics.

Simultaneous two-component measurements of gas velocity and multi-dimensional numerical simulation are employed to characterize the evolution of the in-cylinder turbulent flow structure in a re-entrant bowl-in-piston engine under motored operation. The evolution of the mean flow field, turbulence energy, turbulent length scales, and the various terms contributing to the production of the turbulence energy are correlated and compared, with the objectives of clarifying the physical mechanisms and flow structures that dominate the turbulence production and of identifying the source of discrepancies between the measured and simulated turbulence fields. Additionally, the applicability of the linear turbulent stress modeling hypothesis employed in the k-ε model is assessed using the experimental mean flow gradients, turbulence energy, and length scales.

In-cylinder fluid velocity is measured in an optically accessible, fired HSDI engine at idle. The velocity field is also calculated, including the full induction stroke, using multi-dimensional fluid dynamics and combustion simulation models. A detailed comparison between the measured and calculated velocities is performed to validate the computed results and to gain a physical understanding of the flow evolution. Motored measurements are also presented, to clarify the effects of the fuel injection process and combustion on the velocity field evolution. The calculated mean in-cylinder angular momentum (swirl ratio) and mean flow structures prior to injection agree well with the measurements. Modification of the mean flow by fuel injection and combustion is also well captured.

Managing the injection rate profile is a powerful tool to control engine performance and emission levels. In particular, Common Rail (C.R.) injection systems allow an almost completely flexible fuel injection event in DI-diesel engines by permitting a free mapping of the start of injection, injection pressure, rate of injection and, in the near future, multiple injections. This research deals with the development of a network-based numerical tool for understanding operating condition limits of the Common Rail injector. The models simulate the electro-fluid-mechanical behavior of the injector accounting for cavitation in the nozzle holes. Validation against experiments has been performed. The model has been used to provide insight into the operating conditions of the injector and in order to highlight the application to injection system design.

Cycle-resolved analysis of velocity data obtained in the re-entrant bowl of a fired high-;speed, direct-injection diesel engine, demonstrates an unambiguous, approximately 100% increase in late-cycle turbulence levels over the levels measured during motored operation. Model predictions of the flow field, obtained employing RNG k-ε turbulence modeling in KIVA-3V, do not capture this increased turbulence. A combined experimental and computational approach is taken to identify the source of this turbulence. The results indicate that the dominant source of the increased turbulence is associated with the formation of an unstable distribution of mean angular momentum, characterized by a negative radial gradient. The importance of this source of flow turbulence has not previously been recognized for engine flows. The enhanced late-cycle turbulence is found to be very sensitive to the flow swirl level.

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.

Simultaneous measurements of the radial and the tangential components of velocity are obtained in a high-speed, direct-injection diesel engine typical of automotive applications. Results are presented for engine operation with fuel injection, but without combustion, for three different swirl ratios and four injection pressures. With the mean and fluctuating velocities, the r-θ plane shear stress and the mean flow gradients are obtained. Longitudinal and transverse length scales are also estimated via Taylor's hypothesis. The flow is shown to be sufficiently homogeneous and stationary to obtain meaningful length scale estimates. Concurrently, the flow and injection processes are simulated with KIVA-3V employing a RNG k-ε turbulence model. The measured turbulent kinetic energy k, r-θ plane mean strain rates ( 〈Srθ〉, 〈Srr〉, and 〈Sθθ〉 ), deviatoric turbulent stresses , and the r-θ plane turbulence production terms are compared directly to the simulated results.

The physical in-cylinder processes and ignition during cold starting have been studied using computational models, with particular attention to the influences of blowby, heat transfer during the compression stroke, spray development, vaporization and fuel/air mixture formation and ignition. Two different modeling approaches were used. A thermodynamic zero dimensional cycle analysis program in which the fuel injection effects were not modeled, was used to determine overall and gas exchange effects. The three-dimensional KIVA-II code was used to determine details of the closed cycle events, with modified atomization, blowby and spray/wall impingement models, and a simplified model for ignition. The calculations were used to obtain an understanding of the cold starting process and to identify practical methods for improving cold starting of direct injection diesel engines.