Search Results

A numerical investigation on a series of Diesel spray experiments in constant-volume vessels is proposed. Non reacting conditions were used to assess the spray models and to determine the grid size required to correctly predict the fuel-air mixture formation process. To this end, not only computed liquid and vapor penetrations were compared with experimental data, but also a detailed comparison between computed and experimental mixture fraction distributions was performed at different distances from the injector. Grid dependency was reduced by introducing an Adaptive Local Mesh Refinement technique (ALMR) with an arbitrary level of refinement. Once the capabilities of the current implemented spray models have been assessed, reacting conditions at different ambient densities and temperatures were considered. A Perfectly Stirred Reactor (PSR) combustion model, based on a direct integration of complex chemistry mechanisms over a homogenous cell, was adopted.

Detailed chemistry and turbulence-chemistry interaction need to be properly taken into account for a realistic combustion simulation of IC engines where advanced combustion modes, multiple injections and stratified combustion involve a wide range of combustion regimes and require a proper description of several phenomena such as auto-ignition, flame stabilization, diffusive combustion and lean premixed flame propagation. To this end, different approaches are applied and the most used ones rely on the well-stirred reactor or flamelet assumption. However, well-mixed models do not describe correctly flame structure, while unsteady flamelet models cannot easily predict premixed flame propagation and triple flames. A possible alternative for them is represented by transported probability density functions (PDF) methods, which have been applied widely and effectively for modeling turbulent reacting flows under a wide range of combustion regimes.

A new multi-zone model for the simulation of HCCI engine is here presented. The model includes laminar and turbulent diffusion and conduction exchange between the zones and the last improvements on the numerical aspects. Furthermore, a new strategy for the zone discretization is presented, which allows a better description of the near-wall zones. The aim of the work is to provide a fast and reliable model for carrying out chemical analysis with detailed kinetic schemes. A preliminary sensitivity analysis allows to verify that 10 zones are a convenient number for a good compromise between the computational effort and the description accuracy. The multi-zone predictions are then compared with the CFD ones to find the effective turbulence parameters, with the aim to describe the near-wall phenomena, both in a reactive and non-reactive cases.

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.

This paper presents calculations of engine flows by using the Partially-Averaged Navier Stokes (PANS) method (Girimaji [1]; [2]). The PANS is a scale-resolving turbulence computational approach designed to resolve large scale fluctuations and model the remainder with appropriate closures. Depending upon the prescribed cut-off length (filter width) the method adjusts seamlessly from the Reynolds-Averaged Navier-Stokes (RANS) to the Direct Numerical Solution (DNS) of the Navier-Stokes equations. The PANS method was successfully used for many applications but mainly on static geometries, e.g. Basara et al. [3]; [4]. This is due to the calculation of the cut-off control parameter which requires that the resolved kinetic energy is known and this is usually obtained by suitably averaging of the resolved field. Such averaging process is expensive and impractical for engines as it would require averaging per cycles.

Reynolds-averaged Navier-Stokes (RANS) computations of heat transfer involving wall bounded flows at elevated Prandtl numbers typically suffer from a lack of accuracy and/or increased mesh dependency. This can be often attributed to an improper near-wall turbulence modeling and the deficiency of the wall heat transfer models (based on the so called P-functions) that do not properly account for the variation of the turbulent Prandtl number in the wall proximity (y+< 5). As the conductive sub-layer gets significantly thinner than the viscous velocity sub-layer (for Pr >1), treatment of the thermal buffer layer gains importance as well. Various hybrid strategies utilize blending functions dependent on the molecular Prandtl number, which do not necessarily provide a smooth transition from the viscous/conductive sub-layer to the logarithmic region.

Metallic foams or sponges are materials with a cell structure suitable for many industrial applications, such as reformers, heat catalytic converters, etc. The success of these materials is due to the combination of various characteristics such as mechanical strength, low density, high specific surface, good thermal exchange properties, low flow resistance and sound absorption. Different materials and manufacturing processes produce different type of structure and properties for various applications. In this work a genetic algorithm has been developed and applied to support the design of catalytic devices. In particular, two substrates were considered, namely the traditional honeycomb and an alternative open-cell foam type. CFD simulations of pressure losses and literature based correlations for the heat and mass transfer were used to support the genetic algorithm in finding the best compromise between flow resistance and pollutant abatement.

Many studies have been carried out to optimize the aerodynamic performances of a single car or a single vehicle. In present days the traffic increases and sophisticated technologies are developing to guarantee the drivers safety, to minimize the fuel consumption and be more environmentally friendly. Within this research area a new technique that is being studied is Platooning: this means that different vehicles travel in a configuration that minimizes the aerodynamic drag and therefore the fuel consumption and the longitudinal space. In the present study platoons with different vehicles and configurations are taken into account, to analyze the influence of car shape and relative distance between the vehicles. The research has been carried out using CFD techniques to investigate the different flow fields around different platoons, while wind tunnel tests have been used to validate the results of the CFD simulations.

In the paper, a procedure to assess the quality of the shielding effect of wind barriers to protect large sided vehicles crossing the wake of a bridge pylon under cross wind conditions is proposed. The methodology is based on Multi-Body simulations of the response of the vehicle-driver system when it is subjected to the sudden change of the aerodynamic forces due to the wind-tower interaction. The aerodynamic forces that are instantaneously acting on the vehicle are computed according to a force distribution approach that relies on wind tunnel tests that may be performed on still scaled models. From the knowledge of the aerodynamic force distribution along the vehicle at different yaw angles and of the mean wind profile across the tower wake, the aerodynamic force, acting on the moving vehicle, is reconstructed at each time step taking into consideration the actual vehicle-driver dynamics.

The presented paper deals with a methodology to model cycle-to-cycle variations (CCV) in 0-D/1-D simulation tools. This is achieved by introducing perturbations of combustion model parameters. To enable that, crank angle resolved data of individual cycles (pressure traces) have to be available for a reasonable number of engine cycles. Either experimental data or 3-D CFD results can be applied. In the presented work, experimental data of a single-cylinder research engine were considered while predicted LES 3-D CFD results will be tested in the future. Different engine operating points were selected - both stable ones (low CCV) and unstable ones (high CCV). The proposed methodology consists of two major steps. First, individual cycle data have to be matched with the 0-D/1-D model, i.e., combustion model parameters are varied to achieve the best possible match of pressure traces - an automated optimization approach is applied to achieve that.

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.

Objective of this work is the incorporation of the flame stretch effects in an Eulerian-Lagrangian model for premixed SI combustion in order to describe ignition and flame propagation under highly inhomogeneous flow conditions. To this end, effects of energy transfer from electrical circuit and turbulent flame propagation were fully decoupled. The first ones are taken into account by Lagrangian particles whose main purpose is to generate an initial burned field in the computational domain. Turbulent flame development is instead considered only in the Eulerian gas phase for a better description of the local flow effects. To improve the model predictive capabilities, flame stretch effects were introduced in the turbulent combustion model by using formulations coming from the asymptotic theory and recently verified by means of DNS studies. Experiments carried out at Michigan Tech University in a pressurized, constant-volume vessel were used to validate the proposed approach.

The OpenFOAM® CFD methodology is nowadays employed for simulation in internal combustion engines and a lot of work has been done for an appropriate description of all complex phenomena. At the moment in the RANS turbulence models available in the OpenFOAM® toolbox the turbulence modulation is not yet included, and the present work analyzes the predictive capabilities of the code in simulating high injection pressure fuel sprays after modeling the influence of the dispersed phase on the turbulence structure. Different experiments were employed for the validation. At first, non-evaporating diesel spray was considered in a constant volume and quiescent vessel. The validation was performed via the available experimental spray evolution in terms of penetrations and spatial/temporal fuel distributions. Then the Sandia combustion chamber was chosen for diesel spray simulation in non-reacting conditions.

Wet clutches drag loss simulation is essentially linked to the clutch friction surface patterns in addition to the main geometry and conditions of the interface (relative speed, separation, inner and outer radius, viscosity and boundary pressures). The clutch patterns promote cooling flow and micro-hydrodynamic effects to aid clutch separation but greatly complicate the simulation of drag loss during separation. These drag losses are important in understanding the system losses as well as finding the most effective clutch cooling strategy. Typical clutch models either only consider simple patterns, such as radial grooves, or require significant simulation efforts to evaluate. Additionally, many simple models require calibration to measurement of the actual clutch they try to model before they provide a useful model.

The s.i. combustion process and its corresponding pollutant formation are investigated by means of a quasiD approach and a CFD model. This work has been motivated by the need to better understand the reliability of such models and to assess their accuracies with respect to the prediction of engine performances and emissions. An extended dissertation about the fundamental mechanisms governing the pollutant formation in the turbulent premixed combustion which characterizes the s.i. engines is given. The conclusion of such analysis is the definition of a new reduced chemical scheme, based on the application of partial-equilibrium and steady-state assumptions for the radicals and the solution of a transport equation for each specie which is kinetically controlled. For this purpose the CFD code OpenFOAM [1, 2, 3] and the thermo-fluid dynamic code GASDYN [4, 5] have been applied and enhanced.

Aerodynamic properties of a BMW car model taking over a truck model are studied computationally by applying the scale-resolving PANS (Partially-averaged Navier-Stokes) approach. Both vehicles represent down-scaled (1:2.5), geometrically-similar models of realistic vehicle configurations for which on-road measurements have been performed by Schrefl (2008). The operating conditions of the modelled ‘on-road’ overtaking maneuver are determined by applying the dynamic similarity concept in terms of Reynolds number consistency. The simulated vehicle configuration constitutes of a non-moving truck model and a car model moving against the air flow, the velocity of which corresponds to the car velocity.

Two thread forming processes, rolling and cutting, were studied for their effects on fatigue in cast aluminum 319-T7. Material was excised from cylinder blocks and tested in rotating-bending fatigue in the form of unnotched and notched specimens. The notched specimens were prepared by either rolling or cutting to replicate threads in production-intent parts. Cut threads exhibited conventional notch behavior for notch sensitive materials. In contrast, plastic deformation induced by rolling created residual compressive stresses in the notch root and significantly improved fatigue strength to the point that most of the rolled specimens broke outside the notch. Fractographic and metallographic investigation showed that cracks at the root of rolled notches were deflected upon initiation. This lengthened their incubation period, which effectively increased fatigue resistance.

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 paper presents calculations of external car aerodynamics by using the finite volume (FV) immersed-boundary method. The FV numerical codes primarily employ Reynolds-Averaged Navier-Stokes (RANS) models. In recent years, and due to possibility to run very large computational meshes, these models are usually used in conjunction with the advanced near-wall models. Moreover, it has been often demonstrated that the accuracy of RANS near-wall models relies on the mesh quality near the wall so by the rule, larger number of wall body-fitted cell-layers are employed. An immersed boundary (IB) method becomes an attractive alternative to the ‘standard’ FV approaches especially when applied to low quality CAD data. In general, the IB method is less investigated and validated for the car aerodynamics, particularly in conjunction with advanced near-wall turbulence models and an adaptive mesh refinement (AMR).