Search Results

The evaluation of vehicles real emissions circulating in urban areas is a basic activity for planning and management of implemented traffic measures aiming at emission control and air quality improvement. National, region, and city emission inventories require overall average emission estimation based on modeling technique with a few input parameters such as fleet composition and mission profile, represented by average speed. But in the field of emission modeling an important open issue is the very expensive costs of experimental campaigns needed to obtain driving cycle statistically representative of driving behavior, also if only in a specific link of a network. A possible approach to deal with this problem is represented by the use of traffic microscopic simulation models which are capable to simulate individual car motion on the basis of traffic conditions, road characteristics and management rules.

Modern VVA systems offer new potentialities in improving the fuel consumption for spark-ignition engines at low and medium load, meanwhile they grant a higher volumetric efficiency and performance at high load. Recently introduced systems enhance this concept through the possibility of concurrently modifying the intake valve opening, closing and lift leading to the development of almost "throttle-less" engines. However, at very low loads, the control of the air-flow motion and the turbulence intensity inside the cylinder may require to select a proper combination of the butterfly throttling and the intake valve control, to get the highest BSFC (Brake Specific Fuel Consumption) reduction. Moreover, a low throttling, while improving the fuel consumption, may also produce an increased gas-dynamic noise at the intake mouth. In highly "downsized" engines, the intake valve control is also linked to the turbocharger operating point, which may be changed by acting on the waste-gate valve.

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

During gasoline direct injection (GDI) in spark ignition engines, droplets may hit piston or liner surfaces and be rebounded or deposit in the liquid phase as wallfilm. This may determine slower secondary atomization and local enrichments of the mixture, hence be the reason of increased unburned hydrocarbons and particulate matter emissions at the exhaust. Complex phenomena indeed characterize the in-cylinder turbulent multi-phase system, where heat transfer involves the gaseous mixture (made of air and gasoline vapor), the liquid phase (droplets not yet evaporated and wallfilm) and the solid walls. A reliable 3D CFD modelling of the in-cylinder processes, therefore, necessarily requires also the correct simulation of the cooling effect due to the subtraction of the latent heat of vaporization of gasoline needed for secondary evaporation in the zone where droplets hit the wall. The related conductive heat transfer within the solid is to be taken into account.

One-dimensional (1D) models are commonly employed to study the performances of turbocharged engine. Manufacturers' provided steady turbomachinery maps are usually utilized, although they operate in unsteady conditions as a consequence of pressure pulses propagating into the intake and exhaust systems. This may lead to some inaccuracies in the engine-turbocharger matching calculations, which may be solved through the introduction of proper time-delays (virtual pipe corrections). These drawbacks, however, became more relevant when engine operates under low speed and high load conditions, or during a transient maneuver, because of possibilities of compressor surging.

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.

The present paper deals with a comprehensive analysis of the knocking phenomenon through experiments and numerical simulations. Conventional and non-conventional measurements are performed on a 4-stroke, 4-cylinder, turbocharged GDI engine. The engine exhibits optical accesses to the combustion chamber. Imaging in the UV-visible range is carried out by means of a high spatial and temporal resolution camera through an endoscopic system and a transparent window in the piston head. This last is modified to allow the view of the whole combustion chamber almost until the cylinder walls, to include the so-called eng-gas zones. Optical data are correlated to in-cylinder pressure-based indicated analyses in a cycle resolved approach.

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.

The EU emission standards for new rail Diesel engines are becoming even more stringent. EGR and SCR technologies can both be used to reduce NOx emissions; however, the use of EGR is usually accompanied by an increase in PM emissions and may require a DPF. On the other hand, the use of SCR requires on-board storage of urea. Thus, it is necessary to study these trade-offs in order to understand how these technologies can best be used in rail applications to meet new emission standards. The present study assesses the application of these technologies in Diesel railcars on a quantitative basis using one and three dimensional numerical simulation tools. In particular, the study considers a 560 kW railcar engine with the use of either EGR or SCR based solutions for NOx reduction. The NOx and PM emissions performances are evaluated over the C1 homologation cycle.

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.

Gasoline direct injection (GDI) engines are characterized by complex phenomena involving spray dynamics and possible spray-wall interaction. Control of mixture formation is indeed fundamental to achieve the desired equivalence ratio of the mixture, especially at the spark plug location at the time of ignition. Droplet impact on the piston or liner surfaces has also to be considered, as this may lead to gasoline accumulation in the liquid form as wallfilm. Wallfilms more slowly evaporate than free droplets, thus leading to local enrichment of the charge, hence to a route to diffusive flames, increased unburned hydrocarbons formation and particulate matter emissions at the exhaust. Local heat transfer at the wall obviously changes if a wallfilm is present, and the subtraction of the latent heat of vaporization necessary for secondary phase change is also an issue deserving a special attention.

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.

Ignition and flame inception are well recognised as affecting performance and stable operation of spark ignition engines. The very early stage of combustion is indeed the main source of cycle-to-cycle variability, in particular in gasoline direct injection (GDI) engines, where mixture formation may lead to non-homogenous air-to-fuel distributions, especially under some speed and load conditions. From a numerical perspective, 3D modelling of combustion within Reynolds Averaged Navier Stokes (RANS) approaches is not sufficient to provide reliable information about cyclic variability, unless proper changes in the initial conditions of the flow transport equations are considered. Combustion models based on the flamelet concept prove being particularly suitable for the simulation of the energy conversion process in internal combustion engines, due to their low computational cost. These models include a transport equation for the flame surface density, which needs proper initialization.