Search Results

In this work a detailed model to simulate the transient behavior of catalytic converters is presented. The model is able to predict the unsteady and reacting flows in the exhaust ducts, by solving the system of conservation equations of mass, momentum, energy and transport of reacting chemical species. The en-gine and the intake system have not been included in the simulation, imposing the measured values of mass flow, gas temperature and chemical composition as a boundary condition at the inlet of the exhaust system. A detailed analysis of the diffusion stage triggering is proposed along with simplifications of the physics, finalized to the reduction of the calculation time. Submodels for water condensation and its following evaporation on the monolith surface have been taken into account as well as oxygen storage promoted by ceria oxides.

Nowadays a great attention is paid to the level and quality of noise radiated from the tailpipe end of intake and exhaust systems, to control the gas dynamic noise emitted by the engine as well as the characteristics of the cabin interior sound. The muffler geometry can be optimized consequently, to attenuate or remark certain spectral components of the engine noise, according to the result expected. Evidently the design of complex silencing systems is a time-consuming operation, which must be carried out by means of concurrent experimental measurements and numerical simulations. In particular, 1D and multiD linear/non-linear simulation codes can be applied to predict the silencer behavior in the time and frequency domain. This paper describes the development of a 1D-multiD integrated approach for the simulation of complex muffler configurations such as reverse chambers with inlet and outlet pipe extensions and perforated silencers with the addition of sound absorbing material.

This paper describes the development of a fully 1D and of a 1D-multiD integrated approach for the simulation of complex muffler configurations. The fully 1D approach aims to model the muffler recurring to an equivalent net of 1D pipes. An expansion chamber with offset inlet and outlet pipes was modeled with this preocedure and the resuts compared to CFD simulations, pointing out some critical aspects in the TL prediction. The HLLC Riemann solver and its extension to the second order were implemented both in the 1D and multiD models and exploited to handle the interface between the calculation domains. The integrated 1D-multiD approach was used afterwards to predict the transmission loss of more complex geometries such as series chambers with extended inlet and outlet pipes and with flow reversals. A new procedure has been adopted to calculate the transmission loss, imposing a pressure impulse at the inlet and evaluating the response of the muffler.

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.

Multi-dimensional simulation of hydrodynamics in full-scale wall-flow Diesel Particulate Filters by GpenFQAM®, an open-source C++ object-oriented CFD code, is presented. A new fast and efficient parallel numerical solver has been developed by authors to simulate flows through porous media and it has been tested for the simulation of diesel particulate filters; errors caused by discretization of filter monoliths have been corrected by the formulation of a correction factor, that has been included in the solver. A set of experimental data, available from literature, has been used for code validation.

A new fast and efficient parallel numerical solver for reacting and compressible flows through porous media has been developed in the OpenFOAM® (Open Field Operation and Manipulation) CFD Toolbox. With respect to the macroscopic model for porous media originally available in OpenFOAM®, a different mathematical approach has been followed: the new implemented solver makes use of the physical normal components resulting from the velocity expansion in the unit orthogonal vector basis to compute the Darcy pressure drop across the porous medium. Also, an additional sink term to account for the increased flow friction over the porous wall has been included into the momentum equation. In the new solver, the pressure correction equation is still able to achieve a faster convergency at very low permeability of the medium, also when it is associated with grid non-orthogonality.

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.

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

A new approach for the fluid-dynamic simulation of the Diesel Particulate Filters (DPF) has been developed. A mathematical model has been formulated as a system of nonlinear partial differential equations describing the conservation of mass, momentum and energy for unsteady, compressible and reacting flows, in order to predict the hydrodynamic characteristics of the DPF and to study the soot deposition mechanism. In particular, the mass conservation equations have been solved for each chemical component considered, and the advection of information concerning the chemical composition of the gas has been figured out for each computational mesh. A sub-model for the prediction of the soot cake formation has been developed and predictions of soot deposition profiles have been calculated for different loading conditions. The results of the simulations, namely the calculated pressure drop, have been compared with the experimental data.

This paper describes some recent advances in the field of I.C. engine modeling and simulation, concerning the development and application of a 1D thermo-fluid dynamic research code. An extensive experimental analysis has been concurrently carried out, to support the development and validation of the simulation code. A four-stroke, 10V-cylinder, 5.0 liters automotive S.I. engine has been modeled, in order to predict not only the wave motion in the system and its influence on the cylinder gas exchange process, but also the in-cylinder pressure to get a good prediction of 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 process via a thermodynamic multi-zone model, based on a fractal approach to predict the turbulent combustion.

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.