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.

The objective of this research is a detailed investigation of multiple injections in a highly-dilute diesel low temperature combustion (LTC) regime. This research concentrates on understanding the performance and emissions benefits of multiple injections via experiments and simulations in a 0.48L signal cylinder light-duty engine operating at 2000 r/min and 5.5 bar IMEP. Controlled experiments in the single-cylinder engine are then combined with three computational tools, namely heat release analysis of measured cylinder pressure, a phenomenological spray model using in-cylinder thermodynamics [1], and KIVA-3V Chemkin CFD computations recently tested at LTC conditions [2]. This study examines the effects of fuel split distribution, injection event timing, rail pressure, and boost pressure which are each explored within a defined operation range in LTC.

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 investigates flow and combustion in a full-cycle simulation of a four-stroke, three-valve HCCI engine by visualizing the flow with pathlines. Pathlines trace massless particles in a transient flow field. In addition to visualization, pathlines are used here to trace the history, or evolution, of flow fields and species. In this study evolution is followed from the intake port through combustion. Pathline analysis follows packets of intake charge in time and space from induction through combustion. The local scalar fields traversed by the individual packets in terms of velocity magnitude, turbulence, species concentration and temperatures are extracted from the simulation results. The results show how the intake event establishes local chemical and thermal environments in-cylinder and how the species respond (chemically react) to the local field.

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.

The development of the Diesel Exhaust Filtration Analysis system (DEFA), which utilizes a rectangular wafer of the same substrate material as used in a full-scale Diesel Particulate Filter (DPF), is presented in this paper. Washcoat variations of the wafer substrate (bare, washcoat, and catalyzed washcoat) were available for testing. With this setup, the complications of flow and temperature distribution that arise in the full-scale DPF can be significantly minimized while critical parameters that affect the filtration performance can be fully controlled. Cold flow experiments were performed to test the system's reliability, and determine the permeability of each wafer type. A Computational Fluid Dynamics (CFD) package was utilized to ensure the flow uniformity inside the filter holder during the cold flow test.

An investigation of high speed direct injection (DI) compression ignition (CI) engine combustion fueled with gasoline injected using a triple-pulse strategy in the low temperature combustion (LTC) regime is presented. This work aims to extend the operation ranges for a light-duty diesel engine, operating on gasoline, that have been identified in previous work via extended controllability of the injection process. The single-cylinder engine (SCE) was operated at full load (16 bar IMEP, 2500 rev/min) and computational simulations of the in-cylinder processes were performed using a multi-dimensional CFD code, KIVA-ERC-Chemkin, that features improved sub-models and the Chemkin library. The oxidation chemistry of the fuel was calculated using a reduced mechanism for primary reference fuel combustion chosen to match ignition characteristics of the gasoline fuel used for the SCE experiments.

The paper describes recent advances in the research work concerning the 1-d fluid dynamic modeling of unsteady flows in i.c. engine pipe systems. A comprehensive simulation model has been developed, which is based on different numerical techniques for the solution of the fundamental conservation equations. Classical (MacCormack method plus TVD algorithm) and innovative (the CE-SE method, the discontinuous Galerkin FEM) shock-capturing schemes have been compared, considering the shock-tube problem and the shock-turbulence interaction problem. Moreover, the tracking of the chemical species along the intake and exhaust duct systems has been investigated, introducing the species continuity equations in the numerical model. The engine test case reported in the paper points out the predicted transport of chemical species in the ducts.

There is a significant global effort to study low temperature combustion (LTC) as a tool to achieve stringent emission standards with future light duty diesel engines. LTC utilizes high levels of dilution (i.e., EGR > 60% with <10%O2 in the intake charge) to reduce overall combustion temperatures and to lengthen ignition delay, This increased ignition delay provides time for fuel evaporation and reduces in-homogeneities in the reactant mixture, thus reducing NOx formation from local temperature spikes and soot formation from locally rich mixtures. However, as dilution is increased to the limits, HC and CO can significantly increase. Recent research suggests that CO emissions during LTC result from the incomplete combustion of under-mixed fuel and charge gas occurring after the premixed burn period [1, 2]1. The objective of the present work was to increase understanding of the HC/CO emission mechanisms in LTC at part-load.

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.

Multi-zone CFD simulations with detailed kinetics were used to model iso-octane HCCI experiments performed on a single-cylinder research engine. The modeling goals were to validate the method (multi-zone combustion modeling) and the reaction mechanism (LLNL 857 species iso-octane) by comparing model results to detailed exhaust speciation data, which was obtained with gas chromatography. The model is compared to experiments run at 1200 RPM and 1.35 bar boost pressure over an equivalence ratio range from 0.08 to 0.28. Fuel was introduced far upstream to ensure fuel and air homogeneity prior to entering the 13.8:1 compression ratio, shallow-bowl combustion chamber of this 4-stroke engine. The CFD grid incorporated a very detailed representation of the crevices, including the top-land ring crevice and head-gasket crevice. The ring crevice is resolved all the way into the ring pocket volume. The detailed grid was required to capture regions where emission species are formed and retained.

Predictive models of soot formation during Diesel combustion are of great practical interest, particularly in light of newly proposed strict regulations on particulate emissions. A modified version of the phenomenological model of soot formation developed previously has been implemented in KIVA-II CFD code. The model includes major generic processes involved in soot formation during combustion, i.e., formation of soot precursors, formation of surface growth species, soot particle nucleation, coagulation, surface growth and oxidation. The formulation of the model within the KIVA-II is fully coupled with the mass and energy balances in the system. The model performance has been tested by comparison with the results of optical in-cylinder soot measurements in a single cylinder Cummins NH Diesel engine. The predicted soot volume fraction, number density and particle size agree reasonably well with the experimental data.

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.

Increasing demands on the capabilities of engine simulation and the ability to accurately predict both performance and acoustics has lead to the development of multiple approaches, ranging from fully 3D to simplified 1D models. In this work it will be described the development and application of hybrid 1D-3D approaches and an innovative one based on the 3D cell element. This is designed to model the acoustics of intake and exhaust system components used in internal combustion engines. Models of components are built using a network or grid of 3D cells based primarily on the geometry of the system. This means that these models can be built without fundamental knowledge of acoustically equivalent systems making their range of application larger as well as making them simpler to construct. Due to the 3D nature of these models it is also possible to predict higher order modes and improve the accuracy of models at high frequencies compared to conventional plane wave approaches.

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 objective of this study is to increase fundamental understanding of the effects of fuel composition and properties on low temperature combustion (LTC) and to identify major properties that could enable engine performance and emission improvements, especially under high load conditions. A series of experiments and computational simulations were conducted under LTC conditions using 67% EGR with 9.5% inlet O₂ concentration on a single-cylinder version of the General Motors Corporation 1.9L direct injection diesel engine. This research investigated the effects of Cetane number (CN), volatility and total aromatic content of diesel fuels on LTC operation. The values of CN, volatility, and total aromatic content studied were selected in a DOE (Design of Experiments) fashion with each variable having a base value as well as a lower and higher level. Timing sweeps were performed for all fuels at a lower load condition of 5.5 bar net IMEP at 2000 rpm using a single-pulse injection strategy.

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.

Multi-dimensional computational efforts using comprehensive and skeletal kinetics have been made to investigate the cool flame region in HCCI combustion. The work was done in parallel to an experimental study that showed the impact of the negative temperature coefficient and the cool flame on the start of combustion using different fuels, which is now the focus of the simulation work. Experiments in a single cylinder CFR research engine with n-butane and a primary reference fuel with an octane number of 70 (PRF 70) were modeled. A comparison of the pressure and heat release traces of the experimental and computational results shows the difficulties in predicting the heat release in the cool flame region. The behavior of the driving radicals for two-stage ignition is studied and is compared to the behavior for a single-ignition from the literature. Model results show that PRF 70 exhibits more pronounced cool flame heat release than n-butane.

A modeling study has been conducted in order to characterize the heat transfer in an automotive diesel exhaust system. The exhaust system model, focusing on 2 exhaust pipes, has been created using a transient 1-D engine flow network simulation program. Model results are in excellent agreement with experimental data gathered before commencement of the modeling study. Predicted pipe exit stream temperatures are generally within one percent of experimental values. Sensitivity analysis of the model was the major focus of this study. Four separate variables were chosen for the sensitivity analysis. These being the external convective heat transfer coefficient, external emissivity, mass flow rate of exhaust gases, and amplitude of incoming pressure fluctuations. These variables were independently studied to determine their contribution to changes in exhaust gas stream temperature and system heat flux. There are two primary benefits obtained from conducting this analysis.