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 study examines the impact of slip in aero-thermal conditions of supercooled large droplets (SLD) produced in an Icing Wind Tunnel (IWT) on the ice accretion characteristics. The study identifies potential biases in the SLD model development based on IWT data and numerical predictions that assume the SLD to be in aerothermal equilibrium with the IWT airflow. To obtain realistic temperature and velocity data for each droplet size class in the test section of the Braunschweig Icing Wind Tunnel (BIWT), a Lagrangian droplet tracking solver was used within a Monte Carlo framework. Results showed that SLDs experience considerable slips in velocity and temperature due to their higher inertia and short residence time in the Braunschweig IWT. Large droplets were found to be warmer and slower than the flow in the test section, with larger droplets experiencing larger aerothermal slips.

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.

The continuous pursuit of higher combustion efficiencies, as well as the possible usage of synthetic fuels with different properties than fossil-ones, require reliable and low-cost numerical approaches to support and speed-up engines industrial design. In this context, SI engines operated with homogeneous ultra-lean mixtures both characterized by a classical ignition configuration or equipped with an active prechamber represent the most promising solutions. In this work, for the classical ignition arrangement, a 3DCFD strategy to model the impact of the ignition system type on the CCV is developed using the RANS approach for turbulence modelling. The spark-discharge is modelled through a set of Lagrangian particles, whose velocity is modified with a zero-divergence perturbation at each discharge event, then evolved according to the Simplified Langevin Model (SLM) to simulate stochastic interactions with the surrounding gas flow.

In the recent years, the interest in heavy-duty engines fueled with Compressed Natural Gas (CNG) is increasing due to the necessity to comply with the stringent CO2 limitation imposed by national and international regulations. Indeed, the reduced number of carbon atoms of the NG molecule allows to reduce the CO2 emissions compared to a conventional fuel. The possibility to produce synthetic methane from renewable energy sources, or bio-methane from agricultural biomass and/or animal waste, contributes to support the switch from conventional liquid fuels to CNG. To drive the engine development and reduce the time-to-market, the employment of numerical analysis is mandatory. This requires a continuous improvement of the simulation models toward real predictive analyses able to reduce the experimental R&D efforts. In this framework, 1D numerical codes are fundamental tools for system design, energy management optimization, and so on.

Ducted Fuel Injection (DFI) has the potential to reduce soot emissions in Diesel engines thanks to the enhanced mixing rate resulting from the liquid fuel flow through a small cylindrical pipe located at a certain distance from the nozzle injector hole. A consolidated set of experiments in constant-volume vessel and engine allowed to understand the effects of ambient conditions, duct geometry and shape on fuel-air mixing, combustion and soot formation. However, implementation of this promising technology in compression-ignition engines requires predictive numerical models that can properly support the design of combustion systems in a wide range of operating conditions. This work presents a computational methodology to predict fuel-air mixing and combustion with ducted fuel injection. Attention is mainly focused on turbulence and combustion modelling.

Cycle-To-Cycle Variability (CCV) must be properly considered when modeling the ignition process in SI engines operating with ultra-lean mixtures. In this work, a strategy to model the impact of the ignition type on the CCV was developed using the RANS approach for turbulence modelling, performing multi-cycle simulations for the power-cycle only. The spark-discharge was modelled through a set of Lagrangian particles, introduced along the sparkgap and interacting with the surrounding Eulerian gas flow. Then, at each discharge event, the velocity of each particle was modified with a zero-divergence perturbation of the velocity field with respect to average conditions. Finally, the particles velocity was evolved according to the Simplified Langevin Model (SLM), which keeps memory of the initial perturbation and applies a Wiener process to simulate the stochastic interaction of each channel particle with the surrounding gas flow.

In the next years, the upcoming emission legislations are expected to introduce further restrictions on the admittable level of pollutants from vehicles measured on homologation cycles and real drive tests. In this context, the strict control of pollutant emissions at the cold start will become a crucial point to comply with the new regulation standards. This will necessarily require the implementation of novel strategies to speed-up the light-off of the reactions occurring in the after-treatment system, since the cold start conditions are the most critical one for cumulative emissions. Among the different possible technological solutions, this paper focuses on the evaluation of the potential of a burner system, which is activated before the engine start. The hypothetical burner exploits the lean combustion of an air-gasoline mixture to generate a high temperature gas stream which is directed to the catalyst section promoting a fast heating of the substrate.

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

In this work the 3Dcell method, a quasi3D approach developed by the Internal Combustion Engine Group at Politecnico di Milano, has been extended and applied to the fluid dynamic simulation of turbocharging devices for internal combustion engines, focusing on the compressor side. The 3Dcell is based on a pseudo-staggered leapfrog method applied to the governing equation of a 1D problem arbitrarily oriented in space. The system of equations is solved referring to the relative system in the rotating zone, whereas the absolute reference system has been used elsewhere. The vaneless diffuser has been modelled resorting to the conservation of the angular momentum of the flow stream in the tangential direction, combined with the solution of the momentum equation in the radial direction.

Turbulence modeling for fuel spray simulation plays a prominent role in the understanding of the flow behavior in Internal Combustion Engines (ICEs). Currently, a lot of research work is actively spent on Large Eddy Simulation (LES) turbulence modeling as a replacement option of standard Reynolds averaged approaches in the Eulerian-Lagrangian spray modeling framework, due to its capability to accurately describe flow-induced spray variability and to the lower dependence of the results on the specific turbulence model and/or modeling coefficients. The introduction of LES poses, however, additional questions related to the implementation/adaptation of spray-related turbulence sources and to the rise of conflicting numerics and grid requirements between the Lagrangian and Eulerian parts of the simulated flow.

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.

Metallic open-cell foams have proven to be valuable for many engineering applications. Their success is mainly related to mechanical strength, low density, high specific surface, good thermal exchange, low flow resistance and sound absorption properties. The present work aims to investigate three principal aspects of real foams: the geometrical characterization, the flow regime characterization, the effects of the pore size and the porosity on the pressure drop. The first aspect is very important, since the geometrical properties depend on other parameters, such as porosity, cell/pore size and specific surface. A statistical evaluation of the cell size of a foam sample is necessary to define both its geometrical characteristics and the flow pattern at a given input velocity. To this purpose, a procedure which statistically computes the number of cells and pores with a given size has been implemented in order to obtain the diameter distribution.

The interaction of turbulent premixed methane combustion with the surrounding flow field can be studied using optically accessible test rigs such as a rapid compression expansion machine (RCEM). The high flexibility offered by such a test rig allows its operation at various thermochemical conditions at ignition. However, limitations inherent to such test rigs due to the absence of an intake stroke do not allow turbulence production as found in IC-engines. Hence, means to introduce turbulence need to be implemented and the relevant turbulence quantities have to be identified in order to enable comparability with engine relevant conditions. A dedicated high-pressure direct injection of air at the beginning of the compression phase is considered as a measure to generate adjustable turbulence intensities at spark timing and during the early flame propagation.

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.

This paper reports the studies, design and developments of an electronic electro-actuated gearshifter installed on the DP7, which is Politecnico di Milano car that took part at Formula SAE 2015 competitions in Hockenheim and Varano dè Melegari. The original idea was born to replace the hydraulic gearshift system used until 2011 because of its high weight and cost. After many evaluations about the kind of technology to use, made by previous team members in the electronic department, the final project was a fully electric shifter. This system has proven its qualities among which are lightness and low cost.

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.

The objective of this study is to demonstrate the design and construction of an innovative active gear-shift and clutch for racecars, applied to a Formula Student car, based on the use of DC gear-motors. Racecars require extremely quick gear-shifts and every system to be as light as possible. The proposed solution is designed to reduce energy consumption, weight and improve gear-shift precision compared to traditionally employed electro-hydraulic solutions, although maintaining state of the art performances.

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.