Search Results

The preparation of the air-fuel mixture represents one of the most critical tasks in the definition of a clean and efficient SI engine. Therefore it becomes necessary to consolidate the numerical methods which are able to describe such a complex physical process. Within this context, the authors developed a CFD methodology into an open-source code to investigate the air-fuel mixture formation process in PFI engines. Attention is focused on moving mesh algorithms, Lagrangian spray modeling and spray-wall interaction modeling. Since moving grids are involved and the mesh quality during motion strongly influences the computed in-cylinder flow-field, a FEM-based automatic mesh motion solver combined with topological changes was adopted to preserve the grid quality in presence of high boundary deformations like the interaction between the piston bowl and the valves during the valve-overlap period.

A numerical investigation on a series of Diesel spray experiments in constant-volume vessels is proposed. Non reacting conditions were used to assess the spray models and to determine the grid size required to correctly predict the fuel-air mixture formation process. To this end, not only computed liquid and vapor penetrations were compared with experimental data, but also a detailed comparison between computed and experimental mixture fraction distributions was performed at different distances from the injector. Grid dependency was reduced by introducing an Adaptive Local Mesh Refinement technique (ALMR) with an arbitrary level of refinement. Once the capabilities of the current implemented spray models have been assessed, reacting conditions at different ambient densities and temperatures were considered. A Perfectly Stirred Reactor (PSR) combustion model, based on a direct integration of complex chemistry mechanisms over a homogenous cell, was adopted.

The high complexity of the interplaying physical and chemical phenomena occurring in Diesel combustion has brought an increasing interest versus experimental and computational fundamental studies. A relevant contribution to the availability of suitable and accessible experimental data is represented by the Engine Combustion Network database, which are of great interest for model development and validation because of the well-defined boundary conditions and the wide range of conditions employed [1]. In this paper the authors implemented into a unique open-source CFD code two different combustion models of different complexity and applied them to simulate a selection of significant test cases from this database. The first one is an extension of the Eddy Dissipation Model and it can be used as a development tool because of its limited computational demand.

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.

A 1-D thermo-fluid dynamic simulation code, including a quasi-D combustion model coupled with a detailed kinetic scheme, is used to analyze the combustion process in HCCI engines. The chemical mechanism has previously been validated in comparison with experimental data over a wide range of operating conditions. To explore the impact on model predictions, the cylinder was divided into multiple zones to characterize the conditions of the in-cylinder charge. Particular attention is devoted to the numerical algorithm in order to ensure the robustness and efficiency of the large system solution. This numerical model allows study of the autoignition of the air fuel mixture and determines the chemical evolution of the system. The proposed model was compared with in-cylinder temperature and chemical species profiles. The experimental activity was carried out in the combustion chamber of a single cylinder air cooled engine operating in HCCI mode.

In this work two significant approaches with different level of complexity for Diesel combustion modeling have been analyzed and implemented by the authors in the same open-source code, OpenFOAM: the Eddy Dissipation model and the PaSR (Partially Stirred Reactor) model. Concerning the latter, the potentialities and accuracy of the ISAT (In-Situ Adaptive Tabulation) algorithm, recently implemented by the authors into OpenFOAM, to reduce significantly the required computational time, were evaluated. Similarly, different models to predict pollutant emissions are revised and tested, due to considering the need to fulfill the future emission regulations imposing further reductions in NOx and soot. The mentioned models were applied to simulate a selection of significant Constant-Volume Diesel Combustion test-cases from the Engine Combustion Network database [1].

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.

This work describes the development, application and coupling of two different numerical codes, respectively based on a 1D (Gasdyn) and 3D (OpenFOAM) schematization of the geometrical domain. They have been adopted for the prediction of the wave motion inside the intake and the exhaust systems of internal combustion engines. The HLLC Riemann solver has been implemented both in the CFD and the 1D codes to solve the Euler system of equations, in order to operate with the same solver on the different calculation domains. Moreover, the HLLC solver has been applied to treat the boundary conditions at the interface between the two domains, in such a way to allow the propagation of flow disuniformities through the domain interface, without affecting the solution accuracy. The hybrid approach was used for the simulation of two different test cases: a complex 5 into 1 pipe junction of a high performance V10 engine and a Venturi tube plus a Helmholtz resonator of a single cylinder S.I. engine.

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.

This paper describes the development, application and comparison of two different non-linear numerical codes, respectively based on a 1D and 3D schematization of the geometrical domain, for the prediction of the acoustic behavior of common silencing devices for i.c. engine pulse noise abatement. A white noise approach has been adopted and applied to predict the attenuation curves of silencers in the frequency domain, while a non-reflecting boundary condition was used to represent an anechoic termination. Expansion chambers, Helmholtz and column resonators, Herschel-Quincke tubes have been simulated by both the 1D and the 3D codes and the results compared to the available linear acoustic analytical solutions. Finally, a hybrid approach, in which the CFD code has been integrated with the 1D model, is described and applied to the simulation of a single cylinder engine. The computed results are compared to the measured pressure waves and emitted sound pressure level spectra.

The s.i. combustion process and its corresponding pollutant formation are investigated by means of a quasiD approach and a CFD model. This work has been motivated by the need to better understand the reliability of such models and to assess their accuracies with respect to the prediction of engine performances and emissions. An extended dissertation about the fundamental mechanisms governing the pollutant formation in the turbulent premixed combustion which characterizes the s.i. engines is given. The conclusion of such analysis is the definition of a new reduced chemical scheme, based on the application of partial-equilibrium and steady-state assumptions for the radicals and the solution of a transport equation for each specie which is kinetically controlled. For this purpose the CFD code OpenFOAM [1, 2, 3] and the thermo-fluid dynamic code GASDYN [4, 5] have been applied and enhanced.

Aim of this work is to present and discuss the possibility and the limits of two zone models for spark-ignition engines using a detailed kinetic scheme for the characterization of the evolution of the air-fuel mixture, while an equilibrium approach is used for the burnt zone. Simple experimental measurements of knocking tendency of different fuels in ideal reactors, such as rapid compression machines and shock tube reactors, cannot be directly used for the analysis of octane numbers and sensitivity of hydrocarbon mixtures. Thus a careful investigation is very useful, not only of the combustion chamber behavior, including the modelling of the turbulent flame front propagation, but also of the fluid dynamic behavior of the intake and exhaust system, accounting for the volumetric efficiency of the engine.

The paper describes the results of a parallel 1D thermo-fluid dynamic simulation and experimental investigation of a DI turbocharged Diesel engine. The attention has been focused on the overall engine performances (air flow, torque, power, fuel consumption) as well as on the emissions (NO and particulate) along the after-treatment system, which presents a particulate filter. The 1D research code GASDYN for the simulation of the whole engine system has been enhanced by the introduction of a multi-zone quasi-dimensional combustion model for direct injection Diesel engines. The effect of multiple injections is taken into account (pilot and main injection). The prediction of NO and soot has been carried out respectively by means of a super-extended Zeldovich mechanism and by the Hiroyasu kinetic approach.

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.

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 experimental and simulation work recently carried out to investigate the effects of secondary air injection on the emission conversion in the exhaust after-treatment system of a S.I. automotive engine. The modeling of the 1D unsteady reacting flows in the complete exhaust system of a spark ignition engine, designed to satisfy the Euro IV limits, has been performed including the secondary air injection system, to predict the possible shortening of catalyst light-off time and the speed-up of the after-treatment system warm-up. The transport of chemical species with reactions in gas phase (post-oxidation of unburned HC in the exhaust manifold) and in solid phase (conversion of pollutants in the catalyst) with and without secondary air has been simulated by the 1D thermo-fluid dynamic model GASDYN, developed by the authors.

The paper describes the research work concerning the simulation of 1D unsteady reacting flows in s.i. engine pipe-systems, including pre-catalysts and main catalysts. The numerical model GASDYN has been developed to enable the concurrent prediction of the wave motion in the intake and exhaust ducts, the chemical composition of the gas discharged by the cylinder of a s.i. engine, the chemical and thermal behavior of catalytic converters. The effect of considering the transport of chemical species with reactions in gas phase (post-oxidation of unburned HC in the exhaust manifold) and in solid phase (conversion of pollutants in the catalyst) on the predicted wave motion is reported.

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 calculation of the main pollutant emissions discharged into the atmosphere by means of numerical codes requires the development of integrated models, including either an accurate thermodynamic in-cylinder analysis and the simulation of reacting unsteady flows in the duct system. This paper describes the main features of the numerical model GASDYN developed by the authors, which in the last years has been enhanced in order to achieve this kind of objectives. A multi-zone approach has been adopted to predict the combustion process in s.i. engines, whereas the so called super-extended Zeldovich mechanism has been introduced to perform a more detailed description of all the chemical reactions involved in the NOx production process. The simulation of the reacting flows in the exhaust manifold has been completed by the introduction of further enhancements to predict the chemical behavior of gases inside the catalytic converters.

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.