Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “radically”

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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Frhodo: A program for simulating chemical kinetic measurements and optimizing kinetic mechanisms

Chemical kinetic simulations are frequently used to extract kinetic and mechanistic information from experimental data by fitting simulations to the data. For simple systems of consecutive single channel reactions this is a relatively straightforward task. However, as the complexity of the system increases and particularly for systems involving competitive multi-channel reactions manually optimizing a mechanism to obtain the best fits to a range of experimental data quickly becomes a time-intensive and challenging process, especially when rate coefficients have pressure and temperature dependencies. There is considerable potential for exploiting automated optimization methods to rapidly screen reaction mechanisms and optimize the fit by adjusting rate coefficients within well-defined constraints. A new chemical kinetics simulation program, Frhodo, has been developed for this purpose. The program allows either manual simulation of individual experiments or fully automated optimization against a range of experimental data by adjusting user selected reactions within user defined constraints. Frhodo incorporates machine learning-based optimization routines that work by either minimizing standardized residuals or through Bayesian parameter estimation. Frhodo's optimization capabilities are demonstrated by reexamining two previously studied systems of relevance to combustion. Dissociation of diacetyl followed by recombination of methyl radicals is an example of a sequence of single channel, consecutive reactions. Pyrolysis of 2-methyl furan exemplifies multi-channel, unimolecular reactions that are often encountered and have temperature and pressure dependencies in the rate coefficients of the channels and branching ratios between them. For both systems the optimization routines resulted in solutions similar to the original studies. Finally, comments are made on the strengths and limitations of the approach.

0D simulation↗
DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 2

A wide variety of high temperature experimental data obtained in this study complement the data on the oxidation of the two di-isobutylene isomers presented in Part I and offers a basis for an extensive validation of the kinetic model developed in this study. Due to the increasing importance of unimolecular decomposition reactions in high-temperature combustion, we have investigated the di-isobutylene isomers in high dilution utilizing a pyrolysis microflow reactor and detected radical intermediates and stable products using vacuum ultraviolet (VUV) synchrotron radiation and photoelectron photoion coincidence (PEPICO) spectroscopy. Additional speciation data at oxidative conditions were also recorded utilizing a plug flow reactor at atmospheric pressure in the temperature range 725-1150 K at equivalence ratios of 1.0 and 3.0 and at residence times of 0.35 s and 0.22 s, respectively. Combustion products were analyzed using gas chromatography (GC) and mass spectrometry (MS). Ignition delay time measurements for di-isobutylene were performed at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in 'air' in the temperature range 900-1400 K using a high-pressure shock-tube facility. New measurements of the laminar burning velocities of di-isobutylene/air flames are also presented. The experiments were performed using the heat flux method at atmospheric pressure and initial temperatures of 298-358 K. Moreover, data consistency was assessed with the help of analysis of the temperature and pressure dependencies of laminar burning velocity measurements, which was interpreted using an empirical power-law expression. Electronic structure calculations were performed to compute the energy barriers to the formation of many of the product species formed. The predictions of the present mechanism were found to be in adequate agreement with the wide variety of experimental measurements performed.

09 BIOMASS FUELS↗
DOE OSTI2022

A wide range experimental study and further development of a kinetic model describing propane oxidation

This work presents an experimental and kinetic modeling study of propane oxidation. Ignition delay times of propane were measured in a high-pressure shock tube and in rapid compression machines in the temperature range 689 - 1700 K at equivalence ratios of 0.5, 1.0 and 2.0 in 'air', for a wide range of pressures from 20 to 90 bar. CO and H 2 O mole fraction profiles for propane oxidation were measured in a shock tube behind reflected shock waves in the temperature range 1370-1840 K at equivalence ratios of 0.5, 1.0 and 2.0 and at a pressure of approximately 1.3 atm. Moreover, propane oxidation was studied using a jet-stirred reactor coupled to a synchrotron vacuum ultraviolet photoionization mass spectrometer at low temperatures in the range 565 - 690 K and at a pressure of 1 atm. Further, this wide range of experimental datasets for propane oxidation was used to reoptimize and update our previous kinetic mechanisms, AramcoMech3.0 and NUIGMech1.1. In the current mechanism, NUIGMech1.3, the thermo-chemical parameters of all species relevant to low-temperature propane oxidation chemistry, including propyl-peroxyl, hydroperoxyl-propyl, hydroperoxyl-propyl-peroxyl, and carbonyl-hydroperoxide radicals, are updated based on newly calculated values at the CCSD(T)-F12/TZ-F12//B2PLYPD3/TZ///B2PLYP-D3/TZ level of theory. The improvements made in the thermochemical values and in the kinetic parameters for the low-temperature propane oxidation reactions in NUIGMech1.3 result in better model agreement with the new IDTs and speciation data, including carbon monoxide, formaldehyde, propene, acetaldehyde and various minor products such as ethylene, acetic acid, acrolein as well as various hydroperoxide and cyclic ether species.

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

In-situ two-dimensional temperature measurements using x-ray fluorescence spectroscopy in laminar flames with high silica particle concentrations

X-ray fluorescence spectroscopy (XRF) was used to measure temperatures and study mixing, for the first time in silica particle synthesis flames. Hexamethyldisiloxane (HMDSO) and trimethylsilanol (TMSO) were the particle precursors. A multi-element diffusion burner was used to produce a flat methane flame, and the precursors, dilute in inert gas, were injected via a central jet. Krypton was the fluorescent medium at 3.2 % concentration by volume. Scans with Kr in the central flow and not in the main flow were made to assess the mixing effects between the central and main gas flows. When HMDSO or TMSO were added, a secondary diffusion flame formed between the jet and the main methane flame. The results revealed a dramatic change in the centerline temperature profile of the jet gases when HMDSO or TMSO were added. The main methane flame stoichiometry also affected the temperature profiles. The results show HMDSO and TMSO reactions are initiated in a lowtemperature and low-oxygen concentration region of the jet where thermal decomposition is not expected to be significant. Reaction of the particle precursors is therefore attributed to radical transport from the main methane flame. In the current work, particle number densities of up to 270 g/m 3 locally are estimated. Thus, the study also demonstrates the capability of the XRF technique for high spatial fidelity measurements in flames with high concentrations of condensed-phase particles, leveraging the attribute that the XRF signal is generally not impacted by condensed-phase interferences. In conclusion, the observations and data obtained in this study inform likely reaction pathways for this important class of siloxane compounds.

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

O 2 -Dependence of reactions of 1,2-dimethoxyethanyl and 1,2-dimethoxyethanylperoxy isomers

Reaction mechanisms of Ṙ and ROȮ radicals derived from low-temperature oxidation of 1,2-dimethoxyethane (CH 3 O(CH 2 ) 2 OCH 3 ) were investigated using speciation from multiplexed photoionization mass spectrometry (MPIMS) measurements via Cl-initiated oxidation, in conjunction with electronic structure calculations. The experiments were conducted at 5 bar, from 450 K – 650 K, and O 2 concentrations from 1 · 10 14 cm –3 – 6 · 10 18 cm –3 to probe the effects on competing reaction channels of 1,2-dimethoxyethanyl (Ṙ) and 1,2-dimethoxyethanylperoxy (ROȮ) isomers. Several species were detected with photoionization spectral fitting – ethene, formaldehyde, methyl vinyl ether, and 2-methoxyacetaldehyde – and, as determined by electronic structure calculations, may form via unimolecular decomposition of 1,2-dimethoxyethanyl or 1,2-dimethoxyethanylperoxy. O 2 -dependent yield ratios show that the formation pathways for all species undergo a competition between O 2 -addition and unimolecular decomposition. Here, adiabatic ionization energies were also calculated and utilized along with exact mass determinations to infer contributions for other species derived exclusively from first- and second-O 2 -addition, including 1,2-dimethoxyethene, cyclic ethers, and dicarbonyls.

1,2-dimethoxyethane↗
DOE OSTI2025

Time-resolved measurements of OH during auto-ignition of syngas with trimethylsilanol and hexamethyldisiloxane

The effects of trimethylsilanol (TMSO) and hexamethyldisiloxane (HMDSO) addition on OH time histories during syngas (H 2 and CO) ignition were investigated using the University of Michigan rapid compression facility. Experiments spanned temperatures of 1010–1080 K, at a pressure of approximately 5 atm. Syngas mixtures of 1.2 % H 2 /2.8 % CO/20 % O 2 by volume (balance N 2 and Ar) provided a baseline for comparison with mixtures that included 100, 200, and 1000 ppm of the TMSO and 100 ppm of HMDSO. Narrow-line ultraviolet laser-absorption was used to measure OH mole-fraction during ignition. The addition of TMSO and HMDSO significantly shifted the OH time-histories earlier in time, by up to 51 %, compared with the baseline syngas mixture. The value of the maximum OH mole fraction was consistent between the 100 and 200 ppm TMSO mixtures and the 100 ppm HMDSO mixtures, but the maximum OH increased significantly with the 1000 ppm TMSO mixtures. Here, the OH data indicate TMSO and HMDSO were not direct sources of OH radicals. Analysis further indicates the TMSO and HMDSO decompose rapidly followed by reactions that enhance the production of H atoms, and the increased reactivity observed is via the H + O 2 = OH + O reaction.

Hexamethyldisiloxane↗
DOE OSTI2026

Flame dynamics and kinetic coupling of ammonia and dimethyl-ether in non-premixed cool and warm flames at elevated pressure

Developing advanced low-temperature combustion engines with ammonia-biofuel blends requires a comprehensive understanding of low-temperature flame dynamics and kinetic interactions between ammonia and oxygenated fuels at elevated pressures. This work aims to study the dynamics and kinetics of non-premixed Dimethyl Ether (DME)/Ammonia (NH 3 ) cool and warm flames, and their reignition to hot flames. A counterflow burner is employed to establish DME/NH 3 cool/warm flames at pressures up to 5 atm. The extinction limits of cool flame and the reignition limits of warm flame to hot flame are measured by varying NH 3 concentrations and compared to simulations to quantitatively examine the effects on DME/NH 3 flames. It is found that NH 3 inhibits low-temperature DME oxidation and results in lower cool flame extinction limits. Warm flames in the presence of NH 3 are observed for the first time, revealing a non-monotonic effect of NH 3 addition: a small amount of NH 3 presence enhances warm flame chemistry and promotes reignition to hot flames, while a high NH 3 concentration weakens the warm flame. This trend is further explained by 0-D PSR kinetic simulations and 1-D S-curve flame dynamic calculations. Three flame transition regimes between cool flames (CF), warm flames (WF), and hot flames (HF) by different levels of NH 3 additions at a specific strain rate are identified, namely WFHF reignition, WF-CF transition, and WF extinction. Reaction sensitivity analyses of OH at low temperatures show that NH3 inhibits DME oxidation through OH consumption via H-abstraction and the kinetic couplings of RO 2 /NH 2 , RO 2 /NO x , R/NO x , and O 2 QOOH/NO x further suppress the low-temperature branching. At intermediate-temperatures, NH 2 /NO x /HO 2 coupling promotes warm flames via the pathway NH 2 → H 2 NO → HNO → NO by converting O 2 → HO 2 → OH. At even higher NH₃ concentrations, radical termination reactions of NH 2 + NO/NO 2 and excessive OH consumption via H-abstraction inhibit the flame. The insights into the kinetic coupling between NH 3 and low-temperature chemistry at elevated pressure and its impact on the dynamics of cool-warm-hot flame transitions will contribute to advancing combustion technologies with reduced emissions and improved energy-efficiency.

42 ENGINEERING↗
DOE OSTI2026

An experimental and kinetic modeling study of the ignition of 2-methyl decane

Ignition delay times (IDTs) of 2-methyl decane (C 11 H 24 -2) are measured in a high-pressure shock tube and in a rapid compression machine at equivalence ratios in the range 0.5–2.0 at 90 % dilution, at temperatures in the range 600–1430 K and at pressures of 15 and 30 bar. To clarify the effect of the branched methyl group on fuel oxidation, IDTs of n-undecane (nC 11 H 24 ) are also measured at similar conditions to those measured for C 11 H 24 -2. A new chemical kinetic mechanism, using C3MechV4.0.1 as the core chemistry, is developed and validated against the new experimental data. The thermodynamic properties of the fuel (RH), alkyl (Ṙ), alkyl peroxy (RȮ 2 ), hydroperoxy-alkyl (Q̇OOH), and peroxy hydroperoxy alkyl (Ȯ 2 QOOH) radicals are updated in both the C 11 H 24 -2 and nC 11 H24 models using THERM25. A reaction path flux analysis for C 11 H 24 -2 at different temperatures was conducted. Compared to nC 11 H 24 , C 11 H 24 -2 shows slower reactivity. At low and intermediate temperatures, the chain propagation pathway Q̇OOH ↔ C 11 cyclic ether + ȮH is favored for C 11 H 24 -2, while the chain branching pathway Q̇OOH ↔ Ȯ 2 QOOH ↔ C 11 carbonyl hydroperoxide (KHP) + ȮH is suppressed, leading to lower reactivity compared to nC 11 H 24 . At high temperatures, the presence of the branched methyl group inhibits the direct decomposition of the fuel, resulting in reduced C 2 H 4 formation, which in turn suppresses the reactivity of the fuel.

2-methyl decane↗
DOE OSTI2025

Toucan: A performance portable, scalable implementation of the DECA algorithm

In the field of additive manufacturing (AM), cellular automata (CA) is extensively used to simulate microstructural evolution during solidification. However, while traditional CA approaches are relatively fast, they still require a substantial number of time steps, are limited to moderate volumes, and are relatively difficult to improve through parallelism due to the highly localized nature of the solidification front. Here, to address these issues of time to solution and load balancing, we introduce Toucan, a parallel, performance-portable, and scalable code written in C++ with the Kokkos library that leverages the discrete event inspired cellular automata (DECA) algorithm to perform parallel-in-time (PinT) grain growth simulations. Toucan effectively mitigates load balancing issues by distributing the computational workload more evenly across processors, enhancing scalability and efficiency. We conduct both strong and weak scaling studies on up to 64 GPUs on the Frontier supercomputer, demonstrating that Toucan significantly outperforms the current state-of-the-art, time-stepped CA code, ExaCA, on both single and multi-GPU simulations. Even in AM-specific weak scaling scenarios, Toucan maintains near-ideal scaling, in contrast to the linear increase observed with ExaCA due to the moving laser raster pattern. This study highlights Toucan’s potential to transform microstructural simulations in AM by radically improving both efficiency and scalability over existing methods.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Equation-based and data-driven modeling: Open-source software current state and future directions

Here, a review of current trends in scientific computing reveals a broad shift to open-source and higher-level programming languages such as Python and growing career opportunities over the next decade. Open-source modeling tools accelerate innovation in equation-based and data-driven applications. Significant resources have been deployed to develop data-driven tools (PyTorch, TensorFlow, Scikit-learn) from tech companies that rely on machine learning services to meet business needs while keeping the foundational tools open. Open-source equation-based tools such as Pyomo, CasADi, Gekko, and JuMP are also gaining momentum according to user community and development pace metrics. Integration of data-driven and principles-based tools is emerging. New compute hardware, productivity software, and training resources have the potential to radically accelerate progress. However, long-term support mechanisms are still necessary to sustain the momentum and maintenance of critical foundational packages.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Enhanced boron doping for diamond growth in cylindrical CVD reactors using a Faraday cage

We report on the use of 3D-printed titanium Faraday cages in a cylindrical resonant cavity CVD reactor for the growth of heavily boron-doped single-crystalline superconducting diamond. We show that the cage enhances boron doping uniformity and provides control of the doping efficiency and growth rate by varying the distance between the plasma and the diamond surface. Using finite element modelling, we demonstrate that the Faraday cage focuses the electric field at the top of the cage whilst also suppressing the electric field inside the cage, thereby pinning the plasma at a fixed height above the sample. This eliminates plasma “hot spots” at the edges of the sample and creates a more uniform plasma electron and hydrogen radical distribution over the sample, resulting in a more uniform boron incorporation.

Boron doped diamond↗
DOE OSTI2022

A reaction–diffusion model for grayscale digital light processing 3D printing

We report that digital light processing (DLP) 3D printing is an additive manufacturing process that utilizes light patterns to photopolymerize a liquid resin into a solid. Due to the accuracy of modern digital micromirror devices (DMD) and recent advances in resin chemistry, it is now possible to create functionally graded structures using different light intensity values, also known as grayscale DLP (g-DLP). Different intensities of light lead to differences in the polymer crosslinking density after curing, which ultimately produces a part with gradients of material properties. However, g-DLP is a complicated process. First, the DLP printing is a highly coupled chemical and physical process that involves light propagation, chemical reactions, species diffusion, heat transfer, volume shrinkage, and changes in mechanical behaviors of the curing resin. Second, in g-DLP, light gradients create strong in plane gradients of chemical species concentrations in the curing liquid resin due to the strong dependence of light intensity on the rate of monomer crosslinking. Furthermore, light gradients through the depth create concentration gradients due to the degree of cure dependent light absorption and the use of photoabsorbers. These complex physical features of the printing process must be understood in order to properly control printing parameters such as light exposure time, printing speed, and grayscale variations to achieve accurate mechanical properties. In this paper, a photopolymerization reaction–diffusion model is developed and used in conjunction with experiments to investigate the coupled effects of light propagation, chemical reaction rates, and species diffusion during g-DLP 3D printing. The model is implemented numerically utilizing the finite difference method and simulation results are compared to experimental findings of simple printed structures. The agreement between experimental and model predictions of simple quantities of interest, such as geometric feature sizes, shows that the model can capture the overcure due to free-radical and other species diffusion during printing when grayscale patterns are employed. This model lays the groundwork for future extensions that can incorporate more complex coupled physics such as heat transfer, volume shrinkage, and material property evolution, which are critically important in utilizing g-DLP 3D printing for the fabrication of high-performance parts which excellent geometric and material property tolerances.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Controllable depolymerization of lignin using carbocatalyst graphene oxide under mild conditions

Many efforts have been devoted to depolymerizing lignin to high-value aromatic monomers. However, a controllable and complete lignin inter-unit bond cleavage remains difficult because of the randomness in lignin structure and low selectivity in catalysis. Herein, we demonstrate that a sole carbocatalyst graphene oxide (GO) can selectively oxidize aryl- and alkyl-hydroxyl groups of lignin without the presence of any additional oxidizers, leading to lignin degradation under mild conditions. The results from lignin dimers, lignin oligomers, and organosolv lignin indicate that the depolymerization is catalyzed via an oxidative mechanism through an enol ether intermediate, a stable intermediate that has been chemically isolated here. Quantum chemistry calculations at the density functional theory level indicate that the phenolic –OH group in lignin facilitates the formation of radical cations and thus prohibits ring opening reactions, leading to the complete conversion of lignin dimers to guaiacol through the enol ether intermediate.

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

Non-catalytic gas phase NO oxidation in the presence of decane

In this study, experiments without a catalyst revealed that NO was oxidized in a diesel exhaust gas phase mixture due to the presence of n -C 10 H 22 (decane). This reaction was observed to occur following the low temperature oxidation catalyst test protocol (LTC-D) defined by U.S. DRIVE. The purpose of LTC-D conditions is to simulate an aftertreatment diesel combustion gas mixture in order to test candidate catalyst materials. 100% NO conversion was observed, without a catalyst, after beginning to react at 330 °C accompanied by consumption of decane. After experiments that isolated hydrocarbons, ethylene was also observed to facilitate NO oxidation to a lesser degree (>470 °C). Density functional theory (DFT) calculations were conducted to investigate thermodynamically-possible initiating elementary steps during n -C 10 H 22 consumption and NO oxidation. Two feasible intermediate radicals to oxidize NO to NO 2 are ·C 10 H 21 O 2 and ·HO

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

The effects of turbulence on the flame structure and NO formation of ammonia turbulent premixed combustion at various equivalence ratios

Ammonia is carbon-free and is regarded as a potential fuel to address global warming issues. Here, in this work, three-dimensional direct numerical simulations (DNS) of ammonia/air turbulent premixed flames were performed to explore the influence of turbulence and equivalence ratio on the flame structure and NO formation characteristics. Two equivalence ratios were considered, i.e. Φ=0.9 and Φ=1.1. The general flame structures were presented and species distributions were examined. The NO mass fraction was found to be the highest in the product for the lean case and in the reaction zone for the rich case. The conditional mean values of species mass fractions and reaction rates were compared with those of the unstrained and strained laminar premixed flames to explore how well the laminar flame structures can approach those of the turbulent flame. The budget analysis of the species transport equations showed that turbulent diffusion plays an important role in species transport. The turbulent diffusivity D T was estimated using the gradient hypothesis based on the DNS data. Various laminar flame simulations with different diffusivities were carried out. It was shown that the conditional means of the DNS agree well with those of the laminar flames including D T in the transport property calculation. The global and local NO formation characteristics were investigated. It was found that the mean NO production conditioned on the progress variable is lower compared with the corresponding laminar flame in the rich case. However, the relative contributions from various NO pathways are rarely affected by turbulence. The NO mass fraction is higher in negative curvature regions compared with positive curvature regions of the flame surface for the rich case, which is due to the preferential diffusion of H 2 and other radicals and the enhanced NO pathways in negatively curved regions.

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

Logical approach to tritium vacuum pump selection for fusion applications

Vacuum pumps have been used for tritium processing since the 1940s following the isolation of pure tritium. The largest problem facing tritium vacuum pumps is the tritium radioactive decay. The beta particles, recoil energy, and subsequent radical formations are detrimental to most carbon-based materials. Carbon materials, such as plastics, elastomers, carbon fiber, activated carbon, etc., react with tritium to form methane and in the process strip material from the carbonbased material. Other reactions can be self-catalyzing, as is the case with nitrogen and tritium forming tritiated ammonia at room temperature and pressure. Tritium also readily isotopically exchanges with surface -OH groups, increasing the amount of protium in the gas stream while creating a tritium holdup on the surface. Metals are the preferred material when working with tritium, but tritium will permeate into and through almost all metals at some point. The entrapped tritium can decay inside the metal, creating local helium distortions that can impact the structural properties above and beyond the hydrogen embrittlement issue tritium can cause.

Angelette, Lucas M. [Savannah River National Labor↗
DOE OSTI2021

Lignin-enhanced reduction of structural Fe(III) in nontronite: Dual roles of lignin as electron shuttle and donor

Lignin is a major component of plant-derived soil organic matter (SOM) in soils and sediments. Fe-bearing clay minerals are widely distributed in these environments and often co-exist with lignin. While previous studies have reported the electron shuttling and donating roles of certain redox-active SOM in the dissimilatory reduction of structural Fe(III) in Fe-bearing clay minerals, the role of lignin in this process remains unknown. Here we studied this role by incubating an Fe-rich smectite (nontronite NAu-2) with two types of lignin (soluble and insoluble) in the absence and presence of an Fe(III)-reducing bacterium Shewanella putrefaciens CN32 under anaerobic condition. Labile carbon source (i.e., lactate) was added in some experiments as an extra electron donor. The results demonstrated that both soluble and insoluble lignins abiotically reduced structural Fe(III) in NAu-2. The reduction extent was proportional to lignin concentration. After abiotic reaction, lignin served as either electron shuttle or electron donor in the presence of CN32: (1) When lactate was present, lignin served as an electron shuttle to enhance the rate of Fe(III) reduction; (2) When lactate was absent, lignin served as an electron donor for Fe(III) reduction. Although the biotic Fe(III) reduction extents were similar in the presence of either soluble or insoluble lignin, the reduction rates with soluble lignin were always higher than those with insoluble lignin, likely owing to their different electron transfer mechanisms. After interaction with NAu-2 and/or CN32, soluble lignin structure largely remained intact, but with some decreases of humic/fulvic acid-like and protein-like compounds, aromatic functional groups (e.g., C-H, C=O, COOH), and aliphatic/aromatic compounds. An increase of semiquinone-like organic radicals was observed after lignin interaction with NAu-2. These chemical changes of lignin were likely coupled with reduction of structural Fe(III) in nontronite. The nontronite structure remained largely unchanged without much dissolution and mineral transformation. The findings of this study provide insights into the role of lignin in promoting mineral-microbe interactions and have significant implications for coupled Fe and C biogeochemical cycles in soils and sediments.

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

Experimental determinations of carbon and hydrogen isotope fractionations and methane clumped isotope compositions associated with ethane pyrolysis from 550 to 600°C

Methane clumped isotope compositions signify the relative natural abundances of rare, doubly substituted isotopic species of methane ( 13 CH 3 D and 12 CH 2 D 2 ) and have emerged as a new isotopic tool to trace the sources, sinks, and lifecycles of methane in the environment. Such measurements can identify equilibration (or reequilibration) temperatures if found to be in isotopic equilibrium or non-equilibrium processes (e.g., kinetically controlled reactions or mixing) if not in isotopic equilibrium. Naturally occurring thermogenic methane—formed by the thermally activated breakdown of larger organic molecules—has been found to have clumped isotope compositions consistent with equilibrium at reasonable gas formation temperatures in some settings and non-equilibrium processes occurring during either formation, migration, storage, or extraction in others. To explore the potential controls on the isotopic composition of thermogenic methane, we conducted isothermal time-series ethane pyrolysis experiments at 550 and 600 °C to measure methane and ethane 13 C/ 12 C and D/H fractionations and methane clumped isotope compositions (resolved 13 CH 3 D and 12 CH 2 D 2 ). We explore the effects of modifying the initial clumped isotope composition of ethane and the addition of water vapor to pyrolysis experiments. We observe that ethane and methane 13 C/ 12 C are controlled by kinetic isotope effects and Rayleigh distillation processes. In contrast, ethane and methane D/H and methane clumped isotope compositions appear to be controlled by a combination of these processes and hydrogen isotope exchange. The hydrogen isotope exchange processes lead to isotopic equilibrium as reaction completion is approached for both D/H (ethane/methane) and methane clumped isotope compositions. Here, we develop a chemical model based on a mass balance approach that accounts for inheritance vs. hydrogen-abstraction formation pathways for singly and doubly substituted isotopologues of ethane and methane that is compared to the experimental data. The model allows the determination of carbon and hydrogen kinetic isotope effects associated with ethane cracking and hydrogen abstraction reactions that, where applicable, we compare to prior theoretical constraints. From the comparison of the model to the experimental data, we infer that the kinetically controlled ethane and methane bulk isotope compositions and methane clumped isotope compositions are controlled by kinetic isotope effects (both primary and secondary) associated with both C–C bond and C–H bond cleavage reactions. Specifically, the methane clumped isotope compositions likely result from a combination of clumped isotope effects associated with ethane breakdown and/or assembly of methane isotopologues (expressed in terms of γ-factor parameters ≠ 1) and combinatorial effects that arise probabilistically. We discuss our experimental results in the context of recent pyrolysis experiments and observations of naturally occurring thermogenic methane. We consider a proposal consistent with observations from nature that the hydrogen isotope exchange reactions that promote equilibration of methane isotopic molecules at or near formation temperature may be facilitated by free radicals generated by pyrolysis reactions. In this framework, isotope exchange effectively ceases when pyrolysis effectively ceases locking in compositions that can be consistent with peak formation temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗