Search Results

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.

Hydrogen engines are currently considered as a viable solution to preserve the internal combustion engine as a power unit for vehicle propulsion. In particular, lean-burn gasoline Spark-Ignition (SI) engines have been a major subject of investigations due to the reduced emission levels and high thermodynamic efficiency. This strategy is suitable for the purpose of passenger car applications and cannot be tailored in the field of high performance engine, where the air mass delivered would require oversized turbocharging systems or more complex charging solutions. For this reason, the range of stoichiometric feeding condition is explored in the high performance engine, leading to the consequent issue of abatement of pollutant emissions. In this work a 1D model will be applied to the modeling of a V8 engine fueled with DI of hydrogen. The engine has been derived by a gasoline configuration and adapted to hydrogen in such a way to keep the same performance.

This paper describes some recent advances of the research work concerning the 1D fluid dynamic modeling of unsteady reacting flows in s.i. engine pipe-systems, including pre-catalysts and main catalysts. The numerical model GASDYN developed in previous work has been further enhanced to enable the simulation of the catalyst. The main chemical reactions occurring in the wash-coat have been accounted in the model, considering the mass transfer between gas and solid phase. The oxidation of CO, C3H6, C3H8, H2 and reduction of NO, the steam-reforming reactions of C3H6, C3H8, the water-gas shift reaction of CO have been considered. Moreover, an oxygen-storage sub-model has been introduced, to account for the behavior of Cerium oxides. A detailed thermal model of the converter takes into account the heat released by the exothermic reactions as a source term in the heat transfer equations. The influence of the insulating mat is accounted.

The continuous pursuit of higher combustion efficiencies, as well as the possible usage of synthetic fuels with different properties than fossil-ones, require reliable and low-cost numerical approaches to support and speed-up engines industrial design. In this context, SI engines operated with homogeneous ultra-lean mixtures both characterized by a classical ignition configuration or equipped with an active prechamber represent the most promising solutions. In this work, for the classical ignition arrangement, a 3DCFD strategy to model the impact of the ignition system type on the CCV is developed using the RANS approach for turbulence modelling. The spark-discharge is modelled through a set of Lagrangian particles, whose velocity is modified with a zero-divergence perturbation at each discharge event, then evolved according to the Simplified Langevin Model (SLM) to simulate stochastic interactions with the surrounding gas flow.

Active prechamber combustion systems for SI engines represent a feasible and effective solution in reducing fuel consumption and pollutant emissions for both marine and ground heavy-duty engines. However, reliable and low-cost numerical approaches need to be developed to support and speed-up their industrial design considering their geometry complexity and the involved multiple flow length scales. This work presents a CFD methodology based on the RANS approach for the simulation of active prechamber spark-ignition engines. To reduce the computational time, the gas exchange process is computed only in the prechamber region to correctly describe the flow and mixture distributions, while the whole cylinder geometry is considered only for the power-cycle (compression, combustion and expansion). Outside the prechamber the in-cylinder flow field at IVC is estimated from the measured swirl ratio.

The paper describes the development of a reduced kinetic scheme for the evaluation of the main chemical species (particularly NO and CO) in premixed turbulent flame and its application to a quasi-dimensional combustion model for spark ignition engines. The proposed mechanism is based on the kinetic solution of three transport equations for NO, CO and H, coupled with the partial equilibrium of the so-called water-shuffle equations to derive the OH, O and H2 concentrations. The remaining species are computed applying the element conservation, while the required prompt levels were determined by a separate chemical 1D code for laminar combustion. The proposed chemical scheme was locally validated, considering a turbulent flame inside a premixed flow of air and methane, ignited by a parallel flow of hot gases, by means of a CFD simulation. Successively, it was embedded into a quasi-D thermodynamic combustion model developed by the authors for the simulation of S.I. and C.I. engines.

Modeling plume interaction and collapse for direct-injection gasoline sprays is important because of its impact on fuel-air mixing and engine performance. Nevertheless, the aerodynamic interaction between plumes and the complicated two-phase coupling of the evaporating spray has shown to be notoriously difficult to predict. With the availability of high-speed (100 kHz) Particle Image Velocimetry (PIV) experimental data, we compare velocity field predictions between plumes to observe the full temporal evolution leading up to plume merging and complete spray collapse. The target “Spray G” operating conditions of the Engine Combustion Network (ECN) is the focus of the work, including parametric variations in ambient gas temperature. We apply both LES and RANS spray models in different CFD platforms, outlining features of the spray that are most critical to model in order to predict the correct aerodynamics and fuel-air mixing.

A correct prediction of the initial stages of the combustion process in SI engines is of great importance to understand how local flow conditions, fuel properties, mixture stratification and ignition affect the in-cylinder pressure development and pollutant formation. However, flame kernel growth is governed by many interacting processes including energy transfer from the electrical circuit to the gas phase, interaction between the plasma channel and the flow field, transition between different combustion regimes and gas expansion at very high temperatures. In this work, the authors intend to present a comprehensive, multi-dimensional model that can be used to predict the initial combustion stages in SI engines. In particular, the spark channel is represented by a set of Lagrangian particles where each one of them acts as a single flame kernel.

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.

The present work describes the numerical modeling of medium-speed marine engines, operating in a fumigated dual-fuel mode, i.e. with the second fuel injected in the ports. This engine technology allows reducing engine-out emissions while maintaining the engine efficiency and can be fairly easily retrofitted from current diesel engines. The main premixed fuel that is added can be a low-carbon one and can additionally be of a renewable nature, thereby reducing or even completely removing the global warming impact. To fully optimize the operational parameters of such a large marine engine, computational fluid dynamics can be very helpful. Accurately describing the combustion process in such an engine is key, as the prediction of the heat release and the pollutant formation is crucial. Auto-ignition of the diesel fuel needs to be captured, followed by the combustion and flame propagation of the premixed fuel.

The accurate prediction of pollutant emissions generated by IC engines is a key aspect to guarantee the respect of the emission regulation legislation. This paper describes the approach followed by the authors to achieve a strict numerical coupling of two different 1D modeling tools in a co-simulation environment, aiming at a reliable calculation of engine-out and tailpipe emissions. The main idea is to allow an accurate 1D simulation of the unsteady flows and wave motion inside the intake and exhaust systems, without resorting to an over-simplified geometrical discretization, and to rely on advanced thermodynamic combustion models and kinetic sub-models for the calculation of cylinder-out emissions. A specific fluid dynamic approach is then used to track the chemical composition along the exhaust duct-system, in order to evaluate the conversion efficiency of after-treatment devices, such as TWC, GPF, DPF, DOC, SCR and so on.

In this paper, numerical simulations of an automotive-size optical diesel engine have been conducted employing the Reynolds-Averaged Navier-Stokes (RANS) equations with the standard k-ε turbulence model and a reduced n-heptane chemical mechanism implemented in OpenFOAM. The current paper builds on a previous work where the model has been validated for the same engine using optical diagnostic data. The present study investigates numerically the influence of different operating conditions - relevant for modern diesel engines - on the mixture formation development under non-reactive conditions as well as low- and high-temperature ignition behaviour and flame evolution in the presence of strong jet-wall interactions typically encountered in automotive-size diesel engines. Also, emissions of CO and unburned hydrocarbons (UHC) are considered.

In this paper, computation fluid dynamics (CFD) simulations are employed to describe the effect of flow parameters on the formation of soot and NOx in a heavy duty engine under low load and high load. The complexity of diesel combustion, specially when soot, NOx and other emissions are of interest, requires using a detailed chemical mechanism to have a correct estimation of temperature and species distribution. In this work, Multiple Representative Interactive Flamelets (MRIF) method is employed to describe the chemical reactions, ignition, flame propagation and emissions in the engine. A phenomenological model for soot formation, including soot nucleation, coagulation and oxidation with O2 and OH is incorporated into the flamelet combustion model. Different strategies for modelling NOx are chosen to take into account the longer time scale for NOx formation. The numerical results are compared with experimental data to show the validity of the model for the cases under study.

The 4th Workshop of the Engine Combustion Network (ECN) was held September 5-6, 2015 in Kyoto, Japan. This manuscript presents a summary of the progress in experiments and modeling among ECN contributors leading to a better understanding of soot formation under the ECN “Spray A” configuration and some parametric variants. Relevant published and unpublished work from prior ECN workshops is reviewed. Experiments measuring soot particle size and morphology, soot volume fraction (fv), and transient soot mass have been conducted at various international institutions providing target data for improvements to computational models. Multiple modeling contributions using both the Reynolds Averaged Navier-Stokes (RANS) Equations approach and the Large-Eddy Simulation (LES) approach have been submitted. Among these, various chemical mechanisms, soot models, and turbulence-chemistry interaction (TCI) methodologies have been considered.

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.

This paper exploits the possibilities of achieving an efficient performance optimization methodology to be applied to different spark ignition engine configurations. The objective of the task described here is to determine the combination of parameters which provides the highest volumetric efficiency and effective torque. The definition of general strategy requires first the identification and grouping of the geometric and operating variables to be optimized (duct diameters and lengths, valve timing, spark advance, etc…). The high number of possibilities entails critical choices to reduce, from an engineering design point of view before than from a mathematical point of view, the required computational time. Once proper thermo-fluid dynamic decisions are taken, the most efficient optimization methodology is required. The application of Design of Experiments techniques allows to screen the design space and give a first estimation of the optimal point.

Multiple injections and high EGR rates are now widely adopted for combustion and emissions control in passenger car diesel engines. In a wide range of operating conditions, fuel is provided through one to five separated injection events, and recirculated gas fractions between 0 to 30% are used. Within this context, fast and reliable multi-dimensional models are necessary to define suitable injection strategies for different operating points and reduce both the costs and time required for engine design and development. In this work, the authors have applied a modified version of the characteristic time-scale combustion model (CTC) to predict combustion and pollutant emissions in diesel engines using advanced injection strategies. The Shell auto-ignition model is used to predict auto-ignition, with a suitable set of coefficients that were tuned for diesel fuel.

The definition of a robust methodology to perform a full-cycle CFD simulation of IC engines requires as first step the availability of a reliable grid generation tool, which does not only have to guarantee a high quality mesh but also has to prove to be efficient in terms of required time. In this work the authors discuss a novel approach entirely based on the OpenFOAM technology, in which the available 3D grid generator was employed to automatically create meshes containing hexahedra and split-hexahedra from triangulated surface geometries in Stereolithography (STL) format. The possibility to introduce local refinements and boundary layers makes this tool suitable for IC engine simulations. Grids are sequentially generated at target crank angles which are automatically determined depending on user specified settings such as maximum mesh validity interval and quality parameters like non-orthogonality, skewness and aspect ratio.

Prediction of in-cylinder flows and fuel-air mixing are two fundamental pre-requisites for a successful simulation of direct-injection engines. Over the years, many efforts were carried out in order to improve available turbulence and spray models. However, enhancements in physical modeling can be drastically affected by how the mesh is structured. Grid quality can negatively influence the prediction of organized charge motion structures, turbulence generation and interaction between in-cylinder flows and injected sprays. This is even more relevant for modern direct injection engines, where multiple injections and control of charge motions are employed in a large portion of the operating map. Currently, two different approaches for mesh generation exist: manual and automatic. The first makes generally possible to generate high-quality meshes but, at the same time, it is very time consuming and not completely free from user errors.

Computational fluid dynamics (CFD) codes today represent consolidated tools that cover most physical and chemical processes which occur during operation of internal combustion engines under steady and unsteady conditions. Despite the availability of advanced physical models, the most demanding prerequisite for a CFD engine code is its flexibility in mesh structure and geometry handling capacity to accommodate moving boundaries. In fact, while the motion is solely defined on boundary points, most CFD approaches a-priori specify the position of every mesh vertex in the mesh for every time-step. Alternatives exist, most commonly using mesh generation techniques, like smoothing. In practice this is quite limiting, as it becomes difficult to prescribe solution-dependent motion or perform mesh motion on dynamically adapting meshes. To preserve the mesh quality during extreme boundary deformation due to piston and valve motion, the number of the cells in the mesh needs to be changed.