Search Results

In this work a multilevel CFD analysis have been applied for the design of an intake air-box with improved characteristics of noise reduction and fluid dynamic response. The approaches developed and applied for the optimization process range from the 1D to fully 3D CFD simulation, exploring hybrid approaches based on the integration of a 1D model with quasi-3D and 3D tools. In particular, the quasi-3D strategy is exploited to investigate several configurations, tailoring the best trade-off between noise abatement at frequencies below 1000 Hz and optimization of engine performances. Once the best configuration has been defined, the 1D-3D approach has been adopted to confirm the prediction carried out by means of the simplified approach, studying also the impact of the new configuration on the engine performances.

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.

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.

High-speed OH chemiluminescence imaging is used to measure the lift-off length of diesel sprays in an optical heavy-duty diesel engine of 2 L displacement operated at 1200 rpm and 5 bar IMEP. Stereoscopic images are acquired at two different wavelengths (310 and 330 nm). Subtraction of pairwise images helps reducing the background coming from natural soot incandescence in the OH chemiluminescence images. Intake air temperature (343 to 403 K), motored top dead center density (18 to 22 kg/m3), fuel injection pressure (150 to 250 MPa), intake oxygen concentration (17 to 21 %vol) and nozzle diameter (0.1 and 0.14 mm) are varied and a nonlinear regression model is derived from the experimental results to describe stabilized lift-off length as function of the experimental factors. The lift-off length follows the general trends that are known from spray vessel investigations, but the strength of the dependence on certain variables deviates strongly from those studies.

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 influence of jet-flow and jet-jet interactions on the lift-off length of diesel jets are investigated in an optically accessible heavy-duty diesel engine. High-speed OH chemiluminescence imaging technique is employed to capture the transient evolution of the lift-off length up to its stabilization. The engine is operated at 1200 rpm and at a constant load of 5 bar IMEP. Decreasing the inter-jet spacing shortens the liftoff length of the jet. A strong interaction is also observed between the bulk in-cylinder gas temperature and the inter-jet spacing. The in-cylinder swirl level only has a limited influence on the final lift-off length position. Increasing the inter-jet spacing is found to reduce the magnitude of the cycle-to-cycle variations of the lift-off length.

In the last decades numerical simulations have become reliable tools for the design and the optimization of silencers for internal combustion engines. Different approaches, ranging from simple 1D models to detailed 3D models, are nowadays commonly applied in the engine development process, with the aim to predict the acoustic behavior of intake and exhaust systems. However, the acoustic analysis is usually performed under the hypothesis of infinite stiffness of the silencer walls. This assumption, which can be regarded as reasonable for most of the applications, can lose validity if low wall thickness are considered. This consideration is even more significant if the recent trends in the automotive industry are taken into account: in fact, the increasing attention to the weight of the vehicle has lead to a general reduction of the thickness of the metal sheets, due also to the adoption of high-strength steels, making the vibration of the components a non negligible issue.

In this work an integration between a 1D code (Gasdyn) with a CFD code (OpenFOAM®) has been applied to improve the performance of a Moto3TM engine. The four-stroke, single cylinder S.I. engine was modeled, in order to predict the wave motion in the intake and exhaust systems and study how it affects the cylinder gas exchange process. The engine considered was characterized by having an air induction system with integrated filter cartridge, air-box and intake runner, resulting in a complex air-path form the intake mouth to the intake valves, which presents critical aspects when a 1D modeling is addressed. This paper presents a combined and integrated simulation, in which the intake systems was modeled as a 3D geometry whereas the exhaust system, which presented a simpler geometry, was modeled by means of a 1D approach.

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.

The implementation and the combination of advanced boundary conditions and subgrid scale models for Large Eddy Simulation (LES) in the multi-dimensional open-source CFD code OpenFOAM® are presented. The goal is to perform reliable cold flow LES simulations in complex geometries, such as in the cylinders of internal combustion engines. The implementation of a boundary condition for synthetic turbulence generation upstream of the valve port and of the compressible formulation of the Wall-Adapting Local Eddy-viscosity sgs model (WALE) is described. The WALE model is based on the square of the velocity gradient tensor and it accounts for the effects of both the strain and the rotation rate of the smallest resolved turbulent fluctuations and it recovers the proper y₃ near-wall scaling for the eddy viscosity without requiring dynamic procedure; hence, it is supposed to be a very reliable model for ICE simulation.

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.

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.

The present work presents the concept of a new rotary engine, and provides first investigations for its implementation in the energy sector. The main focus of this work is to provide a theoretical description of the engine and its differences from the state-of-the-art technologies. Its innovative principle consists of concentric operation, with two pistons of different rotation radius and the addition of a third intermediate chamber between the compression and combustion chamber. A description of the engine’s physical model is provided, followed by an analysis of the selected specific geometrical features. Additionally, a thermodynamic analysis clarifies the operational advantage compared to the existing cycles and, finally, a numerical investigation on the engine’s bulk performance is provided to quantify the anticipated results of the theoretical analysis.

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.

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.