Search Results

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.

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.

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 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).

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.

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.

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.

The evaluation of the steady-flow discharge coefficient of ICE port assemble is known to be very sensitive to the capability of the turbulence sub-models in capturing the boundary layer dynamics. Despite the fact that the intrinsically unsteady phenomena related to flow separation claim for LES approach, the present paper aims to demonstrate that RANS simulation can provide reliable design-oriented results by using low-Reynolds cubic k-ε turbulence models. Different engine intake port assemblies and pressure drops have been simulated by using the CFD STAR-CD code and numerical results have been compared versus experiments in terms of both global parameters, i.e. the discharge coefficient, and local parameters, by means of static pressure measurements along the intake port just upstream of the valve seat. Computations have been performed by comparing two turbulence models: Low-Reynolds cubic k-ε and High-Reynolds cubic k-ε.

The paper is focused on the influence of the eddy viscosity turbulence models (EVM) in CFD three-dimensional simulations of steady turbulent engine intake flows in order to assess their reliability in predicting the discharge coefficient. Results have been analyzed by means of the comparison with experimental measurements at the steady flow bench. High Reynolds linear and non-linear and RNG k-ε models have been used for simulation, revealing the strong influence of both the constitutive relation and the ε-equation formulation on the obtained results, while limits in the applicability of more sophisticated near-wall approaches are briefly discussed in the paper. Due to the extreme complexity of typical ICE flows and geometries, the analysis of the behavior of EVM turbulence models has been subsequently applied to a test-case available in literature, i.e. a high-Reynolds compressible flow over a inclined backward facing step (BFS).