Search Results

The reliable prediction of pollutant emissions generated by IC engine powertrains during the WLTP driving cycle is a key aspect to test and optimize different configurations, in order to respect the stringent emission limits. This work describes the application of an integrated modeling tool in a co-simulation environment, coupling a 1D fluid dynamic code for engine simulation with a specific numerical code for aftertreatment modelling by means of a robust numerical approach, to achieve a complete methodology for detailed simulations of driving cycles. The main goal is to allow an accurate 1D simulation of the unsteady flows along the intake and exhaust systems and to apply advanced thermodynamic combustion models for the calculation of cylinder-out emissions.

Nowadays, detailed kinetics is necessary for a proper estimation of both flame structure and pollutant formation in compression ignition engines. However, large mechanisms and the need to include turbulence/chemistry interaction introduce significant computational overheads. For this reason, tabulated kinetics is employed as a possible solution to reduce the CPU time even if table discretization is generally limited by memory occupation. In this work the authors applied tabulated homogeneous reactors (HR) in combination with different turbulent-chemistry interaction approaches to model non-premixed turbulent combustion. The proposed methodologies represent good compromises between accuracy, required memory and computational time. The experimental validation was carried out by considering both constant-volume vessel and Diesel engine experiments.

In the last decades numerical simulations have become reliable tools for the design and the optimization of silencers for internal combustion engines. Different approaches, ranging from simple 1D models to detailed 3D models, are nowadays commonly applied in the engine development process, with the aim to predict the acoustic behavior of intake and exhaust systems. However, the acoustic analysis is usually performed under the hypothesis of infinite stiffness of the silencer walls. This assumption, which can be regarded as reasonable for most of the applications, can lose validity if low wall thickness are considered. This consideration is even more significant if the recent trends in the automotive industry are taken into account: in fact, the increasing attention to the weight of the vehicle has lead to a general reduction of the thickness of the metal sheets, due also to the adoption of high-strength steels, making the vibration of the components a non negligible issue.

Fuel lean combustion and exhaust gas dilution are known to increase the thermal efficiency and reduce NOx emissions. In this study, experiments are performed to understand the effect of equivalence ratio on flame kernel formation and flame propagation around the spark plug for different low turbulent velocities. A series of experiments are carried out for propane-air mixtures to simulate engine-like conditions. For these experiments, equivalence ratios of 0.7 and 0.9 are tested with 20 percent mass-based exhaust gas recirculation (EGR). Turbulence is generated by a shrouded fan design in the vicinity of J-spark plug. A closed loop feedback control system is used for the fan to generate a consistent flow field. The flow profile is characterized by using Particle Image Velocimetry (PIV) technique. High-speed Schlieren visualization is used for the spark formation and flame propagation.

In this paper a three-dimensional time-domain CFD approach has been employed to predict and analyze the acoustic attenuation performance of complex perforated muffler geometries, where strong 3D effects limit the validity of the use of one-dimensional models. A pressure pulse has been imposed at the inlet to excite the wave motion, while unsteady flow computation have been performed to acquire the time histories of the pressures upstream and downstream of the silencer. Pressures in the time domain have been then transformed to acoustic pressures in the frequency domain, to predict the transmission loss.

The implementation and the combination of advanced boundary conditions and subgrid scale models for Large Eddy Simulation (LES) in the multi-dimensional open-source CFD code OpenFOAM® are presented. The goal is to perform reliable cold flow LES simulations in complex geometries, such as in the cylinders of internal combustion engines. The implementation of a boundary condition for synthetic turbulence generation upstream of the valve port and of the compressible formulation of the Wall-Adapting Local Eddy-viscosity sgs model (WALE) is described. The WALE model is based on the square of the velocity gradient tensor and it accounts for the effects of both the strain and the rotation rate of the smallest resolved turbulent fluctuations and it recovers the proper y₃ near-wall scaling for the eddy viscosity without requiring dynamic procedure; hence, it is supposed to be a very reliable model for ICE simulation.

This paper is focused on the development and application of a CFD methodology that can be applied to predict the fuel-air mixing process in stratified charge, sparkignition engines. The Eulerian-Lagrangian approach was used to model the spray evolution together with a liquid film model that properly takes into account its effects on the fuel-air mixing process into account. However, numerical simulation of stratified combustion in SI engines is a very challenging task for CFD modeling, due to the complex interaction of different physical phenomena involving turbulent, reacting and multiphase flows evolving inside a moving geometry. Hence, for a proper assessment of the different sub-models involved a detailed set of experimental optical data is required. To this end, a large experimental database was built by the authors.

Recent advances in the field of numerical modeling of unsteady reacting flows in the exhaust system of s.i. engines are presented in the paper. In particular, it is shown that the integration of a suitable chemical and thermal model for the catalytic converter within a 1D fluid dynamic simulation code has allowed the prediction of the exhaust gas composition from the cylinder to the tailpipe outlet, considering its variation across the catalyst. The composition of the exhaust gas, discharged by the cylinder, is calculated by means of a two- zone combustion model, including emission sub-models. The catalytic converter is simulated by a 1D fluid dynamic and chemical approach, considering laminar flow in each tiny channel of the substrate and chemical reactions in the solid phase, within the wash-coat. The predicted reaction rates are used to determine the specie source terms to be included in the one-dimensional fluid dynamic conservation equations.

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.

Computational fluid dynamics represents a useful tool to support the design and development of Heavy Duty Engines, making possible to test the effects of injection strategies and combustion chamber design for a wide range of operating conditions. Predictive models are required to ensure accurate estimations of heat release and the main pollutant emissions within a limited amount of time. For this reason, both detailed chemistry and turbulence chemistry interaction need to be included. In this work, the authors intend to apply combustion models based on tabulated kinetics for the prediction of Diesel combustion in Heavy Duty Engines. Four different approaches were considered: well-mixed model, presumed PDF, representative interactive flamelets and flamelet progress variable. Tabulated kinetics was also used for the estimation of NOx emissions.

Extensive prior art within the Engine Combustion Network (ECN) using a Bosch single axial-hole injector called ‘Spray A’ in constant-volume vessels has provided a solid foundation from which to evaluate modeling tools relevant to spray combustion. In this paper, a new experiment using a Bosch three-hole nozzle called ‘Spray B’ mounted in a 2.34 L heavy-duty optical engine is compared to sector-mesh engine simulations. Two different approaches are employed to model combustion: the ‘well-mixed model’ considers every cell as a homogeneous reactor and employs multi-zone chemistry to reduce the computational time. The ‘flamelet’ approach represents combustion by an ensemble of laminar diffusion flames evolving in the mixture fraction space and can resolve the influence of mixing, or ‘turbulence-chemistry interactions,’ through the influence of the scalar dissipation rate on combustion.

The acoustic simulation of internal combustion engine exhaust systems is an important aspect to meet customer expectations and legislation targets. One dimensional gas dynamic simulation tools are used for the calculation of the exhaust orifice noise in the early stages of the engine development process. This includes the prediction of the acoustic performance of individual components in the exhaust line. One common element used in exhaust systems to increase the acoustic damping is the plug flow muffler. This study looks at the prediction of acoustic performance of various plug mufflers at different flow velocities. These include a single plug muffler, a double plug muffler and an eccentric plug muffler with different porosities for the perforated sections. To this purpose a generic 3D cell approach was developed and applied.

A detailed understanding of Gasoline Direct Injection (GDI) techniques applied to spark-ignition (SI) engines is necessary as they allow for many technical advantages such as increased power output, higher fuel efficiency and better cold start performances. Within this context, the extensive validation of multi-dimensional models against experimental data is a fundamental task in order to achieve an accurate reproduction of the physical phenomena characterizing the injected fuel spray. In this work, simulations of different Engine Combustion Network (ECN) Spray G conditions were performed with the Lib-ICE code, which is based on the open source OpenFOAM technology, by using a RANS Eulerian-Lagrangian approach to model the ambient gas-fuel spray interaction.

Direct-injection represents a consolidated technology to increase performance and efficiency in spark-ignition engines. It reduces the knock tendency and makes engine downsizing possible through the use of turbocharging. Better control of CO and HC emissions at cold-start is also ensured since there is no wall-impingement in the intake port. However, to take advantages of all the theoretical benefits derived from GDI technology, detailed investigations of both fuel-air mixing and combustion processes are necessary to extend the stratified charge operations in the engine map and to reduce soot emissions, that are now severely regulated by emission standards. In this work, the authors developed a CFD methodology to investigate and optimize the fuel-air mixing process in direct-injection, spark-ignition engines. The Eulerian-Lagrangian approach is used to model the evolution of the fuel spray emerging from a multi-hole injector.

In the most elementary treatment of plane-wave reflection at the open end of a duct system, it is often assumed that the ends are pressure nodes. This implies that pressure is assumed as a constant at the open end termination and that steady flow boundary condition is supposed as instantaneously established. While this simplifying assumption seems reasonable, it does not consider any radiation of acoustic energy from the duct into the surrounding free space; hence, an error in the estimation of the effects of the flow on the acoustical response of an open-end duct occurs. If radiation is accounted, a complicated three-dimensional wave pattern near the duct end is established, which tends to readjust the exit pressure to its steady-flow level. This adjustment process is continually modified by further incident waves, so that the effective instantaneous boundary conditions which determine the reflected waves depend on the flow history.

The current development to set up an automatic procedure for automatic mesh generation and automatic mesh motion for internal combustion engine simulation in OpenFOAM®-2.2.x is here described. In order to automatically generate high-quality meshes of cylinder geometries, some technical issues need to be addressed: 1) automatic mesh generation should be able to control anisotropy and directionality of the grid; 2) during piston and valve motion, cells and faces must be introduced and removed without varying the overall area and volume of the cells, to avoid conservation errors. In particular, interpolation between discrete fields is frequent in computational physics: the use of adaptive and non-conformal meshes necessitates the interpolation of fields between different mesh regions. Interpolation problems also arise in areas such as model coupling, model initialization and visualisation.

Increasing demands on the capabilities of engine simulation and the ability to accurately predict both performance and acoustics has lead to the development of several numerical tools to help engine manufacturers during the prototyping stage. The aid of CFD tools (3D and 1D) can remarkably reduce the duration and the costs of this stage. The need of achieving good accuracy, along with acceptable computational runtime, has given the spur to the development of a geometry based quasi-3D approach. This is designed to model the acoustics and the fluid dynamics of both intake and exhaust system components used in internal combustion engines. Models of components are built using a network of quasi-3D cells based primarily on the geometry of the system. The solution procedure is based on an explicitly time marching staggered grid approach making use of a flux limiter to prevent numerical instabilities.

The paper illustrates and validates a novel predictive combustion model for the estimation of performances and pollutant production in CI engines. The numerical methodology was developed by the authors for near real-time applications, while aiming at an accurate description of the air mixing process by means of a multi-zone approach of the air-fuel mass. Charge stratification is estimated via a 2D representation of the fuel spray distribution that is numerically derived by an axial one-dimensional control-volume description of the direct injection. The radial coordinate of each control volume is reconstructed a posteriori by means of a local distribution function. Fuel mass clustered in each zone is further split in ‘liquid’, ‘unburnt’ and ‘burnt’ sub-zones, given the local properties of the fuel spray control volumes with respect to space-time location of modelled ignition delay, liquid length, and flame lift-off.

This work has the objective to present the extension of a novel quasi-dimensional model, developed to simulate the combustion process in diesel Compression Ignition (CI) engines, to describe this process when Dimethyl ether (DME) is used as fuel. DME is a promising fuel in heavy-duty CI engines application thanks to its high Cetane Number (CN), volatility, high reactivity, almost smokeless combustion, lower CO2 emission and the possibility to be produced with renewable energy sources. In this paper, a brief description of the thermodynamic model will be presented, with particular attention to the implementation of the Tabulated Kinetic Ignition (TKI) model, and how the various models interact to simulate the combustion process. The model has been validated against experimental data derived from constant-volume DME combustion, in this case the most important parameters analyzed and compared were the Ignition Delay (ID) and Flame Lift Off Length (FLOL).

The paper focuses on the development of a mesh moving method based on non-conformal topologically changing grids applied to the simulation of IC engines, where the prescribed motion of piston and valves is accomplished by rigidly translating the sub-domain representing the moving component. With respect to authors previous work, a more robust and efficient algorithm to handle the connectivity of non-conformal interfaces and a mesh-motion solver supporting multiple layer addition/removal of cells, to decouple the time-step constraints of the mesh motion and of the fluid dynamics, has been implemented as a C++ library to extend the already existing classes for dynamic mesh handling of the finite-volume, open-source CFD code OpenFOAM®. Other new features include automatic decomposition of large multiple region domains to preserve processors load balance with topological changes for parallel computations and additional tools for automatic preprocessing and case setup.