Search Results

Although the application of cylinder pressure sensors to gain insight into the combustion process is not a novel topic itself, the recent availability of inexpensive in-cylinder pressure sensors has again prompted an upcoming interest for the utilization of the cylinder pressure signal within engine control and monitoring. Besides the use of the in-cylinder pressure signal for combustion analysis and control the information can also be used to determine related quantities in the exhaust or intake manifold. Within this work two different methods to estimate the pressure inside the exhaust manifold are proposed and compared. In contrary to first principle based approaches, which may require time extensive parameterization, alternative data driven approaches were pursued. In the first method a Principle Component Analysis (PCA) is applied to extract the cylinder pressure information and combined with a polynomial model approach.

Modeling combustion of transportation fuels remains a difficult task due to the extremely large number of species constituting commercial gasoline and diesel. However, for this purpose, multi-component surrogate fuel models with a reduced number of key species and dedicated reaction subsets can be used to reproduce the physical and chemical traits of diesel and gasoline, also allowing to perform CFD calculations. Recently, a detailed surrogate fuel kinetic model, named C3 mechanism, was developed by merging high-fidelity sub-mechanisms from different research groups, i.e. C0-C4 chemistry (NUI Galway), linear C6-C7 and iso-octane chemistry (Lawrence Livermore National Laboratory), and monocyclic aromatic hydrocarbons (MAHs) and polycyclic aromatic hydrocarbons (PAHs) (ITV-RWTH Aachen and CRECK modelling Lab-Politecnico di Milano).

In recent years, the increase of gasoline fuel injection pressure is a way to improve thermal efficiency and lower engine-out emissions in GDI homogenous combustion concept. The challenge of controlling particulate formation as well in mass and number concentrations imposed by emissions regulations can be pursued improving the mixture preparation process and avoiding mixture inhomogeneity with ultra-high injection pressure values up to 100 MPa. The increase of the fuel injection pressure in GDI homogeneous systems meets the demand for increased injector static flow, while simultaneously improves the spray atomization and mixing characteristics with consequent better combustion performance. Few studies quantify the effects of high injection pressure on transient gasoline spray evolution. The aim of this work was to simulate with OpenFOAM the spray morphology of a commercial gasoline injected in a constant volume vessel by a prototypal GDI injector.

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.

Today, Direct-Injection systems are widely used on Spark-Ignition engines in combination with turbo-charging to reduce the fuel-consumption and the knock risks. In particular, the spread of Gasoline Direct Injection (GDI) systems is mainly related to the use of new generations of multi-hole, high-pressure injectors whose characteristics are quite different with respect to the hollow-cone, low-pressure injectors adopted in the last decade. This paper presents the results of an experimental campaign conducted on the spray produced by a GDI six-holes injector into a constant volume vessel with optical access. The vessel was filled with air at atmospheric pressure. Different operating conditions were considered for an injection pressure ranging from 3 to 20 MPa. For each operating condition, spray images were acquired by a CCD camera and then post processed to evaluate the spray penetration and cone angles.

In this paper a three-dimensional time-domain CFD approach has been employed to predict and analyze the acoustic attenuation performance of complex perforated muffler geometries, where strong 3D effects limit the validity of the use of one-dimensional models. A pressure pulse has been imposed at the inlet to excite the wave motion, while unsteady flow computation have been performed to acquire the time histories of the pressures upstream and downstream of the silencer. Pressures in the time domain have been then transformed to acoustic pressures in the frequency domain, to predict the transmission loss.

Partially premixed combustion (PPC) can be applied to decrease emissions and increase fuel efficiency in direct injection, compression ignition (DICI) combustion engines. PPC is strongly influenced by the mixing of fuel and oxidizer, which for a given fuel is controlled mainly by (a) the fuel injection, (b) the in-cylinder flow, and (c) the geometry and dynamics of the engine. As the injection timings can vary over a wide range in PPC combustion, detailed knowledge of the in-cylinder flow over the whole intake and compression strokes can improve our understanding of PPC combustion. In computational fluid dynamics (CFD) the in-cylinder flow is sometimes simplified and modeled as a solid-body rotation profile at some time prior to injection to produce a realistic flow field at the moment of injection. In real engines, the in-cylinder flow motion is governed by the intake manifold, the valve motion, and the engine geometry.

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.

The current development to set up an automatic procedure for automatic mesh generation and automatic mesh motion for internal combustion engine simulation in OpenFOAM®-2.2.x is here described. In order to automatically generate high-quality meshes of cylinder geometries, some technical issues need to be addressed: 1) automatic mesh generation should be able to control anisotropy and directionality of the grid; 2) during piston and valve motion, cells and faces must be introduced and removed without varying the overall area and volume of the cells, to avoid conservation errors. In particular, interpolation between discrete fields is frequent in computational physics: the use of adaptive and non-conformal meshes necessitates the interpolation of fields between different mesh regions. Interpolation problems also arise in areas such as model coupling, model initialization and visualisation.

In this paper, computation fluid dynamics (CFD) simulations are employed to describe the effect of flow parameters on the formation of soot and NOx in a heavy duty engine under low load and high load. The complexity of diesel combustion, specially when soot, NOx and other emissions are of interest, requires using a detailed chemical mechanism to have a correct estimation of temperature and species distribution. In this work, Multiple Representative Interactive Flamelets (MRIF) method is employed to describe the chemical reactions, ignition, flame propagation and emissions in the engine. A phenomenological model for soot formation, including soot nucleation, coagulation and oxidation with O2 and OH is incorporated into the flamelet combustion model. Different strategies for modelling NOx are chosen to take into account the longer time scale for NOx formation. The numerical results are compared with experimental data to show the validity of the model for the cases under study.

The influence of spray-wall interaction on air entrainment in an unsteady non-evaporating diesel spray was studied using laser Doppler anemometry. The spray was injected into confined quiescent air at ambient pressure and temperature and made to impact on a flat wall. The air velocity component normal to a cylindrical surface surrounding the spray was measured during the entire injection period, allowing to evaluate the time history of the entrained air mass flow rate. The influence of wall distance and spray impingement angle on air entrainment characteristics has been investigated and the results indicate that the presence of a wall increases the entrained mass flow rate in the region close to the surface, during the main injection period. Normal impingement appears to produce stronger effects than oblique incidence at 30 and 45 deg. A qualitative explanation of the results is also proposed, based on the drop-gas momentum exchange mechanism.

A comprehensive model for sprays emerging from high pressure swirl injectors for GDI engine application has been developed. The primary and secondary atomization mechanism as well as the evaporation process both in standard and superheated conditions are taken into account. The spray modelling after the injection is based on the Liquid Instability Sheet Atomization (LISA) approach, modified to correctly predict the liquid sheet thickness at the breakup length. The effect of different values of the superheat degree on evaporation and impact on the spray distribution and fuel-air mixing is analyzed. Comparisons with experimental data show good agreements under atmospheric conditions and with different superheated degrees, while some discrepancies occur under higher ambient pressures.

Aim of this work is to present and discuss the possibility and the limits of two zone models for spark-ignition engines using a detailed kinetic scheme for the characterization of the evolution of the air-fuel mixture, while an equilibrium approach is used for the burnt zone. Simple experimental measurements of knocking tendency of different fuels in ideal reactors, such as rapid compression machines and shock tube reactors, cannot be directly used for the analysis of octane numbers and sensitivity of hydrocarbon mixtures. Thus a careful investigation is very useful, not only of the combustion chamber behavior, including the modelling of the turbulent flame front propagation, but also of the fluid dynamic behavior of the intake and exhaust system, accounting for the volumetric efficiency of the engine.

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.

Oxygen enriched air (EA) is a well known industrial mixture in which the content of oxygen is higher respect the atmospheric one, in the range 22-35%. Oxygen EA can be obtained by desorption from water, taking advantage of the higher oxygen solubility in water compared to the nitrogen one, since the Henry constants of this two gases are different. The production of EA by this new approach was already studied by experimental runs and theoretical considerations. New results using salt water are reported. EA promoted combustion is considered as one of the most interesting technologies to improve the performance in diesel engines and to simultaneously control and reduce pollution. This paper explores, by means of 3-dimensional computational fluid dynamics simulations, the effects of EA on the performance and exhaust emissions of a high-speed direct-injection diesel engine.

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 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 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.

In this study two different simulation approaches to large eddy simulation of spark-ignition engines are compared. Additionally, some of the simulation results are compared to experimentally obtained in-cylinder velocity measurements. The first approach applies unstructured grids with an automated meshing procedure, using OpenFoam and Lib-ICE with a mapping approach. The second approach applies the efficient in-house code PsiPhi on equidistant, Cartesian grids, representing walls by immersed boundaries, where the moving piston and valves are described as topologically connected groups of Lagrangian particles. In the experiments, two-dimensional two-component particle image velocimetry is applied in the central tumble plane of the cylinder of an optically accessible engine. Good agreement between numerical results and experiment are obtained by both approaches.

For passenger cars, individual tyre model parameters, used in vehicle models able to simulate vehicle handling behavior, are traditionally derived from expensive component indoor laboratory tests as a result of an identification procedure minimizing the error with respect to force and slip measurements. Indoor experiments on agricultural tyres are instead more challenging and thus generally not performed due to tyre size and applied forces. However, the knowledge of their handling characteristics is becoming more and more important since in the next few years, all agricultural vehicles are expected to run on ordinary asphalt roads at a speed of 80km/h. The present paper presents a methodology to identify agricultural tyres' handling characteristics based only on the measurements carried out on board vehicle (vehicle sideslip angle, yaw rate, lateral acceleration, speed and steer angle) during standard handling maneuvers (step-steers, J-turns, etc.), instead than during indoor tests.