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 quasi-3D approaches are included in many commercial and research 1D numerical codes, in order to increase their simulation accuracy in presence of complex shape 3D volumes, e.g. plenums and silencers. In particular, these are regarded as valuable approaches for application during the design phase of an engine, for their capability of predicting non-planar waves motion and, on the other hand, for their low requirements in terms of computational runtime. However, the generation of a high-quality quasi-3D computational grid is not always straightforward, especially in case of complex elements, and can be a time-consuming operation, making the quasi-3D tool a less attractive option. In this work, a quasi-3D module has been implemented on the basis of the open-source CFD code OpenFOAM and coupled with the 1D code GASDYN.

The paper investigates the interaction between soil and tractor tires through a 2D numerical model. The tire is schematized as a rigid ring presenting a series of rigid tread bars on the external circumference. The outer profile of the tire is divided into a series of elements, each one able to exchange a normal and a tangential contact force with the ground. A 2D soil model was developed to compute the forces at the ground-tire interface: the normal force is determined on the basis of the compression of the soil generated by the sinking of the tire. The soil is modeled through a layer of springs characterized by two different stiffness for the loading (lower stiffness) and unloading (higher stiffness) condition. This scheme allows to introduce a memory effect on the soil which results stiffer and keeps a residual sinking after the passage of the tire. The normal contact force determines the maximum value of tangential force provided before the soil fails.

The use of virtual prototyping early in the design stage of a product has gained popularity due to reduced cost and time to market. The state of the art in vehicle simulation has reached a level where full vehicles are analyzed through simulation but major difficulties continue to be present in interfacing the vehicle model with accurate powertrain models and in developing adequate formulations for the contact between tire and terrain (specifically, scenarios such as tire sliding on ice and rolling on sand or other very deformable surfaces). The proposed work focuses on developing a ground vehicle simulation capability by combining several third party packages for vehicle simulation, tire simulation, and powertrain simulation. The long-term goal of this project consists in promoting the Digital Car idea through the development of a reliable and robust simulation capability that will enhance the understanding and control of off-road vehicle performance.

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.

In this work, the effects of ozone, hydrogen, carbon monoxide, and exhaust gas recirculation (EGR) addition to Haltermann gasoline combustion were investigated. For these additives, laminar and turbulent flame speeds were experimentally determined using spherically propagating premixed flames in a constant volume combustion vessel. Two initial mixture pressures of Po = 1 and 5 bar, two initial mixture temperatures of 358 and 373 K and a range of equivalence ratios (Ф) from 0.5 to 1 were investigated. The additives were added as single, binary and ternary mixtures to Haltermann gasoline over a wide range of concentrations. For the stoichiometric mixture, the addition of 10% H2, 5% CO and 1000 ppm O3 shows remarkable enhancement (80%) in SL0compared to neat Haltermann gasoline. In addition, for this same blend, increasing the mixture initial temperature and pressure results in a significant increase in SL0compared to the neat gasoline.

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.

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.

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.

Striving for sustainable transportation solutions, hydrogen is often identified as a promising energy carrier and internal combustion engines are seen as a cost effective consumer of hydrogen to facilitate the development of a large-scale hydrogen infrastructure. Driven by efficiency and emissions targets defined by the U.S. Department of Energy, a research team at Argonne National Laboratory has worked on optimizing a spark-ignited direct injection engine for hydrogen. Using direct injection improves volumetric efficiency and provides the opportunity to properly stratify the fuel-air mixture in-cylinder. Collaborative 3D-CFD and experimental efforts have focused on optimizing the mixture stratification and have demonstrated the potential for high engine efficiency with low NOx emissions. Performance of the hydrogen engine is evaluated in this paper over a speed range from 1000 to 3000 RPM and a load range from 1.7 to 14.3 bar BMEP.

Pursuing a sustainable energy scenario for transportation requires the blending of renewable oxygenated fuels such as alcohols into commercial hydrocarbon fuels. From a chemical kinetic perspective, this requires the accurate description of both hydrocarbon reference fuels (n-heptane, iso-octane, toluene, etc.) and oxygenated fuels chemistry. A recent systematic investigation of linear C2-C5 alcohols ignition in a rapid compression machine at p = 10-30 bar and T = 650- 900 K has extended the scarcity of fundamental data at such conditions, allowing for a revision of the low temperature chemistry for alcohol fuels in the POLIMI mechanism. Heavier alcohols such as n-butanol and n-pentanol present ignition characteristic of interest for application in HCCI engines, due to the presence of the hydroxyl moiety reducing their low temperature reactivity compared to the parent linear alkanes (i.e. higher octane number).

The dynamics and evolution of turbulent structures inside an engine-like geometry are investigated by means of Large Eddy Simulation. A simplified configuration consisting of a flat-top cylinder head with a fixed, axis-centered valve and low-speed piston has been simulated by the finite volume CFD code OpenFOAM®; the standard version of the software has been extended to include the compressible WALE subgrid-scale model, models for the generation of synthetic turbulence, some improvements to the mesh motion strategy and algorithms for LES data post-processing. In order to study both the initial transient and the quasi- steady operating conditions, ten complete engine cycles have been simulated. Phase and spatial averages have been performed over cycles three to ten in order to extract first and second moment of velocity; these quantities have then been used to validate the numerical procedure by comparison against experimental data.

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.

Highway and roadway surface measurement is a practice that has been ongoing for decades now. This sort of measurement is intended to ensure a safe level of road perturbances. The measurement may be conducted by a slow moving apparatus directly measuring the elevation of the road, at varying distance intervals, to obtain a road profile, with varying degrees of resolution. An alternate means is to measure the surface roughness at highway speeds using accelerometers coupled with high speed distance measurements, such as laser sensors. Vehicles out rigged with such a system are termed inertial profilers. This type of inertial measurement provides a sort of filtered roadway profile. Much research has been conducted on the analysis of highway roughness, and the associated metrics involved. In many instances, it is desirable to maintain an off-road course such that the course will provide sufficient challenges to a vehicle during durability testing.

The development of thermally durable zeolite NH3/Urea-SCR formulations coupled with that of high porosity filters substrates has opened the way to integrate PM and NOx control into a single device, namely an SCR-coated Diesel Particulate Filter (SCRF). A few experimental works are already present in the literature regarding SCRF systems, mainly addressing the DeNOx performances of the system (in both presence and absence of soot) under both steady state and transient conditions. The purpose of the present work is to perform a simulation study focused on phenomena which are expected to play key roles in SCRF systems, such as coupling of reaction and diffusion phenomena, soot effect on DeNOx activity, SCR coating effect on soot regeneration and filtration efficiency and competition between soot oxidation and DeNOx processes involving NO2.

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

There are many challenges in implementing grid aware plug-in electric vehicle (PEV) charging systems with local load control. New opportunities for innovative load control were created as a result of changes to the 2014 National Electric Code (NEC) about automatic load control definitions for EV charging infrastructure. Stakeholders in optimized dispatch of EV charging assets include the end users (EV drivers), site owner/operators, facility managers and utilities. NEC definition changes allow for ‘over subscription’ of more potential EV charging station load than can be continuously supported if the total load at any time is within the supply system safety limit. Local load control can be implemented via compact submeter(s) with locally hosted control algorithms using direct communication to the managed electric vehicle supply equipment (EVSE).

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.

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.