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.

The interaction of turbulent premixed methane combustion with the surrounding flow field can be studied using optically accessible test rigs such as a rapid compression expansion machine (RCEM). The high flexibility offered by such a test rig allows its operation at various thermochemical conditions at ignition. However, limitations inherent to such test rigs due to the absence of an intake stroke do not allow turbulence production as found in IC-engines. Hence, means to introduce turbulence need to be implemented and the relevant turbulence quantities have to be identified in order to enable comparability with engine relevant conditions. A dedicated high-pressure direct injection of air at the beginning of the compression phase is considered as a measure to generate adjustable turbulence intensities at spark timing and during the early flame propagation.

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.

In this paper, an experimental and numerical analysis of combustion process and knock occurrence in a small displacement spark-ignition engine is presented. A wide experimental campaign is preliminarily carried out in order to fully characterize the engine behavior in different operating conditions. In particular, the acquisition of a large number of consecutive pressure cycle is realized to analyze the Cyclic Variability (CV) effects in terms of Indicated Mean Effective Pressure (IMEP) Coefficient of Variation (CoV). The spark advance is also changed up to incipient knocking conditions, basing on a proper definition of a knock index. The latter is estimated through the decomposition and the FFT analysis of the instantaneous pressure cycles. Contemporary, a quasi-dimensional combustion and knock model, included within a whole engine one-dimensional (1D) modeling framework, are developed. Combustion and knock models are extended to include the CV effects, too.

The fluid-mechanic behaviour of straight-sided and re-entrant chamber geometries has been studied using laser doppler anemometry (LDA) technique. Measurements have been carried out during the compression stroke in a direct injection diesel engine, representative of medium size family, operating at 1000 rpm under motored conditions. The mean motion and turbulence intensity have been computed using a filtering procedure on the LDA data. Using the second version of KIVA code, the air flow field evolution during the same crank angle period has been also computed. To perform proper comparisons between measured and computed values of mean velocity and turbulence intensity, a careful choice of the initial conditions for computations has been performed. Reasonable agreement has been found between computed and measured mean swirl velocities for both combustion chamber geometries tested. On the contrary, the computed turbulence intensities underestimate those measured.

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.

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.

In-cylinder cycle-resolved LDV measurements have been made in a diesel engine having a high-squish re-entrant combustion chamber with compression ratio of 21:1. The engine has been motored in the range of 1000 to 3000 rpm thanks to the use of self-lubricating seeding particles. Conventional ensemble-averaging and filtering techniques have been used for analyzing instantaneous velocity data obtained at two points along a diameter located in a horizontal plane at 5 mm below the engine head. The dependence of the mean motion and turbulence on engine speed has been evaluated. The effect of cut-off frequency selection on turbulence values has been also analyzed. Moreover, the Kolmogorov's -5/3 power domain has been investigated in detail by spectral analysis on the instantaneous velocity data.

In-cylinder measurements of turbulent integral length scales, carried out during the last 60 degrees of the compression stroke at 600 and 1,000 rpm by a two-probe volume LDV system, were used to assess the capability of the k-ε model used in KIVA-II code. The objective of the paper is to address the following question: what is the most reasonable definition of turbulent length scale in the k-ε model for engine applications? The answer derived from the comparison between KIVA predictions and experiments that showed a fair agreement between the computed turbulent length scale and the measured lateral integral length scale. The agreement is a result of proper choice of the initial swirl ratio and turbulent kinetic energy at inlet valve closure (IVC) by taking into account the LDV measurements and the value of the constant Cμε in the k-ε model equations that relates the turbulent length scale to k and ε.

Tangential component of velocity and turbulence were measured in three locations in the re-entrant combustion chamber of a motored single-cylinder d.i. Diesel engine (0.435 liter, 21:1 compression ratio) using a Laser Doppler Velocimetry system. Moreover, a modified LDV system with two-probe volume was used to measure directly lateral integral length scales of the velocity tangential component at two engine speeds. The measurements were made on a horizontal plane at 5 mm below the engine head from 100 degrees before TDC to 60 degrees after TDC of both the compression and expansion strokes. The engine was motored at 1,000 and 1,500 rpm respectively. An ensemble-averaging technique was performed to analyze the instantaneous velocity information supplied by two Burst Spectrum Analyzers. The lateral integral length scale was obtained from the integral of the spatial correlation coefficient of the velocity fluctuation for different separation.

The present paper aims at developing a general method to estimate integral and microtime scales of turbulent in-cylinder flow field in reciprocating engines. The ensemble average technique was used to compute the integral time scale from the single point time autocorrelation function, whereas the microtime scale, representative of the most rapid changes that occur in the fluctuation, was computed as the intercept of the parabola that matches the autocorrelation function at the origin. Further, the microtime scale was also estimated by spectral analysis through the energy spectral density function of the ensemble turbulent fluctuation and the results obtained by the two methods were compared. The procedures were applied to the tangential component of the instantaneous velocity data collected, at different engine speeds (1,000, 1,500, 2,000 rpm), within a motored d.i. diesel engine equipped with a re-entrant combustion chamber, using the Laser Doppler Anemometry (LDA) technique.

Ducted Fuel Injection (DFI) has the potential to reduce soot emissions in Diesel engines thanks to the enhanced mixing rate resulting from the liquid fuel flow through a small cylindrical pipe located at a certain distance from the nozzle injector hole. A consolidated set of experiments in constant-volume vessel and engine allowed to understand the effects of ambient conditions, duct geometry and shape on fuel-air mixing, combustion and soot formation. However, implementation of this promising technology in compression-ignition engines requires predictive numerical models that can properly support the design of combustion systems in a wide range of operating conditions. This work presents a computational methodology to predict fuel-air mixing and combustion with ducted fuel injection. Attention is mainly focused on turbulence and combustion modelling.

Turbulence modeling for fuel spray simulation plays a prominent role in the understanding of the flow behavior in Internal Combustion Engines (ICEs). Currently, a lot of research work is actively spent on Large Eddy Simulation (LES) turbulence modeling as a replacement option of standard Reynolds averaged approaches in the Eulerian-Lagrangian spray modeling framework, due to its capability to accurately describe flow-induced spray variability and to the lower dependence of the results on the specific turbulence model and/or modeling coefficients. The introduction of LES poses, however, additional questions related to the implementation/adaptation of spray-related turbulence sources and to the rise of conflicting numerics and grid requirements between the Lagrangian and Eulerian parts of the simulated flow.

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.

The definition of a robust methodology to perform a full-cycle CFD simulation of IC engines requires as first step the availability of a reliable grid generation tool, which does not only have to guarantee a high quality mesh but also has to prove to be efficient in terms of required time. In this work the authors discuss a novel approach entirely based on the OpenFOAM technology, in which the available 3D grid generator was employed to automatically create meshes containing hexahedra and split-hexahedra from triangulated surface geometries in Stereolithography (STL) format. The possibility to introduce local refinements and boundary layers makes this tool suitable for IC engine simulations. Grids are sequentially generated at target crank angles which are automatically determined depending on user specified settings such as maximum mesh validity interval and quality parameters like non-orthogonality, skewness and aspect ratio.

Today, multi-hole Diesel injectors can be mainly characterized by three different nozzle hole shapes: cylindrical, k-hole, and ks-hole. The nozzle hole layout plays a direct influence on the injector internal flow field characteristics and, in particular, on the cavitation and turbulence evolution over the hole length. In turn, the changes on the injector internal flow correlated to the nozzle shape produce immediate effects on the emerging spray. In the present paper, the fluid dynamic performance of three different Diesel nozzle hole shapes are evaluated: cylindrical, k-hole, and ks-hole. The ks-hole geometry was experimentally characterized in order to find out its real internal shape. First, the three nozzle shapes were studied by a fully transient CFD multiphase simulation to understand their differences in the internal flow field evolutions. In detail, the attention was focused on the turbulence and cavitation levels at hole exit.