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.

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.

As the number of actuators and sensors increases in modern combustion engines, the task of optimizing engine performance becomes increasingly complex. Efficient information processing techniques are therefore important, both for off-line calibration of engine maps, and on-line adjustments based on sensor data. In-cylinder pressure sensors are slowly spreading from laboratory use to production engines, thus making data with high temporal resolution of the combustion process available. The standard way of using the cylinder pressure data for control and diagnostics is to focus on a few important physical features extracted from the pressure trace, such as the combustion phasing CA50, the indicated mean effective pressure IMEP, and the ignition delay. These features give important information on the combustion process, but much information is lost as the information from the high-resolution pressure trace is condensed into a few key parameters.

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 offers a two-zone NOx model suitable for vehicle on-board, on-line implementation. Similar NOx modeling attempts have previously been undertaken. The hereby suggested method does however offer clear and important benefits over the previously methods, utilizing a significantly different method to handle temperature calculations within the (two) different zones avoiding iterative computation. The new method significantly improves calculation speed and, most important of all, reduces implementation complexity while still maintaining reasonable accuracy and the physical interpretation of earlier suggested methods. The equations commonly used to compute NOx emissions is also rewritten in order to suit a two-zone NOx model. An algorithm which can be used to compute NOx emissions is presented and the intended contribution of the paper is a NOx model, implementation feasible for an embedded system, e.g. embedded processor or embedded electronic hardware (FPGA).

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.

The emission control in heavy-duty vehicles today is based on predefined injection strategies and after-treatment systems such as SCR (selective catalytic reduction) and DPF (diesel particulate filter). State-of-the-art engine control is presently based on cycle-to-cycle resolution. The introduction of the crank angle resolved pressure measurement, from a piezo-based pressure sensor, enables the possibility to control the fuel injection based on combustion feedback while the combustion is occurring. In this paper a study is presented on the possibility to control NOx (nitrogen oxides) formation with a crank angle resolved NOx estimator as feedback. The estimator and the injection control are implemented on an FPGA (Field-Programmable Gate Array) to manage the inherent time constraints. The FPGA is integrated with the rest of the engine control system for injection control and measurement.

In recent years, stricter regulations on emissions and higher demands for more fuel efficient vehicles have led to a greater focus on increasing the efficiency of the internal combustion engine. Nowadays, there is increasing interest in the recovery of waste heat from different engine sources such as the coolant and exhaust gases using, for example, a Rankine cycle. In diesel engines 15% to 30% of the energy from the fuel can be lost to the coolant and hence, does not contribute to producing work on the piston. This paper looks at reducing the heat losses to the coolant by increasing coolant temperatures within a single cylinder Scania D13 engine and studying the effects of this on the energy balance within the engine as well as the combustion characteristics. To do this, a GT Power model was first validated against experimental data from the engine.

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.

In the last couple of decades, countries have enacted new laws concerning environmental pollution caused by heavy-duty commercial and passenger vehicles. This is done mainly in an effort to reduce smog and health impacts caused by the different pollutions. One of the legislated pollutions, among a wide range of regulated pollutions, is nitrogen oxides (commonly abbreviated as NOx). The SCR (Selective Catalytic Reduction) was introduced in the automotive industry to reduce NOx emissions leaving the vehicle. The basic idea is to inject a urea solution (AdBlue™) in the exhaust gas before the gas enters the catalyst. The optimal working temperature for the catalyst is somewhere in the range of 300 to 400 °C. For the reactions to occur without a catalyst, the gas temperature has to be at least 800 °C. These temperatures only occur in the engine cylinder itself, during and after the combustion.

The novel wave-shaped bowl piston geometry design with protrusions has been proved in previous studies to enhance late-cycle mixing and therefore significantly reduce soot emissions and increase engine thermodynamic efficiency. The wave-shaped piston is characterized by the introduction of evenly spaced protrusions around the inner wall of the bowl, with a matching number with the number of injection holes, i.e., flames. The interactions between adjacent flames strongly affect the in-cylinder flow and the wave shape is designed to guide the near-wall flow. The flow re-circulation produces a radial mixing zone (RMZ) that extends towards the center of the piston bowl, where unused air is available for oxidation promotion. The waves enhance the flow re-circulation and thus increase the mixing intensity of the RMZ.

This article deals with study of ionization current sensing technique's signal characteristics while operating with pre-chamber spark plug to achieve plasma jet ignition in a 6 cylinder 9 liter turbo-charged natural gas engine under EGR and excess air dilution. Unlike the signal with conventional spark plug which can be divided into distinct chemical and thermal ionization peaks, the signal with pre-chamber spark plug shows a much larger first peak and a negligible second peak thereafter. Many studies in past have found the time of second peak coinciding with the time of maximum cylinder pressure and this correlation has been used as an input to combustion control systems but the absence of second peak makes application of this concept difficult with pre-chamber spark plug.