Search Results

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

The paper presents a method for the indoor testing of road vehicle suspension systems. A suspension is positioned on a rotating drum which is located in the Laboratory for the Safety of Transport at Politecnico di Milano. Special six-axis load cells have been designed and used for measuring the forces/moments acting at each suspension-chassis joints. The forces/moments, wheel accelerations, displacements are measured up to 100 Hz. Two different types of test can be performed. The tire/wheel unbalance effect on the suspension system behavior (Vibration and Harshness, VH) has been analyzed by testing the suspension system from zero to the vehicle maximum speed on a flat surface and by monitoring the forces transmitted to the chassis. In the second kind of test, the suspension system has been excited as the wheel passes over different cleats fixed on the drum.

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.

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.

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.

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.

The paper focuses on the development of a mesh moving method based on non-conformal topologically changing grids applied to the simulation of IC engines, where the prescribed motion of piston and valves is accomplished by rigidly translating the sub-domain representing the moving component. With respect to authors previous work, a more robust and efficient algorithm to handle the connectivity of non-conformal interfaces and a mesh-motion solver supporting multiple layer addition/removal of cells, to decouple the time-step constraints of the mesh motion and of the fluid dynamics, has been implemented as a C++ library to extend the already existing classes for dynamic mesh handling of the finite-volume, open-source CFD code OpenFOAM®. Other new features include automatic decomposition of large multiple region domains to preserve processors load balance with topological changes for parallel computations and additional tools for automatic preprocessing and case setup.

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.

The tightening trend of regulations on the levels of admitted pollutant emissions has given a great spur to the research work in the field of combustion and after-treatment devices. Despite the improvements that can be applied to the development of the combustion process, pollutant emissions cannot be reduced to zero; for this reason, the aftertreatment system will become a key component in the path to achieving near-zero emission levels. This study focuses on the numerical analysis and optimization of different metallic substrates, specifically developed for three-way catalyst (TWC) and Diesel oxidation catalyst (DOC) applications, to improve their thermal efficiency by reducing radial thermal losses through the outer mantle. The optimization process relies on computational fluid dynamics (CFD) simulations supported by experimental measurements to validate the numerical models carried out under uncoated conditions, where chemical reactions do not occur.

Electronic Stability Program (ESP) and Antilock Braking System (ABS) are nowadays a standard equipment for passenger cars. ESP increases vehicle safety by applying differential braking torque to the wheels while cornering, thus it extends the area of intervention of ABS which prevents the wheels from being locked up in emergency braking, especially on low friction road surfaces, allowing the driver to maintain steering control of the vehicle, to avoid obstacles and to reduce vehicle stopping distance on most road surfaces. This paper describes a flexible mechatronic test bench for ESP/ABS Electronic Control Unit (ECU) based on Hardware-In-the-Loop (HIL) simulation technique. It consists of a passenger car hydraulic braking system (from master cylinder to brake calipers), with the ESP/ABS ECU integrated and a flexible real-time platform, which simulates vehicle dynamics.

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.

The paper deals with the bifurcation analysis of a simple mathematical model describing an automobile running on an even surface. Bifurcation analysis is adopted as the proper procedure for an in-depth understanding of the stability of steady-state motion of cars (either cornering or running straight ahead). The aim of the paper is providing the fundamental information for inspiring further studies on vehicle dynamics with or without a human driver. The considered mechanical model of the car has two degrees of freedom, nonlinear tire characteristics are included. A simple driver model is introduced. Experimental validations of the model are produced. As a first step, bifurcation analysis is performed without driver (fixed control). Ten different combinations of front and rear tire characteristics (featuring understeer or oversteer automobiles) are considered. Steering angle and speed are varied. Many different dynamical behaviors of the model are found.

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.

The employment of hydrogen as energy carrier for transportation sector represents a significant challenge for powertrains. Spark-ignition (SI) engines are feasible and low-cost devices to convert the hydrogen chemical energy into mechanical work. However, significant efforts are needed to successfully retrofit the available configurations. The computational fluid dynamics (CFD) modelling represents a useful tool to support experiments, clarifying the impact of the engine characteristics on both the mixture preparation and the combustion development. In this work, a CFD investigation is carried out on typical light-duty SI engine configurations, exploring the two main strategies of hydrogen addition: port fuel injection (PFI) and direct injection (DI). The purpose is to assess the behaviour of widely-used numerical models and methodologies when hydrogen is employed instead of traditional carbon-based fuels.

NH₃/urea SCR is a very effective and widely used technology for the abatement of NOx from diesel exhaust. The SCR mechanism is well understood and the catalyst behavior can be predicted by mathematical models - as long as operation above the temperature limit for AdBlue® injection is considered. The behavior below this level is less understood. During the first seconds up to minutes after cold start, complete NOx abatement can be observed over an SCR catalyst in test bench experiments, together with a significant increase in temperature after the converter (ca. 100 K). In this work these effects have been investigated over a monolith Cu-zeolite SCR catalyst. Concentration step experiments varying NO, NO₂ and H₂O have been carried out in lab scale, starting from room temperature. Further, the interaction of C₃H₆ and CO with NOx over the SCR has been investigated.

Metallic foams or sponges are materials with a cell structure suitable for many industrial applications, such as reformers, heat catalytic converters, etc. The success of these materials is due to the combination of various characteristics such as mechanical strength, low density, high specific surface, good thermal exchange properties, low flow resistance and sound absorption. Different materials and manufacturing processes produce different type of structure and properties for various applications. In this work a genetic algorithm has been developed and applied to support the design of catalytic devices. In particular, two substrates were considered, namely the traditional honeycomb and an alternative open-cell foam type. CFD simulations of pressure losses and literature based correlations for the heat and mass transfer were used to support the genetic algorithm in finding the best compromise between flow resistance and pollutant abatement.

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.

A multi-objective optimal design of a brushless DC electric motor for a brake system application is presented. Fifteen design variables are considered for the definition of the stator and rotor geometry, pole pieces and permanent magnets included. Target performance indices (peak torque, efficiency, rotor mass and inertia) are defined together with design constraints that refer to components stress levels and temperature thresholds, not to be surpassed after heavy duty cycles. The mathematical models used for optimization refer to electromagnetic field and related currents computation, to thermo-fluid dynamic simulation, to local stress and vibration assessment. An Artificial Neural Network model, trained with an iterative procedure, is employed for global approximation purposes. This allows to reduce the number of simulation runs needed to find the optimal configurations. Some of the Pareto-optimal solutions resulting from the optimal design process are analysed.

The ride comfort of three Alfa Romeo cars, namely Giulietta (1955), Alfetta (1972) and 159 (2005) has been assessed both objectively and subjectively. The three cars belong to the same market segment. The aim is to let young engineers or graduate students understand how technology has evolved and eventually learn a lesson from the assessed trend. A number of cleats have been fixed at the ground and the three cars have traversed such uneven surface. The objective assessment of the ride comfort has been performed by means of accelerometers fixed at the seat rails, additionally a special dummy developed at Politecnico di Milano has been employed. The subjective assessment has been performed by a panel of passengers. The match between objective and subjective ratings is very good. Simple mathematical models have been employed to establish a (successful) comparison between experimental and computational results. The ride comfort differs substantially among the cars.