Simulation of the autoignition and combustion of <Emphasis Type="Italic">n</Emphasis>-heptane droplets using a detailed kinetic mechanism

The autoignition and combustion of n-heptane droplets are simulated using a detailed kinetic mechanism. A mathematical model, based on first principles, contains no adjustable parameters. The burning rate constants for the combustion of droplets are calculated over a wide range of pressures, temperature, fuel-to-oxidizer equivalence ratios of the gas-droplet suspension, and droplet diameters. The calculated and measured delay times of autoignition of droplets are compared. The calculation results agree well with the available experimental data. The detonability of gas-droplet suspensions with partial pre-evaporation of fuel is estimated.     

Ignition limit and shock-to-detonation transition mode of n-heptane/air mixture in high-speed wedge flows

In this work, oblique detonation of n-heptane/air mixture in high-speed wedge flows is simulated by solving the reactive Euler equations with a two-dimensional (2D) configuration. This is a first attempt to model complicated hydrocarbon fuel oblique detonation waves (ODWs) with a detailed chemistry (44 species and 112 reactions). Effects of freestream equivalence ratios and velocities are considered, and the abrupt and smooth transition from oblique shock to detonation are predicted. Ignition limit, ODW characteristics, and predictability of the transition mode are discussed. Firstly, homogeneous constant-volume ignition calculations are performed for both fuel-lean and stoichiometric mixtures. The results show that the ignition delay generally increases with the wedge angle. However, a negative wedge angle dependence is observed, due to the negative temperature coefficient effects. The wedge angle range for successful ignition of n-heptane/air mixtures decreases when the wedge length is reduced. From two-dimensional simulations of stationary ODWs, the initiation length generally decreases with the freestream equivalence ratio, but the transition length exhibits weakly non-monotonic dependence. Smooth ODW typically occurs for lean conditions (equivalence ratio < 0.4). The interactions between shock/compression waves and chemical reaction inside the induction zone are also studied with the chemical explosive mode analysis. Moreover, the predictability of the shock-to-detonation transition mode is explored through quantifying the relation between ignition delay and chemical excitation time. It is demonstrated that the ignition delay (the elapsed time of the heat release rate, HRR, reaches the maximum) increases, but the excitation time (the time duration from the instant of 5% maximum HRR to that of the maximum) decreases with the freestream equivalence ratio for the three studied oncoming flow velocities. Smaller excitation time corresponds to stronger pressure waves from the ignition location behind the oblique shock. When the ratio of excitation time to ignition delay is high (e.g., > 0.5 for n-C₇H₁₆, > 0.3 for C₂H₂ and > 0.2 for H₂, based on the existing data compilation in this work), smooth transition is more likely to occur.     

An experimental and modeling study of the shock tube ignition of a mixture of n-heptane and n-propylbenzene as a surrogate for a large alkyl benzene

Alkyl aromatics are an important chemical class in gasoline, jet and diesel fuels. In the present work, an n-propylbenzene and n-heptane mixture is studied as a possible surrogate for large alkyl benzenes contained in diesel fuels. To evaluate it as a surrogate, ignition delay times have been measured in a heated high pressure shock tube (HPST) for a mixture of 57% n-propylbenzene/43% n-heptane in air (≈21% O₂, ≈79% N₂) at equivalence ratios of 0.29, 0.49, 0.98 and 1.95 and compressed pressures of 1, 10 and 30 atm over a temperature range of 1000–1600 K. The effects of reflected-shock pressure and equivalence ratio on ignition delay time were determined and common trends highlighted. A combined n-propylbenzene and n-heptane reaction mechanism was assembled and simulations of the shock tube experiments were carried out. The simulation results showed very good agreement with the experimental data for ignition delay times. Sensitivity and reaction pathway analyses have been performed to reveal the important reactions responsible for fuel oxidation under the shock tube conditions studied. It was found that at 1000 K, the main consumption pathways for n-propylbenzene are abstraction reactions on the alkyl chain, with particular selectivity to the allylic site. In comparison at 1500 K, the unimolecular decomposition of the fuel is the main consumption pathway.     

Experimental and modeling study of the autoignition characteristics of commercial diesel under engine-relevant conditions

Knowledge of the autoignition characteristics of diesel fuels is of great importance for understanding the combustion performance in engines and developing surrogate fuels. Here ignition delays of China's stage 6 diesel, a commercial fuel, were measured in a heated rapid compression machine (RCM) under engine-relevant conditions. Gas-phase autoignition experiments were carried out at equivalence ratios ranging from 0.37 to 1.0, under compressed pressures of 10, 15, and 20?bar, and within a temperature range of 685–865?K. In all investigated conditions, negative temperature coefficient (NTC) behavior of the total ignition delays is observed. The autoignition of the diesel fuel exhibits pronounced two-stage characteristics with strong low-temperature reactivity. Experimental results indicate that the total ignition delays shorten with increasing compressed pressure, oxygen mole fraction and fuel mole fraction. The first-stage ignition delays are mainly controlled by compressed temperature and also affected by oxygen mole fraction and compressed pressure but show a very weak dependence on fuel mole fraction. Correlations describing the first-stage ignition delay and the total ignition delay were proposed to further clarify the ignition delay dependence on the multiple factors. Additionally, it is found that the newly measured ignition delays well coincide with and complement the diesel ignition data in the literature. A recently developed diesel mechanism was used to simulate the diesel autoignition on the RCM. The simulation results are found to agree well the experimental measurements over the whole temperature ranges. Species concentration analysis and brute force sensitivity analysis were also conducted to identify the crucial species and reactions controlling the autoignition of the diesel fuel.     

Autoignition and combustion of hydrocarbon-hydrogen-air homogeneous and heterogeneous ternary mixtures

A numerical simulation of the ignition and combustion of hydrocarbon-hydrogen-air homogeneous and heterogeneous (gas-drop) ternary mixtures for three hydrocarbon fuels (n-heptane, n-decane, and n-dodecane) is for the first time performed. The simulation is carried out based on a fully validated detailed kinetic mechanism of the oxidation of n-dodecane, which includes the mechanisms of the oxidation of n-decane, n-heptane, and hydrogen as constituent parts. It is demonstrated that the addition of hydrogen to a homogeneous or heterogeneous hydrocarbon-air mixture increases the total ignition delay time at temperatures below 1050 K, i.e., hydrogen acts as an ignition inhibitor. At low temperatures, even ternary mixtures with a very high hydrogen concentration show multistage ignition, with the temperature dependence of the ignition delay time exhibiting a negative temperature coefficient region. Conversely, the addition of hydrogen to homogeneous and heterogeneous hydrocarbon-air mixtures at temperatures above 1050 K reduces the total ignition delay time, i.e., hydrogen acts as an autoignition promoter. These effects should be kept in mind when discussing the prospects for the practical use of hydrogen-containing fuel mixtures, as well as in solving the problems of fire and explosion safety.     

Investigation of the Livengood–Wu integral for modelling autoignition in a high-pressure bomb

The reaction progress variable, which is widely used in premixed and diffusion combustion studies, comprises a set of pre-selected intermediate species to denote reaction progress. Progress towards autoignition can also be traced by the Livengood–Wu (LW) integral. Autoignition occurs when the LW integral attains a value of unity. This concept is further explored by applying it to an inhomogeneous mixture scenario, to determine the time and place of autoignition occurrence. A semidetailed mechanism (137 species and 633 reactions) for n-heptane/iso-octane/toluene is used in this study. Two numerical schemes based on the LW integral are proposed and incorporated into a computational fluid dynamics platform, to model autoignition in a 3D configuration, when a spray is injected into a constant volume bomb under diesel engine conditions. Tabulated chemistry, a traditional method of modelling autoignition using information from pre-calculated igniting diffusion flames, is also used for comparison purposes. The associated predicted pressure profiles are compared with experimental measurements.     

Simulation of spray development and turbulent combustion processes in low and high speed diesel engines by the CMC-ISR model

Simulation is performed to analyse the characteristics of turbulent spray combustion in conventional low and high speed diesel engine conditions. Turbulence–chemistry interaction is resolved by the Conditional Moment Closure (CMC) model in the spatially integrated form of an Incompletely Stirred Reactor (ISR). After validation against measured pressure traces, characteristic length and time scales and dimensionless numbers are estimated at the locations of sequentially injected fuel groups. Conditional flame structures are calculated for sequentially evaporated fuel groups to consider different available periods for ignition chemistry. Injection overlaps the combustion period in the high rpm engine, while most combustion occurs after injection and evaporation are complete in the low rpm engine. Ignition occurs in rich premixture with the initial peak temperature at the equivalence ratio around 2–4 as observed in Dec [2]. It corresponds to the most reactive mixture fraction of the minimum ignition delay for the given mixture states. Combustion proceeds to lean and rich sides in the mixture fraction space as a diffusion process by turbulence. The mean scalar dissipation rates (SDRs) are lower than the extinction limit to show stability of diffusion flames throughout the combustion period.     

Multidimensional chemistry coordinate mapping approach for combustion modelling with finite-rate chemistry

A multidimensional chemistry coordinate mapping (CCM) approach is presented for efficient integration of chemical kinetics in numerical simulations of turbulent reactive flows. In CCM the flow transport is integrated in the computational cells in physical space, whereas the integration chemical reactions are carried out in a phase space made up of a few principal variables. Each cell in the phase space corresponds to several computational cells in the physical space, resulting in a speedup of the numerical integration. In reactive flows with small hydrocarbon fuels two principal variables have been shown to be satisfactory to construct the phase space. The two principal variables are the temperature (T) and the specific element mass ratio of the H atom (J _H). A third principal variable, σ=?J _H·?J _H, which is related to the dissipation rate of J _H, is required to construct the phase space for combustion processes with an initially non-premixed mixture. For complex higher hydrocarbon fuels, e.g. n-heptane, care has to be taken in selecting the phase space in order to model the low-temperature chemistry and ignition process. In this article, a multidimensional CCM algorithm is described for a systematic selection of the principal variables. The method is evaluated by simulating a laminar partially remixed pre-vaporised n-heptane jet ignition process. The CCM approach is then extended to simulate n-heptane spray combustion by coupling the CCM and Reynolds averaged Navier–Stokes (RANS) code. It is shown that the computational time for the integration of chemical reactions can be reduced to only 3–7%, while the result from the CCM method is identical to that of direct integration of the chemistry in the computational cells.     

Influence of pilot-fuel mixing on the spatio-temporal progression of two-stage autoignition of diesel-sprays in low-reactivity ambient fuel-air mixture

The spatial and temporal locations of autoignition for direct-injection compression-ignition engines depend on fuel chemistry, temperature, pressure, and mixing trajectories in the fuel jets. Dual-fuel systems can provide insight into both fuel-chemistry and physical effects by varying fuel reactivities and engine operating conditions. In this context, the spatial and temporal progression of two-stage autoignition of a diesel-fuel surrogate, n-heptane, in a lean-premixed charge of synthetic natural-gas (NG) and air is imaged in an optically accessible heavy-duty diesel engine. The lean-premixed charge of NG is prepared by fumigation upstream of the engine intake manifold. Optical diagnostics include high-speed (15kfps) cool-flame chemiluminescence-imaging as an indicator of low-temperature heat-release (LTHR) and OH* chemiluminescence-imaging as an indicator high-temperature heat-release (HTHR). NG prolongs the ignition delay of the pilot fuel and increases the combustion duration. Zero-dimensional chemical-kinetics simulations provide further understanding by replicating a Lagrangian perspective for mixtures evolving along streamlines originating either at the fuel nozzle or in the ambient gas, for which the pilot-fuel concentration is either decreasing or increasing, respectively. The zero-dimensional simulations predict that LTHR initiates most likely on the air streamlines before transitioning to HTHR, either on fuel-streamlines or on air-streamlines in regions of near-constant ?. Due to the relatively short pilot-fuel injection-durations, the transient increase in entrainment near the end of injection (entrainment wave) is important for quickly creating auto-ignitable mixtures. To achieve desired combustion characteristics, e.g., multiple ignition-kernels and favorable combustion phasing and location (e.g., for reducing wall heat-transfer or optimizing charge stratification), adjusting injection parameters could tailor mixing trajectories to offset changes in fuel ignition chemistry.     

Effects of pre-spark heat release on engine knock limit

It is well known that spark ignited engine efficiency is limited by end gas autoignition, commonly known as knock. This study focuses on a recently discovered phenomena, pre-spark heat release (PSHR) due to low-temperature chemistry, and its impact on knock behavior. Boosted operating conditions are more common as engines are downsizing and downspeeding in efforts to increase fuel economy and prone to PSHR. Experiments were prone at fixed fueling and air fuel ratio for a range of intake temperature that spanned the threshold for PSHR. It was found that when PSHR occurred, the knock-limited combustion phasing was insensitive to intake temperature; higher intake temperatures did not require retarded timings as it is usual. Inspection of the temperature–pressure history overlaid on ignition delay contours allow the results to be explained. The temperature rise from the low-temperature reactions moves the end gas state into the negative temperature coefficient (NTC) region, which terminates the heat release reactions. The end gas then resides in the long ignition delay peninsula, which inhibits knock.     

Autoignition of propane–air mixtures behind reflected shock waves

Ignition times and autoignition modes for propane–air mixtures have been studied behind reflected shock waves. Experiments were performed over temperatures between 1000 and 1750 K, pressures between 2 and 20 atm, and equivalence ratios of = 0.5, 1.0, and 2.0. Ignition delay times were determined using pressure measurements, C₂ emission profiles, and luminosity measurements in the visible spectrum (380–680 nm). Empirical correlations for ignition time for low temperature (1000–1300 K) and high temperature (1300–1800 K) ranges have been deduced from the experimental data. Different autoignition modes of the mixture (strong, transient, and weak) were identified by comparing velocities of reflected shock wave at different distances from the reflecting wall.     

Influence of intermediate temperature heat release on autoignition reactivity of single-stage ignition fuels with varying octane sensitivity

This study investigates the effects of intermediate temperature heat release (ITHR) on autoignition reactivity of full boiling range gasolines with different octane sensitivity through intake temperature and simulated exhaust gas recirculation (EGR) sweeps in a homogenous charge compression ignition (HCCI) engine. To isolate the ITHR effects, low temperature reactivity was suppressed through the use of high intake temperature and low intake oxygen mole fraction. For quantification of ITHR, a new method was applied to the engine data by examining the maximum value of the second derivative of heat release rate. Combustion phasing comparisons of fuels with octane sensitivity showed that fuel with less octane sensitivity became more reactive as intake temperature and simulated EGR ratio decreased, while fuel with higher octane sensitivity had a reverse trend. For all of the fuels that were tested, the amount of ITHR increased as the intake temperature and oxygen mole fraction increased. These ITHR trends, depending on octane sensitivity, were almost identical with the trends of combustion phasing, showing that ITHR significantly affects fuel autoignition reactivity and determines octane sensitivity.     

Conditional statistics of nonreacting and reacting sprays in turbulent flows by direct numerical simulation

Conditional statistics concerning evaporation and combustion of a spray are investigated in homogeneous, isotropic, and decaying two-dimensional (2D) turbulence. Randomly distributed, polydisperse droplets of n-heptane go through single-step combustion chemistry. Attention is focused on parametric effects of initial Sauter mean radius (SMR), turbulence level and droplet velocity in both reacting and nonreacting cases. A simple linear model for the conditional evaporation rate is proposed and validated against DNS data. A conventional β-probability density function (pdf) is shown to be valid with no peak occurring on the fuel side. The amplitude mapping closure (AMC) model works well for the conditional scalar dissipation rate with evaporating and reacting sprays. Parametric study shows that initial SMR and droplet velocity are major factors affecting conditional flame structures, whereas the effect of reaction is not significant except during autoignition.     

Two-stage autoignition and combustion mode evolution in boundary layer flows above a cold flat plate

Boundary layers are omnipresent in fundamental kinetic experimental facilities and practical combustion engines, which can cause ambiguity and misleading results in kinetic target acquisition and even abnormal engine combustion. In this paper, using n-heptane as a representative large hydrocarbon fuel exhibiting pronounced low-temperature chemistry (LTC), two-dimensional numerical simulation is conducted to resolve the transient autoignition phenomena affected by a boundary layer. We focus on the ignition characteristics and the subsequent combustion mode evolution of a hot combustible mixture flowing over a colder flat plate in an isobaric environment. For cases with autoignition occurring within the boundary layer, similarity is observed in the first-stage ignition as manifested by a constant temperature at all locations. The first-stage ignition is found to be rarely affected by heat and radical loss within the boundary layer. While for the main ignition event, an obvious dependence of ignition process on boundary layer thickness is identified, where the thermal-chemical process exhibits similarity at locations with similar boundary layer thickness, and the main ignition tends to first occur within the boundary layer at the domain end and generates a C-shape reaction front. It is found that sequential spontaneous autoignition is the dominant subsequent combustion mode at high-pressure conditions. At low to intermediate pressures, auto-ignition assisted flame propagation is nevertheless the dominant mode for combustion evolution. This research identifies novel features of autoignition and the subsequent combustion mode evolution affected by a cold, fully developed boundary layer, and provides useful guidance to the interpretation of abnormal combustion and combustion mode evolution in boundary layer flows.     

Autoignition of reacting mixtures at engine-relevant conditions using confined spherically expanding flames

Propagation of a confined spherically expanding flame induces isentropic compression that can culminate in autoignition and/or detonation under conducive thermodynamic conditions. This relatively simple technique measures a distinct ‘characteristic ignition delay time’ and complements other established approaches such as the rapid compression machine and shock tube. The present study details this methodology by examining the autoignition characteristics of dimethyl-ether/oxygen/nitrogen/helium reactive mixtures for equivalence ratios of 0.6 and 0.9, an initial temperature of 468 K, and initial pressures of 3 to 6 atm. The experimental results display the classic two-stage ignition typical of dimethyl-ether oxidation at low-temperatures with first-stage ignition occurring at approximately 3.6 times the initial pressure. To aid in the interpretation of the experimental results, two numerical models were used: a zero-dimensional batch reactor model, which accepts experimental pressure-time history and calculates the sensitivities of characteristic ignition delay times to kinetics, and a low Mach number, Lagrangian one-dimensional code that was developed to model both flame propagation and end-gas autoignition. Simulation results were shown to adequately capture the physics of unsteady flame propagation, end-gas autoignition, and the controlling reactions of the latter. It was found also that under certain conditions the behavior of first and second ignition stages could be modified by unsteady pressure effects.     

Skeletal mechanism generation with CSP and validation for premixed n-heptane flames

An automated procedure has been previously developed to generate simplified skeletal reaction mechanisms for the combustion of n-heptane/air mixtures at equivalence ratios between 0.5 and 2.0 and different pressures. The algorithm is based on a Computational Singular Perturbation (CSP)-generated database of importance indices computed from homogeneous n-heptane/air ignition solutions. In this paper, we examine the accuracy of these simplified mechanisms when they are used for modeling laminar n-heptane/air premixed flames. The objective is to evaluate the accuracy of the simplified models when transport processes lead to local mixture compositions that are not necessarily part of the comprehensive homogeneous ignition databases. The detailed mechanism was developed by Curran et al. and involves 560 species and 2538 reactions. The smallest skeletal mechanism considered consists of 66 species and 326 reactions. We show that these skeletal mechanisms yield good agreement with the detailed model for premixed n-heptane flames, over a wide range of equivalence ratios and pressures, for global flame properties. They also exhibit good accuracy in predicting certain elements of internal flame structure, especially the profiles of temperature and major chemical species. On the other hand, we find larger errors in the concentrations of many minor/radical species, particularly in the region where low-temperature chemistry plays a significant role. We also observe that the low-temperature chemistry of n-heptane can play an important role at very lean or very rich mixtures, reaching these limits first at high pressure. This has implications to numerical simulations of non-premixed flames where these lean and rich regions occur naturally.     

Skeletal mechanism generation and analysis for n-heptane with CSP

We use a procedure based on the decomposition into fast and slow dynamical components offered by the Computational Singular Perturbation (CSP) method to generate automatically skeletal kinetic mechanisms for the simplification of the kinetics of n-heptane oxidation. The detailed mechanism of the n-heptane oxidation here considered has been proposed by Curran et al. and involves 561 species and 2538 reactions. After carrying out a critical assessment of important aspects of this procedure, we show that the comprehensive skeletal kinetic mechanisms so generated are able to reproduce the main features of n-heptane ignition at various initial pressures and temperatures and equivalence ratios. A by-product of the algorithm that generates the skeletal mechanisms is the identification of the network of important species and reactions at a given state of the kinetic system. The analysis of this network is carried out by resorting to a visual representation of the pathways at selected time instants of the ignition process. Visual inspection of the pathways enables the identification and comparison of the relevant kinetic processes as obtained at different ignition regimes. The graphs are generated by interfacing the model reduction procedure with the open-source package graphviz.     

Numerical simulations of the ignition of n-heptane droplets in the transition diameter range from heterogeneous to homogeneous ignition

Spontaneous ignition of single n-heptane droplets in a constant volume filled with air is numerically simulated with the spherical symmetry. The volume is closed against mass, species, and energy transfer. The numerical model is fully transient. It continues calculation even after the droplet has completely vaporized, and therefore can predict pre-vaporized ignition. Initial pressure and initial air temperature are fixed at 3 MPa and 773 K, respectively. The droplet is initially at room temperature, and its diameter is between 1 and 100 μm. When the overall equivalence ratio is fixed to be sufficiently large, there exists no ignition limit in terms of initial droplet diameter d₀, and the ignition delay takes a minimum value at certain d₀. In such a case, transition from the heterogeneous ignition to the homogeneous ignition with decreasing d₀ is observed. When d₀ is fixed to be so small that the ignition would not occur in an infinite volume of air, the ignition delay takes a minimum value at certain , which is less than unity. Two-stage ignition behavior is investigated with this model. Ignition delay of a cool flame has the dependence on d₀ that is similar to that of ignition delay of a hot flame when is unity. When is almost zero, the ignition limit for cool flame in terms of d₀ is not identified unlike that for hot flame.     

Microgravity experiments of flame spreading along a fuel droplet array in fuel vapor-air mixture

Combustion experiments of fuel droplet array in fuel vapor-air mixture were performed at microgravities to investigate growth mechanism of group combustion of fuel droplets. A 10-droplet array was inserted into the test section filled with a saturated fuel vapor-air mixture as a simple model of prevaporized sprays. Gas equivalence ratio of the fuel vapor-air mixture was regulated by the test section temperature. n-Decane droplets of 0.8 mm in the initial diameter were suspended at the crossing points of 10 sets of X-shaped suspenders. The first droplet was ignited by a hot wire to initiate flame spread along a fuel droplet array. Flame spread speed was obtained from the history of the leading edge position of a spreading flame. Effects of droplet spacing and gas equivalence ratio on the flame spreading behavior and the flame spread speed were examined. The droplet spacing and the gas equivalence ratio were varied from 1.6 to 10.2 mm and from 0.2 to 0.7, respectively. The gas equivalence ratio has little effect on the relationship between the flame spreading behavior and the droplet spacing. The flame spread speed increases as the increase in the gas equivalence ratio at all droplet spacings. The influence of the gas equivalence ratio on the flame spread speed becomes strong as the increase in the droplet spacings. The flame spread speed increases as the increase in the droplet spacing, and then decreases. The maximum flame spread speed appears in the range from 2.4 to 3 mm at all gas equivalence ratios.     

Ignition delay of fatty acid methyl ester fuel droplets: Microgravity experiments and detailed numerical modeling

Recent optical engine studies have linked increases in NO_x emissions from fatty acid methyl ester combustion to differences in the premixed autoignition zone of the diesel fuel jet. In this study, ignition of single, isolated liquid droplets in quiescent, high temperature air was considered as a means of gaining insight into the transient, partially premixed ignition conditions that exist in the autoignition zone of a fatty acid methyl ester fuel jet. Normal gravity and microgravity (10⁻⁴ m/s²) droplet ignition delay experiments were conducted by use of a variety of neat methyl esters and commercial soy methyl ester. Droplet ignition experiments were chosen because spherically symmetric droplet combustion represents the simplest two-phase, time-dependent chemically reacting flow system permitting a numerical solution with complex physical submodels. To create spherically symmetric conditions for direct comparison with a detailed numerical model, experiments were conducted in microgravity by use of a 1.1 s drop tower. In the experiments, droplets were grown and deployed onto 14 μm silicon carbide fibers and injected into a tube furnace containing atmospheric pressure air at temperatures up to 1300 K. The ignition event was characterized by measurement of UV emission from hydroxyl radical (OH*) chemiluminescence. The experimental results were compared against predictions from a time-dependent, spherically symmetric droplet combustion simulation with detailed gas phase chemical kinetics, spectrally resolved radiative heat transfer and multi-component transport. By use of a skeletal chemical kinetic mechanism (125 species, 713 reactions), the computed ignition delay period for methyl decanoate (C₁₁H₂₂O₂) showed excellent agreement with experimental results at furnace temperatures greater than 1200 K.