Search Results

Detailed chemistry and turbulence-chemistry interaction need to be properly taken into account for a realistic combustion simulation of IC engines where advanced combustion modes, multiple injections and stratified combustion involve a wide range of combustion regimes and require a proper description of several phenomena such as auto-ignition, flame stabilization, diffusive combustion and lean premixed flame propagation. To this end, different approaches are applied and the most used ones rely on the well-stirred reactor or flamelet assumption. However, well-mixed models do not describe correctly flame structure, while unsteady flamelet models cannot easily predict premixed flame propagation and triple flames. A possible alternative for them is represented by transported probability density functions (PDF) methods, which have been applied widely and effectively for modeling turbulent reacting flows under a wide range of combustion regimes.

A new multi-zone model for the simulation of HCCI engine is here presented. The model includes laminar and turbulent diffusion and conduction exchange between the zones and the last improvements on the numerical aspects. Furthermore, a new strategy for the zone discretization is presented, which allows a better description of the near-wall zones. The aim of the work is to provide a fast and reliable model for carrying out chemical analysis with detailed kinetic schemes. A preliminary sensitivity analysis allows to verify that 10 zones are a convenient number for a good compromise between the computational effort and the description accuracy. The multi-zone predictions are then compared with the CFD ones to find the effective turbulence parameters, with the aim to describe the near-wall phenomena, both in a reactive and non-reactive cases.

In recent years, the increase of gasoline fuel injection pressure is a way to improve thermal efficiency and lower engine-out emissions in GDI homogenous combustion concept. The challenge of controlling particulate formation as well in mass and number concentrations imposed by emissions regulations can be pursued improving the mixture preparation process and avoiding mixture inhomogeneity with ultra-high injection pressure values up to 100 MPa. The increase of the fuel injection pressure in GDI homogeneous systems meets the demand for increased injector static flow, while simultaneously improves the spray atomization and mixing characteristics with consequent better combustion performance. Few studies quantify the effects of high injection pressure on transient gasoline spray evolution. The aim of this work was to simulate with OpenFOAM the spray morphology of a commercial gasoline injected in a constant volume vessel by a prototypal GDI injector.

The interaction of fuel sprays and in-cylinder flow in direct-injection engines is expected to alter kinetic energy and integral length scales at least during some portions of the engine cycle. High-speed particle image velocimetry was implemented in an optical four-valve, pent-roof spark-ignition direct-injection single-cylinder engine to quantify this effect. Non-firing motored engine tests were performed at 1300 RPM with and without fuel injection. Two fuel injection timings were investigated: injection in early intake stroke represents quasi-homogenous engine condition; and injection in mid compression stroke mimics the stratified combustion strategy. Two-dimensional crank angle resolved velocity fields were measured to examine the kinetic energy and integral length scale through critical portions of the engine cycle. Reynolds decomposition was applied on the obtained engine flow fields to extract the fluctuations as an indicator for the turbulent flow.

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.

Today, Direct-Injection systems are widely used on Spark-Ignition engines in combination with turbo-charging to reduce the fuel-consumption and the knock risks. In particular, the spread of Gasoline Direct Injection (GDI) systems is mainly related to the use of new generations of multi-hole, high-pressure injectors whose characteristics are quite different with respect to the hollow-cone, low-pressure injectors adopted in the last decade. This paper presents the results of an experimental campaign conducted on the spray produced by a GDI six-holes injector into a constant volume vessel with optical access. The vessel was filled with air at atmospheric pressure. Different operating conditions were considered for an injection pressure ranging from 3 to 20 MPa. For each operating condition, spray images were acquired by a CCD camera and then post processed to evaluate the spray penetration and cone angles.

Modeling plume interaction and collapse for direct-injection gasoline sprays is important because of its impact on fuel-air mixing and engine performance. Nevertheless, the aerodynamic interaction between plumes and the complicated two-phase coupling of the evaporating spray has shown to be notoriously difficult to predict. With the availability of high-speed (100 kHz) Particle Image Velocimetry (PIV) experimental data, we compare velocity field predictions between plumes to observe the full temporal evolution leading up to plume merging and complete spray collapse. The target “Spray G” operating conditions of the Engine Combustion Network (ECN) is the focus of the work, including parametric variations in ambient gas temperature. We apply both LES and RANS spray models in different CFD platforms, outlining features of the spray that are most critical to model in order to predict the correct aerodynamics and fuel-air mixing.

A numerical methodology to simulate the high pressure spray evolution and the fuel-air mixing in diesel engines is presented. Attention is focused on the employed atomization model, a modified version of the Huh and Gosman, on the definition of a turbulence length scale limiter and of an adaptive local mesh refinement technique to minimize the result grid dependency. All the discussed models were implemented into Lib-ICE, which is a set of libraries and solvers, specifically tailored for engine simulations, which runs under the open-source CFD technology OpenFOAM®. To provide a comprehensive assessment of the proposed methodology, the validation procedure consisted into simulating, with a unique and coherent setup of all models, two different sets of experiments: a non-evaporating diesel fuel spray in a constant-volume vessel with optical access and an evaporating non-reacting diesel fuel spray in an optical engine.

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.

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.

Developing a complete understanding of the structure and behavior of the near-wall region (NWR) in reciprocating, internal combustion (IC) engines and of its interaction with the core flow is needed to support the implementation of advanced combustion and engine operation strategies, as well as predictive computational models. The NWR in IC engines is fundamentally different from the canonical steady-state turbulent boundary layers (BL), whose structure, similarity and dynamics have been thoroughly documented in the technical literature. Motivated by this need, this paper presents results from the analysis of two-component velocity data measured with particle image velocimetry near the head of a single-cylinder, optical engine. The interaction between the NWR and the core flow was quantified via statistical moments and two-point velocity correlations, determined at multiple distances from the wall and piston positions.

Objective of this work is the incorporation of the flame stretch effects in an Eulerian-Lagrangian model for premixed SI combustion in order to describe ignition and flame propagation under highly inhomogeneous flow conditions. To this end, effects of energy transfer from electrical circuit and turbulent flame propagation were fully decoupled. The first ones are taken into account by Lagrangian particles whose main purpose is to generate an initial burned field in the computational domain. Turbulent flame development is instead considered only in the Eulerian gas phase for a better description of the local flow effects. To improve the model predictive capabilities, flame stretch effects were introduced in the turbulent combustion model by using formulations coming from the asymptotic theory and recently verified by means of DNS studies. Experiments carried out at Michigan Tech University in a pressurized, constant-volume vessel were used to validate the proposed approach.

RNG k-ε closure turbulence dissipation equations are evaluated employing the CFD code KIVA-3V Release 2. The numerical evaluations start by considering simple jet flows, including incompressible air jets and compressible helium jets. The results show that the RNG closure turbulence model predicts lower jet tip penetration than the "standard" k-ε model, as well as being lower than experimental data. The reason is found to be that the turbulence kinetic energy is dissipated too slowly in the downstream region near the jet nozzle exit. In this case, the over-predicted R term in RNG model becomes a sink of dissipation in the ε-equation. As a second step, the RNG turbulence closure dissipation models are further tested in complex engine flows to compare against the measured evolution of turbulence kinetic energy, and an estimate of its dissipation rate, during both the compression and expansion processes.

Premixed combustion modes in compression ignition engines are studied as a promising solution to meet fuel economy and increasingly stringent emissions regulations. Nevertheless, PCCI combustion systems are not yet consolidated enough for practical applications. The high complexity of such combustion systems in terms of both air-fuel charge preparation and combustion process control requires the employment of robust and reliable numerical tools to provide adequate comprehension of the phenomena. Object of this work is the development and validation of suitable models to evaluate the effects of charge premixing levels in diesel combustion. This activity was performed using the Lib-ICE code, which is a set of applications and libraries for IC engine simulations developed using the OpenFOAM® technology.

Detailed chemistry represents a fundamental pre-requisite for a realistic simulation of combustion process in Diesel engines to properly reproduce ignition delay and flame structure (lift-off and soot precursors) in a wide range of operating conditions. In this work, the authors developed reduced mechanisms for n-dodecane starting from the comprehensive kinetic mechanism developed at Politecnico di Milano, well validated and tested in a wide range of operating conditions [1]. An algorithm combining Sensitivity and Flux Analysis was employed for the present skeletal reduction. The size of the mechanisms can be limited to less than 100 species and incorporates the most important details of low-temperature kinetics for a proper prediction of the ignition delay. Furthermore, the high-temperature chemistry is also properly described both in terms of reactivity and species formation, including unsaturated compounds such as acetylene, whose concentration controls soot formation.

The turbulent in-cylinder air flow and the unsteady high-pressure fuel injection lead to a highly transient air fuel mixing process in spark-ignition direct-injection (SIDI) engines, which is the leading cause for combustion cycle-to-cycle variation (CCV) and requires further investigation. In this study, crank-angle resolution particle image velocimetry (PIV) was employed to simultaneously measure the air flow and fuel spray structure at 1300 rpm in an optically accessible single-cylinder SIDI engine. The measurement was conducted at the center tumble plane of the four-valve pent-roof engine, bisecting the spark plug and fuel injector. 84 consecutive cycles were recorded for three engine conditions, i.e. (1) none-fueled motored condition, (2) homogeneous-charge mode with start of injection (SOI) during intake (50 crank-angle degree (CAD) after top dead center exhaust, aTDCexh), and (3) stratified-charge mode with SOI during mid compression (270 aTDCexh).

A primary goal of large eddy simulation, LES, is to capture in-cylinder cycle-to-cycle variability, CCV. This is a first step to assess the efficacy of 35 consecutive computed motored cycles to capture the kinetic energy in the TCC-III engine. This includes both the intra-cycle production and dissipation as well as the kinetic energy CCV. The approach is to sample and compare the simulated three-dimensional velocity equivalently to the available two-component two-dimensional PIV velocity measurements. The volume-averaged scale-resolved kinetic energy from the LES is sampled in three slabs, which are volumes equal to the two axial and one azimuthal PIV fields-of-view and laser sheet thickness. Prior to the comparison, the effects of sampling a cutting plane versus a slab and slabs of different thicknesses are assessed. The effects of sampling only two components and three discrete planar regions is assessed.

The influence of the bulk in-cylinder flow on the spray evolution, evaporation, fuel-air mixing and subsequent flame propagation has been studied in an optical SIDI engine. Quantitative LIF imaging of equivalence ratios was used to characterize the mixture formation and the influence of the local equivalence ratio at the time of spark on the flame propagation. Two extreme bulk flow conditions - high and low swirl - were investigated and pronounced differences in mixture homogeneity and flame propagation were found and characterized.