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.

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.

Nowadays development of automotive HVAC is a challenging task wherein thermal comfort and safety are very critical factors to be met. HVAC system is responsible for the demisting and defrosting of the vehicle’s windshield and for creating/maintaining a pleasing environment inside the cabin by controlling airflow, velocity, temperature and purity of air. Fog or ice which forms on the windshield is the main reason for invisibility and leads to major safety issues to the customers while driving. It has been shown that proper clear visibility for the windshield could be obtained with a better flow pattern and uniform flow distribution in the defrost mode of the HVAC system and defrost duct. Defroster performance has received significant attention from OEMs to meet the specific global performance standards of FMVSS103 and SAE J902. Therefore, defroster performance is seriously taken into consideration during the design of HVAC system and defroster duct.

This paper focuses on detailed numerical simulations of direct injection diesel and gasoline sprays from production grade, multi-hole injectors. In a dual-fuel engine the direct injection of both the fuels can facilitate appropriate mixture preparation prior to ignition and combustion. Diesel and gasoline sprays were simulated using high-fidelity Large Eddy Simulations (LES) with the dynamic structure sub-grid scale model. Numerical predictions of liquid penetration, fuel density distribution as well as transverse integrated mass (TIM) at different axial locations versus time were compared against x-ray radiography data obtained from Argonne National Laboratory. A necessary, but often overlooked, criterion of grid-convergence is ensured by using Adaptive Mesh Refinement (AMR) for both diesel and gasoline. Nine different realizations were performed and the effects of random seeds on spray behavior were investigated.

In recent years, there has been renewed attention focused on open jet correction methods, in particular on the two-measurement method of E. Mercker, K. Cooper, and co-workers. This method accounts for blockage and static pressure gradient effects in automotive wind tunnels and has been shown by both computations and experiments to appropriately adjust drag coefficients towards an on-road condition, thus allowing results from different wind tunnels to be compared on a more equitable basis. However, most wind tunnels have yet to adopt the method as standard practice due to difficulties in practical application. In particular, it is necessary to measure the aerodynamic forces on every vehicle configuration in two different static pressure gradients to capture that portion of the correction. Building on earlier proof-of-concept work, this paper demonstrates a practical method for implementing the two-measurement procedure and demonstrates how it can be used for production testing.

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.

Typical production vehicle development includes road testing of a vehicle towing a trailer to evaluate powertrain thermal performance. In order to correlate tests with simulations, the aerodynamic effects of pulling a trailer behind a vehicle must be estimated. During real world operation a vehicle often encounters cross winds. Therefore, the effects of cross winds on the drag of a vehicle–trailer combination should be taken into account. Improving the accuracy of aerodynamic load prediction for a vehicle-trailer combination should in turn lead to improved simulations and better thermal performance. In order to best simulate conditions for real world trailer towing, a study was performed using reduced scale models of a Sport Utility Vehicle (SUV) and a Pickup Truck (PT) towing a medium size cargo trailer. The scale model vehicle and trailer combinations were tested in a full scale wind tunnel.

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

A vehicle driving down the road naturally pitches, rolls and heaves due to road inputs (for example, bumps, potholes, driving dynamics, etc.) and also due to the influence of aerodynamic forces. The vehicle attitude changes directly as a result of aerodynamic forces that can be seen during wind tunnel testing of production level vehicles, with some measurements possible in order to evaluate the aerodynamics effects. This naturally occurring phenomenon is not always represented in aerodynamics simulations, either for reduced scale models or computational fluid dynamics (CFD) simulations or even rigid body full scale testing. It can be shown through visual techniques how much deflection is typically occurring, including both vehicle attitude changes as well as vehicle body distortions. From the analysis, an adjustment to the CFD models can be made to compensate for the aerodynamics effects.

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.

Due to stringent environmental requirements, the vehicle under-hood and underbody temperatures have been steadily increasing. The increased temperatures affect components life and therefore, more thermal protection measures may be necessary. In this paper, we present an algorithm for estimation of automotive component life due to thermal effects through the vehicle life. Traditional approaches consider only the maximum temperature that a component will experience during severe driving maneuvers. However, that approach does not consider the time duration or frequency of exposure to temperature. We have envisioned a more realistic and science based approach to estimate component life based on vehicle duty cycles, component temperature profile, frequency and characteristics of material thermal degradation. In the proposed algorithm, a transient thermal analysis model provides the exhaust gas and exhaust surface temperatures for all exhaust system segments, and for any driving scenario.

The application of Computational Fluid Dynamics (CFD) for three-dimensional (3D) combustion analysis coupled with detailed chemistry in engine development is hindered by its expensive computational cost. Chemistry computation may occupy as much as 90% of the total computational cost. In the present paper, a new two-step iso-octane combustion model was developed for spark-ignited (SI) engine to maximize computational efficiency while maintaining acceptable accuracy. Starting from the model constants of an existing global combustion model, the new model was developed using an approach based on sensitivity analysis to approximate the results of a reference skeletal mechanism. The present model involves only five species and two reactions and utilizes only one uniform set of model constants. The validation of the new model was performed using shock tube and real SI engine cases.

In recent years reducing the automobile HVAC (Heating Ventilation and automobile conditioning) noise inside the vehicle cabin is one of the main criterions for all OEMs to provide comfort level to the passengers. The primary function of the HVAC is to deliver more air to the cabin with less noise generation for various blower speeds. Designing the optimum HVAC with less noise is one of the major challenges for all automotive manufacturers and HVAC suppliers. During the design stage, physical parts are not available and hence the simulation technique helps to evaluate the noise level of HVAC. In this study, a computational 1D (one dimensional) analysis is carried out to compute the airflow noise originated from the HVAC unit and propagated to the passenger cabin. Modeling has been done using unigraphics and the analysis is carried out using the commercial 1D software GT suite.