Search Results

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.

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.

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.

Nowadays a great attention is paid to the level and quality of noise radiated from the tailpipe end of intake and exhaust systems, to control the gas dynamic noise emitted by the engine as well as the characteristics of the cabin interior sound. The muffler geometry can be optimized consequently, to attenuate or remark certain spectral components of the engine noise, according to the result expected. Evidently the design of complex silencing systems is a time-consuming operation, which must be carried out by means of concurrent experimental measurements and numerical simulations. In particular, 1D and multiD linear/non-linear simulation codes can be applied to predict the silencer behavior in the time and frequency domain. This paper describes the development of a 1D-multiD integrated approach for the simulation of complex muffler configurations such as reverse chambers with inlet and outlet pipe extensions and perforated silencers with the addition of sound absorbing material.

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.

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.

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.

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 durability performance of brake hoses is a crucial issue for such components. Accelerated fatigue testing of brake hoses is necessary for understanding achievable lifetime, actually computation of durability is quite cumbersome due to the many different materials the hoses are made from. Despite SAE standards are available, accelerated testing of brake hoses subject to actual torsional and bending stresses seem important to provide relevant feedback to designers. In this paper, an innovative methodology for assessing the fatigue behavior of brake hoses of road vehicles is proposed. A dynamic testbed is specifically designed and realized, able to reproduce the actual assembly conditions of the hoses fitted into a vehicle suspension. The designed testbed allows to replicate actual loading conditions on the brake hoses by simulating the vertical dynamics and steering of the suspension system together with brake pressure.

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.

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.

Modeling the combustion process of a diesel-biodiesel fuel spray in a 3-dimensional (3D) computational fluid dynamics (CFD) domain remains challenging and time-consuming despite the recent advancement in computing technologies. Accurate representation of the in-cylinder processes is essential for CFD studies to provide invaluable insights into these events, which are typically limited when using conventional experimental measurement techniques. This is especially true for emerging new fuels such as biodiesels since fundamental understanding of these fuels under combusting environment is still largely unknown. The reported work here is dedicated to evaluating the Adaptive Local Mesh Refinement (ALMR) approach in OpenFOAM® for improved simulation of reacting biodiesel fuel spray. An in-house model for thermo-physical and transport properties is integrated to the code, along with a chemical mechanism comprising 113 species and 399 reactions.

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 definition of a robust methodology to perform a full-cycle CFD simulation of IC engines requires as first step the availability of a reliable grid generation tool, which does not only have to guarantee a high quality mesh but also has to prove to be efficient in terms of required time. In this work the authors discuss a novel approach entirely based on the OpenFOAM technology, in which the available 3D grid generator was employed to automatically create meshes containing hexahedra and split-hexahedra from triangulated surface geometries in Stereolithography (STL) format. The possibility to introduce local refinements and boundary layers makes this tool suitable for IC engine simulations. Grids are sequentially generated at target crank angles which are automatically determined depending on user specified settings such as maximum mesh validity interval and quality parameters like non-orthogonality, skewness and aspect ratio.

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 most common automotive drivelines transmit the engine torque to the driven axle through a differential. Semi-active versions of this device ([4], [5], [6]) have been recently conceived to improve vehicle handling at limit and under particular conditions; these differentials are based on the structural scheme of the passive one but they try to manipulate the vehicle dynamics by controlling the distribution of the driving torque on the wheels of the same axle thus generating a yaw moment. Unfortunately a semi-active differential is not able to perform a complete yaw control since the torque can only be transferred from the faster wheel to the slower one; on the other hand, active differentials ([11], [12], [13]) allow to generate the most appropriate yaw moment controlling both the amount of transferred torque and its direction.