Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “REACTION”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Mode specificity in the OH + CHD{sub 3} reaction: Reduced-dimensional quantum and quasi-classical studies on an ab initio based full-dimensional potential energy surface

An initial state selected time-dependent wave packet method is applied to study the dynamics of the OH + CHD{sub 3} reaction with a six-dimensional model on a newly developed full-dimensional ab initio potential energy surface (PES). This quantum dynamical (QD) study is complemented by full-dimensional quasi-classical trajectory (QCT) calculations on the same PES. The QD results indicate that both translational energy and the excitation of the CH stretching mode significantly promote the reaction while the excitation of the umbrella mode has a negligible effect on the reactivity. For this early barrier reaction, interestingly, the CH stretching mode is more effective than translational energy in promoting the reaction except at very low collision energies. These QD observations are supported by QCT results. The higher efficacy of the CH stretching model in promoting this early barrier reaction is inconsistent with the prediction of the naively extended Polanyi’s rules, but can be rationalized by the recently proposed sudden vector projection model.

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

Ring-polymer molecular dynamics: Rate coefficient calculations for energetically symmetric (near thermoneutral) insertion reactions (X + H{sub 2}) → HX + H(X = C({sup 1}D), S({sup 1}D))

Following our previous study of prototypical insertion reactions of energetically asymmetric type with the RPMD (Ring-Polymer Molecular Dynamics) method [Y. Li, Y. Suleimanov, and H. Guo, J. Phys. Chem. Lett. 5, 700 (2014)], we extend it to two other prototypical insertion reactions with much less exothermicity (near thermoneutral), namely, X + H{sub 2} → HX + H where X = C({sup 1}D), S({sup 1}D), in order to assess the accuracy of this method for calculating thermal rate coefficients for this class of reactions. For both chemical reactions, RPMD displays remarkable accuracy and agreement with the previous quantum dynamic results that make it encouraging for the future application of the RPMD to other barrier-less, complex-forming reactions involving polyatomic reactants with any exothermicity.

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

THE C({sup 3}P) + NH{sub 3} REACTION IN INTERSTELLAR CHEMISTRY. I. INVESTIGATION OF THE PRODUCT FORMATION CHANNELS

The product formation channels of ground state carbon atoms, C({sup 3}P), reacting with ammonia, NH{sub 3}, have been investigated using two complementary experiments and electronic structure calculations. Reaction products are detected in a gas flow tube experiment (330 K, 4 Torr) using tunable vacuum-ultraviolet (VUV) photoionization coupled with time of flight mass spectrometry. Temporal profiles of the species formed and photoionization spectra are used to identify primary products of the C + NH{sub 3} reaction. In addition, H-atom formation is monitored by VUV laser induced fluorescence (LIF) from room temperature to 50 K in a supersonic gas flow generated by the Laval nozzle technique. Electronic structure calculations are performed to derive intermediates, transition states, and complexes formed along the reaction coordinate. The combination of photoionization and LIF experiments supported by theoretical calculations indicate that in the temperature and pressure range investigated, the H + H{sub 2}CN production channel represents 100% of the product yield for this reaction. Kinetics measurements of the title reaction down to 50 K and the effect of the new rate constants on interstellar nitrogen hydride abundances using a model of dense interstellar clouds are reported in Paper II.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Oxygen Storage Properties of La 1– x Sr x FeO 3-δ for Chemical-Looping Reactions—An In Situ Neutron and Synchrotron X-ray Study

Oxygen storage materials (OSMs) provide lattice oxygen for a number of chemical-looping reactions including natural gas combustion and methane reforming. La 1–x Sr x FeO 3-δ has shown promise for use as an OSM in methane reforming reactions due to its high product selectivity, fast oxide diffusion, and cycle stability. Here, we investigate the structural evolution of the series La 1–x Sr x FeO 3-δ for x = 0, 1/3, 1/2, 2/3, and 1, using in situ synchrotron X-ray and neutron diffraction, as it is cycled under the conditions of a chemical-looping reactor (methane and oxygen atmospheres). In the compositions x = 1/3, 1/2, 2/3, and 1, we discover an envelope , or temperature range, of oxygen storage capacity (OSC), where oxygen can easily and reversibly be inserted and removed from the OSM. Our in situ X-ray and neutron diffraction results reveal that while samples with higher Sr contents had a higher OSC, those same samples suffered from slower reaction kinetics and some, such as the x = 1/2 and x = 2/3 compositions, had local variations in Sr content, which led to inhomogeneous regions with varying reaction rates. Therefore, we highlight the importance of in situ diffraction studies, and we propose that these measurements are required for the thorough evaluation of future candidate OSMs. Furthermore, we recommend La 2/3 Sr 1/3 FeO 3-δ as the optimal OSM in the series because its structure remains homogeneous throughout the reaction, and its OSC envelope is similar to that of the higher doped materials.

Taylor, Daniel D.↗
DOE OSTI2016

Bond energies of ThO{sup +} and ThC{sup +}: A guided ion beam and quantum chemical investigation of the reactions of thorium cation with O{sub 2} and CO

Kinetic energy dependent reactions of Th{sup +} with O{sub 2} and CO are studied using a guided ion beam tandem mass spectrometer. The formation of ThO{sup +} in the reaction of Th{sup +} with O{sub 2} is observed to be exothermic and barrierless with a reaction efficiency at low energies of k/k{sub LGS} = 1.21 ± 0.24 similar to the efficiency observed in ion cyclotron resonance experiments. Formation of ThO{sup +} and ThC{sup +} in the reaction of Th{sup +} with CO is endothermic in both cases. The kinetic energy dependent cross sections for formation of these product ions were evaluated to determine 0 K bond dissociation energies (BDEs) of D{sub 0}(Th{sup +}–O) = 8.57 ± 0.14 eV and D{sub 0}(Th{sup +}–C) = 4.82 ± 0.29 eV. The present value of D{sub 0} (Th{sup +}–O) is within experimental uncertainty of previously reported experimental values, whereas this is the first report of D{sub 0} (Th{sup +}–C). Both BDEs are observed to be larger than those of their transition metal congeners, TiL{sup +}, ZrL{sup +}, and HfL{sup +} (L = O and C), believed to be a result of lanthanide contraction. Additionally, the reactions were explored by quantum chemical calculations, including a full Feller-Peterson-Dixon composite approach with correlation contributions up to coupled-cluster singles and doubles with iterative triples and quadruples (CCSDTQ) for ThC, ThC{sup +}, ThO, and ThO{sup +}, as well as more approximate CCSD with perturbative (triples) [CCSD(T)] calculations where a semi-empirical model was used to estimate spin-orbit energy contributions. Finally, the ThO{sup +} BDE is compared to other actinide (An) oxide cation BDEs and a simple model utilizing An{sup +} promotion energies to the reactive state is used to estimate AnO{sup +} and AnC{sup +} BDEs. For AnO{sup +}, this model yields predictions that are typically within experimental uncertainty and performs better than density functional theory calculations presented previously.

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

Experimental study of the Si 30 ( He 3 , d ) P 31 reaction and thermonuclear reaction rate of Si 30 ( p , γ ) P 31

Background: Abundance anomalies in some globular clusters, such as the enhancement of potassium and the depletion of magnesium, can be explained in terms of an earlier generation of stars polluting the presently observed ones. It was shown that the potential range of temperatures and densities of the polluting sites depends on the strength of a few critical reaction rates. The 30 Si(p,γ) 31 P reaction has been identified as one of these important reactions. Purpose: The key ingredient for evaluating the thermonuclear 30 Si(p,γ) 31 P reaction rate is the strength of the resonances which, at low energy, are proportional to their proton width. Therefore, the goal of this work is to determine the proton widths of unbound 31 P states. Method: States in 31 P were studied at the Maier-Leibnitz-Laboratorium using the one-proton 30 Si( 3 He,d) 31 P transfer reaction. Deuterons were detected with the high resolution one quadrupole three dipoles (Q3D) magnetic spectrometer. Angular distribution and spectroscopic factors were extracted for 27 states, and proton widths and resonance strengths were calculated for the unbound states. Results: Several 31 P unbound states have been observed for the first time in a one-proton transfer reaction. Above 20 MK, the 30 Si(p,γ) 31 P reaction rate is now entirely estimated from the observed properties of 31 P states. The reaction rate uncertainty from all resonances other than the E$^{c.m.}_{r}$ = 149 keV resonance has been reduced down to less than a factor of 2 above that temperature. The unknown spin and parity of the E$^{c.m.}_{r}$ = 149 keV resonance dominates the uncertainty in the rate in the relevant temperature range. Conclusion: Here, the remaining source of uncertainty on the 30 Si(p,γ) 31 P reaction rate comes from the unknown spin and parity of the E$^{c.m.}_{r}$ = 149 keV resonance which can change the reaction rate by a factor of 10 in the temperature range of interest.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Developing the S 32 ( p , d ) S * 31 ( p ) ( γ ) reaction to probe the P 30 ( p , γ ) S 31 reaction rate in classical novae

Background: The 30 P(p,γ) 31 S reaction rate is one of the largest remaining sources of uncertainty in the final abundances of nuclei created in a classical nova involving a ONe white dwarf. The reaction rate directly influences silicon isotopic ratios, which are used as identifiers of presolar grains with nova origins. Additionally, the uncertainty in the 30 P(p,γ) 31 S reaction rate has been found to limit the use of nova nuclear thermometers based on observations of elemental ratios in nova ejecta. Purpose: Reduce uncertainties in the nuclear data for proton-unbound states in 31 S, which act as resonances for the 30 P(p,γ) 31 S reaction at classical nova temperatures, and develop a technique for high efficiency, high-resolution reaction-decay coincidence measurements. Methods: Here, the 32 S(p,d) 31 S reaction was used to populate the states of interest in 31 S. The experiment was performed at the Texas A&M Cyclotron Institute using the LLNL Hyperion array for the detection of charged particles and γ rays. A downstream silicon telescope was used to select reaction deuterons, and a single upstream silicon detector was used to measure protons emitted in the decay of unbound 31 S levels. Results: Several states in 31 S above the proton separation energy were observed to have been populated. Decay protons from the resonant states in 31 S were identified as events in the upstream silicon detectors that came in coincidence with deuterons in the downstream telescope. Protons emitted from these states were measured and branching ratios extracted. Conclusions: While no new reaction rate is derived, spin-parity assignments for several higher-lying proton unbound states have been confirmed. Measured p 0 branching ratios for these levels have been compared to previous measurements with good agreement, and in some cases provided a reduction in uncertainty. The previously identified T = 3/2 state may have been incorrectly assigned a large p 0 branching ratio in a previous measurement. The technique of measuring reaction-decay coincidences with a particle-gamma setup appears promising.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Sensitivity to Thermonuclear Reaction Rates in Modeling the Abundance Anomalies of NGC 2419

Abundance anomalies in globular clusters provide strong evidence for multiple stellar populations within each cluster. These populations are usually interpreted as distinct generations, with the currently observed second-generation stars having formed in part from the ejecta of massive, first-generation “polluter” stars, giving rise to the anomalous abundance patterns. The precise nature of the polluters and their enrichment mechanism are still unclear. Even so, the chemical abundances measured in second-generation stars within the globular cluster NGC 2419 provide insight into this puzzling process. Previous work used Monte Carlo nuclear reaction network calculations to constrain the temperature–density conditions that could reproduce the observed abundances, thereby placing robust limits on the origins of the polluter material. The effect of individual reaction rates on these conditions has not been studied, however. Thus, we perform an exhaustive sensitivity study on the nuclear physics input to determine which reactions have the greatest impact on these predictions. We find that the {sup 30}Si(p,γ){sup 31}P, {sup 37}Ar(p,γ){sup 38}K, {sup 38}Ar(p,γ){sup 39}K, and {sup 39}K(p,γ){sup 40}Ca reactions are all critical in determining the temperature–density conditions, and ultimately, the origins of the polluter material. We conclude with recommendations for future experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The highly exothermic hydrogen abstraction reaction H 2 Te + OH → H 2 O + TeH: comparison with analogous reactions for H 2 Se and H 2 S

The “gold standard” CCSD(T) method is adopted along with the correlation consistent basis sets up to aug-cc-pV5Z-PP to study the mechanism of the hydrogen abstraction reaction H 2 Te + OH. Here, the predicted geometries and vibrational frequencies for reactants and products are in good agreement with the available experimental results. With the ZPVE corrections, the transition state in the favorable pathway of this reaction energetically lies 1.2 kcal mol -1 below the reactants, which is lower than the analogous relative energies for the H 2 Se + OH reaction (-0.7 kcal mol -1 ), the H 2 S + OH reaction (+0.8 kcal mol -1 ) and the H 2 O + OH reaction (+9.0 kcal mol -1 ). Accordingly, the exothermic reaction energies for these related reactions are predicted to be 47.8 (H 2 Te), 37.7 (H 2 Se), 27.1 (H 2 S), and 0.0 (H 2 O) kcal mol -1 , respectively. Geometrically, the low-lying reactant complexes for H 2 Te + OH and H 2 Se + OH are two-center three-electron hemibonded structures, whereas those for H 2 S + OH and H 2 O + OH are hydrogen-bonded. With ZPVE and spin–orbit coupling corrections, the relative energies for the reactant complex, transition state, product complex, and the products for the H 2 Te + OH reaction are estimated to be -13.1, -1.0, -52.0, and -52.6 kcal mol -1 , respectively. Finally, twenty-eight DFT functionals have been tested systematically to assess their ability in describing the potential energy surface of the H 2 Te + OH reaction. The best of these functionals for the corresponding energetics are -9.9, -1.4, -46.4, and -45.4 kcal mol -1 (MPWB1K), or -13.1, -2.4, -57.1, and -54.6 kcal mol -1 (M06-2X), respectively.

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

Active learning of chemical reaction networks via probabilistic graphical models and Boolean reaction circuits

Discerning networks of many reactions among multiple interconverting species is challenging. Here, we present a reaction network identification methodology. Our methodology enumerates all stoichiometrically and chemically feasible reactions and requires statistical evidence from effluent concentrations for the inclusion or exclusion of each from the reaction network, contrasting with the commonly seen incremental approach and other work of relying heavily upon chemical intuition and assuming the reactions occurring. Using graph theory alongside an active learning design of experiments that propose maximally informative feeds, we identify the underlying reaction network with minimal laboratory runs. Here, we introduce chemistry-probabilistic graphical modeling and Boolean reaction circuits to statistically quantify which reactions occur from effluent concentrations. Our methodology accurately discerns active reactions, as showcased upon a laboratory network of cross-ketonization of furoic and lauric acid and validated upon simulated networks of thermal and CO 2 -assisted ethane dehydrogenation.

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

Tracking the reaction networks of acetaldehyde oxide and glyoxal oxide Criegee intermediates in the ozone-assisted oxidation reaction of crotonaldehyde

The reaction of unsaturated compounds with ozone (O 3 ) is recognized to lead to the formation of Criegee intermediates (CIs), which play a key role in controlling the atmospheric budget of hydroxyl radicals and secondary organic aerosols. The reaction network of two CIs with different functionality, i.e. acetaldehyde oxide (CH 3 CHOO) and glyoxal oxide (CHOCHOO) formed in the ozone-assisted oxidation reaction of crotanaldehyde (CA), is investigated over a temperature range between 390 K and 840 K in an atmospheric pressure jet-stirred reactor (JSR) at a residence time of 1.3 s, stoichiometry of 0.5 with a mixture of 1% crotonaldehyde, 10% O 2 , at an fixed ozone concentration of 1000 ppm and 89% Ar dilution. Molecular-beam mass spectrometry in conjunction with single photon tunable synchrotron vacuum-ultraviolet (VUV) radiation is used to identify elusive intermediates by means of experimental photoionization energy scans and ab initio threshold energy calculations for isomer identification. Addition of ozone (1000 ppm) is observed to trigger the oxidation of CA already at 390 K, which is below the temperature where the oxidation reaction of CA was observed in the absence of ozone. The observed CA + O 3 product, C 4 H 6 O 4 , is found to be linked to a ketohydroperoxide (2-hydroperoxy-3-oxobutanal) resulting from the isomerization of the primary ozonide. Products corresponding to the CIs uni- and bi-molecular reactions were observed and identified. A network of CI reactions is identified in the temperature region below 600 K, characterized by CIs bimolecular reactions with species like aldehydes, i.e., formaldehyde, acetaldehyde, and crotonaldehyde and alkenes, i.e., ethene and propene. The region below 600 K is also characterized by the formation of important amounts of typical low-temperature oxidation products, such as hydrogen peroxide (H 2 O 2 ), methyl hydroperoxide (CH 3 OOH), and ethyl hydroperoxide (C 2 H 5 OOH). Detection of additional oxygenated species such as alcohols, ketene, and aldehydes are indicative of multiple active oxidation routes. This study provides important information about the initial step involved in the CIs assisted oligomerization reactions in complex reactive environments where CIs with different functionalities are reacting simultaneously. It provides new mechanistic insights into ozone-assisted oxidation reactions of unsaturated aldehydes, which is critical for the development of improved atmospheric and combustion kinetics models.

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

Reaction Rate Enhancement for Cu(In,Ga)Se 2 Absorber Materials Using Ag-Alloying

The addition of Ag to Cu-Ga-In precursors for synthesizing (Ag,Cu)(In,Ga)Se 2 (ACIGS) thin films has shown benefits including improved adhesion, greater process tolerance, and potential for improved device performance. In this study, reaction pathways to form Cu(In,Ga)Se 2 (CIGS) and ACIGS were studied by time-progressive reactions at 450 °C in a 5% Ar/H 2 Se atmosphere followed by ex situ characterization. Results indicated that the addition of 25% Ag/(Ag+Cu) to the CIGS film reduces the reaction time by 50%. X-ray diffraction (XRD) analysis of CIGS films showed that the CuInSe 2 phase initially formed after 3.5 min. The slow reaction of the stable γ-Cu 9 (In,Ga) 4 phase, however, required more than 20 min to complete. Importantly, the addition of Ag to the CIGS film accelerated the reaction. Energy-dispersive X-ray spectroscopy shows that Ga/(Ga+In) grading occurs in the first 10 min of the reaction. XRD analysis showed that the chalcopyrite phase fully forms after 10 min and no significant changes were observed in samples selenized from 10-45 min. Reaction pathways of Ag-alloyed films were further characterized using in situ high temperature XRD analysis. The onset temperature of Se reaction was detected at 230 °C and a AgIn 2 phase transformation to (Ag,Cu)In 2 occurred during the early stage of the reaction.

36 MATERIALS SCIENCE↗
NASA NTRS2006

Ab initio Quantum Chemical Studies of Reactions in Astrophysical Ices. Reactions Involving CH3OH, CO2, CO, HNCO in H2CO/NH3/H2O Ices

While reactions between closed shell molecules generally involve prohibitive barriers in the gas phase, prior experimental and theoretical studies have demonstrated that some of these reactions are significantly enhanced when confined within an icy grain mantle and can occur efficiently at temperatures below 100 K with no additional energy processing. The archetypal case is the reaction of formaldehyde (H2CO) and ammonia (NH3) to yield hydroxymethylamine (NH2CH2OH). In the present work we have characterized reactions involving methanol (CH3OH), carbon dioxide (CO2), carbon monoxide (CO), and isocyanic acid (HNCO) in search of other favorable cases. Most of the emphasis is on CH3OH, which was investigated in the two-body reaction with one H2CO and the three-body reaction with two H2CO molecules. The addition of a second H2CO to the product of the reaction between CH3OH and H2CO was also considered as an alternative route to longer polyoxymethylene polymers of the -CH2O- form. The reaction between HNCO and NH3 was studied to determine if it can compete against the barrierless charge transfer process that yields OCN(-) and NH4(+). Finally, the H2CO + NH3 reaction was revisited with additional benchmark calculations that confirm that little or no barrier is present when it occurs in ice.

Woon, David E.↗
NASA NTRS2013

Estimating the Backup Reaction Wheel Orientation Using Reaction Wheel Spin Rates Flight Telemetry from a Spacecraft

A report describes a model that estimates the orientation of the backup reaction wheel using the reaction wheel spin rates telemetry from a spacecraft. Attitude control via the reaction wheel assembly (RWA) onboard a spacecraft uses three reaction wheels (one wheel per axis) and a backup to accommodate any wheel degradation throughout the course of the mission. The spacecraft dynamics prediction depends upon the correct knowledge of the reaction wheel orientations. Thus, it is vital to determine the actual orientation of the reaction wheels such that the correct spacecraft dynamics can be predicted. The conservation of angular momentum is used to estimate the orientation of the backup reaction wheel from the prime and backup reaction wheel spin rates data. The method is applied in estimating the orientation of the backup wheel onboard the Cassini spacecraft. The flight telemetry from the March 2011 prime and backup RWA swap activity on Cassini is used to obtain the best estimate for the backup reaction wheel orientation.

Rizvi, Farheen↗
DOE OSTI2011

High-precision {sup 28}Si(p,t){sup 26}Si reaction to determine {sup 22}Mg({alpha},p){sup 25}Al reaction rates

The rise time of stellar x-ray bursts is a signature of thermonuclear runaway processes in the atmosphere of neutron stars and is highly sensitive to a series of ({alpha},p) reactions via high-lying resonances in sd-shell nuclei. Lacking data for the relevant resonance levels, the stellar reaction rates have been calculated using statistical, Hauser-Feshbach models, assuming a high-level density. This assumption may not be correct in view of the selectivity of the ({alpha},p) reaction to natural parity states. We measured the {sup 28}Si(p,t){sup 26}Si reaction with a high-resolution spectrometer to identify resonance levels in {sup 26}Si above the {alpha}-emission threshold at 9.164 MeV excitation energy. These resonance levels are used to calculate the stellar reaction rate of the {sup 22}Mg({alpha},p){sup 25}Al reaction and to test the validity of the statistical assumption.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

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↗