Search Results

This paper presents a nonlinear control approach to achieve good performances in vehicle path following and collision avoidance when the vehicle is driving under cruise highway conditions. Nonlinear model predictive control (NLMPC) is adopted to achieve online trajectory control based on a simplified vehicle model. GMRES/Continuation algorithm is used to solve the online optimization problem. Simulations show that the proposed controller is capable of tracking the desired path as well as avoiding the obstacles.

The friction estimation at the tire-ground contact is crucial for the active safety of vehicles. Friction estimation is a key problem of vehicle dynamics and the ultimate solution is still unknown. However the proposed approach, based on a simple idea and on a simple hardware, provides an actual solution. The idea is to compare the tire characteristic at a given friction (nominal characteristic) with the actual characteristic that the tire has while running. The comparison among these two characteristics (the nominal one and the actual one) gives the desired friction coefficient. The friction coefficient is expressed in vector form and a number of running parameters are identified. The mentioned comparison is an efficient but complex algorithm based on a mathematical formulation of the tire characteristic. The actual tire characteristic is somehow measured in real time by a relatively simple smart wheel which is able to detect the three forces and the three moments acting at the hub.

The modeling of fuel sprays under well-characterized conditions relevant for heavy-duty Diesel engine applications, allows for detailed analyses of individual phenomena aimed at improving emission formation and fuel consumption. However, the complexity of a reacting fuel spray under heavy-duty conditions currently prohibits direct simulation. Using a systematic approach, we extrapolate available spray models to the desired conditions without inclusion of chemical reactions. For validation, experimental techniques are utilized to characterize inert sprays of n-dodecane in a high-pressure, high-temperature (900 K) constant volume vessel with full optical access. The liquid fuel spray is studied using high-speed diffused back-illumination for conditions with different densities (22.8 and 40 kg/m3) and injection pressures (150, 80 and 160 MPa), using a 0.205-mm orifice diameter nozzle.

Large Eddy Simulation (LES) represents nowadays one of the most promising techniques for the evaluation of the dynamics and evolution of turbulent structures characterizing internal combustion engines (ICE). In the present paper, subdivided into two parts, the capabilities of the open-source CFD code OpenFOAM® v2.3.0 are assessed in order to evaluate its suitability for engine cold flow LES analyses. Firstly, the code dissipative attitude is evaluated through an inviscid vortex convection test to ensure that the levels of numerical dissipation are compatible with LES needs. Quality and completeness estimators for LES simulations are then proposed. In particular the Pope M parameter is used as a LES completeness indicator while the LSR parameter provides useful insights far calibrating the grid density. Other parameters such as the two-grid LESIQk index are also discussed.

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 flow field resulting from injecting a gas jet into a crossflow confined in a narrow square duct has been studied under steady regime using schlieren imaging and laser Doppler velocimetry (LDV). This transparent duct is intended to simulate the intake port of an internal combustion engine fueled by gaseous mixture, and the jet is issued from a round nozzle. The schlieren images show that the relative small size of the duct would confine the development of the transverse jet, and the interaction among jet and sidewalls strongly influences the mixing process between jet and crossflow. The mean velocity and turbulence fields have been studied in detail through LDV measurements, at both center plane and several cross sections. The well-known flow feature formed by a counter rotating vortex pair (CVP) has been observed, which starts to appear at the jet exit section and persists far downstream contributing to enhancing mixing process.

This paper is focused on the development and application of a CFD methodology that can be applied to predict the fuel-air mixing process in stratified charge, sparkignition engines. The Eulerian-Lagrangian approach was used to model the spray evolution together with a liquid film model that properly takes into account its effects on the fuel-air mixing process into account. However, numerical simulation of stratified combustion in SI engines is a very challenging task for CFD modeling, due to the complex interaction of different physical phenomena involving turbulent, reacting and multiphase flows evolving inside a moving geometry. Hence, for a proper assessment of the different sub-models involved a detailed set of experimental optical data is required. To this end, a large experimental database was built by the authors.

Dynamic tests have been performed on carbon fiber racing seats following the FIA regulations. The tests have shown, in rear impact tests, a relatively strong rebound leading to large forward bending of neck, and, in side impact tests, very large lateral displacement of the head, the latter protruding dangerously towards hard portions of the car structure. Stiffening the seat back by steel struts results in reducing strongly both the motion and the acceleration of the head. Simulations of the dynamics of the tests have been done with multi-body models, including the Hybrid III dummy and seat deflection, by means of the program VEDYAC. It has been found that computer simulation can predict very accurately the result of a test, provided the numerical models have been carefully calibrated to match the dummy tolerance bands. Once they have been calibrated and validated with a number of tests, the computer models can be very useful to extend the test results to different test conditions.

Shock absorbers and damper systems are important parts of automobiles and motorcycles because they have effects on safety, ride comfort, and handling. In particular, for vehicle safety, shock absorber system plays a fundamental role in maintaining the contact between tire and road. Generally, to assure the best trade-off between safety and ride comfort, a fine experimental tuning on all shock absorber components is necessary. Inside a common damper system the presence of several conjugated actions made by springs, oil and pressurized air requires a significant experimental support and a great number of prototypes and test. Aimed to reduce the design and tuning phases of a damper system, it is necessary to join these phases together with a numerical modelling phase. The aim of this paper is to present the development of a mono-dimensional (1D) model for simulating dynamic behaviour of damper system.

The paper presents a new test rig (named RuotaVia) composed basically by a drum (2,6 m diameter), providing a running contact surface for vehicle wheels. A number of measurements on either full vehicles or vehicle sub-systems (single suspension system or single tyre) can be performed. Tire characteristics influencing rollover can be assessed. The steady-state maximum loads are as follows: Radial: 100kN, tangential: 100kN, lateral (axial with respect to the drum): 100kN. The superstructure carrying a measuring hub can excite the wheel under test up to 20 Hz in lateral and vertical directions. The steer angle range is ± 25 deg, the camber range is ± 80 deg. The minimum eigenfrequency of the drum is higher than 90 Hz and its maximum tangential speed is 440 km/h.

A three dimensional numerical simulation of the ECN “Spray A” is presented. Both primary and secondary breakup of the spray are included. The fuel is n-Dodecane. The n-Dodecane kinetic mechanism is modeled using a skeletal mechanism that consists of 103 species and 370 reactions [9]. The kinetic mechanism is computationally heavy when coupled with three dimensional numerical simulations. Multidimensional chemistry coordinate mapping (CCM) approach is used to speedup the simulation. CCM involves two-way mapping between CFD cells and a discretized multidimensional thermodynamic space, the so called multidimensional chemistry coordinate space. In the text, the cells in the discretized multidimensional thermodynamic space are called zone to discriminate them from the CFD cells. In this way, the CFD cells which are at the similar thermodynamic state are identified and grouped into a unique zone. The stiff ODEs operates only on the zones containing at least one CFD cell.

In the most elementary treatment of plane-wave reflection at the open end of a duct system, it is often assumed that the ends are pressure nodes. This implies that pressure is assumed as a constant at the open end termination and that steady flow boundary condition is supposed as instantaneously established. While this simplifying assumption seems reasonable, it does not consider any radiation of acoustic energy from the duct into the surrounding free space; hence, an error in the estimation of the effects of the flow on the acoustical response of an open-end duct occurs. If radiation is accounted, a complicated three-dimensional wave pattern near the duct end is established, which tends to readjust the exit pressure to its steady-flow level. This adjustment process is continually modified by further incident waves, so that the effective instantaneous boundary conditions which determine the reflected waves depend on the flow history.

This paper describes the design and development of a virtual environment conceived to support flight operations of an Unmanned Air Vehicle (UAV) used for wind mapping in the proximity of existing or planned wind farms. The virtual environment can be used in pre-flight briefings aiming to define a trajectory from a list of waypoints, to change and eventually re-plan the mission in case of intersection with no fly zones, to simulate the mission, and to preview images/videos taken from the UAV on-board cameras. During flight, the tool can be used to compute the wind speed along the trajectory by analyzing the data streaming from the UAV. The integration of Augmented Reality (AR) techniques in the flight environment provides assistance in remotely piloted landings, and allows visualizing flight and environmental information that are critical to the mission.

Wind turbines in cold climates are likely to suffer from icing events, deteriorating the aerodynamic performances of the blades and decreasing their power output. Continuous ice accretion causes an increase in the ice mass and, consequently, in the centrifugal force to which the ice shape is subjected. This can result in the shedding of chunks of ice, which can jeopardize the aeroelastic properties of the blade and, most importantly, the safety of the surrounding people and of the wind turbine structure itself. In this work, ice shedding analysis is performed on a quasi-3D, multi-step ice geometry accreted on the NREL 5MW reference wind turbine. A preliminary investigation is performed by including the presence of an ice protection system to decrease the adhesion surface of the ice on the blade. A reference test case with a simple geometry is used as verification for the correct implementation of the procedure.

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.

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.

To improve the combustion and emission characteristics of the large-bore marine engines, the spray is usually designed as an inter-spray impingement to promote the fuel-air mixing process, which implies frequent droplet collisions. Properly describing the collision dynamics of liquid droplets has been of interest in the field of spray modeling for marine engine applications. In this context, this work attempts to develop an accurate and efficient methodology for modeling impinging sprays under engine-like conditions. Experimental validations in terms of spray penetration and morphology are initially carried out at different operating conditions considering the parametric variations of ambient temperature and pressure, where the measurements are performed on a large-scale constant volume chamber with two symmetrical injectors.

The effects of fuel temperature on both the geometry and the droplet size and velocity of a GDI swirled injector spray were investigated by means of visualizations and PDA measurements. Isooctane was used as model fuel and was injected in a quiescent bomb at injection pressure of 7 MPa. Bomb pressure ranged from 40 kPa to 800 kPa with injector nozzle temperature ranging from 293 K to 393 K. A drastic change in spray geometry was observed when conditions above the vaporization curve were reached. The temperature increase has two macroscopic effects on the spray geometry: at the nozzle exit the liquid flash boiling strongly enlarges the spray angle, at a certain distance from the nozzle the air entrainment collapses the spray. Raising the fuel temperature up to flash boiling conditions causes a significant decrease of the average droplet size.

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.