Search Results

The accurate representation of Internal Combustion Engine (ICE) flows via CFD is an extremely complex task: it strongly depends on a combination of highly impacting factors, such as grid resolution (both local and global), choice of the turbulence model, numeric schemes and mesh motion technique. A well-founded choice must be made in order to avoid excessive computational cost and numerical difficulties arising from the combination of fine computational grids, high-order numeric schemes and geometrical complexity typical of ICEs. The paper focuses on the comparison between different mesh motion technologies, namely layer addition and removal, morphing/remapping and overset grids. Different grid strategies for a chosen mesh motion technology are also discussed. The performance of each mesh technology and grid strategy is evaluated in terms of accuracy and computational efficiency (stability, scalability, robustness).

In-cylinder large scale and small scale structures are widely recognized to strongly influence the mixing process in HSDI Diesel engines, and therefore combustion and pollutant emissions. In particular, swirl motion intensity and temporal evolution during the intake and compression strokes must be correctly estimated to properly target the spray jets. The experimental characterization of the attitude of a valve/port assembly to promote swirl is traditionally limited to the steady flow bench, in which the analysis is carried out for fixed valve positions / fixed pressure drops and with no piston. Since flow bench analyses cannot reproduce the highly complex instantaneous flow conditions typical of actual engine operations, the use of fully-transient in-cylinder numerical simulations can become extremely useful to correctly address the engine ability to promote adequate flow structures and patterns.

CFD and FE tools are intensively adopted by engine manufacturers in order to prevent thermo-mechanical failures reducing time- and cost-to market. The capability to predict correctly the physical factors leading to damages is hence essential for their application in the industrial practice. This is even more important for last generation SI engines, where the more and more stringent need to lower fuel consumption and pollutant emissions is pushing designers to reduce engine displacement in favor of higher specific power, usually obtained by means of turbocharging. This brings to a new generation of SI engines characterized by higher and higher adiabatic efficiency and thermo-mechanical loads. A recent research highlighted the different behavior of the thermal boundary layer of such engines operated at high revving speeds and high loads if compared to the same engines operated at low loads and revving speeds or even engines with a lower specific power.

It is widely recognized that spatial and temporal evolution of both macro- and micro- turbulent scales inside internal combustion engines affect air-fuel mixing, combustion and pollutants formation. Particularly, in spark ignition engines, tumbling macro-structure induces the generation of a proper turbulence level to sustain the development and propagation of the flame front. As known, 3D-CFD codes are able to describe the evolution of the in-cylinder flow and turbulence fields with good accuracy, although a high computational effort is required. For this reason, only a limited set of operating conditions is usually investigated. On the other hand, thanks to a lower computational burden, 1D codes can be employed to study engine performance in the whole operating domain, despite of a less detailed description of in-cylinder processes. The integration of 1D and 3D approaches appears hence a promising path to combine the advantages of both.

The paper aims at providing information about the in-cylinder flow structure and its evolution of a high speed direct injection (HSDI) four valve per cylinder engine for off-highway applications. Fully transient CFD analyses by means of state-of-the-art tools and methodologies are carried out for the whole intake and compression strokes, in order to evaluate port effects on both engine permability and in-cylinder flow field evolution. Organized mean motions (i.e., swirl, tumble and squish) are investigated, trying to establish general rules in the port design optimization process, addressing relationships between the relative port orientation and the in-cylinder flow structure. Different port configurations are compared, each deriving from the rotation of the BASE port configuration on two different planes, the former being perpendicular to the cylinder axis, while the latter being parallel to the cylinder axis.

Hybrid URANS/LES turbulence modeling is rapidly emerging as a valuable complement to standard LES for full-engine multi-cycle simulation. Among the available approaches, zonal hybrids are potentially attractive due to the possibility of clearly identify URANS and LES zones, eventually introducing further zone types with dynamically switching behavior. The present work aims at evaluating the impact of different zonal configurations on the simulated flow statistics using the well-assessed TCC-III experimental engine setup. More specifically, different methods (URANS, LES or seamless DES) are applied inside the cylinder volume, as well as into the intake/exhaust ports and plenums. For each of the five tested configurations, in-cylinder flow features are compared against the reference TCC-III experimental measurements, in terms of velocity field statistics and quality indices.

The use of hydrogen in internal combustion engines is an effective approach to significantly support the reduction of CO2 emissions from the transportation sector using technically affordable solutions. The use of direct injection is the most promising approach to fully exploit hydrogen potential as a clean fuel, while preserving targets in terms of power density and emissions. In this frame, the development of an effective combustion system largely relies on the hydrogen-air mixture formation process, so to adequately control the charge stratification to mitigate pre-ignitions and knock and to minimize NOx formation. Hence, improving capabilities of designing a correct gas jet-air interaction is of paramount importance. In this paper the analysis of the evolution of a high-pressure gas jet produced by a single-hole prototype injector operated with different pressure ratios is presented.

The paper reports a wide numerical analysis of the well-known “Darmstadt engine” operated under motored condition. The engine, which features multiple optical accesses and is representative of currently made four-valve pentroof GDI production engines, is simulated using computational grids of increasing density and two widely adopted approaches to model turbulence, Reynolds Averaged Navier Stokes (RANS) and Large Eddy Simulation (LES). In the first part of the paper, attention is focused on the increase of grid density within the RANS modelling framework: both bulk-flow grid density and near-wall grid density are varied in order to analyse potentials and limitations of the different grid strategies and evaluate the trade-off between accuracy and computational cost.

Large-eddy simulation (LES) applications for internal combustion engine (ICE) flows are constantly growing due to the increase of computing resources and the availability of suitable CFD codes, methods and practices. The LES superior capability for modeling spatial and temporal evolution of turbulent flow structures with reference to RANS makes it a promising tool for describing, and possibly motivating, ICE cycle-to-cycle variability (CCV) and cycle-resolved events such as knock and misfire. Despite the growing interest towards LES in the academic community, applications to ICE flows are still limited. One of the reasons for such discrepancy is the uncertainty in the estimation of the LES computational cost. This in turn is mainly dependent on grid density, the CFD domain extent, the time step size and the overall number of cycles to be run. Grid density is directly linked to the possibility of reducing modeling assumptions for sub-grid scales.

The increasing interest in the application of Large Eddy Simulation (LES) to Internal Combustion Engines (hereafter ICEs) flows is motivated by its capability to capture spatial and temporal evolution of turbulent flow structures. Furthermore, LES is universally recognized as capable of simulating highly unsteady and random phenomena driving cycle-to-cycle variability (CCV) and cycle-resolved events such as knock and misfire. Several quality criteria were proposed in the recent past to estimate LES uncertainty: however, definitive conclusions on LES quality criteria for ICEs are still far to be found. This paper describes the application of LES quality criteria to the TCC-III single-cylinder optical engine from University of Michigan and GM Global R&D; the analyses are carried out under motored condition.

The paper critically discusses Large-Eddy Simulation (LES) potential to investigate cycle-to-cycle variability (CCV) in internal combustion engines. Particularly, the full load/peak power engine speed operation of a high-performance turbocharged GDI unit, for which ample cycle-to-cycle fluctuations were observed during experimental investigations at the engine test bed, is analyzed through a multi-cycle approach covering 25 subsequent engine cycles. In order to assess the applicability of LES within the research and development industrial practice, a modeling framework with a limited impact on the computational cost of the simulations is set up, with particular reference to the extent of the computational domain, the computational grid size, the choice of boundary conditions and numerical sub-models [1, 2, 3].

After-treatment sensors are used in the ECU feedback control to calibrate the engine operating parameters. Due to their contact with exhaust gases, especially NOx sensors are prone to soot deposition with a consequent decay of their performance. Several phenomena occur at the same time leading to sensor contamination: thermophoresis, unburnt hydrocarbons condensation and eddy diffusion of submicron particles. Conversely, soot combustion and shear forces may act in reducing soot deposition. This study proposes a predictive 3D-CFD model for the analysis of the development of soot deposition layer on the sensor surfaces. Alongside with the implementation of deposit and removal mechanisms, the effects on both thermal properties and shape of the surfaces are taken in account. The latter leads to obtain a more accurate and complete modelling of the phenomenon influencing the sensor overall performance.

The interest towards hydrogen fueling in internal combustion engines (ICEs) is rapidly growing, due to its potential impact on the reduction of the carbon footprint of the road transportation sector in a short-term scenario. While the conversion of the existing fleet to a battery-electric counterpart is highly debated in terms of both technical feasibility and life-cycle-based environmental impact, automotive researchers and technicians are exploring other solutions to reduce, if not to nullify, the carbon footprint of the existing ICE fleet. Indeed, ICE conversion to “green” fuels is seen as a promising short-term solution which does not require massive changes in powertrain production and end-of-life waste management. To better evaluate potentials and challenges of hydrogen fueling, a clear understanding of fuel injection and mixture formation prior to combustion is mandatory.

Recently, a growing interest in the development of more accurate phenomenological turbulence models is observed, since this is a key pre-requisite to properly describe the burn rate in quasi-dimensional combustion models. The latter are increasingly utilized to predict engine performance in very different operating conditions, also including unconventional valve control strategies, such as EIVC or LIVC. Therefore, a reliable phenomenological turbulence model should be able to physically relate the actuated valve strategy to turbulence level during the engine cycle, with particular care in the angular phase when the combustion takes place.

As known, reliable information about underlying turbulence intensity is a mandatory pre-requisite to predict the burning rate in quasi-dimensional combustion models. Based on 3D results reported in the companion part I paper, a quasi-dimensional turbulence model, embedded under the form of “user routine” in the GT-Power™ software, is here presented in detail. A deep discussion on the model concept is reported, compared to the alternative approaches available in the current literature. The model has the potential to estimate the impact of some geometrical parameters, such as the intake runner orientation, the compression ratio, or the bore-to-stroke ratio, thus opening the possibility to relate the burning rate to the engine architecture. Preliminarily, a well-assessed approach, embedded in GT-Power commercial software v.2016, is utilized to reproduce turbulence characteristics of a VVA engine.