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.

Cycle-to-cycle variation in combustion phasing and combustion rate cause knock to occur differently in every cycle. This is found to be true even if the end gas thermo-chemical time history is the same. Three cycles are shown that have matched combustion phasing, combustion rate, and time of knock onset, but have knock intensity that differs by a factor of six. Thus, the prediction of knock intensity must include a stochastic component. It is shown that there is a relationship between the maximum possible knock intensity and the unburned fuel energy at the time of knock onset. Further, for a small window of unburned energy at knock onset, the probability density function of knock intensity is self similar when scaled by the 95th percentile of the cumulative distribution, and log-normal in shape.

The long-term goal of this work is to develop a conceptual model for multiple injections of diesel jets. The current work contributes to that effort by performing a detailed modeling investigation into mechanisms that are predicted to control 1st and 2nd stage ignition in single-pulse diesel (n-dodecane) jets under different conditions. One condition produces a jet with negative ignition dwell that is dominated by mixing-controlled heat release, and the other, a jet with positive ignition dwell and dominated by premixed heat release. During 1st stage ignition, fuel is predicted to burn similarly under both conditions; far upstream, gases at the radial-edge of the jet, where gas temperatures are hotter, partially react and reactions continue as gases flow downstream. Once beyond the point of complete fuel evaporation, near-axis gases are no longer cooled by the evaporation process and 1st stage ignition transitions to 2nd stage ignition.

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.

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 method for automatically generating hex-dominant meshes for Computational Fluid Dynamics (CFD) applications is presented in this article. Two important regions of the mesh for any CFD simulation are the interior mesh and the boundary layer mesh. The interior mesh needs to be fine in the critical flow regions to ensure accurate solutions. The proposed method uses Bubble Mesh algorithm which packs bubbles inside the geometry to generate the mesh nodes. Algorithm was tested for sample flow problems and improvements were made to interior and boundary layer mesh generation methods. The interior mesh is generated using directionality and sizing control functions specified on the points of a 3D grid generated over the entire geometry. This offers a flexible control over mesh sizing and local mesh refinement. Boundary layer mesh is important to accurately model the physics of boundary layer near the geometry walls.

Diesel piston bowl geometry can affect turbulent mixing and therefore it impacts heat-release rates, thermal efficiency, and soot emissions. The focus of this work is on the effects of bowl geometry and injection timing on turbulent flow structure. This computational study compares engine behavior with two pistons representing competing approaches to combustion chamber design: a conventional, re-entrant piston bowl and a stepped-lip piston bowl. Three-dimensional computational fluid dynamics (CFD) simulations are performed for a part-load, conventional diesel combustion operating point with a pilot-main injection strategy under non-combusting conditions. Two injection timings are simulated based on experimental findings: an injection timing for which the stepped-lip piston enables significant efficiency and emissions benefits, and an injection timing with diminished benefits compared to the conventional, re-entrant piston.

Life histories of droplets evaporating on a hot plate under pressure were obtained. The curves are similar to those obtained by one investigator at atmospheric pressure but are displaced to higher temperatures at higher pressures. Similarities between boiling heat transfer and the life history curves are pointed out. Also, that the liquid will most probably reach critical pressure and temperature at temperatures existing inside an engine. The effects of reaching the critical temperature on heat transfer and on vaporization and diffusion are discussed.

This research investigates how the handling of mixing and heat transfer in a multi-zone kinetic solver affects the prediction of carbon monoxide and hydrocarbon emissions for simulations of HCCI engine combustion. A detailed kinetics multi-zone model is now more closely coordinated with the KIVA3V computational fluid dynamics code for simulation of the compression and expansion processes. The fluid mechanics is solved with high spatial and temporal resolution (40,000 cells). The chemistry is simulated with high temporal resolution, but low spatial resolution (20 computational zones). This paper presents comparison of simulation results using this enhanced multi-zone model to experimental data from an isooctane HCCI engine.

The effects of numerical methodology in defining the initial conditions and simulating the compression stroke in D.I. diesel engine CFD computations are studied. Lumped and pointwise approaches were adopted in assigning the initial conditions at IVC. The lumped approach was coupled with a two-dimensional calculation of the compression stroke. The pointwise methodology was based on the results of an unsteady calculation of the intake stroke performed by using the STAR-CD code in the realistic engine and port geometry. Full engine and 60 deg. sector meshes were used in the compression stroke calculations in order to check the accuracy of the commonly applied axi-symmetric fluid dynamics assumption. Analysis of the evolution of the main fluid dynamics parameters revealed that local conditions at the time of injection strongly depend on the numerical procedure adopted.

Surface film formation is an important phenomenon during spray impingement in a combustion chamber. The film that forms on the chamber walls and piston bowl produces soot post-combustion. While some droplets stick to the wall surface, others splash and interact with the gas present inside the combustion chamber. Accurate prediction of both the film thickness and splashed mass is crucial for surface film model development since it leads to a precise estimation of the amount of soot and other exhaust gases formed. This information could guide future studies aimed at a comprehensive understanding of the combustion process and might enable development of engines with reduced emissions. Dynamic structure Large Eddy Simulation (LES) turbulence model implemented for in-cylinder sprays [1] has shown to predict the flow structure of a spray more accurately than the Reynolds-averaged Navier-Stokes turbulence model.

The incorporation of detailed chemistry models in internal combustion engine simulations is becoming mandatory as local, globally lean, low-temperature combustion strategies are setting the path towards a more efficient and environmentally sustainable use of energy resources in transportation. In this paper, we assessed the computational efficiency of a recently developed sparse analytical Jacobian chemistry solver, namely ‘SpeedCHEM’, that features both direct and Krylov-subspace solution methods for maximum efficiency for both small and large mechanism sizes. The code was coupled with a high-dimensional clustering algorithm for grouping homogeneous reactors into clusters with similar states and reactivities, to speed-up the chemical kinetics solution in multi-dimensional combustion simulations.

Dual-fuel reactivity controlled compression ignition (RCCI) engine experiments were conducted with port fuel injection of isooctane and direct injection of n-heptane. The experiments were conducted at a nominal load of 4.75 bar IMEPg, with low isooctane equivalence ratios. Two sets of experiments explored the effects of direct injection timing with single and double injections, and multi-dimensional CFD modeling was used to explore mixture preparation and timing effects. The findings were that if fuel-liner impingement is to be avoided, double injections provide a 40% reduction in CO and HC emissions, resulting in a 1% increase in thermal efficiency. The second engine experiment showed that there is a linear relationship between reactivity (PRF number) and intake temperature. It was also found that if the premixed fuel fraction is above a certain limit, the high-temperature heat release (HTHR) can be manipulated by changing the global PRF number of the in-cylinder fuel blend.

Engine experiments have revealed the importance of fuel condensation on the emission characteristics of low temperature combustion. However, direct in-cylinder experimental evidence has not been reported in the literature. In this paper, the in-cylinder condensation processes observed in optically accessible engine experiments are first illustrated. The observed condensation processes are then simulated using state-of-the-art multidimensional engine CFD simulations with a phase transition model that incorporates a well-validated phase equilibrium numerical solver, in which a thermodynamically consistent phase equilibrium analysis is applied to determine when mixtures become unstable and a new phase is formed. The model utilizes fundamental thermodynamics principles to judge the occurrence of phase separation or combination by minimizing the system Gibbs free energy.

The tradeoff between NOx emissions and combustion instability in an engine operating in the dual-fuel Reactivity Controlled Compression Ignition (RCCI) combustion mode was investigated using a combination of engine experiments and detailed CFD modeling. Experiments were performed on a single cylinder version of a General Motors/Fiat JTD 1.9L four-cylinder diesel engine. Gasoline was injected far upstream of the intake valve using an air assisted injector and fuel vaporization system and diesel was injected directly into the cylinder using a common rail injector. The timing of the diesel injection was swept from −70° ATDC to −20° ATDC while the gasoline percentage was adjusted to hold the average combustion phasing (CA50) and load (IMEPg) constant at 0.5° ATDC and 7 bar, respectively. At each operating point the variation in IMEP, peak PRR, and CA50 was calculated from the measured cylinder pressure trace and NOx, CO, soot and UHC were recorded.

Sector mesh modeling is the dominant computational approach for combustion system design optimization. The aim of this work is to quantify the errors descending from the sector mesh approach through three geometric modeling approaches to an optical diesel engine. A full engine geometry mesh is created, including valves and intake and exhaust ports and runners, and a full-cycle flow simulation is performed until fired TDC. Next, an axisymmetric sector cylinder mesh is initialized with homogeneous bulk in-cylinder initial conditions initialized from the full-cycle simulation. Finally, a 360-degree azimuthal mesh of the cylinder is initialized with flow and thermodynamics fields at IVC mapped from the full engine geometry using a conservative interpolation approach. A study of the in-cylinder flow features until TDC showed that the geometric features on the cylinder head (valve tilt and protrusion into the combustion chamber, valve recesses) have a large impact on flow complexity.

The performance of a structural design significantly depends upon the assumptions made on input load. In order to estimate the input load, during the design and development stage of the suspension assembly of a BAJA car, designers and analysts invest immense amount of time and effort to formulate the mathematical model of the design. These theoretical formulations may include idealization errors which can affect the performance of the car as a final product. Due to the errors associated with the assumption of design load, several components might have more weight or may have less strength than needed. This discrepancy between the assumed input load (lab or theoretical studies) and the actual load from the environment can be eliminated by performing a real life testing process using load recovery methodology. Commercial load cells exist in industry to give engineers insight to understanding the complex real world loading of their structures.

Knock in a Rotax-914 engine was modeled and investigated using an improved version of the KIVA-3V code with a G-equation combustion model, together with a reduced chemical kinetics model. The ERC-PRF mechanism with 47 species and 132 reactions [1] was adopted to model the end gas auto-ignition in front of the flame front. The model was validated by a Caterpillar SI engine and a Rotax-914 engine in different operating conditions. The simulation results agree well with available experimental results. A new engineering quantified knock criterion based on chemical mechanism was then proposed. Hydroperoxyl radical (HO2) shows obvious accumulation before auto-ignition and a sudden decrease after auto-ignition. These properties are considered to be a good capability for HO2 to investigate engine knock problems.

3-D Computational Fluid Dynamics (CFD) simulations have been performed to study particulate formation in a Spark-Ignition (SI) engine under premixed conditions. A semi-detailed soot model and a chemical kinetic model, including poly-aromatic hydrocarbon (PAH) formation, were coupled with a spark ignition model and the G equation flame propagation model for SI engine simulations and for predictions of soot mass and particulate number density. The simulation results for in-cylinder pressure and particle size distribution (PSDs) are compared to available experimental studies of equivalence ratio effects during premixed operation. Good predictions are observed with regard to cylinder pressure, combustion phasing and engine load. Qualitative agreements of in-cylinder particle distributions were also obtained and the results are helpful to understand particulate formation processes.

The three-dimensional CFD codes KIVA-II and KIVA-3 have been used together to study the effects of intake generated in-cylinder flow structure on fuel-air mixing and combustion in a direct injected (DI) Diesel engine. In order to more accurately account for the effect of intake flow on in-cylinder processes, the KIVA-II code has been modified to allow for the use of data from other CFD codes as initial conditions. Simulation of the intake and compression strokes in a heavy-duty four-stroke DI Diesel engine has been carried out using KIVA-3. Flow quantities and thermodynamic field information were then mapped into a computational grid in KIVA-II for use in the study of mixing and combustion. A laminar and turbulent timescale combustion model, as well as advanced spray models, including wave breakup atomization, dynamic drop drag, and spray-wall interaction has been used in KIVA-II.