SEARCH · Search NASA

Results for “COMBUSTION - EFFECTS OF FUEL ATOMIZATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2021

High-speed phase contrast imaging of spray breakup of jet fuels under combusting conditions

X-ray phase contrast imaging is conducted at 90,517 Hz on a fuel spray to characterize the breakup of the bulk fuel into ligaments and individual droplets in a single-sector, swirl-stabilized combustor under reacting conditions. Phase contrast imaging is used to qualitatively assess spray breakup for two fuels: Jet-A, which represents a normal viscosity fuel, and a blend of JP-5 (64%) and farnesane (36% by volume), which represents a high-viscosity fuel, at three fuel-flow rates and two air inlet temperatures. Time-series images of the fuels show that atomization occurs much more rapidly at higher fuel pressure with thinner ligaments being rapidly shredded. The high-viscosity fuel is qualitatively observed to have markedly longer, thicker ligaments and larger droplet diameters after the primary breakup than the normal viscosity fuel. Using an image processing and droplet quantification strategy, the mean droplet velocity magnitudes were found to decrease with decreasing fuel pressure and fuel-flow rate. The increase in preheat temperature from 323 to 370 K results in an increase in the mean droplet velocity magnitude and a decrease in numerical average and Sauter mean diameter, indicating improved atomization. The high-viscosity fuel was found to have larger average droplet diameters at each condition than the normal viscosity fuel. The increased preheat is found to enhance the combustion efficiency, which is hypothesized to couple with the improved atomization due to the higher combustion temperatures. This indicates the importance of examining spray atomization and breakup under realistic combusting conditions for relevant fuels as the heat release directly couples with the atomization process.

42 ENGINEERING↗

DOE OSTI2020

SYSTEMATIC MULTI-INDEX VALIDATIONS OF SIDI ENGINE FLOW FIELD LES COMPUTATIONS USING CRANK ANGLE-RESOLVED PIV MEASUREMENTS

In-cylinder flow fields make significant impacts on the fuel atomization, fuel mixture formation, and combustion process in spark ignition direct injection (SIDI) engines. In recent years, model-based simulation approaches are preferred in regard to investigating the transient in-cylinder flow field characteristics. Most commonly, the simulation models are validated using single representative flow field at a typical crank angle measured by particle image velocimetry (PIV). However, it provides only limited knowledge about the flow field which is highly three-dimensional and of transient nature. In this study, crank angle-resolved Ply measurements are conducted on three distinct planes inside the cylinder to capture the transient process of flow field characteristics which vary with the crank angle. These three planes consist of one tumble plane through the spark plug tip, one tumble plane along two intake ports, and one swirl plane at 30 mm below the cylinder head. Large eddy simulation (LES) is employed for the numerical computations using the CONVERGE codes. On the basis of large datasets for both temporal and spatial domains, a multi-index systematic validation approach is conducted. Crank angle-resolved calculations of global indices and local indices are implemented using the flow fields velocity data obtained from both PIV and LES on select planes. Global indices reveal the trends in similarities of different crank angle degrees and locations, while local indices give the detail comparison results. In summary, with the systematic multi-index validation approach, the crucial crank angle degrees and locations for model verification will be detected. Furthermore, the corresponding critical flow features are analyzed. Practical guideline of flow field validation is proposed.

Crank angle-resolved PIV↗

DOE OSTI2021

Multitude Characterization and Prediction of DOE Advanced Biofuels Properties

Advanced multitude of experiments ranging from the liquid fuel to combustion are conducted on Co-OPTIMA fuels to aid the characterization of the fuels. The series of targeted experiments characterized Co-OPTIMA fuel spray atomization, flame topology, flame speed, autoignition, volatility, viscosity, soot/coking, and compatibility. The fuels are selected and prioritized based on input from national lab members. The research characterized and predicted biomass-based, low greenhouse gas fuels and blends combustion, autoignition, and physical properties of mixtures of identified compounds at engine-relevant conditions, in particular those properties that blend non-linearly. The main challenge of Co-OPTIMA is the evaluation of a variety of biofuels and blends in all the reaction conditions that might be encountered in new high-efficiency engines. Despite the improved high-throughput experimental techniques, it seems unlikely that all of the performance metrics could be measured for all relevant petroleum derived, bio-derived molecules and mixtures, and reaction conditions. A series of targeted experiments ranging from the liquid fuel to the combustion process is required, and to extract the maximum information from each experiment. These targeted experiments evaluated how a specific fuel will perform in an engine. The series of targeted experiments are as follows: (a) Spray Atomization, Vaporization and Droplet Formation (b) Combustion Flame and Local Fuel/Air Image-Based Measurements (c) Laminar Flame Speed Measurements (d) Autoignition and Soot Measurements (e) Synchrotron Coupled Fundamental Autoignition Experiments (f) Fuel Coking and Hot Surface Deposit g) Fuel Volatility Measurements (h) Viscosity Measurements (i) Seal Flexible Fuel Compatibility These experimental processes provide an essential pathway for the prediction of fuel behaviors in engines and systematic process for fuel down select.

09 BIOMASS FUELS↗

DOE OSTI2026

The Reaction of Atomic Hydrogen with an Imine: Direct Kinetic Measurements and Computational Studies

Imines are of interest in the combustion modeling of nitrogen-based fuels such as ammonia, as intermediates in atmospheric oxidation of amines, and may offer a new class of environmentally benign substitutes for perhalogenated reagents. Elementary gas-phase rate constants for the consumption of H atoms by 1,1,1,3,3,3-hexafluoro-2-propanimine were measured over 294–736 K using the laser flash photolysis/atomic resonance fluorescence technique. The results are summarized as 3.0 × 10 –11 exp(−14.5 kJ mol –1 /RT) cm 3 molecule –1 s –1 with a confidence interval of ±13%. These are the first temperature-dependent kinetic experiments on an imine of the form R 2 C═NH and thereby provide the only validation of transition-state modeling for such systems. Several pathways are feasible; comparison of theory and experiment suggests the dominant channel is addition of H atoms to the nitrogen atom. On this basis, the observations are rationalized quantitatively. Computations also yield the thermochemistry of (CF 3 ) 2 CNH, with Δ f H 298 = −1232 kJ mol –1 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Experimental and Kinetic Modeling Study of 3-Methyl-2-butenol (Prenol) Oxidation

Longer chain alcohols with 4–5 carbon atoms are attractive alternative fuels as they can be derived from biological sources and since their combustion leads to lower exhaust gas levels of NO x and soot compared to commercial fossil fuels. The auto-ignition behavior of fuels that contain both a hydroxyl group and a C=C double bond in their molecular structure is not well established in the literature. Understanding the influence of these functional groups on the ignition behavior of fuels is critical in the development of tailor-made fuels for advanced combustion engines. In this study, ignition delay times of an unsaturated alcohol, 3-methyl-2-butenol (prenol), are measured using a high-pressure shock tube and a rapid compression machine at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 in “air” in the temperature range 600–1400 K. Furthermore, a detailed kinetic model is developed and validated using the new experimental data in this study. In addition, speciation data in a jet-stirred reactor, ignition delay times, and laminar burning velocities available in the literature were also used to validate the new kinetic model. Fuel flux and sensitivity analyses are performed using this new model to determine the important fuel consumption pathways and critical reactions that affect the reactivity of prenol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Combined Experimental/Numerical Study of the Soot Formation Process in a Gasoline Direct-Injection Spray in the Presence of Laser-Induced Plasma Ignition

Here, combustion issued from an eight-hole, direct-injection spray was experimentally studied in a constant-volume pre-burn combustion vessel using simultaneous high-speed diffused back-illumination extinction imaging (DBIEI) and OH* chemiluminescence. DBIEI has been employed to observe the liquid-phase of the spray and to quantitatively investigate the soot formation and oxidation taking place during combustion. The fuel-air mixture was ignited with a plasma induced by a single-shot Nd:YAG laser, permitting precise control of the ignition location in space and time. OH* chemiluminescence was used to track the high-temperature ignition and flame. The study showed that increasing the delay between the end of injection and ignition drastically reduces soot formation without necessarily compromising combustion efficiency. For long delays between the end of injection and ignition (1.9 ms) soot formation was eliminated in the main downstream charge of the fuel spray. However, poorly atomized and large droplets formed at the end of injection (dribble) eventually do form soot near the injector even when none is formed in the main charge. The quantitative soot measurements for these spray and ignition scenarios, resolved in time and space, represents a significant new achievement. Reynolds-averaged Navier-Stokes (RANS) simulations were performed to assess spray mixing and combustion. An analysis of the predicted fuel-air mixture in key regions, defined based upon experimental observations, was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing. The mixture analysis indicates that soot production can be avoided if the flame propagates into regions where the equivalence ratio (Φ) is already below 2. Reactive RANS simulations have also been performed, but with a poor match against the experiment, as the flame speed and heat-release rate are largely over estimated. This modeling weakness appears related to a very high level of turbulent viscosity predicted for the high-momentum spray in the RANS simulations, which is an important consideration for modeling ignition and flame propagation in mixtures immediately created by the spray.

42 ENGINEERING↗

DOE OSTI2023

A physics-driven Σ-Y atomization model for heavy-duty engine simulations

The atomization of a liquid jet is a multi-scale and multi-physics problem of interest for many engineering applications. Particularly, it drives the fuel-oxidizer mixing and dictates the efficiency of combustion engines. High-fidelity multi-phase simulations remain challenging due to the excessive computational cost required to capture all the atomization scales. Thus, atomization models are necessary to represent sub-grid liquid structures. Here, in this work, a modified Σ-Y model in the context of Eulerian-Lagrangian Spray Atomization (ELSA) is used to transport the surface area density of the spray. The model's predictive performance is assessed under various operating conditions relevant to heavy-duty engines using the Engine Combustion Network (ECN) Spray C and Spray D research-grade injectors. The classic droplet collision formulation of the Σ-Y model alone does not replicate the response of the measured spray surface area to changes in injector, ambient pressure, and injection pressure, requiring individual tuning of the model's parameters. Instead, a transition between dense and dilute spray breakup mechanisms is proposed in terms of the average droplet spacing. The collision breakup mechanism represents the dilute spray, whereas the droplet size in the dense spray is driven by a competition between the integral scale of turbulence and the balance between the turbulent kinetic energy and the surface energy of the droplets. Such an approach minimizes the requirements for model tuning. Moreover, the role of the constants of a compressible standard k-ε RANS model is assessed in the injector's internal flow and external spray simulation framework, and an updated set is proposed. The results are validated against x-ray radiography and Ultra-Small Angle X-ray Scattering (USAXS) data and highlight the predictive capabilities of the proposed physics-driven Σ-Y model, which is compatible with engine simulation turnaround times.

42 ENGINEERING↗

DOE OSTI2021

Fatty acid derived alkyl ether fuels for compression ignition

A fuel for an internal combustion engine includes a fatty alkyl ether having a formula corresponding to formula (II): wherein x is 1-8, and y is 0 to 3; and the alkyl is an alkyl group having a number of carbon atoms that is less than the number of carbon atoms in the alkyl chain on the opposite side of the oxygen atom. The fatty alky ether can be used as a neat fuel or blend with biodiesel, diesel, ethanol or other fuels. The fatty alkyl ethers are improved in cetane number and cold flow properties over a biodiesel with fatty acid methyl acid methyl ester compounds. This is particularly valuable for compression ignition engines. A method of combustion in several types of engines is also disclosed.

Monroe, Eric↗

DOE OSTI2020

Comparative analysis of characteristics of a supersonic cw chemical HF laser on molecular fluorine and nitrogen trifluoride

We have experimentally studied the characteristics of a supersonic cw chemical HF laser with a flat nozzle array based on a nozzle – nozzle reagent mixing scheme operating with the use of F{sub 2} – D{sub 2} – He and NF{sub 3} – D{sub 2} – He fuel compositions in an atomic fluorine generator. The active medium flow field pattern, its gas-dynamic characteristics, and the laser radiation spectrum are considered. The comparison of these characteristics has made it possible to establish a strong impact on the flow field (lasing zone length) of the chemical composition of the combustion products from the atomic fluorine generator. The gas-dynamic parameters of the active medium formed by using a NF{sub 3} – D{sub 2} – He fuel have turned out to be less acceptable in terms of the gas flow pressure recovery than in the case of F{sub 2} – D{sub 2} – He fuel, while the radiation spectrum was more preferable from the viewpoint of passing through the atmosphere. (paper)

36 MATERIALS SCIENCE↗

DOE OSTI2022

Numerical and Experimental Investigations on the Ignition Behavior of OME

On the path towards climate-neutral future mobility, the usage of synthetic fuels derived from renewable power sources, so-called e-fuels, will be necessary. Oxygenated e-fuels, which contain oxygen in their chemical structure, not only have the potential to realize a climate-neutral powertrain, but also to burn more cleanly in terms of soot formation. Polyoxymethylene dimethyl ethers (PODE or OMEs) are a frequently discussed representative of such combustibles. However, to operate compression ignition engines with these fuels achieving maximum efficiency and minimum emissions, the physical-chemical behavior of OMEs needs to be understood and quantified. Especially the detailed characterization of physical and chemical properties of the spray is of utmost importance for the optimization of the injection and the mixture formation process. The presented work aimed to develop a comprehensive CFD model to specify the differences between OMEs and dodecane, which served as a reference diesel-like fuel, with regards to spray atomization, mixing and auto-ignition for single- and multi-injection patterns. The simulation results were validated against experimental data from a high-temperature and high-pressure combustion vessel. The sprays’ liquid and vapor phase penetration were measured with Mie-scattering and schlieren-imaging as well as diffuse back illumination and Rayleigh-scattering for both fuels. To characterize the ignition process and the flame propagation, measurements of the OH* chemiluminescence of the flame were carried out. Significant differences in the ignition behavior between OMEs and dodecane could be identified in both experiments and CFD simulations. Liquid penetration as well as flame lift-off length are shown to be consistently longer for OMEs. Zones of high reaction activity differ substantially for the two fuels: Along the spray center axis for OMEs and at the shear boundary layers of fuel and ambient air for dodecane. Additionally, the transient behavior of high temperature reactions for OME is predicted to be much faster.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Alkyl dialkoxyalkanoates as bioderived, high cetane diesel fuels

A fuel for an internal combustion engine includes a C 5 to C 30 dialkoxyalkanoate corresponding to formula (I), wherein the R 1 group is —H or a —CH 3 group, the R 2 and R 2 ′ groups are alkyl groups independently selected to have 1 to 9 carbon atoms; and the R 3 group is selected to have 1 to 9 carbon atoms. The compounds described herein may be used as neat fuels or mixed fuels (with diesel, biodiesel, jet fuel, marine fuel or other fuel compounds) in autoignition or spark ignition engines, such as diesel engines, gasoline (spark ignition) engines, MCCI, Homogeneous Charge Compression Ignition (HCCI) engines, or more generally in Low-Temperature Gasoline Combustion (LTGC) engines (using gasoline-like fuels), that have the high-efficiency advantages of HCCI but can operate with some level of charge inhomogeneities. Methods of making these compounds are environmentally friendly and can be done in a continuous manner.

Myllenbeck, Nicholas↗

DOE OSTI2025

Enabling Low-Temperature (LTP) Ignition Technologies for Multi-Mode Engines through the Development of a Validated High-Fidelity LTP Model for Predicative Simulations Tools

The goal of multi-mode engine architectures is to extend current lean-burn dilution limits with renewable fuels, which requires spark plugs to deposit high energies (hundreds of mJ) in order to initiate ignition and complete combustion. At elevated energy deposition rates, spark plugs experience increased electrode erosion and thermal losses, which ultimately shortens the spark-plug lifetime and lowers ignition efficiency. As such, in order to safeguard the efficiency gains of multi-mode concepts, new and improved ignition technologies are required. Recently, non-equilibrium low-temperature plasmas (LTP) have been shown to promote energy-efficient ignition via quenching and transport of electronically excited atoms and molecules, selective radical production and fast heating of hydrocarbon/air mixtures [1-2]. Thus, LTP is seen as a technology that can potentially improve the energy extraction efficiency of fuels, while enabling kinetically controlled combustion modes towards fuel leaner conditions to realize current DOE VTO goals of improving the sustainability of future mobility [3]. Although many previous studies have demonstrated the efficacy of plasma-assisted ignition to enhance combustion, the detailed enhancement mechanisms remain largely unknown, especially for oxygenated fuels and at elevated pressures that are most relevant to practical engine conditions. These barriers hinder the development of accurate and comprehensive numerical models that seek to describe LTP-based ignition in existing engine design software tools and methods. Current state-of-the-art simulation capabilities for LTP ignition systems are in need of improvements since they deliver qualitative results only due to important limitations of existing approaches. Firstly, validated kinetic models with elementary steps for plasma discharges in oxygenated fuel/air mixtures of relevance to the transportation sector are required. Such kinetic models do not exist at present and will be developed and validated within this project. Secondly, plasma discharges and reactive mixture ignition are multi-scale, unsteady processes requiring high-performance numerical methods and software that execute efficiently on DOE supercomputers. Such software does not exist at present and will be developed and applied to practical LTP ignition scenarios as part of this project. Thirdly, experimental databases that are tailored to serve as benchmark in support of the development of predictive computational models of LTP ignition do not exist and will be part of this project.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2023

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstraction by $CH_3O\dot{O}$ radicals

Hydrogen atom abstraction by methyl peroxy ($CH_3O\dot{O}$) radicals can play an important role in gasoline/ethanol interacting chemistry for fuels that produce high concentrations of methyl radicals. Detailed kinetic reactions for hydrogen atom abstraction by $CH_3O\dot{O}$ radicals from the components of FGF-LLNL (a gasoline surrogate) including cyclopentane, toluene, 1-hexene, n -heptane, and isooctane have been systematically studied in this work. Here, the M06-2X/6-311 ++ G(d,p) level of theory was used to obtain the optimized structure and vibrational frequency for all stationary points and the low-frequency torsional modes. The 1-D hindered rotor treatment for low-frequency torsional modes was treated at M06-2X/6-31G level of theory. The UCCSD(T)-F12a/cc-pVDZ-F12 and QCISD(T)/CBS level of theory were used to calculate single point energies for all species. High pressure limiting rate constants for all hydrogen atom abstraction channels were performed using conventional transition state theory with unsymmetric tunneling corrections. Individual rate constants are reported in the temperature range from 298.15 to 2000 K. Our computed results show that the abstraction of allylic hydrogen atoms from 1-hexene is the fastest at low temperatures. When the temperature increases, the hydrogen atom abstraction reaction channel at the primary alkyl site gradually becomes dominant. Thermodynamics properties for all stable species and high-pressure limiting rate constants for each reaction pathway obtained in this work were incorporated into the latest gasoline surrogate/ethanol model to investigate the influence of the rate constants calculated here on model predicted ignition delay times.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

New reactions of diazene and related species for modelling combustion of amine fuels

Potential energy surfaces for reactions involving N 2 H 2 isomers of diazene (diimide) have been explored using density functional theory, with energies based on coupled-cluster theory. A focus is on processes that create or consume these species, and isomerisation between the E (trans) and Z (cis) forms of HNNH. Furthermore, these include isomerisation and dissociation pathways for HNNH, addition of H atoms to form N 2 H 3 , abstraction by H atoms yielding short-lived NNH, and abstraction reactions of H with N 2 H 3 . Transition state and capture theories are applied for high-pressure-limiting behaviour, while low-pressure and falloff regions are characterised via the methods of Troe and coworkers. Rate constants and thermochemistry are provided to improve models of diamine chemistry, relevant to the combustion of NH 3 especially at high concentrations, high pressures or under reducing conditions. Results indicate that amine radical recombination mainly yields the E HNNH isomer, while H-abstraction from N 2 H 3 results in E HNNH and H 2 NN. However, at elevated temperature E → Z isomerisation becomes competitive, and Z HNNH, being more reactive, acts to enhance the diazene consumption rate.

10 SYNTHETIC FUELS↗

DOE OSTI2022

Room-Temperature Optical Detection of 14 CO 2 below the Natural Abundance with Two-Color Cavity Ring-Down Spectroscopy

Radiocarbon’s natural production, radiative decay, and isotopic rarity make it a unique tool to probe carbonaceous systems in the life and earth sciences. However, the difficulty of current radiocarbon ( 14 C) detection methods limits scientific adoption. Here, two-color cavity ring-down spectroscopy detects 14 CO 2 in room-temperature samples with an accuracy of one-tenth the natural abundance in 3 min. The intracavity pump–probe measurement uses two cavity-enhanced lasers to cancel out cavity ring-down rate fluctuations and strong one-photon absorption interference (>10 000 1/s) from hot-band transitions of CO 2 isotopologues. Selective, room-temperature detection of small 14 CO 2 absorption signals (<1 1/s) reduces the technical and operational burdens for cavity-enhanced measurements of radiocarbon, which can benefit a wide range of applications like biomedical research and field-detection of combusted fossil fuels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Characterizing Impacts of Dry Coal Feeding in High Pressure Oxy-Coal Combustion Systems

Reaction Engineering International (REI) has managed a team of experts from the University of Utah, Southeast University (SEU) in Nanjing, China, Electric Power Research Institute (EPRI), Corrosion Management Ltd. (C-M), Praxair, and Brigham Young University (BYU) to investigate dry pulverized coal feeding for pressurized oxy-coal combustion. Dry feed firing systems for entrained flow, pressurized, oxy-coal combustors have not been well developed, although related technologies have been used in the Shell Gasification Process and for pressurized fluid bed combustion. DOE-funded research recently completed at REI and the University of Utah focused on characterizing impacts of high temperatures and pressures in oxy-coal combustion systems. For high pressure combustion, that research used a coal slurry feed into a 17 bar pressurized combustor. As a consequence of that research, it was identified that fuel feeding and firing system flexibility are challenges that require attention. Based on that experience, the approach of using a coal slurry feed system leads to challenges in producing consistent atomization of the slurry, which causes burnout problems, especially at high pressures. In addition, slurry atomization processes may be difficult to scale to sizes appropriate for practical commercial use. Dry pressurized coal burner systems, on the other hand, have the potential to yield efficiency gains, improve flexibility and facilitate applications at larger scales. Experimental work was conducted at the University of Utah Industrial Combustion and Gasification Research Facility as well as the 100 kW pressurized oxy-coal combustor (POC) facility at Brigham Young University. Mechanism development and CFD-based combustion and dense-phase flow modeling were performed at REI. Successful completion of the project objectives has resulted in the following key deliverables: 1) Design and prototype of a pressurized pulverized coal feeding and oxy-firing system 2) Data from a 100kW, pressurized (15 bar) entrained flow reactor with a dry feeding delivery and burner system that describes flame characteristics, radiative heat flux profiles, carbon burnout, along with characteristics of ash aerosols, fouling, and slagging. 3) Validated and transportable models that describe the relevant conditions in pressurized oxy-combustion systems and that can be used for scale-up and optimization. 4) Principles to guide design of high pressure, pilot-scale and full-scale coal oxy-firing systems. 5) Assessment of pressurized oxy-combustion impacts on key parameters relevant to oxy-coal fired utility boilers such as coal devolatilization, char oxidation, mineral matter transformation, deposition, and corrosion. The experimental data, pressurized oxy-firing system principles, and process mechanisms provided by this work can be used by electric utilities, boiler OEMs, equipment suppliers, design firms, software vendors, consultants and government agencies to assess the use of high temperature and high pressure oxy-combustion in current research and to guide development of new oxy-coal boiler designs.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2022

Characterizing Impacts of Dry Coal Feeding in High Pressure Oxy-Coal Combustion Systems

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2023

High-Accuracy Heats of Formation for Alkane Oxidation: From Small to Large via the Automated CBH-ANL Method

It is generally challenging to obtain high-accuracy predictions for the heat of formation for species with more than a handful of heavy atoms, such as those of importance in standard combustion mechanisms. To this end, we construct the CBH-ANL approach and illustrate that, for a set of 194 alkane oxidation species, it can be used to produce ΔH f (o K) values with 2σ uncertainties of 0.2-0.5 kcal mol -1 . This set includes the alkanes, hydroperoxides, and alkyl, peroxy, and hydroperoxyalkyl radicals for 17 representative hydrocarbon fuels containing up to 10 heavy atoms with various degrees of branching in the alkane backbone. The CBH-ANL approach, automated in the QTC and AutoMech software suites, builds balanced chemical equations for the calculation of ΔH f (o K), in which the reference species may be up to five heavy atoms. The high-level ANLo and ANL1 reference ΔH f (o K) values are further refined for even the largest of these reference species with a novel laddering approach. We perform a comprehensive quantification of the uncertainties for both the individual reference species (the largest of which is 0.15 kcal mol -1 ) and the propagation of those uncertainties when used in the calculation of ΔH f (o K) for the 194 target species. Here we examine the sensitivity of the predicted ΔH f (o K) values to (i) electronic energies from various methods, including ωB97X-D/cc-pVTZ, B2PLYP-D3/cc-pVTZ, CCSD(T)-F12b/cc-pVDZ-F12//B2PLYP-D3/cc-pVTZ, and CCSD(T)-F12b/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ; (ii) the zero-point vibrational energies (ZPVEs), where we consider harmonic ZPVEs as well as two scaling-based estimates of the anharmonic ZPVEs, all implemented for both ωB97X-D/cc-pVTZ and B2PLYP-D3/cc-pVTZ calculations; (iii) the particular CBH-ANL scheme employed; and (iv) the procedure for choosing the reference conformer for the analyses. The discussion concludes with a summary of the estimated overall uncertainty in the predictions and a validation of the predictions for the alkane subset.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗