Search Results

Hydrogen engines are currently considered as a viable solution to preserve the internal combustion engine as a power unit for vehicle propulsion. In particular, lean-burn gasoline Spark-Ignition (SI) engines have been a major subject of investigations due to the reduced emission levels and high thermodynamic efficiency. This strategy is suitable for the purpose of passenger car applications and cannot be tailored in the field of high performance engine, where the air mass delivered would require oversized turbocharging systems or more complex charging solutions. For this reason, the range of stoichiometric feeding condition is explored in the high performance engine, leading to the consequent issue of abatement of pollutant emissions. In this work a 1D model will be applied to the modeling of a V8 engine fueled with DI of hydrogen. The engine has been derived by a gasoline configuration and adapted to hydrogen in such a way to keep the same performance.

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.

Nowadays quasi-3D approaches are included in many commercial and research 1D numerical codes, in order to increase their simulation accuracy in presence of complex shape 3D volumes, e.g. plenums and silencers. In particular, these are regarded as valuable approaches for application during the design phase of an engine, for their capability of predicting non-planar waves motion and, on the other hand, for their low requirements in terms of computational runtime. However, the generation of a high-quality quasi-3D computational grid is not always straightforward, especially in case of complex elements, and can be a time-consuming operation, making the quasi-3D tool a less attractive option. In this work, a quasi-3D module has been implemented on the basis of the open-source CFD code OpenFOAM and coupled with the 1D code GASDYN.

The paper investigates the interaction between soil and tractor tires through a 2D numerical model. The tire is schematized as a rigid ring presenting a series of rigid tread bars on the external circumference. The outer profile of the tire is divided into a series of elements, each one able to exchange a normal and a tangential contact force with the ground. A 2D soil model was developed to compute the forces at the ground-tire interface: the normal force is determined on the basis of the compression of the soil generated by the sinking of the tire. The soil is modeled through a layer of springs characterized by two different stiffness for the loading (lower stiffness) and unloading (higher stiffness) condition. This scheme allows to introduce a memory effect on the soil which results stiffer and keeps a residual sinking after the passage of the tire. The normal contact force determines the maximum value of tangential force provided before the soil fails.

Active prechamber combustion systems for SI engines represent a feasible and effective solution in reducing fuel consumption and pollutant emissions for both marine and ground heavy-duty engines. However, reliable and low-cost numerical approaches need to be developed to support and speed-up their industrial design considering their geometry complexity and the involved multiple flow length scales. This work presents a CFD methodology based on the RANS approach for the simulation of active prechamber spark-ignition engines. To reduce the computational time, the gas exchange process is computed only in the prechamber region to correctly describe the flow and mixture distributions, while the whole cylinder geometry is considered only for the power-cycle (compression, combustion and expansion). Outside the prechamber the in-cylinder flow field at IVC is estimated from the measured swirl ratio.

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.

A correct prediction of the initial stages of the combustion process in SI engines is of great importance to understand how local flow conditions, fuel properties, mixture stratification and ignition affect the in-cylinder pressure development and pollutant formation. However, flame kernel growth is governed by many interacting processes including energy transfer from the electrical circuit to the gas phase, interaction between the plasma channel and the flow field, transition between different combustion regimes and gas expansion at very high temperatures. In this work, the authors intend to present a comprehensive, multi-dimensional model that can be used to predict the initial combustion stages in SI engines. In particular, the spark channel is represented by a set of Lagrangian particles where each one of them acts as a single flame kernel.

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

SI engines fueled with hydrogen represent a promising powertrain solution to meet the ambitious target of carbon-free emissions at the tailpipe. Therefore, fast and reliable numerical tools can significantly support the automotive industry in the optimization of such technology. In this work, a 1D-3D methodology is presented to simulate in detail the combustion process with minimal computational effort. First, a 1D analysis of the complete engine cycle is carried out on the user-defined powertrain configuration. The purpose is to achieve reliable boundary conditions for the combustion chamber, based on realistic engine parameters. Then, a 3D simulation of the power-cycle is performed to mimic the combustion process. The flow velocity and turbulence distributions are initialized without the need of simulating the gas exchange process, according to a validated technique.

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 deals with the design and manufacturing of a McPherson suspension arm made from short glass fiber reinforced polyamide (PA66). The design of the arm and the design of the molds have been made jointly. According to Industry 4.0 paradigms, a full digitalization of both the product and process has been performed. Since the mechanical behavior of the suspension arm strongly depends on constraints which are difficult to be modelled, a simpler structure with well-defined mechanical constraints has been developed. By means of such simple structure, the model for the behavior of the material has been validated. Since the suspension arm is a hybrid structure, the associated simple structure is hybrid as well, featuring a metal sheet with over-molded polymer. The issues referring to material flow, material to material contact, weld lines, fatigue strength, high and low temperature behavior, creep, dynamic strength have been investigated on the simple structure.

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.

A new electric powertrain and axle for light/medium trucks is presented. The indoor testing and the simulation of the dynamic behavior are performed. The powertrain and axle has been produced by Streparava and tested at the Laboratory for the Safety of Transport of the Politecnico di Milano. The tests were aimed at defining the multi-physics perfomance of the powertrain and axle (efficiency, acceleration and braking, temperature and NVH). The whole system for indoor tests was composed by the powertrain and axle (electric motor, driveline, suspensions, wheels) and by the test rig (drums, driveline and electric motor). The (driving) axle was positioned on a couple of drums, and the drums provided the proper torques to the wheels to reproduce acceleration and braking. Additionally a cleat fixed on one drum excited the vibration of the suspensions and allowed assessing NVH performance. The simulations were based on a special co-simulation between 1D-AMESIM and VIRTUAL.LAB.

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.

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.

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.

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

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.