Search Results

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.

In the last few years, the effect of diabatic test conditions on compressor performance maps has been widely investigated, leading some Authors to propose different correction models. The accuracy of turbocharger performance map constitute the basis for the tuning and validation of a numerical method, usually adopted for the prediction of engine-turbocharger matching. Actually, it is common practice in automotive applications to use simulation codes, which can either require measured compression ratio and efficiency maps as input values or calculate them “on the fly” throughout specific sub-models integrated in the numerical procedures. Therefore, the ability to correct the measured performance maps taking into account internal heat transfer would allow the implementation of commercial simulation codes used for engine-turbocharger matching calculations.

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 flow field resulting from injecting a gas jet into a crossflow confined in a narrow square duct has been studied under steady regime using schlieren imaging and laser Doppler velocimetry (LDV). This transparent duct is intended to simulate the intake port of an internal combustion engine fueled by gaseous mixture, and the jet is issued from a round nozzle. The schlieren images show that the relative small size of the duct would confine the development of the transverse jet, and the interaction among jet and sidewalls strongly influences the mixing process between jet and crossflow. The mean velocity and turbulence fields have been studied in detail through LDV measurements, at both center plane and several cross sections. The well-known flow feature formed by a counter rotating vortex pair (CVP) has been observed, which starts to appear at the jet exit section and persists far downstream contributing to enhancing mixing process.

The most common automotive drivelines transmit the engine torque to the driven axle through a differential. Semi-active versions of this device ([4], [5], [6]) have been recently conceived to improve vehicle handling at limit and under particular conditions; these differentials are based on the structural scheme of the passive one but they try to manipulate the vehicle dynamics by controlling the distribution of the driving torque on the wheels of the same axle thus generating a yaw moment. Unfortunately a semi-active differential is not able to perform a complete yaw control since the torque can only be transferred from the faster wheel to the slower one; on the other hand, active differentials ([11], [12], [13]) allow to generate the most appropriate yaw moment controlling both the amount of transferred torque and its direction.

The potentialities shown by controlled differentials is making the automotive industry to explore this field. While VDC systems can only guarantee a safe behaviour at limit, a controlled differential can also increase the handling performance. The system derives from a RWD driveline with a semi-active differential, to which has been added a controlled wet clutch that directly connects the engine to the front axle. This device allows to distribute the drive torque between the two axles. It can be easily understood that in this device the torque distribution doesn't depend only from the central clutch action, but also from the engaged gear. Because of this particular layout this system can't work in the whole gear because thermal problems due to kinematical reasons. So the central clutch controller has to consider the gear position too.

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.

The scope of the work presented in this paper was to apply the latest open source CFD achievements to design a state of the art, direct-injection (DI), heavy-duty, natural gas-fueled engine. Within this context, an initial steady-state analysis of the in-cylinder flow was performed by simulating three different intake ducts geometries, each one with seven different valve lift values, chosen according to an estabilished methodology proposed by AVL. The discharge coefficient (Cd) and the Tumble Ratio (TR) were calculated in each case, and an optimal intake ports geometry configuration was assessed in terms of a compromise between the desired intensity of tumble in the chamber and the satisfaction of an adequate value of Cd. Subsequently, full-cycle, cold-flow simulations were performed for three different engine operating points, in order to evaluate the in-cylinder development of TR and turbulent kinetic energy (TKE) under transient conditions.

The paper presents results of an experimental investigation of the fluid dynamic processes during the air/fuel mixture formation period between an evaporating diesel spray and swirl air flow under realistic engine conditions. Particle Image Velocimetry (PIV) experiments have been carried out using an optically accessible prototype 2-stroke diesel engine equipped with a swirled combustion chamber. The flow within the chamber assumes a well structured swirl motion, similar to that developing in a real diesel engine, operating at high swirl ratio. The engine has been equipped with a common rail injection system and a solenoid-controlled injector, in use on automotive engines for the European market, able to manage multiple injection strategies. Two injector nozzles have been tested: a micro-sac 5-hole nozzle, 0.13 mm diameter, 150° spray angle and a 7-hole, 0.141 mm diameter, 148° spray angle.

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 current development to set up an automatic procedure for automatic mesh generation and automatic mesh motion for internal combustion engine simulation in OpenFOAM®-2.2.x is here described. In order to automatically generate high-quality meshes of cylinder geometries, some technical issues need to be addressed: 1) automatic mesh generation should be able to control anisotropy and directionality of the grid; 2) during piston and valve motion, cells and faces must be introduced and removed without varying the overall area and volume of the cells, to avoid conservation errors. In particular, interpolation between discrete fields is frequent in computational physics: the use of adaptive and non-conformal meshes necessitates the interpolation of fields between different mesh regions. Interpolation problems also arise in areas such as model coupling, model initialization and visualisation.

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

The air entrainment in a transient diesel spray was studied using laser Doppler anemometry to provide information on the effect of gas density and temperature. The spray was injected vertically into a confined quiescent atmosphere and the entrained mass flow rate was evaluated by measuring the air velocity component normal to a cylindrical geometric surface surrounding the spray, and extending to about 200 nozzle diameters (50 mm). The experimental results, relative to a density range from 0.84 to 7.02 kg/m3 and a temperature range from 293 to 473 K, indicate that the non dimensional entrainment rate, averaged in time over the main injection period, depends on the distance from the nozzle and both gas density and temperature. A first analysis, based on the available data, allowed to quantify the dependence and provided a correlation with such variables.

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.

The effects of fuel temperature on both the geometry and the droplet size and velocity of a GDI swirled injector spray were investigated by means of visualizations and PDA measurements. Isooctane was used as model fuel and was injected in a quiescent bomb at injection pressure of 7 MPa. Bomb pressure ranged from 40 kPa to 800 kPa with injector nozzle temperature ranging from 293 K to 393 K. A drastic change in spray geometry was observed when conditions above the vaporization curve were reached. The temperature increase has two macroscopic effects on the spray geometry: at the nozzle exit the liquid flash boiling strongly enlarges the spray angle, at a certain distance from the nozzle the air entrainment collapses the spray. Raising the fuel temperature up to flash boiling conditions causes a significant decrease of the average droplet size.