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.

Methane abatement in the exhaust gas of natural gas engines is much more challenging in respect to the oxidation of other higher order hydrocarbons. Under steady state λ sweep, the methane conversion efficiency is high at exact stoichiometric, and decreases steeply under both slightly rich and slightly lean conditions. Transient lean to rich transitions can improve methane conversion at the rich side. Previous experimental work has attributed the enhanced methane conversion to activation of methane steam reforming. The steam reforming rate, however, attenuates over time and the methane conversion rate gradually converges to the low steady state values. In this work, a reactor model is established to predict steady state and transient transition characteristics of a three-way catalyst (TWC) mounted in the exhaust of a natural gas heavy-duty engine.

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.

The sustainability of Martian outposts development is strongly based on the capability of achieving a high level of autonomy both in terms of operations management and of resources availability. In situ production of consumables is a key point to allow humans to work and live on Mars avoiding or limiting the need for re-supplies of materials from Earth. Required consumables can be produced in situ exploiting the locally available resources, but also by means of green-houses and waste recycle systems. Dedicated robotic missions for in situ demonstration of this type of technologies are a fundamental step of the Martian In Situ Resources Utilization (ISRU) development roadmap. This paper is focused on the extraction of oxygen and fuels (e.g. methane) from the Martian atmosphere, and presents a feasibility study for a small technological demonstration plant.

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 scope of the work presented in this paper was to apply the latest open source CFD achievements to design a state of the art, direct-injection (DI), heavy-duty, natural gas-fueled engine. Within this context, an initial steady-state analysis of the in-cylinder flow was performed by simulating three different intake ducts geometries, each one with seven different valve lift values, chosen according to an estabilished methodology proposed by AVL. The discharge coefficient (Cd) and the Tumble Ratio (TR) were calculated in each case, and an optimal intake ports geometry configuration was assessed in terms of a compromise between the desired intensity of tumble in the chamber and the satisfaction of an adequate value of Cd. Subsequently, full-cycle, cold-flow simulations were performed for three different engine operating points, in order to evaluate the in-cylinder development of TR and turbulent kinetic energy (TKE) under transient conditions.

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.

This paper presents results of an experimental investigation on a flexible port dual fuel injection using different ethanol to gasoline mass fractions. A four stroke, two cylinder turbocharged SI engine was used for the experiments. The engine speed was set at 3000 rpm, tests were carried out at medium-high load and two air-fuel-ratio. The initial reference conditions were set running the engine, fueled with full gasoline at the KLSA boundary, in accordance with the standard ECU engine map. This engine point was representative of a rich mixture (λ=0.9) in order to control the knock and the temperature at turbine inlet. The investigated fuels included different ethanol-gasoline mass fractions (E10, E20, E30 and E85), supplied by dual injection within the intake manifold. A spark timing sweep, both at stoichiometric and lean (λ=1.1) conditions, up to the most advanced one without knock was carried out.

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.

In the recent years, the interest in heavy-duty engines fueled with Compressed Natural Gas (CNG) is increasing due to the necessity to comply with the stringent CO2 limitation imposed by national and international regulations. Indeed, the reduced number of carbon atoms of the NG molecule allows to reduce the CO2 emissions compared to a conventional fuel. The possibility to produce synthetic methane from renewable energy sources, or bio-methane from agricultural biomass and/or animal waste, contributes to support the switch from conventional liquid fuels to CNG. To drive the engine development and reduce the time-to-market, the employment of numerical analysis is mandatory. This requires a continuous improvement of the simulation models toward real predictive analyses able to reduce the experimental R&D efforts. In this framework, 1D numerical codes are fundamental tools for system design, energy management optimization, and so on.