As powertrains and IC engines continue to grow in complexity, many vehicle manufacturers (OEMs) are turning to simulation in an effort to reduce design validation and calibration costs. Ultimately, their aim is to complete this process entirely within the virtual domain, without the need for any physical testing. Practical simulation techniques for the prediction of knock in spark ignition (SI) engines rely on empirical ignition delay correlations (IDCs). These IDCs are used to approximate the complex ignition delay characteristics of real and surrogate fuel compositions with respect to temperature, pressure and mixture composition. Over the last 40 years, a large number of IDCs have been put forward in the literature, spanning a broad range of fuels, operating conditions and calibration methods. However, the applicability of these tools has yet to be verified at the high brake mean effective pressure (BMEP) operating conditions relevant to highly boosted, downsized engines. Here, the applicability of 16 gasoline-relevant IDCs for predicting knock onset at high loads (BMEP > 30bar) has been investigated by comparing the knock predictions from each IDC against experimentally measured knock onset times. Firstly, a detailed investigation into cylinder pressure data processing techniques was performed to determine which knock detection and angle of knock onset (aKO) measurement methods were most appropriate at high loads. A method based on the maximum amplitude pressure oscillation (MAPO) during knock-free operation best estimated cycle classifications, whilst Shahlari’s Signal Energy Ratio technique [1] most accurately predicted knock onset. To the author’s knowledge, this is the first time that such a comprehensive study on the accuracy of these techniques at such high loads has been conducted. Importantly, these findings represent a valuable framework to inform other researchers in the field of knocking combustion on which techniques are needed to extract accurate and relevant information from measured cylinder pressure records. Secondly, the data processing techniques derived were applied to experimental data collected across a wide range of high BMEP operating conditions (up to a maximum of 32 bar) using a 1.6 litre, 4-cylinder SI engine. Trapped charge composition and temperature were predicted using a calibrated 1D model of the engine, whilst the temperature of a hypothetical hotspot in the unburned zone was estimated separately by assuming adiabatic compression from a point after intake valve closing and by mapping γ (the ratio of specific heat capacities) as a function of temperature. This revealed that none of the IDCs tested performed well at conditions relevant to modern, downsized engines. The IDC that achieved the best overall balance between aKO accuracy and cycle-classification agreement was the “cool-flame” correlation for iso-octane proposed by Ma [2]. However, this had an unacceptably high average aKO error of ±3.5° compared to the ±2°CA limit observed within the literature, and its average cycle-classification accuracy was below 60%. The main reason for this relatively modest accuracy was a large number of false-positive cycle classifications, which mainly occurred in slow or late burning cycles. Further work should therefore focus on methods to reduce the number of false positive classifications obtained with this correlation, which could be achieved using empirical correlations to describe the latest point in the cycle for which knock would be permitted to occur in terms other measureable combustion parameters. Overall, this research has generated a unique insight into combustion at very high loads, as well as an extensive dataset that can be used for future research to improve the accuracy of empirical knock modelling techniques. Furthermore, this work has demonstrated that for the purposes of virtual spark timing calibration and the avoidance of knock, the current crop of practical simulation tools is not accurate enough at the conditions relevant to modern SI engines and has provided a better understanding of their limitations. These findings represent a major contribution to the field from both a research perspective and for industrial applications.

NOx produced during combustion can persist in the exhaust gases of a gas turbine engine in quantities significant to induce regulatory concerns. There has been much research which has led to important insights into NOx chemistry. One method of NOx reduction is exhaust gas recirculation. In exhaust gas recirculation, a portion of the exhaust gases that exit are redirected to the inlet air stream that enters the combustion chamber, along with fuel. Due to the presence of NOx in the exhaust gases which are subsequently introduced into the burner, knowledge of the effects of NOx on combustion is advantageous. Contrary to general NOx research, little has been conducted to investigate the sensitizing effects of NO2 and N2O addition to methane/oxygen combustion. Experiments were made with dilute and real fuel air mixtures of CH4/O2/Ar with the addition of NO2 and N2O. The real fuel air concentrations were made with the addition of NO2 only. The equivalence ratios of mixtures made were 0.5, 1 and 2. The experimental pressure range was 1 - 44 atm and the temperature range tested was 1177 – 2095 K. The additives NO2 and N2O were added in concentrations from 831 ppm to 3539 ppm. The results of the mixtures with NO2 have a reduction in ignition delay time across the pressure ranges tested, and the mixtures with N2O show a similar trend. At 1.3 atm, the NO2 831 ppm mixture shows a 65% reduction and shows a 75% reduction at 30 atm. The NO2 mixtures showed a higher decrease in ignition time than the N2O mixtures. The real fuel air mixture also showed a reduction. Sensitivity Analyses were performed. The two most dominant reactions in the NO2 mixtures are the reaction O+H2 = O+OH and the reaction CH3+NO2 = CH3O+NO. The presence of this second reaction is the means by which NO2 decreases ignition delay time, which is indicated in the experimental results. The reaction produces CH3O which is reactive and can participate in chain propagating reactions, speeding up ignition. The two dominant reactions for the N2O mixture are the reaction O+H2 = O+OH and, interestingly, the other dominant reaction is the reverse of the initiation reaction in the N2O-mechanism: O+N2+M = N2O+M. The reverse of this reaction is the direct oxidation of nitrous oxide. The O produced in this reaction can then speed up ignition by partaking in propagation reactions, which was experimentally observed.

No description available.

A new shock tube targeting low temperature, high pressure, and long test times was designed and installed at the Turbomachinery Laboratory in December of 2008. The single-pulse shock tube uses either lexan diaphragms or die-scored aluminum disks of up to 4 mm in thickness. The modular design of the tube allows for optimum operation over a large range of thermodynamic conditions from 1 to 100 atm and between 600-4000 K behind the reflected shock wave. The new facility allows for ignition delay time, chemical kinetics, high-temperature spectroscopy, vaporization, atomization, and solid particulate experiments. An example series of ignition delay time experiments was made on mixtures of CH4/C2H6/O2/Ar at pressures from 1 to 30.7 atm, intermediate temperatures from 1082 to 2248 K, varying dilutions (between 75 and 98% diluent), and equivalence ratios ranging from fuel lean (0.5) to fuel rich (2.0) in this new facility. The percentage by volume variation and equivalence ratios for the mixtures studied were chosen to cover a wide parameter space not previously well studied. Results are then used to validate and improve a detailed kinetics mechanism which models the oxidation and ignition of methane and other higher order hydrocarbons, through C4, with interest in further developing reactions important to methane- and ethane-related chemistry.

Future thrusts for gasoline engine development can be broadly summarized into two categories: (i) efficiency improvements in conventional spark ignition engines, and (ii) development of advance compression ignition (ACI) concepts. Efficiency improvements in conventional spark ignition engines requires downsizing (and turbocharging) which may be achieved by using high octane gasolines, whereas, low octane gasolines fuels are anticipated for ACI concepts. The current work provides the essential combustion kinetic data, targeting both thrusts, that is needed to develop high fidelity gasoline surrogate mechanisms and surrogate complexity guidelines. Ignition delay times of a wide range of certified gasolines and surrogates are reported here. These measurements were performed in shock tubes and rapid compression machines over a wide range of experimental conditions (650 – 1250 K, 10 – 40 bar) relevant to internal combustion engines. Using the measured the data and chemical kinetic analyses, the surrogate complexity requirements for these gasolines in homogeneous environments are specified. For the discussions presented here, gasolines are classified into three categories: (i) Low octane gasolines including Saudi Aramco’s light naphtha fuel (anti-knock index, AKI = (RON + MON)/2 = 64; Sensitivity (S) = RON – MON = 1), certified FACE (Fuels for Advanced Combustion Engines) gasoline I and J (AKI ~ 70, S = 0.7 and 3 respectively), and their Primary Reference Fuels (PRF, mixtures of n-heptane and iso-octane) and multi-component surrogates. (ii) Mid octane gasolines including FACE A and C (AKI ~ 84, S ~ 0 and 1 respectively) and their PRF surrogates. Laser absorption measurements of intermediate and product species formed during gasoline/surrogate oxidation are also reported. (iii) A wide range of n-heptane/iso-octane/toluene (TPRF) blends to adequately represent the octane and sensitivity requirements of high octane gasolines including FACE gasoline F and G (AKI ~ 91, S = 5.6 and 11 respectively) and certified Haltermann (AKI ~ 87, S = 7.6) and Coryton (AKI ~ 92, S = 10.9) gasolines. To assess conditions where shock tubes may not be ideal devices for ignition delay measurements, this work also presents a detailed discussion on shock tube pre-ignition affected ignition data and the ignition regimes in homogeneous environments. The shock tube studies on pre-ignition and associated bulk ignition advance may help engines research community understand and control super-knock events.

Indiana University-Purdue University Indianapolis (IUPUI) / This study contributes to a better understanding of ignition by hot combustion gases which finds application in internal combustion chambers with pre-chamber ignition as well as in wave rotor engine applications. The experimental apparatus consists of two combustion chambers: a pre chamber that generates the transient hot jet of gas and a main chamber which contains the main fuel air blend under study. Variables considered are three fuel mixtures (Hydrogen, Methane, 50\% Hydrogen-Methane), initial pressure in the pre-chamber ranging from 1 to 2 atm, equivalence ratio of the fuel air mixture in the main combustion chamber ranging from 0.4 to 1.5, and initial temperature of the main combustion chamber mixture ranging from 297 K to 500 K. Experimental data makes use of 4 pressure sensors with a recorded sampling rate up to 300 kHz, as well as high speed Schlieren imaging with a recorded frame rate up to 20,833 frame per seconds. Results shows an overall increase in ignition delay with increasing equivalence ratio. High temperature of the main chamber blend was found not to affect hot jet ignition delay considerably. Physical mixing effects, and density of the main chamber mixture have a greater effect on hot jet ignition delay.

The combustion behavior of methane and ethane is important to the study of natural gas and other alternative fuels that are comprised primarily of these two basic hydrocarbons. Understanding the transition from methane-dominated ignition kinetics to ethane-dominated kinetics for increasing levels of ethane is also of fundamental interest toward the understanding of hydrocarbon chemical kinetics. Much research has been conducted on the two fuels individually, but experimental data of the combustion of blends of methane and ethane is limited to ratios that recreate typical natural gas compositions (up to ~20% ethane molar concentration). The goal of this study was to provide a comprehensive data set of ignition delay times of the combustion of blends of methane and ethane at near atmospheric pressure. A group of ten diluted CH4/C2H6/O2/Ar mixtures of varying concentrations, fuel blend ratios, and equivalence ratios (0.5 and 1.0) were studied over the temperature range 1223 to 2248 K and over the pressure range 0.65 to 1.42 atm using a new shock tube at the University of Central Florida Gas Dynamics Laboratory. Mixtures were diluted with either 75 or 98% argon by volume. The fuel blend ratio was varied between 100% CH4 and 100% C2H6. Reaction progress was monitored by observing chemiluminescence emission from CH* at 431 nm and the pressure. Experimental data were compared against three detailed chemical kinetics mechanisms. Model predictions of CH* emission profiles and derived ignition delay times were plotted against the experimental data. The models agree well with the experimental data for mixtures with low levels of ethane, up to 25% molar concentration, but show increasing error as the relative ethane fuel concentration increases. The predictions of the separate models also diverge from each other with increasing relative ethane fuel concentration. Therefore, the data set obtained from the present work provides valuable information for the future improvement of chemical kinetics models for ethane combustion.

Shock-tube experiments were performed with syngas mixtures near atmospheric pressure with varying equivalence ratios behind reflected shock waves. Pressure and hydroxyl radical (OH*) emission traces were recorded and used to calculate ignition delay time for a single mixture at equivalence ratios of [phi ]=0.4, 0.7, 1.0, and 2.0 over a range of temperatures from 913-1803 K. The syngas mixture was tested at full concentration as well as with 98% dilution in Argon. The full concentration mixtures were used to compare ignition delay time measurements with the theoretical calculations obtained through the use of chemical kinetics modeling using the Davis et al. mechanism. The dilute mixtures were used to study the OH* emission profiles compared to those of the kinetics model. The model was in poor agreement with the experimental data especially at lower temperatures with an ignition delay difference of more than an order of magnitude. These ignition delay time data supplement the few existing data and are in relative agreement. The species profile comparison of OH* compared to the model also showed poor agreement, with the worst agreement at the highest temperatures. While the disagreements with ignition delay time and profile comparisons cannot be explained at this time, the data presented support other findings. The data provide additional information towards understanding this disagreement relative to syngas mixtures despite the relatively well known kinetics of the primary constituents Hydrogen and Carbon Monoxide.

Autoignition of non-premixed methane-air mixtures is investigated using first-order Conditional Moment closure (CMC). In CMC, scalar quantities are conditionally averaged with respect to a conserved scalar, usually the mixture fraction. The conditional fluctuations are often of small order, allowing the chemical source term to be modeled as a function of the conditional species concentrations and the conditional enthalpy (temperature). The first-order CMC derivation leaves many terms unclosed such as the conditional scalar dissipation rate, velocity and turbulent fluxes, and the probability density function. Submodels for these quantities are discussed and validated against Direct Numerical Simulations (DNS). The CMC and the turbulent velocity and mixing fields calculations are decoupled based on the frozen mixing assumption, and the CMC equations are cross-stream averaged across the flow following the shear flow approximation. Finite differences are used to discretize the equations, and a two-step fractional method is implemented to treat separately the stiff chemical source term. The stiff ODE solver LSODE is used to solve the resulting system of equations. The recently developed detailed chemical kinetics mechanism UBC-Mech 1.0 is employed throughout this study, and preexisting mechanisms are visited. Several ignition criteria are also investigated. Homogeneous and inhomogeneous CMC calculations are performed in order to investigate the role of physical transport in autoignition. Furthermore, the results of the perfectly homogeneous reactor calculations are presented and the critical value of the scalar dissipation rate for ignition is determined. The results are compared to the shock tube experimental data of Sullivan et al. The current results show good agreement with the experiments in terms of both ignition delay and ignition kernel location, and the trends obtained in the experiments are successfully reproduced. The results were shown to be sensitive to the scalar dissipation model, the chemical kinetics, and the ignition criterion.

Autoignition of non-premixed methane-air mixtures is investigated using first-order Conditional Moment closure (CMC). In CMC, scalar quantities are conditionally averaged with respect to a conserved scalar, usually the mixture fraction. The conditional fluctuations are often of small order, allowing the chemical source term to be modeled as a function of the conditional species concentrations and the conditional enthalpy (temperature). The first-order CMC derivation leaves many terms unclosed such as the conditional scalar dissipation rate, velocity and turbulent fluxes, and the probability density function. Submodels for these quantities are discussed and validated against Direct Numerical Simulations (DNS). The CMC and the turbulent velocity and mixing fields calculations are decoupled based on the frozen mixing assumption, and the CMC equations are cross-stream averaged across the flow following the shear flow approximation. Finite differences are used to discretize the equations, and a two-step fractional method is implemented to treat separately the stiff chemical source term. The stiff ODE solver LSODE is used to solve the resulting system of equations. The recently developed detailed chemical kinetics mechanism UBC-Mech 1.0 is employed throughout this study, and preexisting mechanisms are visited. Several ignition criteria are also investigated. Homogeneous and inhomogeneous CMC calculations are performed in order to investigate the role of physical transport in autoignition. Furthermore, the results of the perfectly homogeneous reactor calculations are presented and the critical value of the scalar dissipation rate for ignition is determined. The results are compared to the shock tube experimental data of Sullivan et al. The current results show good agreement with the experiments in terms of both ignition delay and ignition kernel location, and the trends obtained in the experiments are successfully reproduced. The results were shown to be sensitive to the scalar dissipation model, the chemical kinetics, and the ignition criterion.