Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy intermediate”

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 883 records · Page 49

Advancing technologies for lignin-based jet fuel production in aqueous phase

Integrating lignin into a cellulosic ethanol plant for the co-production of lignin-based jet fuel (LJF) in aqueous phase offers a significant opportunity to boost operational efficiency, economic viability, carbon conversion, and the overall sustainability of biofuel and chemical production. LJF is lignin-structure-based jet fuel blendstocks primarily composed of alkyl-substituted mono-, bi-, and tri-cyclohexanes. It exhibits high energy density, potential for low emissions, and favourable blend characteristics that comply with drop-in specifications. An overview of lignin feedstock, catalytic processes, LJF chemical compositions, fuel properties tests, and techno-economic analysis (TEA) and life cycle assessment (LCA) indicate that (1) the reactivity of lignin plays a crucial role in its structure transformation to LJF molecules; (2) catalytic processing of lignin to LJF can occur through a simultaneous depolymerization and hydrodeoxygenation process, bypassing the intermediate step of producing and upgrading lignin-derived oil; (3) LJF's uniqueness molecules making it more promising for high energy content and low emission jet fuel properties for next generation sustainable aviation fuel (SAF); and (4) TEA and LCA demonstrate that LJF is not only potentially cost-effective but also offers favourable carbon footprint compared to other SAFs. In conclusion, this review highlights the most recent advancements in LJF technology, along with the challenges and opportunities that lie ahead in fulfilling its potential.

09 BIOMASS FUELS↗

Single-Step Mechanism for Regioselective Nitration of 9,10-BN-Naphthalene with Acetyl Nitrate in the Gas Phase

Here, the energetics of the regioselective mononitration of 9,10-BN-naphthalene with acetyl nitrate (H 3 C 2 NO 4 ) were modeled with ab initio simulations in the gas phase and an acetonitrile solvent. The single-electron-transfer (SET) nitration mechanism leading to a σ-complex and a single-step nitration mechanism were modeled. The energy barrier for the single-step mechanism was lower than that for the SET mechanism in the gas phase. However, the two are much more energetically competitive in the solvent. The σ-complex was found to be unstable in the gas phase owing to the interaction with the counterion. Using the single-step mechanism, the carbon site 1 nearest boron had the lowest activation energy for nitration of 22.6 kcal/mol, while site 3 had the second lowest barrier of 24.6 kcal/mol. Finally, details on the molecular structures at intermediate and transition states as well as charges in different configurations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Charge transfer as a mechanism for chlorophyll fluorescence concentration quenching

Highly concentrated solutions of chlorophyll display rapid fluorescence quenching. The same devastating energy loss is not seen in photosynthetic light-harvesting antenna complexes, despite the need for chromophores to be in close proximity to facilitate energy transfer. A promising, though unconfirmed mechanism for the observed quenching is energy transfer from an excited chlorophyll monomer to a closely associated chlorophyll pair that subsequently undergoes rapid nonradiative decay to the ground state via a short-lived intermediate charge-transfer state. In this work, we make use of newly emerging fast methods in quantum chemistry to assess the feasibility of this proposed mechanism. We calculate rate constants for the initial charge separation, based on Marcus free-energy surfaces extracted from molecular dynamics simulations of solvated chlorophyll pairs, demonstrating that this pathway will compete with fluorescence (i.e., drive quenching) at experimentally measured quenching concentrations. We show that the rate of charge separation is highly sensitive to interchlorophyll distance and the relative orientations of chromophores within a quenching pair. We discuss possible solvent effects on the rate of charge separation (and consequently the degree of quenching), using the light-harvesting complex II (LH2) protein from rps. acidophila as a specific example of how this process might be controlled in a protein environment. Crucially, we reveal that the LH2 antenna protein prevents quenching, even at the high chlorophyll concentrations required for efficient energy transfer, by restricting the range of orientations that neighboring chlorophyll pairs can adopt.

59 BASIC BIOLOGICAL SCIENCES↗

Electrolyzers for Hydrogen Production: Solid Oxide, Alkaline, and Proton Exchange Membrane

The effects of climate change have led to a push for cleaner energy sources across multiple sectors. Hydrogen has emerged as a promising energy carrier in this context, as it can be produced using various water electrolysis technologies capable of utilizing clean power (e.g., nuclear and renewable electricity). However, a comprehensive environmental assessment of these technologies requires an understanding of the environmental impacts during their complete life cycle, including the embodied emissions in their material composition and during their manufacturing stages; these emissions are often neglected. This report provides a brief overview of major water electrolysis technologies, viz., proton exchange membrane electrolyzer cell (PEMEC) or polymer electrolyte membrane (PEM), alkaline electrolysis cell (AEC) and solid oxide electrolysis cell (SOEC), along with a detailed bill of materials for these technologies that has been incorporated in the updated G reenhouse gases, R egulated E missions, and E nergy use in T ransportation (GREET ® ) 2022 model. We also provide an inventory for the intermediate materials used to produce these electrolyzers, which has not been covered in prior releases. Using these material and energy flows, the GREET model provides a detailed life cycle inventory (energy use and emissions) from raw material extraction through complete production of electrolyzers for major electrolysis technologies.

08 HYDROGEN↗

Effect of reactive intermediates on the propagation of turbulent premixed flames

The most simple version of the Bray and Libby model of turbulent flame propagation has been extended to include a two-step reaction mechanism in which reactants decompose into reactive intermediates which then form products. A transport equation for the joint pdf of intermediate and product concentrations was derived, the pdf parameterized, and transport equations for the parameters obtained. The one-dimensional concentration equations were then solved for the parameters of the pdf and the turbulent kinetic energy. Flame speed and angle are recovered as eigenvalues of the problem. The results show the essential features observed in experiments. The product pdf is composed of peaks at zero and unity product concentration, with a continuous structure representing partially reacted gas. The intermediates pdf has a peak at zero concentration and a broad bell-shaped distribution at nonzero concentration, representing partially reacted gas. The results also show the effect of partitioning of energy release on the flame thickness and speed.

Azzazy, M.↗

Resonant inelastic x-ray scattering investigation of Hund's and spin-orbit coupling in 5 d 2 double perovskites

𝐵 site ordered 5⁢𝑑 2 double perovskites (𝐴 2 ⁢𝐵⁢𝐵′⁢O 6 , 𝐵′ = 5⁢𝑑 2 ) display a remarkable range of physical properties upon variation of the chosen 𝐵 and 𝐵′ site ions. This sensitivity to chemical substitution reflects the delicate balance and profound impact of strong electronic correlation and spin-orbit coupling in such systems. Here, we present rhenium 𝐿 2 and 𝐿 3 resonant inelastic x-ray scattering (RIXS) measurements of two such physically dissimilar materials, Mott-insultating Ba 2⁢ Y⁢Re⁢O 6 and semiconducting Sr 2⁢ Cr⁢Re⁢O 6 . Despite these differences, our RIXS results reveal similar energy scales of Hund's (𝐽 𝐻 ) and spin-orbit coupling (𝜁) in the two materials, with both systems firmly in the intermediate Hund's coupling regime where 𝐽 𝐻 /𝜁∼1. However, there are clear differences in their RIXS spectra. The conductive character of Sr 2⁢ Cr⁢Re⁢O 6 broadens and obfuscates the atomic transitions within an electron-hole continuum, while the insulating character of Ba 2 ⁢Y⁢Re⁢O 6 results in sharp atomic excitations. This contrast in their RIXS spectra despite their similar energy scales reflects a difference in the itinerancy-promoting hopping integral and illustrates the impact of the local crystal environment in double perovskites. Finally, 𝐿 2 and 𝐿 3 edge analyses of the atomic excitations in Ba 2 ⁢Y⁢Re⁢O 6 reveal that the ordering of the low lying excited states is inverted compared to previous reports, such that the appropriate energy scales of Hund's and spin-orbit coupling are significantly modified. We present exact diagonalization calculations of the RIXS spectra at both edges which show good agreement with our results for new energy scales of 𝜁 = 0.290⁢(5)⁢eV and 𝐽 𝐻 = 0.38⁢(2)⁢eV [𝐽 𝐻 /𝜁 = 1.30⁢(5)].

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Designing electrolytes with high solubility of sulfides/disulfides for high-energy-density and low-cost K-Na/S batteries

Alkaline metal sulfur (AMS) batteries offer a promising solution for grid-level energy storage due to their low cost and long cycle life. However, the formation of solid compounds such as M 2 S 2 and M 2 S (M = Na, K) during cycling limits their performance. Here we unveil intermediate-temperature K-Na/S batteries utilizing advanced electrolytes that dissolve all polysulfides and sulfides (K 2 S x , x = 1–8), significantly enhancing reaction kinetics, specific capacity, and energy density. These batteries achieve near-theoretical capacity (1655 mAh g -1 sulfur) at 75 °C with a 1 M sulfur concentration. At a 4 M sulfur concentration, they deliver 830 mAh g -1 at 2 mA cm -2 , retaining 71% capacity after 1000 cycles. This new K-Na/S battery with specific energy of 150-250 Wh kg -1 only employs earth-abundant elements, making it attractive for long-duration energy storage.

36 MATERIALS SCIENCE↗

Creating Frustrated Lewis Pairs in Defective Boron Carbon Nitride for Electrocatalytic Nitrogen Reduction to Ammonia

Abstract The electrocatalytic nitrogen reduction reaction (NRR) on metal‐free catalysts is an attractive alternative to the industrial Haber–Bosch process. However, the state‐of‐the‐art metal‐free electrocatalysts still suffer from low Faraday efficiencies and low ammonia yields. Herein, we present a molecular design strategy to develop a defective boron carbon nitride (BCN) catalyst with the abundant unsaturated B and N atoms as Lewis acid and base sites, which upgrades the catalyst from a single “Lewis acid catalysis” to “frustrated Lewis pairs (FLPs) catalysis.” 14 N 2 / 15 N 2 exchange experiments and density functional theory (DFT) calculations reveal that FLPs can adsorb an N 2 molecule to form a six‐membered ring intermediate, which enables the cleavage of N 2 via a pull–pull effect, thereby significantly reducing the energy barrier to −0.28 eV. Impressively, BCN achieves a high Faraday efficiency of 18.9 %, an ammonia yield of 20.9 μg h −1 mg −1 cat. , and long‐term durability.

Lin, Wenwen↗

Creating Frustrated Lewis Pairs in Defective Boron Carbon Nitride for Electrocatalytic Nitrogen Reduction to Ammonia

Abstract The electrocatalytic nitrogen reduction reaction (NRR) on metal‐free catalysts is an attractive alternative to the industrial Haber–Bosch process. However, the state‐of‐the‐art metal‐free electrocatalysts still suffer from low Faraday efficiencies and low ammonia yields. Herein, we present a molecular design strategy to develop a defective boron carbon nitride (BCN) catalyst with the abundant unsaturated B and N atoms as Lewis acid and base sites, which upgrades the catalyst from a single “Lewis acid catalysis” to “frustrated Lewis pairs (FLPs) catalysis.” 14 N 2 / 15 N 2 exchange experiments and density functional theory (DFT) calculations reveal that FLPs can adsorb an N 2 molecule to form a six‐membered ring intermediate, which enables the cleavage of N 2 via a pull–pull effect, thereby significantly reducing the energy barrier to −0.28 eV. Impressively, BCN achieves a high Faraday efficiency of 18.9 %, an ammonia yield of 20.9 μg h −1 mg −1 cat. , and long‐term durability.

Lin, Wenwen↗

Charge-Transfer States Enable Spin-Selective Formation of Quartet State Qudits in Luminescent Radical-Chromophore Dyads

Luminescent tris(2,4,6-trichlorophenyl)methyl (TTM) radicals are promising doublet emitters that have been used to generate spin-optical interfaces in molecular electron spin qudits. In particular, photoexcitation of covalent TTM-chromophore dyads can generate quartet spin states that can be detected via photoluminescence. However, the mechanism of quartet spin initialization is complicated by competing charge-transfer (CT) and energy-transfer (EnT) pathways between the radical and the chromophore. Here, we demonstrate the role of CT intermediates by engineering a covalent dyad of TTM and naphthalene-(1,4:5,8)-bis(dicarboximide) (NDI), a well-known electron acceptor chromophore. Transient absorption spectroscopy indicates that photoexcitation of TTM results in ultrafast electron transfer from 2* TTM to NDI, which kinetically outcompetes EnT. Electron paramagnetic resonance spectroscopy reveals that charge recombination proceeds via spin–orbit charge-transfer intersystem crossing to generate 3* NDI, followed by spin mixing with 2TTM to form the quartet state. Furthermore, this process uniquely spin-polarizes the quartet state, providing fundamental design principles to purify the initial wave function for quantum information processing.

Charge transfer↗

Observation of partonic flow in proton—proton and proton—nucleus collisions

Quantum Chromodynamics predicts a phase transition from hadronic matter to quark–gluon plasma (QGP) at high temperatures and energy densities, where quarks and gluons (partons) are no longer confined within hadrons. The QGP forms in ultrarelativistic heavy-ion collisions. Anisotropic flow coefficients, quantifying the azimuthal expansion of produced matter, probe QGP properties. Flow measurements in high-energy heavy-ion collisions show a distinctive grouping of anisotropic flow for baryons and mesons at intermediate transverse momentum – a feature associated with flow imparted at the quark level, confirming QGP existence. The observation of QGP-like features in proton–proton and proton–ion collisions has sparked debate about QGP formation in smaller systems. For the first time, we demonstrate the distinctive grouping of anisotropic flow for baryons and mesons in high-multiplicity proton–lead and proton–proton collisions at the Large Hadron Collider (LHC). These results are described by a model including hydrodynamic flow followed by hadron formation via quark coalescence, consistent with the formation of partonic flowing systems in these collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Photochemical phosphorus-enabled scaffold remodeling of carboxylic acids

The excitation of carbonyl compounds by light to generate radical intermediates has historically been restricted to ketones and aldehydes; carboxylic acids have been overlooked because of high energy requirements and low quantum efficiency. A successful activation strategy would necessitate a bathochromic shift in the absorbance profile, an increase in triplet diradical lifetime, and ease of further functionalization. We present a single-flask transformation of carboxylic acids to acyl phosphonates that can access synthetically useful triplet diradicals under visible light or near-ultraviolet irradiation. The use of phosphorus circumvents unproductive Norrish type I processes, promoting selectivity that enables hydrogen-atom transfer reactivity. Use of this strategy promotes the efficient scaffold remodeling of carboxylic acids through various annulation, contraction, and expansion manifolds.

Science & Technology - Other Topics↗

Predissociation in N2(C'4, 1 Sigma u +) and other N2 states and its importance in the atmospheres of Titan and Triton

The objectives of this program are to further the understanding of the upper atmospheres of Titan, Triton, and the Earth in terms of the observed emissions of the 13-14 eV states of N2. These states are generated at quite high rates, yet very little emission is observed from them. The reasons are complex, involving resonance trapping and predissociation, and it is desired to quantify the effects of predissociation, particularly on the c(sub 4)' 1 Sigma(sub u),(sup +) state of N2. Earlier experiments had indicated that predissociation of the c(sub 4)' state was of little importance, yet over the last two years a growing body of evidence has shown that for levels above v = 2, predissociation is in fact a major process. It is the v = 0 level for which production by electron bombardment and photoexcitation is highest, and so it has been most important to evaluate the effects of predissociation on this particular level. The goal has been to target c(sub 4)' (v = 0) for a thorough analysis, in which both the extent of predissociation as a function of rotational level and the atomic product branching ratio, where the only possible products are N(4S) + N(4S) and N(2D) + N(4S), are determined. For the first year of funding, the intention was to demonstrate two-photon excitation of the intermediate N2(a(sup 1) Pi(sub g)) state, so that the gap to the 13 eV energy region could be bridged, and then use a second laser to reach the c(sub 4)' state itself.

Slanger, Tom G.↗

Overview of the Predictive Simulation Capability Element of the Plume Surface Interaction Project

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The GCD PSI Project, planned to be completed over a four year period, contains a Predictive Simulation Capability (PSC) Element focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the PSC Element, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper describes the PSC Element of the PSI Project as well as providing descriptions of recent accomplishments and remaining work. The overall structure of the PSC Element is broken down into four areas of focus. The first area is the Prediction of Plume Flow in low pressure environments. The development approach taken is the augmentation of the existing production-mode computational fluid dynamics (CFD) tool Loci/Chem, with targeted extensions necessary to accurately model rarefied conditions found in both Martian and Lunar applications. Production readiness and validation of predictive capability are the major objectives of this task. The second area is the effect of mixed continuum/rarefied flow on crater development and ejecta sheets. A new CFD application, Loci/GGFS (Gas Granular Flow Solver), is being developed which implements an Eulerian/Eulerian two-phase model of gas- and soil-phases in order to simulate the soil erosion, crater formation, and soil ejecta transport in a fully coupled simulation. This task seeks to verify Loci/GGFS is production-ready as well as perform validation studies to determine the degree of predictive capability achieved by Loci/GGFS. The third area is focused on the details and extension of particle phase modeling of soil. In this task, Discrete Element Modeling (DEM) techniques are used to perform direct simulations of complex soil particles under the action of forcing similar to that to be cause by PSI. The simulation results are then used to construct closures to the Eulerian model of the soil phase used by Loci/GGFS. The fourth area is gas-particle interaction modeling. In this task, experiments are being conducted as well as detailed simulation results are being studied to further understand the complexities of gas-particle interactions in dilute, intermediate, and high soil volume fraction regimes. Improved models of particle drag and the particle turbulent kinetic energy (PTKE) resulting from the interaction of gas flows within particle clouds are the objective of this task.

Jeff West↗

Cobalt-Group 13 Complexes Catalyze CO 2 Hydrogenation via a Co(-I)/Co(I) Redox Cycle

The Co(-I) dihydrogen complexes, [(η 2 -H 2 )CoML]-, where ML is the group 13 metalloligand, N(o-(NCH 2 PiPr 2 )C 6 H 4 ) 3 M, and M is Al, Ga, or In, were previously reported ( J. Am. Chem. Soc. 2017, 139, 6570-6573). In this work, the related Co(-I) end-on dinitrogen adducts, [(N 2 )CoML]-, were isolated and investigated as precatalysts for CO 2 hydrogenation. The Co–Ga catalyst was highly active, achieving 19,200 formate turnovers with an initial turnover frequency of 27,000 h –1 under 34 atm of 1:1 CO 2 /H 2 and using Verkade’s proazaphosphatrane as a base at ambient temperature. The Co–Al catalyst was moderately active, while the Co–In complex was inactive. Hence, tuning the group 13 ion greatly influences the catalytic activity at the Co site. To elucidate the role of the group 13 support, experimental and theoretical mechanistic studies of the Co–Ga and Co–Al catalysts were conducted. The Co(-I) H 2 species are potent hydride donors with estimated thermodynamic hydricities (ΔG° H– ) of 32.0(1) and 37.4(1) kcal/mol in CH 3 CN for M = Al and Ga, respectively. By acting as masked Co(I) dihydrides, the Co(-I) H 2 species operate via an unusual Co(-I)/Co(I) redox cycle. After hydride transfer to CO 2 , the resulting intermediate is the Co(I) hydride complex, HCoML, which was independently synthesized and structurally characterized for M = Al and Ga. The Gibbs free energy for H 2 binding, ΔG° bind (1 atm), to generate (η 2 -H 2 )HCoML was slightly more favorable for HCoGaL (-4.2(1) kcal/mol) than for HCoAlL (-2.7(1) kcal/mol). In the subsequent step, the deprotonation reaction to regenerate the initial catalyst was much more favorable for (η 2 -H 2 )HCoGaL (pK a of 31.4, CH 3 CN) than for (η 2 -H 2 )HCoAlL (pK a of 34.3). Finally, the straightforward substitution of Al with Ga perturbs the energy profile of the catalytic reaction (|ΔΔG° H– | = 5.4 kcal/mol, |ΔΔG° bind | = 1.5 kcal/mol, and |ΔΔG° K a | = 4.0 kcal/mol) and thus provides a thermodynamic rationale for the higher catalytic efficiency of Co–Ga over Co–Al.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reactive Molecular Dynamics Simulations and Quantum Chemistry Calculations To Investigate Soot-Relevant Reaction Pathways for Hexylamine Isomers

Sooting tendencies of a series of nitrogen-containing hydrocarbons (NHCs) have been recently characterized experimentally using the yield sooting index (YSI) methodology. This work aims to identify soot-relevant reaction pathways for three selected C6H15N amines, namely, dipropylamine (DPA), diisopropylamine (DIPA), and 3,3-dimethylbutylamine (DMBA) using ReaxFF molecular dynamics (MD) simulations and quantum mechanical (QM) calculations and to interpret the experimentally observed trends. ReaxFF MD simulations are performed to determine the important intermediate species and radicals involved in the fuel decomposition and soot formation processes. QM calculations are employed to extensively search for chemical reactions involving these species and radicals based on the ReaxFF MD results and also to quantitatively characterize the potential energy surfaces. Specifically, ReaxFF simulations are carried out in the NVT ensemble at 1400, 1600, and 1800 K, where soot has been identified to form in the YSI experiment. These simulations account for the interactions among test fuel molecules and pre-existing radicals and intermediate species generated from rich methane combustion, using a recently proposed simulation framework. ReaxFF simulations predict that the reactivity of the amines decrease in the order DIPA > DPA > DMBA, independent of temperature. Both QM calculations and ReaxFF simulations predict that C2H4, C3H6, and C4H8 are the main nonaromatic soot precursors formed during the decomposition of DPA, DIPA, and DMBA, respectively, and the associated reaction pathways are identified for each amine. Both theoretical methods predict that sooting tendency increases in the order DPA, DIPA, and DMBA, consistent with the experimentally measured trend in YSI. This work demonstrates that sooting tendencies and soot-relevant reaction pathways of fuels with unknown chemical kinetics can be identified efficiently through combined ReaxFF and QM simulations. Overall, predictions from ReaxFF simulations and QM calculations are consistent, in terms of fuel reactivity, major intermediates, and major nonaromatic soot precursors.

fuel decomposition↗