Search Results

Proton Exchange Fuel Cells (PEMFCs) are considered one of the most prominent technologies to decarbonize the transportation sector, with emphasis on long-haul/long-range trucks, off-highway, maritime and railway. The flow field of reactants is dictated by the layout of machined channels in the bipolar plates, and several established designs (e.g., parallel channels, single/multi-pass serpentine) coexist both in research and industry. In this context, the flow behavior at cathode embodies multiple complexities, namely an accurate control of the inlet/outlet humidity for optimal membrane hydration, pressure losses, water removal at high current density, and the limitation of laminar regime. However, a robust methodology is missing to compare and quantify such aspects among the candidate designs, resulting in a variety of configurations in use with no justification of the specific choice.

With growing concern about global warming, alternatives to fossil fuels in internal combustion engines are searched. In this context, hydrogen is one of the most interesting fuels as it shows excellent combustion properties such as laminar flame speed and energy density. In this work a CFD methodology for 3D-CFD in-cylinder simulations of engine combustion is proposed and its predictive capabilities are validated against test-bench data from a direct injection spark-ignition (DISI) prototype. The original engine is a naturally aspirated, single cylinder compression ignition (Diesel fueled) unit. It is modified substituting the Diesel injector with a spark plug, adding two direct gas injectors, and lowering the compression ratio to run with hydrogen fuel. A 3D-CFD model is built, embedding in-house developed ignition and heat transfer models besides G-equation one for combustion.

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.

Knock is one of the main limitations on Spark-Ignited (SI) Internal Combustion Engine (ICE) performance and efficiency and so has been the object of study for over one hundred years. Great strides have been made in terms of understanding in that time, but certain rather elementary practical problems remain. One of these is how to interpret if a running engine is knocking and how likely this is to result in damage. Knocking in a development environment is typically quantified based on numerical descriptions of the high frequency content of a cylinder pressure signal. Certain key frequencies are observed, which Draper [1] explained with fundamental acoustic theory back in 1935. Since then, a number of approaches of varying complexity have been employed to correlate what is happening within the chamber with what is measured by a pressure transducer.

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 paper aims at defining a methodology for the prediction and understanding of knock tendency in internal combustion engine piston crevices by means of CFD simulations. The motivation for the analysis comes from a real design requirement which appeared during the development of a new high performance SI unit: it is in fact widely known that, in high performance engines (especially the turbocharged ones), the high values of pressure and temperature inside the combustion chamber during the engine cycle may cause knocking phenomena. “Standard” knock can be easily recognized by direct observation of the in-cylinder measured pressure trace; it is then possible to undertake proper actions and implement design and control improvements to prevent it with relatively standard 3D-CFD analyses.

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

The detailed study of part-load conditions is essential to characterize engine-out emissions in key operating conditions. The relevance of part-load operations is further emphasized by the recent regulations such as the new WLTP standard. Combustion development at part-load operations depends on a complex interplay between moderate turbulence levels (low engine speed and tumble ratio), low in-cylinder pressure and temperature, and stoichiometric-to-lean mixture quality (to maximize fuel efficiency). From a modelling standpoint, the reduced turbulence intensity compared to full-load operations complicates the interaction between different sub-models (e.g., reconsideration of the flamelet hypothesis adopted by common combustion models). In this article, the authors focus on chemistry-based simulations for laminar flame speed of gasoline surrogates at conditions typical of part-load operations. The analysis is an extension of a previous study focused on full-load operations.

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.

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.

High power-density Diesel engines are characterized by remarkable thermo-mechanical loads. Therefore, compared to spark ignition engines, designers are forced to increase component strength in order to avoid failures. 3D-CFD simulations represent a powerful tool for the evaluation of the engine thermal field and may be used by designers, along with FE analyses, to ensure thermo-mechanical reliability. The present work aims at providing an integrated in-cylinder/CHT methodology for the estimation of a Diesel engine thermal field. On one hand, in-cylinder simulations are fundamental to evaluate not only the integral amount of heat transfer to the combustion chamber walls, but also its point-wise distribution. To this specific aim, an improved heat transfer model based on a modified thermal wall function is adopted to estimate correctly wall heat fluxes due to combustion.

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.

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.

New SI engine generations are characterized by a simultaneous reduction of the engine displacement and an increase of the brake power; such targets are achieved through the adoption of several techniques such as turbocharging, direct fuel injection, variable valve timing and variable port lengths. This design approach, called “downsizing”, leads to a marked increase in the thermal loads acting on the engine components, in particular on those facing the combustion chamber. Hence, an accurate evaluation of the thermal field is of primary importance in order to avoid mechanical failures. Moreover, the correct evaluation of the temperature distribution improves the prediction of pointwise abnormal combustion onset.

In spark-ignited direct-injected engines, the formation of fuel pools on the piston is one of the major promoters of unburnt hydrocarbons and soot: in order to comply with the increasingly stringent emission regulations (EU6 and forthcoming), it is therefore necessary to limit fuel deposit formation. The combined use of advanced experimental techniques and detailed 3D-CFD simulations can help to understand the mechanisms driving fuel pool formation. In the paper, a combined experimental and numerical characterization of pool formation in a GDI engine is carried out to investigate and understand the complex interplay of all the mentioned factors. In particular, a low-load low-rpm engine operation is investigated for different ignition phasing, and the impact of both fuel formulation and instantaneous piston temperature variations in the CFD analyses are evaluated.

A new generation of highly downsized SI engines with specific power output around or above 150 HP/liter is emerging in the sport car market sector. Technologies such as high-boosting, direct injection and downsizing are adopted to increase power density and reduce fuel consumption. To counterbalance the increased risks of pre-ignition, knock or mega-knock, currently made turbocharged SI engines usually operate with high fuel enrichments and delayed (sometimes negative) spark advances. The former is responsible for high fuel consumption levels, while the latter induce an even lower A/F ratio (below 11), to limit the turbine inlet temperature, with huge negative effects on BSFC. A possible solution to increase knock resistance is investigated in the paper by means of 3D-CFD analyses: water/methanol emulsion is port-fuel injected to replace mixture enrichment while preserving, if not improving, indicated mean effective pressure and knock safety margins.

Engine downsizing is gaining popularity in the high performance engine market sector, where a new generation of highly downsized engines with specific power outputs around or above 150 HP/litre is emerging. High-boost and downsizing, adopted to increase power density and reduce fuel consumption, have to face the increased risks of pre-ignition, knock or mega-knock. To counterbalance autoignition of fuel/air mixture, such engines usually operate with high fuel enrichments and delayed (sometimes negative) spark advances. The former is responsible for high fuel consumption levels, while the latter reduces performance and induces an even lower A/F ratio (below 11), to limit the turbine inlet temperature, with huge negative effects on BSFC.

In the present paper, two different methodologies are adopted and critically integrated to analyze the knock behavior of a last generation small size spark ignition (SI) turbocharged VVA engine. Particularly, two full load operating points are selected, exhibiting relevant differences in terms of knock proximity. On one side, a knock investigation is carried out by means of an Auto-Regressive technique (AR model) to process experimental in-cylinder pressure signals. This mathematical procedure is used to estimate the statistical distribution of knocking cycles and provide a validation of the following 1D-3D knock investigations. On the other side, an integrated numerical approach is set up, based on the synergic use of 1D and 3D simulation tools. The 1D engine model is developed within the commercial software GT-Power™. It is used to provide time-varying boundary conditions (BCs) for the 3D code, Star-CD™.