Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Soot”

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.

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Co-optimization of fuel properties, combustion system geometry, and injection strategy for conventional diesel fuel

Here, studies have shown that fuel properties can impact an engine’s operation in several ways, including ignition delay, sooting tendency, mixture formation, and combustion temperature. In mixing-controlled compression ignition (MCCI) engines, the fuel system design and piston bowl geometry significantly affect combustion performance and emissions. Based on current information, it is difficult to draw conclusions about fuel property effects and sensitivities. The central fuel hypothesis approach used in the US Department of Energy Co-Optima program has worked well for spark ignition fuels: identifying critical fuel property ranges is sufficient to screen fuel blends that are expected to maximize efficiency and reduce pollutant emissions. However, for MCCI-relevant fuels, the information gained from past studies is not sufficient to build such a merit function or to allow for performing a similar screening of fuel blends. It is hypothesized that a co-optimization of a fuel’s physical and chemical properties, combustion system geometry, and injection strategy could leverage synergies between the effects of the fuel properties and geometries, resulting in improved performance over state-of-the-art. A machine learning–assisted unconstrained global optimization algorithm was used to explore a design space comprising 23 independent variables. The results show that physical property effects were minimal even for large variations in fuel properties, and the only interaction effect that was observed was the effect of varied fuel density parameters on fuel/air mixture formation. Nevertheless, these interactions were not sufficient in magnitude to significantly affect optimization results. Therefore, analysis of the results suggests that fuel physical properties cannot be leveraged in a co-optimization context to increase engine efficiency.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

CFD modeling of non-catalytic, partial-oxidation engine reformer for flare mitigation

Flaring associated natural gas is commonly employed in the oil and gas industry to reduce methane (CH 4 ) emissions but generates carbon dioxide (CO 2 ) and harmful pollutants, significantly contributing to air pollution and posing risks to public health. To mitigate this impact, M2X Energy Inc. has developed a small-scale, modular gas-to-methanol system. This system features an engine reformer that performs fuel-rich partial oxidation of wellhead gas to produce syngas—a mixture of carbon monoxide (CO) and hydrogen (H 2 )—followed by a downstream reactor for methanol synthesis. This study focused on computational fluid dynamics (CFD) modeling of the engine reformer to simulate partial oxidation chemistry, predict the rich-burn operating limit, and assess syngas quality, ultimately aiding in design and operational optimization. The CFD model, developed within a Reynolds-Averaged Navier-Stokes (RANS) turbulence framework, incorporated sub-models for turbulent combustion, a chemical mechanism with polycyclic aromatic hydrocarbon (PAH) pathways, and soot emissions to accurately capture the fuel-rich, turbulent jet ignition and combustion processes. Model validation against experimental data showed good agreement across pre- and main-chamber pressures, apparent heat release rates, and exhaust gas concentrations of key species (H 2 , CO, CO 2 , CH 4 ) for varying intake equivalence ratios. Here, the model identified a rich-burn operating limit near a fuel-air equivalence ratio of 2.35, consistent with experimental observations. Furthermore, syngas quality analysis revealed that extending the rich-burn limit through engine reformer optimization could enhance syngas production, contributing to higher methanol synthesis efficiency.

Computational Fluid Dynamics↗
DOE OSTI2024

Compact and accurate chemical mechanism for methane pyrolysis with PAH growth

In this study, a reliable and compact chemical mechanism of gas-phase methane pyrolysis leading to formation of large polycyclic aromatic hydrocarbon (PAH) molecules has been developed. This model is designed for studies of carbon nanostructure synthesis such as carbon black and graphene flakes, including soot growth kinetics. Methane pyrolysis with carbon nanostructure synthesis is a two-stage process where conversion of CH 4 to C 2 H 2 precedes the growth of PAH molecules from acetylene. We present a single chemical mechanism that accurately describes both stages. We have constructed a compact and accurate chemical mechanism capable of modeling both stages of methane pyrolysis based on the ABF mechanism which was expanded with the most prominent reaction pathways from the mechanism by Tao for small PAH molecules and HACA pathways for larger PAH molecules, up to 37 aromatic rings. The resulting mechanism was validated through comparison to multiple available sets of experimental data. Good agreement with the experimental data for both processes was obtained. Performance of the mechanism was tested for pyrolysis of methane-rich mixtures under long residence times leading to abundant formation of PAH molecules. It is shown that the inclusion of larger PAH species (up to A37) in the chemical mechanism is important for accurate prediction of the fraction of carbon converted to PAH molecules and, correspondingly, the residual fraction of acetylene in the mixture. The mechanism file is available upon request.

08 HYDROGEN↗
DOE OSTI2023

Uncertainty quantification of a deep learning fuel property prediction model

Deep learning models are being widely used in the field of combustion. Given the black-box nature of typical neural network based models, uncertainty quantification (UQ) is critical to ensure the reliability of predictions as well as the training datasets, and for a principled quantification of noise and its various sources. Deep learning surrogate models for predicting properties of chemical compounds and mixtures have been recently shown to be promising for enabling data-driven fuel design and optimization, with the ultimate goal of improving efficiency and lowering emissions from combustion engines. In this study, UQ is performed for a multi-task deep learning model that simultaneously predicts the research octane number (RON), Motor Octane Number (MON), and Yield Sooting Index (YSI) of pure components and multicomponent blends. The deep learning model is comprised of three smaller networks: Extractor 1, Extractor 2, and Predictor, and a mixing operator. The molecular fingerprints of individual components are encoded via Extractor 1 and Extractor 2, the mixing operator generates fingerprints for mixtures/blends based on linear mixing operation, and the predictor maps the fingerprint to the target properties. Two different classes of UQ methods, Monte Carlo ensemble methods and Bayesian neural networks (BNNs), are employed for quantifying the epistemic uncertainty. Combinations of Bernoulli and Gaussian distributions with DropConnect and DropOut techniques are explored as ensemble methods. All the DropConnect, DropOut and Bayesian layers are applied to the predictor network. Aleatoric uncertainty is modeled by assuming that each data point has an independent uncertainty associated with it. The results of the UQ study are further analyzed to compare the performance of BNN and ensemble methods. Although this study is confined to UQ of fuel property prediction, the methodologies are applicable to other deep learning frameworks that are being widely used in the combustion community.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Probing the low-temperature chemistry of methyl hexanoate: Insights from oxygenate intermediates

Understanding the combustion of methyl esters is crucial to elucidate kinetic pathways and predict combustion parameters, soot yields, and fuel performance of biodiesel, however most kinetic studies of methyl esters have focused on smaller, surrogate model esters. Methyl hexanoate is a larger methyl ester approaching the chain length of methyl esters found in biodiesel and has not received as much research attention as other smaller esters. The purpose of this work is to present the first atmospheric pressure combustion data of methyl hexanoate, CH 3 CH 2 CH 2 CH 2 CH 2 COOCH 3 . Mixtures of 2% methyl hexanoate in O 2 and N 2 are studied using a plug flow reactor at atmospheric pressure, wall temperatures from 573 to 973 K, residence times from roughly 1-2 s., and fuel equivalence ratios of 1, 1.5, and 2. Exhaust gases are analyzed by a gas chromatograph-mass spectrometer system and species mole fractions are presented. Here, the literature model shows satisfactory agreement with the experimental species profiles and improvements for future mechanistic studies are suggested. In particular, this work proposes new unimolecular decomposition pathways of methyl hexanoate to form methanol or methyl acetate. Furthermore, the experiment detected three unsaturated esters that are direct products of the low temperature oxidation chemistry and it provides more insight into branching ratios for the formation of methyl hexanoate radicals and for the decomposition of hydroperoxyalkyl radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The role of cool-flame fluctuations in high-pressure spray flames, studied using high-speed optical diagnostics and Large-Eddy Simulations

Here this work investigates the low- and high-temperature ignition and combustion processes, applied to the Engine Combustion Network Spray A flame, combining advanced optical diagnostics and large-eddy simulations (LES). Simultaneous high-speed (50 kHz) formaldehyde (CH 2 0) planar laser-induced fluorescence (PLIF) and line-of-sight OH* chemiluminescence imaging were used to measure the low- and high-temperature flame, during ignition as well as during quasi-steady combustion. While tracking the cool flame at the laser sheet plane, the present experimental setup allows detection of distinct ignition spots and dynamic fluctuations of the lift-off length over time, which overcomes limitations for flame tracking when using schlieren imaging. After significant development to improve LES prediction of the low-and high-temperature flame position, both during the ignition processes and quasi-steady combustion, the simulations were analyzed to gain understanding of the mixture variance and how this variance affects formation/consumption of CH 2 0. Analysis of the high-temperature ignition period shows that a key improvement in the LES is the ability to predict heterogeneous ignition sites, not only in the head of the jet, but in shear layers at the jet edge close to the position where flame lift-off eventually stabilizes. The LES analysis also shows concentrated pockets of CH 2 0, in the center of jet and at 20 mm downstream of the injector (in regions where the equivalence ratio is greater than 6), that are of similar length scale and frequency as the experiment (approximately 5–6 kHz). The periodic oscillation of CH 2 0 match the frequency of pressure waves generated during auto-ignition and reflected within the constant-volume vessel throughout injection. The ability of LES to capture the periodic appearance and destruction of CH 2 0 is particularly important because these structures travel downstream and become rich premixed flames that affect soot production.

42 ENGINEERING↗
DOE OSTI2022

The role of radical-radical chain-propagating pathways in the phenyl + propargyl reaction

Well-skipping radical-radical reactions can provide a chain-propagating pathway for formation of polycyclic radicals implicated in soot inception. Here we use controlled pyrolysis in a microreactor to isolate and examine the role of well-skipping channels in the phenyl (C 6 H 5 ) + propargyl (C 3 H 3 ) radical-radical reaction at temperatures of 800-1600 K and pressures near 25 Torr. The temperature and concentration dependence of the closed-shell (C 9 H 8 ) and radical (C 9 H 7 ) products are observed using electron-ionization mass spectrometry. The flow in the reactor is simulated using a boundary layer model employing a chemical mechanism based on recent rate coefficient calculations. Further, comparison between simulation and experiment shows reasonable agreement, within a factor of 3, while suggesting possible improvements to the model. In contrast, eliminating the well-skipping reactions from the chemistry mechanism causes a much larger discrepancy between simulation and experiment in the temperature dependence of the radical concentration, revealing that the well-skipping pathways, especially to form indenyl radical, are significant at temperatures of 1200 K and higher. While most C 9 H 7 forms by well-skipping at 25 Torr, an additional simulation indicates that the wellskipping channels only contribute around 3% of the C 9 H x yield at atmospheric pressure, thus indicating a negligible role of the well-skipping pathways at atmospheric and higher pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Visualization of post-detonation fireball flowfields and comparison to CFD modeling

Visualization of flow structures within post-detonation fireballs has been performed for benchmark validation of numerical simulations. Custom pressed PETN explosives with a 12-mm diameter hemispherical form factor were used to produce a spherically symmetric post-detonation flow with low soot yield. Hydroxyl-radical planar laser induce fluorescence (OH-PLIF) was employed to visualize the structure ranging from approximately 10 μs to 35 μs after shock breakout from the explosive pellet. Fireball simulations were performed using the HyBurn Computational Fluid Dynamics (CFD) package. Experimental OH-PLIF results were compared to synthetic OH-PLIF from post-processing of CFD simulations. From the comparison of experimental and synthetic OH-PLIF images, CFD is shown to replicate much of the flow structure observed in the experiments, revealing potential differences in turbulent length scales and OH kinetics. Furthermore, results provide significant advancement in experimental resolution of these harsh turbulent combustion environments and validate physical models thereof.

Detonator↗
DOE OSTI2024

An improved semi-global intrinsic kinetics model for high temperature carbon oxidation

Measurements of the oxidation rates of various forms of carbon (soot, graphite, coal char) have often shown an unexplained attenuation with increasing temperatures in the vicinity of 2000 K, even when accounting for diffusional transport limitations and gas-phase chemical effects (e.g. CO 2 dissociation). With the development of oxy-fuel combustion approaches for pulverized coal utilization with carbon capture, high particle temperatures are readily achieved in sufficiently oxygen-enriched environments. Here, in this work, a new semi-global intrinsic kinetics model for high temperature carbon oxidation is created by starting with a previously developed 5-step mechanism that was shown to reproduce all major known trends in carbon oxidation, except for its high temperature kinetic falloff, and incorporating a recently discovered surface oxide decomposition step. The predictions of this new model are benchmarked by deploying the kinetic model in a steady-state reacting particle code (SKIPPY) and comparing the simulated results against a carefully measured set of pulverized coal char combustion temperature measurements over a wide range of oxygen concentrations in N 2 and CO 2 environments. The results show that the inclusion of the spontaneous surface oxide decomposition reaction step significantly improves predictions at high particle temperatures. Furthermore, the simulations reveal that O atoms released from the oxide decomposition step enhance the radical pool in the near-surface region and within the particle interior itself. Incorporation of literature rates for O and OH reactions with the carbon surface results in a reduction in the predicted radical pool concentrations and a very minor enhancement of the overall carbon oxidation rate.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Flame stabilization in DME spray flames under engine-relevant conditions characterized by OH* chemiluminescence and formaldehyde laser-induced fluorescence

The transient and quasi-steady flame structures of Dimethyl Ether (DME) fuel sprays, produced by a single-hole injector (Spray D), were investigated using Planar Laser-Induced Fluorescence (PLIF) and chemiluminescence imaging in a constant-volume chamber under Engine Combustion Network (ECN) Spray A conditions (900 K ambient temperature, 60 bar ambient pressure, 1500 bar injection pressure, and 22.8 kg/m 3 ambient density). Low-temperature chemical reaction zones were visualized using formaldehyde (CH 2 O) PLIF with 355 nm excitation, while high-temperature flame regions were captured via chemiluminescence imaging of excited-state hydroxyl radicals (OH*). Both transient and quasi-steady flame structures clearly show the transition from CH 2 O to OH*, highlighting the progression from low- to high-temperature combustion, while the position of the flame is displaced for DME compared to reference hydrocarbon n-dodecane. Homogeneous reactor calculations with detailed chemistry and using adiabatic mixing for initial temperature show that CH 2 O peaks are significantly higher for DME at the same equivalence ratio, with a higher heat-release during the cool-flame regime with respect to the fuel heating value. Thus, the cool-flame dynamic as a precursor to high-temperature combustion and flame stabilization exhibit distinct behavior for DME relative to conventional hydrocarbons, and these phenomena are effectively resolved through the soot-free nature of DME and the high-speed, time-resolved diagnostics.

CH2O laser-induced fluorescence↗
DOE OSTI2025

Molecular Mass Growth Processes to Polycyclic Aromatic Hydrocarbons through Radical–Radical Reactions Exploiting Photoionization Reflectron Time-of-Flight Mass Spectrometry

Polycyclic aromatic hydrocarbons (PAHs) represent critical building blocks in molecular mass growth processes to carbonaceous nanoparticles, referred to as interstellar and circumstellar grains along with soot particles in astrophysical environments and combustion systems, respectively. Recent advancements on elucidating elementary steps to PAHs have utilized reactions of aromatic radicals, resonantly stabilized free radicals, and aliphatic radicals with closed shell hydrocarbons. However, the role of radical–radical reactions (RRRs) leading to PAHs has remained largely unexplored on the molecular level due to preceding experimental challenges in producing sufficiently high number densities of radical reactants for isomer-selective detection of products from bimolecular and termolecular reactions. This Account offers the latest developments in our knowledge on the mechanisms and pathways to PAHs via RRRs probed in a chemical microreactor at temperatures as high as 1600 K. Product preservation in a molecular beam coupled with synchrotron vacuum ultraviolet photoionization reflectron time-of-flight mass spectrometry and photoelectron photoion coincidence spectroscopy enabled isomer-selective detection of PAHs of up to three rings by their photoionization efficiency curves, which were fit with a linear combination of reference curves for identification. Experiments were combined with computational fluid dynamics modeling of the physicochemical processes in the microreactor, as well as high-level electronic structure calculations to reveal the reaction pathways of each system. Six distinct reaction mechanisms were discovered in this work: propargyl addition─benzannulation (PABA), methyl addition─ring expansion (MARE), cyclopentadienyl addition─naphthylization (CPAN), fulvenallenyl addition─cyclization─aromatization (FACA), benzyl addition─aromatization (BAA), and phenyl addition─pentacyclization (PAP). By systematically varying the number of carbon atoms in the radical reactants, molecular mass growth processes involving reactions between radicals with odd numbers of carbon atoms access aromatics carrying one, two, or three six-membered rings, whereas reactions between even- and odd-carbon-numbered radicals produce aromatics combining five- and six-membered rings. Our investigations reveal unconventional cycloadditions on excited state triplet surfaces, additions of radicals to low spin density carbon-centered radicals, spiroaromatic and fulvene-type intermediates, and highly strained bicyclic reaction intermediates, challenging current perceptions of PAH molecular mass growth processes. All of the listed mechanisms, except for FACA, feature endoergic reactions or barriers which lie above the separated reactants and therefore might be central to circumstellar environments of carbon-rich stars and planetary nebulae as their descendants, but they play no role in the gas phase of cold molecular clouds where temperatures as low as 10 K dominate. Altogether, this work provides detailed reaction mechanisms of PAH growth processes, advancing our knowledge of the chemistry of carbonaceous matter in the universe.

Addition reactions↗
DOE OSTI2021

Use of Transmission Electron Microscopy for Analysis of Aerosol Particles and Strategies for Imaging Fragile Particles

For over 25 years, transmission electron microscopy (TEM) has provided a method for the study of aerosol particles with sizes from below the optical diffraction limit to several microns, resolving the particles as well as smaller features. The wide use of this technique to study aerosol particles has contributed important insights about environmental aerosol particle samples and model atmospheric systems. TEM produces an image that is a 2D projection of aerosol particles that have been impacted onto grids and, through associated techniques and spectroscopies, can contribute additional information such as the determination of elemental composition, crystal structure, and 3D particle structures. Soot, mineral dust, and organic/inorganic particles have all been analyzed using TEM and spectroscopic techniques. TEM, however, has limitations that are important to understand when interpreting data including the ability of the electron beam to damage and thereby change the structure and shape of particles, especially in the case of particles composed of organic compounds and salts. In this paper, we concentrate on the breadth of studies that have used TEM as the primary analysis technique. Another focus is on common issues with TEM and cryogenic-TEM. Insights for new users on best practices for fragile particles, that is, particles that are easily susceptible to damage from the electron beam, with this technique are discussed. Tips for readers on interpreting and evaluating the quality and accuracy of TEM data in the literature are also provided and explained.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Biodiesel Ethers: Fatty Acid-Derived Alkyl Ether Fuels as Improved Bioblendstocks for Mixing-Controlled Compression Ignition Engines

In the last 20 years, biodiesel consumption in the United States has swiftly increased to ~2 billion gallons per year as a renewable supplement to fossil fuel. However, further expansion of biodiesel use is currently limited in part by poor cold weather performance, which prevents year-round blending and necessitates blend walls ≤5% v/v. In order to provide a diesel fuel blendstock with improved cold weather performance (cloud point, pour point, and cold filter plug point), while at the same time maintaining other required fuel performance specifications, several biodiesel redox analogues were synthesized and tested. The best performing candidate fuels from this class showed improvement in the derived cetane number (29.3% shorter ignition delay), lower heating value (+4.7 MJ/kg), relative sooting tendency (–7.4 YSI/MJ), and cloud point (15 °C lower) when compared to a B100 biodiesel composed of an identical fatty acid profile. It was observed as a general trend that the reduced form of biodiesel, fatty alkyl ethers (FAEs), shows performance improvements in all fuel property metrics. The suite of improved properties provided by FAEs gives biodiesel producers the opportunity to diversify their portfolio of products derived from lipid and alcohol feedstocks to include long-chain alkyl ethers, a biodiesel alternative with particular applicability for winter weather conditions across the US.

09 BIOMASS FUELS↗
DOE OSTI2021

Probing the Chemical Kinetics of Minimalist Functional Group Gasoline Surrogates

We report surrogate mixtures are routinely used for understanding gasoline fuel combustion in engine simulations. The general trend in surrogate formulation has been to increase the number of fuel components in a mixture to better emulate real fuel properties. Recently, a new surrogate design strategy based on functional group analysis of real gasolines was proposed using a minimal number of species [minimalist functional group (MFG)-approach]. MFG surrogates (having just one or two components) could experimentally capture the ignition delay time (IDT), threshold sooting index, and smoke point of different gasoline fuels with hundreds of components. However, other combustion characteristics were not explored, and kinetic modeling of MFG surrogates was not reported. These aspects are addressed in this paper, where the combustion behavior of MFG surrogates for various gasolines was assessed by simulating IDT, jet-stirred reactor oxidation, and premixed laminar flame speeds using chemical kinetic modeling. MFG simulations were compared with experimental data of the real gasolines as well as with the more complex multicomponent (five to nine species) surrogates. This study reveals that binary MFG surrogate mixtures are capable of accurately simulating the combustion behavior of more complex gasoline fuels with hundreds of components.

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 OSTI2021

Measurements of Intermediate Species in Fuel-Rich Oxidation of Ethylene, Toluene, and n -Decane

Polycyclic aromatic hydrocarbons (PAHs) are important precursors to the formation of soot and are also pollutant emissions from combustion devices, including internal combustion engines. To understand the formation of PAHs, the low-molecular-weight hydrocarbon chemistry leading to their formation needs to be understood. Toward this goal, fuel-rich oxidation of three different hydrocarbons, i.e., ethylene, toluene, and n-decane, has been investigated in an atmospheric pressure flow reactor at varying temperatures (1000–1350 K), equivalence ratios (φ = 3.0–12.0), and residence times of 0.25–1.5 s. The major C 1 to C 7 intermediates, such as methane, acetylene, ethylene, allene, propyne, propylene, vinylacetylene, 1,3-butadiene, cyclopentadiene, benzene, and toluene, were quantified using a gas chromatograph equipped with a flame ionization detector. The experimental tendency of intermediate species formation, e.g., the dependence of temperature, equivalence ratio, and residence time, was similar in ethylene oxidation and n-decane oxidation. The concentrations of intermediate species up to C 4 were higher in ethylene/n-decane oxidation than toluene oxidation, while a nearly equal or larger amount of C5–C7 species was produced in toluene oxidation. The experimental data were compared with modeling results using a detailed chemical kinetic mechanism. The calculated data using the kinetic model were in agreement with the experimental results. Additionally, a comprehensive kinetic analysis on the reaction pathways of each species was conducted to assess the differences in the oxidation chemistry with the change in the structure of hydrocarbons. In particular, in-depth analysis of benzene formation was performed, elucidating that two benzene formation pathways were important in ethylene and n-decane oxidation: (1) H elimination of the 2,4-cyclohexadienyl radical produced from the isomerization of the 2,4-cyclopentadienylmethyl radical, and (2) dehydrogenation of cyclohexadiene produced from reactions of vinyl radical with 1,3-butadiene. In toluene oxidation, it was found that benzene was primarily produced from toluene through the replacement of the methyl group with hydrogen.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

A Conical Intersection Influences the Ground State Rearrangement of Fulvene to Benzene

The rearrangement of fulvene to benzene is believed to play an important role in the formation of soot during hydrocarbon combustion. Previous work has identified two possible mechanisms for the rearrangement-a unimolecular path and a hydrogen-atom-assisted, bimolecular path. Computational results to date have suggested that the unimolecular mechanism faces a barrier of about 74 kcal/mol, which makes it unable to compete with the bimolecular mechanism under typical combustion conditions. This computed barrier is about 10 kcal/mol higher than the experimental value, which is an unusually large discrepancy for modern electronic structure theory. In the present work, we have reinvestigated the unimolecular mechanism computationally, and we have found a second transition state that is approximately 10 kcal/mol lower in energy than the previously identified one and, therefore, in excellent agreement with the experimental value. The existence of two transition states for the same rearrangement arises because there is a conical intersection between the two lowest singlet states which occurs in the vicinity of the reaction coordinates. The two possible paths around the cone on the lower adiabatic surface give rise to the two distinct saddle points. The lower barrier for the unimolecular mechanism now makes it competitive with the bimolecular one, according to our calculations. In support of this conclusion, we have reanalyzed some previous experimental results on anisole pyrolysis, which leads to benzene as a significant product and have shown that the unimolecular and bimolecular mechanisms for fulvene ? benzene must be occurring competitively in that system. Finally, we have identified that similar conical intersections arise during the isomerizations of benzofulvene and isobenzofulvene to naphthalene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Comprehensive Kinetics on the C 7 H 7 Potential Energy Surface under Combustion Conditions

The automated kinetics workflow code, KinBot, was used to explore and characterize the regions of the C 7 H 7 potential energy surface that are relevant to combustion environments and especially soot inception. We first explored the lowest-energy region, which includes the benzyl, fulvenallene + H, and cyclopentadienyl + acetylene entry points. We then expanded the model to include two higher-energy entry points, vinylpropargyl + acetylene and vinylacetylene + propargyl. The automated search was able to uncover the pathways from the literature. In addition, three important new routes were discovered: a lower-energy pathway connecting benzyl with vinylcyclopentadienyl, a decomposition mechanism from benzyl that results in side-chain hydrogen atom loss to produce fulvenallene + H, and shorter and lower energy routes to the dimethylene-cyclopentenyl intermediates. We systematically reduced the extended model to a chemically relevant domain composed of 63 wells, 10 bimolecular products, 87 barriers, and 1 barrierless channel and constructed a master equation using the CCSD(T)-F12a/cc-pVTZ//ωB97X-D/6-311++G(d,p) level of theory to provide rate coefficients for chemical modeling. Our calculated rate coefficients show excellent agreement with measured ones. We also simulated concentration profiles and calculated branching fractions from the important entry points to provide an interpretation of this important chemical landscape.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗