Search Results

Ammonia/urea-SCR is a mature technology, applied worldwide for the control of NOx emissions in combustion exhausts from thermal power plants, cogeneration units, incinerators and stationary diesel engines and more recently also from mobile sources. However a greater DeNOx activity at low temperatures is desired in order to meet more and more restrictive legislations. In this paper we report transient and steady state data collected over commercial Fe-ZSM-5 and V₂O₅-WO₃/TiO₂ catalysts showing high NOx reduction efficiencies in the 200 - 350°C T-range when NO and ammonia react with nitrates, e.g., in the form of an aqueous solution of ammonium nitrate. Under such conditions a new reaction occurs, the so-called "Enhanced SCR" reaction, 2 NH₃ + 2 NO + NH₄NO₃ → 3 N₂ + 5 H₂O.

In this work a multilevel CFD analysis have been applied for the design of an intake air-box with improved characteristics of noise reduction and fluid dynamic response. The approaches developed and applied for the optimization process range from the 1D to fully 3D CFD simulation, exploring hybrid approaches based on the integration of a 1D model with quasi-3D and 3D tools. In particular, the quasi-3D strategy is exploited to investigate several configurations, tailoring the best trade-off between noise abatement at frequencies below 1000 Hz and optimization of engine performances. Once the best configuration has been defined, the 1D-3D approach has been adopted to confirm the prediction carried out by means of the simplified approach, studying also the impact of the new configuration on the engine performances.

Swirling flows are very dominant in applied technical problems, especially in IC engines, and their prediction requires rather sophisticated modeling. An adaptive low-pass filtering procedure for the modeled turbulent length and time scales is derived and applied to Menter' original k - ω SST turbulence model. The modeled length and time scales are compared to what can potentially be resolved by the computational grid and time step. If the modeled scales are larger than the resolvable scales, the resolvable scales will replace the modeled scales in the formulation of the eddy viscosity; therefore, the filtering technique helps the turbulence model to adapt in accordance with the mesh resolution and the scales to capture.

In the next years, the upcoming emission legislations are expected to introduce further restrictions on the admittable level of pollutants from vehicles measured on homologation cycles and real drive tests. In this context, the strict control of pollutant emissions at the cold start will become a crucial point to comply with the new regulation standards. This will necessarily require the implementation of novel strategies to speed-up the light-off of the reactions occurring in the after-treatment system, since the cold start conditions are the most critical one for cumulative emissions. Among the different possible technological solutions, this paper focuses on the evaluation of the potential of a burner system, which is activated before the engine start. The hypothetical burner exploits the lean combustion of an air-gasoline mixture to generate a high temperature gas stream which is directed to the catalyst section promoting a fast heating of the substrate.

Methane abatement in the exhaust gas of natural gas engines is much more challenging in respect to the oxidation of other higher order hydrocarbons. Under steady state λ sweep, the methane conversion efficiency is high at exact stoichiometric, and decreases steeply under both slightly rich and slightly lean conditions. Transient lean to rich transitions can improve methane conversion at the rich side. Previous experimental work has attributed the enhanced methane conversion to activation of methane steam reforming. The steam reforming rate, however, attenuates over time and the methane conversion rate gradually converges to the low steady state values. In this work, a reactor model is established to predict steady state and transient transition characteristics of a three-way catalyst (TWC) mounted in the exhaust of a natural gas heavy-duty engine.

In the last decades numerical simulations have become reliable tools for the design and the optimization of silencers for internal combustion engines. Different approaches, ranging from simple 1D models to detailed 3D models, are nowadays commonly applied in the engine development process, with the aim to predict the acoustic behavior of intake and exhaust systems. However, the acoustic analysis is usually performed under the hypothesis of infinite stiffness of the silencer walls. This assumption, which can be regarded as reasonable for most of the applications, can lose validity if low wall thickness are considered. This consideration is even more significant if the recent trends in the automotive industry are taken into account: in fact, the increasing attention to the weight of the vehicle has lead to a general reduction of the thickness of the metal sheets, due also to the adoption of high-strength steels, making the vibration of the components a non negligible issue.

In this work an integration between a 1D code (Gasdyn) with a CFD code (OpenFOAM®) has been applied to improve the performance of a Moto3TM engine. The four-stroke, single cylinder S.I. engine was modeled, in order to predict the wave motion in the intake and exhaust systems and study how it affects the cylinder gas exchange process. The engine considered was characterized by having an air induction system with integrated filter cartridge, air-box and intake runner, resulting in a complex air-path form the intake mouth to the intake valves, which presents critical aspects when a 1D modeling is addressed. This paper presents a combined and integrated simulation, in which the intake systems was modeled as a 3D geometry whereas the exhaust system, which presented a simpler geometry, was modeled by means of a 1D approach.

Selective Catalytic Reduction (SCR) systems are nowadays widely applied for the reduction of NOx emitted from Diesel engines. The typical process is based on the injection of aqueous urea in the exhaust gases before the SCR catalyst, which determines the production of the ammonia needed for the catalytic reduction of NOx. However, this technology is affected by two main limitations: a) the evaporation of the urea water solution (UWS) requires a sufficiently high temperature of the exhaust gases and b) the formation of solid deposits during the UWS evaporation is a frequent phenomenon which compromise the correct operation of the system. In this context, to overcome these issues, a technology based on the injection of gaseous ammonia has been recently proposed: in this case, ammonia is stored at the solid state in a cartridge containing a Strontium Chloride salt and it is desorbed by means of electrical heating.

Transient and steady-state kinetic data are herein presented to analyze the inhibiting effect of ammonia on the NH₃-SCR of NO at low temperatures over a Fe-zeolite commercial catalyst for vehicles. It is shown that in SCR converter models a rate expression accounting for NH₃ inhibition of the Standard SCR reaction is needed in order to predict the specific dynamics observed both in lab-scale and in engine test bench runs upon switching on and off the ammonia feed. Two redox, dual site kinetic models are developed which ascribe such inhibition to the spill-over of ammonia from its adsorption sites, associated with the zeolite, to the redox sites, associated with the Fe promoter. Better agreement both with lab-scale intrinsic kinetic runs and with engine test-bench data, particularly during transients associated with dosing of ammonia to the SCR catalyst, is obtained assuming slow migration of NH₃ between the two sites.

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.

Transport activities contribute significantly to the air pollution and its impact on emissions is a key element in the evaluation of any transport policy or plan. Calculation of emissions has therefore gained institutional importance in the European Community. To obtain emission factors several methods make use of only vehicle mean velocity, which can be easily obtained by vehicle flow and density in the road. Recently in ARTEMIS project by Rapone et al. (2005–2007) a meso scale emission model, named KEM (Kinematic Emission Model), able to calculate emission factor has been developed. This model is based on a new statistical methodology, capable to consider more attributes than the simple mean speed to characterize driving behaviour. An interesting approach to determine the exact mix of driving cycles is represented by the use of microscopic traffic simulation models that could be used to avoid the very expensive costs of experimental campaigns needed to obtain real driving cycle.

Transport activities contribute significantly to air pollution. For this reason any policy or plan, carried out by administration or institution, requires the assessment of its impact on the emissions. To assess the overall pollutant production from transport, it is necessary to calculate emission factors. For this aim several methods exist which only use the average speed of the traffic stream, which can be theoretically obtained by vehicles flow and density on the road. Recently, a new statistical approach has been developed capable to consider more attributes than the simple mean speed to characterize driving behaviour, not only in the determination of driving cycles but also in the emission modelling. In this context, a meso scale emission model, named KEM, Kinematic Emission Model, able to calculate emission factor was developed. However, it is necessary to consider that the input to this model is, in any case, the driving cycle.

In the last years, as a response to the more and more restrictive emission legislation, new devices (SRC, DOC, NOx-trap, DPF) have been progressively introduced as standard components of modern after-treatment system for Diesel engines. In addition, the adoption of electrical heating is nowadays regarded with interest as an effective solution to promote the light-off of the catalyst at low temperature, especially at the start-up of the engine and during the low load operation of the engine typical of the urban drive. In this work, a state-of-the-art 48 V electrical heated catalyst is considered, in order to investigate its effect in increasing the abatement efficiency of a standard DOC. The electrical heating device considered is based on a metallic support, arranged in a spiral layout, and it is heated by the Joule effect due to the passage of the electrical current.

In this paper the low temperature reduction process of nitrates stored at high temperatures over model Pt-Ba/Al2O3 LNT catalysts using both H2 and C3H6 is analyzed. The results indicate that over the Pt-Ba/Al2O3 catalyst the reduction of stored NOx with both H2 and C3H6 occurs at temperature below those corresponding to their thermal stability. Accordingly, the reduction process occurs through a Pt-catalyzed surface reaction, which does not involve, as a preliminary step, the thermal decomposition of the adsorbed NOx species. The occurrence of such a pathway also requires the co-presence of the storage element and of the noble metal on the same support.

In-cylinder measurements of turbulent integral length scales, carried out during the last 60 degrees of the compression stroke at 600 and 1,000 rpm by a two-probe volume LDV system, were used to assess the capability of the k-ε model used in KIVA-II code. The objective of the paper is to address the following question: what is the most reasonable definition of turbulent length scale in the k-ε model for engine applications? The answer derived from the comparison between KIVA predictions and experiments that showed a fair agreement between the computed turbulent length scale and the measured lateral integral length scale. The agreement is a result of proper choice of the initial swirl ratio and turbulent kinetic energy at inlet valve closure (IVC) by taking into account the LDV measurements and the value of the constant Cμε in the k-ε model equations that relates the turbulent length scale to k and ε.

This work presents a comprehensive numerical model for ice accretion and Ice Protection System (IPS) simulation over a 2D component, such as an airfoil. The model is based on the Myers model for ice accretion and extended to include the possibility of a heated substratum. Six different icing conditions that can occur during in-flight ice accretion with an Electro-Thermal Ice Protection System (ETIPS) activated are identified. Each condition presents one or more layers with a different water phase. Depending on the heat fluxes, there could be only liquid water, ice, or a combination of both on the substratum. The possible layers are the ice layer on the substratum, the running liquid film over ice or substratum, and the static liquid film between ice and substratum caused by ice melting. The last layer, which is always present, is the substratum. The physical model that describes the evolution of these layers is based on the Stefan problem. For each layer, one heat equation is solved.

This paper presents a novel fully-automatic remeshing procedure, based on the level-set method and Delaunay triangulation, to model three-dimensional boundary problems and generate a new conformal body-fitted mesh. The proposed methodology is applied to long-term in-flight ice accretion, which is characterized by the formation of extremely irregular ice shapes. Since ice accretion is coupled with the aerodynamic flow field, a multi-step procedure is implemented. The total icing exposure time is subdivided into smaller time steps, and at each time step a three-dimensional body-fitted mesh, suitable for the computation of the aerodynamic flow field around the updated geometry, is generated automatically. The methodology proposed can effectively deal with front intersections, as shown with a manufactured example.

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 introduces and discusses a new 2D physical model which has been developed and validated in order to study and optimize the handling behavior of the tire. It can be divided into two parts, the structural model and the contact area model. The parameters, that are function of the vertical load, are identified or calculated by comparison with the results provided by 3D finite element models. The input data for the identification procedure consist of a set of frequency responses performed on the finite element model. A second set of simulations on the 3D model of the tread pattern gives the characteristics of the contact model. Once built the 2D model it is easy to carry out both steady state and transient analysis. The steady state analysis returns the cornering carpet, in terms of lateral force and self-aligning moment as function of the cornering angle. The transient analysis allows the evaluation of the relaxation length and dynamic characteristics.

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.