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.

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.

Particle image velocimetry (PIV) is used to investigate the structure and evolution of the mean velocity field in the swirl (r-θ) plane of a motored, optically accessible diesel engine with a typical production combustion chamber geometry under motoring conditions (no fuel injection). Instantaneous velocities were measured were made at three swirl-plane heights (3 mm, 10 mm and 18 mm below the firedeck) and three swirl ratios (2.2, 3.5 and 4.5) over a range of crank angles in the compression and expansion strokes. The data allow for a direct analysis of the structures within the ensemble mean flow field, the in-cylinder swirl ratio, and the radial profile of the tangential velocity. At all three swirl ratios, the ensemble mean velocity field contains a single dominant swirl flow structure that is tilted with respect to the cylinder axis. The axis of this structure precesses about the cylinder axis in a manner that is largely insensitive to swirl ratio.

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.

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.

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.

A novel algebraic similarity model for subgrid scalar dissipation rate has been developed as part of the Large Eddy Simulation (LES) package KIVA3V-LES for diesel engine study. The model is proposed from an a priori study using Direct Numerical Simulation (DNS) of forced isotropic turbulence. In the a posteriori test, fully resolved turbulent passive scalar field measurements are used to validate the model in actual engine flows. For reason of the length limit by SAE and the specific interest in engine applications, only a prior test and a posteriori test in engine flows are included in this paper. A posteriori tests in isotropic cube flow, turbulent round jet and flame cases will be presented in separate papers. An engine LES simulation of multi consecutive cycles was performed in this study.

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.

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.

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.

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.

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.

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.

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.