Search Results

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.

Aim of this work is to develop a general purpose model for combustion and knocking prediction in SI engines, by coupling a thermo-fluid dynamic model for engine simulation with a general detailed kinetic scheme, including the low-temperature oxidation mechanism, for the prediction of the auto-ignition behavior of hydrocarbons. A quasi-D approach is used to describe the in-cylinder thermodynamic processes, applying the conservation of mass and energy over the cylinder volume, modeled as a single open system. The complex chemistry model has been embedded into the code, by using the same integration algorithm for the conservation equations and the reacting species, and taking into account their mutual interaction in the energy balance. A flame area evolution predictive approach is used to evaluate the turbulent flame front propagation as function of the engine operating parameters.

In order to fulfill the more and more stringent emission levels (Euro V, SULEV…), catalytic converter light-off time has to be reduced as much as possible. Consequently, all the parts upstream of the catalytic converter have to be designed in order to minimize the gas heat loss. As a matter of fact, considering the emission performance, all components of the hot end contribute to a better after-treatment. In this study, we focus on the exhaust manifold, that has a major contribution to the thermal mass upstream of the catalyst. The study carried out aims at highlighting the impact of fabricated manifold length and thickness on emissions and engine performance. Several manifold designs, dedicated to different naturally aspirated gasoline engine applications, have been tested on a dynamic engine bench or chassis dyno. Emission results were also supported by temperature measurements.

This paper describes a detailed analysis of the unsteady flows in the intake and exhaust systems of a modern four-cylinder, turbocharged Diesel engine by means of advanced numerical tools and experimental measurements. In particular, a 1D-3D integrated fluid dynamic model, based on the GASDYN (1D) and Lib-ICE (3D) codes, has been developed and applied for the schematization of the geometrical domain and the prediction of the wave motion in the whole intake and the exhaust systems, including the air cleaner, the intercooler, the after-treatment devices and the silencers. Firstly, a detailed 1D simulation has been carried out to predict the pressure pulses, average pressures and temperatures in several cross-sections of the pipe systems for different speeds and loads, considering the complex geometry of the air filter, the intake manifold, the intercooler and the exhaust manifold.

The paper describes the development and validation of a quasi-dimensional multi-zone combustion model for Gasoline Direct Injection engines. The model has been embedded in the 1D thermo-fluid-dynamic code for the simulation of the whole engine system named GASDYN and developed by the authors [1, 2 and 3]. The GDI engine combustion model solves mass, energy and species equations using a 4th order Runge-Kutta integration method; the fuel spray is initially divided into a number of zones fixed regardless of the injected amount and the time step, considering the following break-up, droplet evaporation and air entrainment in each single zone. Experimental correlations have been used for the spray penetration and spatial information. Once the ignition begins it is assumed that the flame propagates spherically, evaluating its velocity by means of a fractal combustion approach and considering the local air-fuel ratio, which is the result of the spray evolution within the combustion chamber.

The paper describes the research work carried out on the thermo-fluid dynamic modeling of an S.I. engine coupled to the vehicle in order to predict the engine and tailpipe emissions during the ECE European driving cycle. The numerical code GASDYN has been extended to simulate the engine + vehicle operation during the first 90 seconds of the NEDC driving cycle, taking account of the engine and exhaust system warm-up after the cold start. The chemical composition of the engine exhaust gas is calculated by means of a thermodynamic multi-zone combustion model, augmented by kinetic emission sub-models for the prediction of pollutant emissions. A simple procedure has been implemented to model the vehicle dynamic behavior (one degree of freedom model). A closed-loop control strategy (proportional-derivative) has been introduced to determine the throttle opening angle, corresponding to the engine operating point when the vehicle is following the ECE cycle.

Nowadays 1D fluid dynamic models are widely used by engine designers, since they can give sufficiently accurate predictions in short times, allowing to support the optimization and development work of any prototype. According to the last requirements in terms of pollutant emission control, some enhancements have been introduced in the 1D code GASDYN, to improve its ability in predicting the composition of the exhaust gas discharged by the cylinders and the transport of the chemical species along the exhaust system. The main aspects of the methods adopted to model the combustion process and the related formation of pollutants are described in the paper. To account for the burnt gas stratification, two different approaches have been proposed, depending on the expected turbulence levels inside the combustion chamber. The reliability of the simulation of the pollutant formation process has been enhanced by the integration of the thermodynamic module with the Chemkin code.

The scope of the work presented in this paper was to apply the latest open source CFD achievements to design a state of the art, direct-injection (DI), heavy-duty, natural gas-fueled engine. Within this context, an initial steady-state analysis of the in-cylinder flow was performed by simulating three different intake ducts geometries, each one with seven different valve lift values, chosen according to an estabilished methodology proposed by AVL. The discharge coefficient (Cd) and the Tumble Ratio (TR) were calculated in each case, and an optimal intake ports geometry configuration was assessed in terms of a compromise between the desired intensity of tumble in the chamber and the satisfaction of an adequate value of Cd. Subsequently, full-cycle, cold-flow simulations were performed for three different engine operating points, in order to evaluate the in-cylinder development of TR and turbulent kinetic energy (TKE) under transient conditions.

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 automotive industry is striving to adopt model-based engine design and optimization procedures to reduce development time and costs. In this scenario, first-principles gas dynamic models predicting the mass, energy and momentum transport in the engine air path system with high accuracy and low computation effort are extremely important today for performance prediction, optimization and cylinder charge estimation and control. This paper presents a comparative study of two different modeling approaches to predict the one-dimensional unsteady compressible flow in the engine air path system. The first approach is based on a quasi-3D finite volume method, which relies on a geometrical reconstruction of the calculation domain using networks of zero-dimensional elements. The second approach is based on a model-order reduction procedure that projects the nonlinear hyperbolic partial differential equations describing the 1D unsteady flow in engine manifolds onto a predefined basis.

This paper describes the development and application of the 1D thermo-fluid dynamic research code GASDYN to the simulation of a Lamborghini 12 cylinder, V 60°, 6.2 L automotive S.I. engine. The model has been adopted to carry out an integrated simulation (thermodynamic, fluid dynamic and chemical) of the engine coupled to its intake and exhaust manifolds, in order to predict not only the wave motion in the ducts and its influence on the cylinder gas exchange process, but also the in-cylinder combustion process and the pollutant emission concentration along the exhaust system. The gas composition in the exhaust pipe system is dictated by the cylinder discharge process, after the calculation of the combustion via a thermodynamic multi-zone model, based on a “fractal geometry” approach.

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.

This paper deals with some recent advances in the field of 1D fluid dynamic modeling of unsteady reacting flows in complex s.i. engine pipe-systems, involving a catalytic converter. In particular, a numerical simulation code has been developed to allow the simulation of chemical reactions occurring in the catalyst, in order to predict the chemical specie concentration in the exhaust gas from the cylinder to the tailpipe outlet, passing through the catalytic converter. The composition of the exhaust gas, discharged by the cylinder and then flowing towards the converter, is calculated by means of a thermodynamic two-zone combustion model, including emission sub-models. The catalytic converter can be simulated by means of a 1D fluid dynamic and chemical approach, considering the laminar flow in each tiny channel of the substrate.

The paper describes the 1-D fluid dynamic modeling of unsteady flows with chemical specie tracking in the ducts of a four-cylinder s.i. automotive engine, to predict the composition of the exhaust gas reaching the catalyst inlet. A comprehensive simulation model, based on classical and innovative numerical techniques for the solution of the governing equations, has been developed. The non-traditional shock-capturing CE-SE (Conservation Element-Solution Element) method has been extended to deal with the propagation of chemical species. A comparison of the MacCormack method plus FCT or TVD algorithms with the CE-SE method has pointed out the superiority of the latter scheme in the propagation of contact discontinuities. A realistic composition of the exhaust products in the cylinder, evaluated by a two-zone combustion model including emission sub-models, has been imposed at the opening of the exhaust valve, considering the effect of short-circuit of air during valve overlap.

Turbocharging is playing today a fundamental role not only to improve automotive engine performance, but also to reduce fuel consumption and exhaust emissions for both Spark Ignition and Diesel engines. Dedicated experimental investigations on turbochargers are therefore necessary to assess a better understanding of its performance. The availability of experimental information on turbocharger steady flow performance is an essential requirement to optimize the engine-turbocharger matching, which is usually achieved by means of simulation models. This aspect is even more important when referred to the turbine efficiency, since its swallowing capacity can be accurately evaluated through the measurement of mass flow rate, inlet temperature and pressure ratio across the machine.

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.

The dynamics and evolution of turbulent structures inside an engine-like geometry are investigated by means of Large Eddy Simulation. A simplified configuration consisting of a flat-top cylinder head with a fixed, axis-centered valve and low-speed piston has been simulated by the finite volume CFD code OpenFOAM®; the standard version of the software has been extended to include the compressible WALE subgrid-scale model, models for the generation of synthetic turbulence, some improvements to the mesh motion strategy and algorithms for LES data post-processing. In order to study both the initial transient and the quasi- steady operating conditions, ten complete engine cycles have been simulated. Phase and spatial averages have been performed over cycles three to ten in order to extract first and second moment of velocity; these quantities have then been used to validate the numerical procedure by comparison against experimental data.