Search Results

A simulation method is presented for the analysis of combustion in spark ignition (SI) engines operated at elevated exhaust gas recirculation (EGR) level and employing multiple spark plug technology. The modeling is based on a zero-dimensional (0D) stochastic reactor model for SI engines (SI-SRM). The model is built on a probability density function (PDF) approach for turbulent reactive flows that enables for detailed chemistry consideration. Calculations were carried out for one, two, and three spark plugs. Capability of the SI-SRM to simulate engines with multiple spark plug (multiple ignitions) systems has been verified by comparison to the results from a three-dimensional (3D) computational fluid dynamics (CFD) model. Numerical simulations were carried for part load operating points with 12.5%, 20%, and 25% of EGR. At high load, the engine was operated at knock limit with 0%, and 20% of EGR and different inlet valve closure timing.

In scooter/motorbike engines coherent and stable tumble motion generation is still considered an effective mean in order to both reduce engine emissions and promote higher levels of combustion efficiency. The scientific research also assessed that squish motion is an effective mean for speeding up the combustion in a combustion process already fast. In a previous technical paper the authors demonstrated that for an engine having a high C/D ratio the squish motion is not only not necessary but also detrimental for the stability of the tumble motion itself, because there is a strong interaction between these two motions with the consequent formation of secondary vortices, which in turn penalizes the tumble breakdown and the turbulent kinetic energy production.

A computer program which rapidly calculates the equilibrium mole fractions and the partial derivatives of the mole fractions with respect to temperature, pressure and equivalence ratio for the products of combustion of any hydrocarbon fuel and air is described. A subroutine is also given which calculates the gas constant, enthalpy, internal energy and the partial derivatives of these with respect to temperature, pressure and equivalence ratio. Some examples of the uses of the programs are also given.

This paper presents a methodology for the 3D CFD simulation of the intake and compression processes of four stroke internal combustion engines.The main feature of this approach is to provide very accurate initial conditions by means of a cost-effective initialization step. Calculations are applied to a low stroke-to-bore SI engine, operated at full load and maximum engine speed. It is demonstrated that initial conditions for this kind of engines have an important influence on flow field development, particularly in terms of mean velocities close to the firing TDC. Simulation results are used to discuss the choice of a set of parameters for the flow field characterization of low stroke-to-bore engines, as well as to provide an insight into the flow patterns during the overlapping period.

Nowadays, an automotive DI Diesel engine is demanded to provide an adequate power output together with limit-complying NOx and soot emissions so that the development of a specific combustion concept is the result of a trade-off between conflicting objectives. In other words, the development of a low-emission DI diesel combustion concept could be mathematically represented as a multi-objective optimization problem. In recent years, genetic algorithm and CFD simulations were successfully applied to this kind of problem. However, combining GA optimization with actual CFD-3D combustion simulations can be too onerous since a large number of simulations is usually required, resulting in a high computational cost and, thus, limiting the suitability of this method for industrial processes.

Fuels are formulated by a variety of different components characterized by chemical and physical properties spanning a wide range of values. Changing the ratio between the mixture component molar fractions, it is possible to fulfill different requirements. One of the main properties that can be strongly affected by mixture composition is the volatility that represents the fuel tendency to vaporize. For example, changing the mixture ratio between alcohols and hydrocarbons, it is possible to vary the mixture saturation pressure, therefore the fuel vaporization ratio during the injection process. This paper presents a 1D numerical model to simulate the superheated injection process of a gasoline-ethanol mixture through real nozzle geometries. In order to test the influence of the mixture properties on flash atomization and flash evaporation, the simulation is repeated for different mixtures characterized by different gasoline-ethanol ratio.

As a consequence of reduced fuel consumption, direct injection gasoline engines have already prevailed against port fuel injection. However, in-cylinder fuel homogenization strongly depends on charge motion and injection strategies and can be challenging due to the reduced available time for mixture formation. An insufficient homogenization has generally a negative impact on the combustion and therefore also on efficiency and emissions. In order to reach the targets of the intensified CO2 emission reduction, further increase in efficiency of SI engines is essential. In this connection, 0D/1D simulation is a fundamental tool due to its application area in an early stage of development and its relatively low computational costs. Certainly, inhomogeneities are still not considered in quasi dimensional combustion models because the prediction of mixture formation is not included in the state of the art 0D/1D simulation.

Knocking combustions heavily limits the efficiency of Spark Ignition engines. The compression ratio is limited in the design stage of the engine development, letting to Spark Advance control the task of reducing the odds of abnormal combustions. A detailed analysis of knocking events can help improving engine performance and diagnosis strategies. An effective way is to use advanced 3D CFD (Computational Fluid Dynamics) simulation for the analysis and prediction of combustion performance. Standard 3D CFD approach is based on RANS (Reynolds Averaged Navier Stokes) equations and allows the analysis of the mean engine cycle. However knocking phenomenon is not deterministic and it is heavily affected by the cycle to cycle variation of engine combustions. A methodology for the evaluation of the effects of CCV (Cycle by Cycle Variability) on knocking combustions is here presented, based on both the use of Computation Fluid Dynamics (CFD) tools and experimental information.

Fine-wire resistance thermometers were used to measure compression gas temperatures in a motoring (nonfiring) cycle CFR engine. Temperature versus crankangle curves were obtained for the compression and expansion strokes by means of tungsten wires ranging in diameter from 0.15–1.00 mils and at speeds from 600–1800 rpm. The results were compared with the infrared pyrometer; the peak temperature and peak crankangle lags were determined as a function of the wire diameter and engine speed. Attempts to evaluate the instantaneous energy balance around the wire resulted in a negative heat transfer coefficient, for which no current satisfactory explanation is available, although other observers have reported similar phenomena. The tungsten resistance thermometer is simple to build, easy to install, and requires no modification of the engine block for use during motoring. Thus, it is suitable for comparing the compression temperatures of different design engines.

Cycle-to-cycle variation in combustion phasing and combustion rate cause knock to occur differently in every cycle. This is found to be true even if the end gas thermo-chemical time history is the same. Three cycles are shown that have matched combustion phasing, combustion rate, and time of knock onset, but have knock intensity that differs by a factor of six. Thus, the prediction of knock intensity must include a stochastic component. It is shown that there is a relationship between the maximum possible knock intensity and the unburned fuel energy at the time of knock onset. Further, for a small window of unburned energy at knock onset, the probability density function of knock intensity is self similar when scaled by the 95th percentile of the cumulative distribution, and log-normal in shape.

Modeling ignition kernel development in spark ignition engines is crucial to capturing the sources of cyclic variability, both with RANS and LES simulations. Appropriate kernel modeling must ensure that energy transfer from the electrodes to the gas phase has the correct timing, rate and locations, until the flame surface is large enough to be represented on the mesh by the G-Equation level-set method. However, in most kernel models, geometric details driving kernel growth are missing: either because it is described as Lagrangian particles, or because its development is simplified, i.e., down to multiple spherical flames. This paper covers the geometric aspects of kernel development, which makes up the core of a Triangulated Lagrangian Ignition Kernel model. One (or multiple, if it restrikes) spark channel is initialized as a one-dimensional Lagrangian particle thread.

The long-term goal of this work is to develop a conceptual model for multiple injections of diesel jets. The current work contributes to that effort by performing a detailed modeling investigation into mechanisms that are predicted to control 1st and 2nd stage ignition in single-pulse diesel (n-dodecane) jets under different conditions. One condition produces a jet with negative ignition dwell that is dominated by mixing-controlled heat release, and the other, a jet with positive ignition dwell and dominated by premixed heat release. During 1st stage ignition, fuel is predicted to burn similarly under both conditions; far upstream, gases at the radial-edge of the jet, where gas temperatures are hotter, partially react and reactions continue as gases flow downstream. Once beyond the point of complete fuel evaporation, near-axis gases are no longer cooled by the evaporation process and 1st stage ignition transitions to 2nd stage ignition.

In the near future, emissions legislation will become more and more restrictive for direct injection SI engines by adopting a stringent limitation of particulate number emissions in late 2017. In order to cope with the combustion system related challenges coming along with the introduction of this new standard, Hitachi Automotive Systems Ltd., Hitachi Europe GmbH and IAV GmbH work collaboratively on demonstrating technology that allows to satisfy EU6c emissions limitations by application of Hitachi components dedicated to high pressure injection (1). This paper sets out to describe both the capabilities of a new high pressure fuel system improving droplet atomization and consequently mixture homogeneity as well as the process of utilizing the technology during the development of a demonstrator vehicle called DemoCar. The Hitachi system consists of a fuel pump and injectors operating under a fuel pressure of 30 MPa.

The acoustic attenuation performance of a single-pass, perforated concentric silencer filled with continuous strand fibers is investigated theoretically and experimentally. One-dimensional analytical and three-dimensional boundary element methods are employed to predict the acoustic attenuation in the absence of mean flow. Measured complex characteristic impedance and wave number are used to account for the wave propagation through absorbing fiber. The perforation impedance facing the fiber is also presented in terms of the complex characteristic impedance and wave number. The effects of perforate duct porosity and the fiber density are examined. Comparisons of predictions with the experiments illustrate the need for multidimensional analysis at higher frequencies, while the one-dimensional treatment provides a reasonable accuracy at lower frequencies, as expected. The study also shows a significant improvement in the acoustic attenuation of the silencer due to fiber absorption.

Computational fluid dynamics of gas-fueled large-bore spark ignition engines with pre-chamber ignition can speed up the design process of these engines provided that 1) the reliability of the results is not affected by poor meshing and 2) the time cost of the meshing process does not negatively compensate for the advantages of running a computer simulation. In this work a flame propagation model that runs with arbitrary hybrid meshes was developed and coupled with the KIVA4-MHI CFD solver, in order to address these aims. The solver follows the G-Equation level-set method for turbulent flame propagation by Tan and Reitz, and employs improved numerics to handle meshes featuring different cell types such as hexahedra, tetrahedra, square pyramids and triangular prisms. Detailed reaction kinetics from the SpeedCHEM solver are used to compute the non-equilibrium composition evolution downstream and upstream of the flame surface, where chemical equilibrium is instead assumed.

In the next incoming future the necessity of reducing the raw emissions leads to the challenge of an increment of the thermal engine efficiency. In particular it is necessary to increase the engine efficiency not only at full load but also at partial load conditions. In the open literature very few technical papers are available on the partial load conditions analysis. In the present paper the analysis of the effect of the throttle valve rotational direction on the mixture formation is analyzed. The engine was a PFI 4-valves motorcycle engine. The throttle valve opening angle was 17.2°, which lays between the very partial load and the partial load condition. The CFD code adopted for the analysis was the FIRE AVL code v. 2013.2. The exhaust, intake and compression phases till TDC were simulated: inlet/outlet boundary conditions from 1D simulations were imposed.

The component response synthesis approach utilizing frequency response function (FRF) has been used to analyze the dynamic interaction of two or more vehicle components coupled at discrete interface points. This method is somewhat suitable for computing higher frequency response because experimental component FRFs can be incorporated into the formulation directly. However its calculations are quite sensitive to measurement errors in the FRFs due to the several matrix inversion steps involved. In the past, researchers have essentially used a combined direct inverse and truncated singular valued decomposition (TSVD) technique to ensure a stable calculation, which is typically applied semi-empirically due to the lack of understanding of the influence of measurement error.

Recent studies in sheet metal forming, conducted at universities world wide, emphasize the development of computer aided techniques for process simulation. To be practical and acceptable in a production environment, these codes must be easy to use and allow relatively quick solutions. Often, it is not necessary to make exact predictions but rather to establish the influence of process variables upon part quality, tool stresses, material flow, and material thickness variation. In cooperation with its industrial partners, the ERC for Net Shape Manufacturing of the Ohio State University has applied a number of computer codes for analysis and design of sheet metal forming operations. This paper gives a few selected examples taken from automotive applications and illustrates practical uses of computer simulations to improve productivity and reduce tool development and manufacturing costs.

Large Eddy Simulation (LES) represents nowadays one of the most promising techniques for the evaluation of the dynamics and evolution of turbulent structures characterizing internal combustion engines (ICE). In the present paper, subdivided into two parts, the capabilities of the open-source CFD code OpenFOAM® v2.3.0 are assessed in order to evaluate its suitability for engine cold flow LES analyses. Firstly, the code dissipative attitude is evaluated through an inviscid vortex convection test to ensure that the levels of numerical dissipation are compatible with LES needs. Quality and completeness estimators for LES simulations are then proposed. In particular the Pope M parameter is used as a LES completeness indicator while the LSR parameter provides useful insights far calibrating the grid density. Other parameters such as the two-grid LESIQk index are also discussed.

The increasing of the overall engine performance requires the investigation of the unsteady engine phenomena affecting intake air flow and the air-fuel mixing process. The “standard” RANS methodology often doesn't allow one to achieve a qualitative and quantitative accurate prediction of these phenomena. The aim of this paper is to show the potential and the limits of LES numerical technique in the simulation of actual IC engine flows and to assess the influence of some basic parameters on the LES simulation results. The paper introduces the use of a merit parameter suggested by Pope for evaluating the quality of the LES solution. The CFD code used is Fluent v6.2 and two basic test cases have been simulated. The first one is the flow over a backward facing step in order to perform a preliminary parametric numerical analysis. A one-equation dynamic subgrid-scales turbulence model is used.