Search Results

Polymer electrolyte membrane (PEM) fuel cells will play a crucial role in the decarbonization of the transport sector, in particular for heavy duty applications. However, performance and durability of PEMFC stacks is still a concern especially when operated under high power density conditions, as required in order to improve the compactness and to reduce the cost of the system. In this context, the optimization of the geometry of hydrogen and air distributors represents a key factor to improve the distribution of the reactants on the active surface, in order to guarantee a proper water management and avoiding membrane dehydration.

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.

We present a framework for the robust optimization of the heat flux distribution for an anti-ice electro-thermal ice protection system (AI-ETIPS) and iced airfoil performance analysis under uncertain conditions. The considered uncertainty regards a lack of knowledge concerning the characteristics of the cloud i.e. the liquid water content and the median volume diameter of water droplets, and the accuracy of measuring devices i.e., the static temperature probe, uncertain parameters are modeled as uniform random variables. A forward uncertainty propagation analysis is carried out using a Monte Carlo approach. The optimization framework relies on a gradient-free algorithm (Mesh Adaptive Direct Search) and three different problem formulations are considered in this work. Two bi-objective deterministic optimizations aim to minimize power consumption and either minimize ice formations or the iced airfoil drag coefficient.

Multi-purpose agricultural tractors are vehicles that are usually used in rough paths and on off-road situations characterized by strong slope variations. The main feature of this kind of vehicles is the stability in working conditions to avoid overturning while it is on duty. This characteristic is given by the interaction between the suspension system and the vehicle frame. Due to the limited size of this kind of vehicle, the stability feature could be given by chassis deformation or using a two-piece frame connected by a spherical joint. This paper presents the validation of a numerical lumped-parameters model able to reproduce the vertical dynamics of a multi-purpose tractor featured by a yielding chassis. The unknown model parameters have been estimated firstly with static tests to study the vertical tire and suspension stiffnesses. The dynamic tests using a four-post-test rig have been performed to tune the unknown dynamic parameters.

Wheel and wheelhouses contribute up to 20-30% of the aerodynamic drag of passenger cars. Simulating the flow field around wheels is challenging due to the complexity of the flow structures generated by tires and rims, wheel rotation, tire deformation and contact with the ground. High accuracy is usually obtained with transient simulations that treat rim rotation with the Sliding Mesh (SM) approach, which is also computationally expensive. Previous studies have confirmed that the application of a tangential velocity component to the rim surface is unphysical for open rims, while a Moving Reference Frame (MRF) is lacking accuracy and the averaged results depend on the initial spokes position. These methods do not consider the dynamic nature of the problem. This work proposes the use of the Actuator Line (AL) and Rotor Disk (RD) approaches as alternatives for simulating open rims with much lower computational cost.

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.

The paper investigates the dynamics of the tire and brake during hard braking or wheel locking, from the view point of a brake manufacturer. A new test rig, named BRAD (BRembo Automotive Dynamometer) is presented which measures the forces acting both at the brake and at the tire-ground interface. Lateral forces are not measured. In the test rig, the ground is represented by a drum. The features of the test rig are presented. The measurement accuracy is declared. The first result is that, near wheel locking, a substantial part of the braking power is generated by the tire and not by the brake. The test rig quantifies such a partitioning of brake power, which is important for current and future electric motorsport activities. Some 30% of the braking power is due to tire during hard braking. The second result is that, due to such important braking power at the tire, the tire is heated up, which increases considerably the maximum friction.

In recent years, the increase of gasoline fuel injection pressure is a way to improve thermal efficiency and lower engine-out emissions in GDI homogenous combustion concept. The challenge of controlling particulate formation as well in mass and number concentrations imposed by emissions regulations can be pursued improving the mixture preparation process and avoiding mixture inhomogeneity with ultra-high injection pressure values up to 100 MPa. The increase of the fuel injection pressure in GDI homogeneous systems meets the demand for increased injector static flow, while simultaneously improves the spray atomization and mixing characteristics with consequent better combustion performance. Few studies quantify the effects of high injection pressure on transient gasoline spray evolution. The aim of this work was to simulate with OpenFOAM the spray morphology of a commercial gasoline injected in a constant volume vessel by a prototypal GDI injector.

The increasing diffusion of gasoline direct injection (GDI) engines requires a more detailed and reliable description of the phenomena occurring during the fuel injection process. As well known the thermal and fluid-dynamic conditions present in the combustion chamber greatly influence the air-fuel mixture process deriving from GDI injectors. GDI fuel sprays typically evolve in wide range of ambient pressure and temperatures depending on the engine load. In some particular injection conditions, when in-cylinder pressure is relatively low, flash evaporation might occur significantly affecting the fuel-air mixing process. In some other particular injection conditions spray impingement on the piston wall might occur, causing high unburned hydrocarbons and soot emissions, so currently representing one of the main drawbacks of GDI engines.

Aqueous Urea is a non-toxic and stable ammonia carrier and its injection and mixing represent the basis for the most common de-NOx technology for mobile applications. The reactant feed preparation process is defined by evaporation, thermolysis and hydrolysis of the liquid mixture upstream the Selective Catalytic Reduction reactor, and it is strongly dependent on the interaction between spray and gaseous flow. Low-pressure driven injectors are the common industrial standard for these applications, and their behavior in almost-ambient pressure cross flows is significantly different from any in-cylinder application. For this reason, two substantially different injectors in terms of geometry and design are experimentally studied, characterizing drop sizes and velocities through Phase Doppler Anemometry (PDA) and liquid mass spatial distribution through Shadow Imaging (SI).

Very little is currently known of the aerodynamic interaction between neighboring cycloidal rotors. Such knowledge is, however, of crucial importance to tune the controller and rotor disposition of a cyclogiro aircraft. Thus, a three-dimensional computational fluid dynamics (CFD) model is developed, validated, and used to analyze the D-Dalus L1 four-rotor unmanned aircraft operating under several configurations. The model solves the Euler equations using the OpenFOAM toolbox in order to provide fast results on a desktop computer. Validation is performed against thrust forces and flow streamlines obtained during wind tunnel experiments at various flight velocities. Numerical results from CFD match the trends of the experimental data. Flow behavior matches the video footage of the wind tunnel tests. Although boundary layer effects are neglected, satisfactory results are obtained both qualitatively and quantitatively.

This paper provides a review on state-of-art modern cooling systems employed for thermal cooling of electric motors for vehicle applications. In recent years, the pursue of a more sustainable and ecofriendly mobility has pushed the research towards the development of electric vehicle powertrain systems. Besides the evident advantages of the adoption of electric traction systems in terms of pollution and efficiency, the need of an effective cooling system for the electric machine components gained more and more importance in order to maintain high efficiency and ensure high durability. In fact, it is known that high temperatures can be harmful for the electric motor: besides the evident damages for mechanical parts, the influence on the permanent magnet properties is not negligible [1] [2]. In this fast-evolving environment, different solutions for the thermal problem have been researched and adopted, each one with its own pros and cons.

In the last years, as a response to the more and more restrictive emission legislation, new devices (SRC, DOC, NOx-trap, DPF) have been progressively introduced as standard components of modern after-treatment system for Diesel engines. In addition, the adoption of electrical heating is nowadays regarded with interest as an effective solution to promote the light-off of the catalyst at low temperature, especially at the start-up of the engine and during the low load operation of the engine typical of the urban drive. In this work, a state-of-the-art 48 V electrical heated catalyst is considered, in order to investigate its effect in increasing the abatement efficiency of a standard DOC. The electrical heating device considered is based on a metallic support, arranged in a spiral layout, and it is heated by the Joule effect due to the passage of the electrical current.

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.

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.

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.

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.

Today, multi-hole Diesel injectors can be mainly characterized by three different nozzle hole shapes: cylindrical, k-hole, and ks-hole. The nozzle hole layout plays a direct influence on the injector internal flow field characteristics and, in particular, on the cavitation and turbulence evolution over the hole length. In turn, the changes on the injector internal flow correlated to the nozzle shape produce immediate effects on the emerging spray. In the present paper, the fluid dynamic performance of three different Diesel nozzle hole shapes are evaluated: cylindrical, k-hole, and ks-hole. The ks-hole geometry was experimentally characterized in order to find out its real internal shape. First, the three nozzle shapes were studied by a fully transient CFD multiphase simulation to understand their differences in the internal flow field evolutions. In detail, the attention was focused on the turbulence and cavitation levels at hole exit.

Detailed chemistry represents a fundamental pre-requisite for a realistic simulation of combustion process in Diesel engines to properly reproduce ignition delay and flame structure (lift-off and soot precursors) in a wide range of operating conditions. In this work, the authors developed reduced mechanisms for n-dodecane starting from the comprehensive kinetic mechanism developed at Politecnico di Milano, well validated and tested in a wide range of operating conditions [1]. An algorithm combining Sensitivity and Flux Analysis was employed for the present skeletal reduction. The size of the mechanisms can be limited to less than 100 species and incorporates the most important details of low-temperature kinetics for a proper prediction of the ignition delay. Furthermore, the high-temperature chemistry is also properly described both in terms of reactivity and species formation, including unsaturated compounds such as acetylene, whose concentration controls soot formation.