Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical kinetics”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Performance of U 3 Si 2 in an LWR following a cladding breach during normal operation

High-density fuels have been proposed as a possible replacement for uranium-dioxide as a fuel for light water reactors (LWRs) due to their increased loading of fissionable material. The objectives of this proposal are (1) to increase reactor cycle length and reactor power, and (2) to offset any neutronic penalty associated with advanced cladding systems. Of the high-density fuels under consideration, there is particular interest in triuranium disilicide (U 3 Si 2 ) due to its increased metal density and favorable thermal properties as compared to UO 2 . However, there are concerns regarding the chemical compatibility of U 3 Si 2 with water and steam as used for LWR coolant.This paper summarizes research on fuel-coolant chemical compatibility for UO 2 LWR fuel during a cladding breach and highlights that generally, because of its chemical inertness, UO 2 -coolant reactions are of little consequence to reactor operation. However, the volumetric expansion associated with the reaction of UO 2 and oxygen is a concern for possible conditions encountered during air ingress of dry storage. These same concerns arise for U 3 Si 2 , which exhibits greater volumetric expansion than UO 2 when exposed to water or steam. These reactions ultimately result in increased fuel volume that the cladding must accommodate, as well as additional heat generated as the fuel reacts.The BISON fuel performance code was used to perform a comparative analysis on the behavior of UO 2 and U 3 Si 2 under normal operation. Silicide fuel simulations were then extended to demonstrate how varying thermodynamic and chemical kinetics influence fuel expansion and subsequent cladding performance during a cladding breach. These simulations were further extended to a 3D subsection of a fuel rod to demonstrate the characteristics of the resulting cladding crack.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Effects of isoalcohol blending with gasoline on autoignition behavior in a rapid compression machine: Isopropanol and isobutanol

Alcohols, and particularly isoalcohols, are potentially advantageous blendstocks towards achieving efficient, low-carbon intensity internal combustion engines. Their use in advanced configurations, such as boosted spark-ignition or spark-assisted compression ignition, requires a comprehensive understanding of their blending effects on the low- and intermediate-temperature autoignition behavior of petroleum-derived gasoline. This work reports an experimental and modeling study of such autoignition characteristics quantified in a twin-piston rapid compression machine. Isopropanol and isobutanol are blended into a research-grade gasoline (FACE-F) at oxygenate blend levels of 0 to 30% vol/vol, with tests conducted at pressures of 20 and 40 bar, temperatures from 700 to 1000 K, and dilute stoichiometric fuel loadings. Changes to overall reactivity, including first-stage and main ignition times, and preliminary exothermicity are established, with comparisons made to previous measurements with ethanol-blended FACE-F gasoline. Furthermore, it is found that at low-temperature/NTC conditions (700–860 K) the isoalcohols suppress first-stage reactivity and associated heat release while main ignition times are extended. At NTC/intermediate-temperature (860–1000 K) conditions changes to fuel reactivity are less significant with isopropanol slightly suppressing reactivity and isobutanol promoting ignition. Detailed chemical kinetic modeling is used to interpret the experimental measurements. Overall trends of suppression or promotion in the blending behavior are reasonably captured by the model. Sensitivity and rate of production analyses indicate that at lower temperatures H-atom abstraction reactions from the surrogate fuel molecules (e.g., cyclopentane, isooctane) and the isoalcohols via OH are important leading to TC 3 H 6 OH and IC 4 H 8 OH–C radicals, for isopropanol and isobutanol respectively, which act as scavengers in the system. At higher temperatures, similar chemistries are dominant, but there is an increasing importance of abstraction by HO 2 . The kinetic modeling also indicates that the promoting effect of isobutanol at higher temperatures is due to the increased abstractions at the γ-sites, while at lower temperatures abstraction at the α-site leads to greater reactivity suppression.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimental and modeling study of water time histories during H 2 S-N 2 O combustion in a shock tube

Here, with the interest in directly burning sour gas in gas turbines, and the fact that even small amounts of H 2 S or its combustion products can alter combustion characteristics, many research studies have been performed to better understand the combustion chemistry of H 2 S. In the present study, the water formation was followed by laser absorption with N 2 O as an oxidant, instead of O 2 . Nitrous Oxide being essentially consumed via N 2 O (+M) ⇌ N 2 + O (+M), the water formation via the H 2 S + O route can then be probed to further validate the models. Three H 2 S/N 2 O mixtures diluted in 98% Ar were studied to cover the following range of equivalence ratios: 0.5, 1.0, and 2.0, over a wide range of temperatures (1580–1940 K) around atmospheric pressure. A chemical kinetic model was then developed first by validating the base N 2 O kinetics mechanism, then by investigating the experimental results presented herein. The N 2 O kinetics mechanism was updated with recent work on the low- and high-pressure limits for N 2 O decomposition (N 2 O (+M) ⇌ N 2 + O (+M)) as well as a review of rate coefficients for N 2 O + H ⇌ N 2 + OH from the literature. Good agreement is shown for H 2 /N 2 O mixtures. Updates were then made to the H 2 S kinetics mechanism, specifically, an update from the literature on SO 2 + H (+M) ⇌ SO + OH (+M) and an adjustment to SO + SH ⇌ S 2 + OH. Additionally, reactions between SH and N 2 O were determined using W1BD data and transition state theory which required the addition of an NNS sub-mechanism. With these updates, the mechanism provides good agreement with the H 2 S/N 2 O experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Probing O 2 dependence of hydroperoxy-butyl reactions via isomer-resolved speciation

Degenerate chain-branching mechanisms of n-alkanes are centered on the formation of hydroperoxy-alkyl radicals ($\dot{Q}$OOH), formed via $\dot{R}$ + O 2 reactions, and the ensuing competition between unimolecular decomposition and second-O 2 -addition. Quantitative measurements of partially oxidized intermediates formed via reactions of $\dot{Q}$OOH provide critical constraints that are required for accurate modeling of combustion chemistry. To examine the influence of temperature and oxygen concentration on intermediates from unimolecular decomposition of $\dot{Q}$OOH, isomer-resolved speciation measurements were conducted on n-butane oxidation at 835 Torr in a jet-stirred reactor (JSR) from 500 – 900 K. Resulting from negative-temperature coefficient behavior, cyclic ether formation peaked at two temperatures, 650 K and 800 K, which were selected for separate experiments to quantify the O 2 -dependence of species profiles using O 2 concentrations of 4.2 · 10 17 – 1.1 · 10 19 molecules cm –3 . Utilizing vacuum-ultraviolet absorption spectroscopy and electron-impact mass spectrometry, cyclic ether isomers were quantified separately, including explicit resolution of cis– and trans– isomers of 2,3-dimethyloxirane. Stereoisomers of 2-butene were also quantified explicitly. For all cyclic ethers, a common trend in O 2 -dependence emerged: species concentrations reach a maximum near 3.0 · 10 18 molecules cm –3 (equivalence ratio of 0.5). Although quantitative disparities are evident, chemical kinetics modeling qualitatively reproduces the O 2 dependence of species at 650 K. However, at 800 K, weak dependence on O 2 is predicted, which is in contrast with the measurements. Two carbonyls, diacetyl and methyl vinyl ketone, were also quantified and follow similar dependence on [O 2 ] and temperature as the cyclic ethers, which indicates some fraction forms via $\dot{Q}$OOH-mediated reactions. The discrepancies between the measured and model-predicted species profiles indicate that sub-mechanisms for important intermediates may require additional elementary reactions, including stereochemical-specific reactions, to improve the fidelity of n-alkane combustion modeling.

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

Thermal diffusion, exhaust gas recirculation and blending effects on lean premixed hydrogen flames

Thermodiffusively-unstable lean premixed hydrogen flames are investigated using two-dimensional direct numerical simulation employing finite-rate chemical kinetics. Three databases are generated focussing on the inclusion of the Soret effect, the recirculation of exhaust gas, and blending with methane. A simple rescaling of a classic thermal diffusion model is presented and shown to mimic multicomponent diffusion with very low computational cost and little-to-no loss in accuracy. It is also shown that a previously developed model for mean local flame speeds in lean premixed hydrogen flames can still be used provided Soret effects are taken into account in one-dimensional calculations. The addition of exhaust gas to the unburned mixture is found to enhance thermodiffusive instability; the primary mechanism for this was shown to be the highly-efficient third-body nature of water, with the reduction of adiabatic flame temperature a second-order effect. Again, the existing mean local flame speed model proved sufficient. Finally, blending with methane was found to reduce the thermodiffusive response of the flame, more so than the existing model suggests, despite adjustment of the fuel Lewis number; an adapted model is presented to account for this.

08 HYDROGEN↗
DOE OSTI2026

Pulse radiolysis and transient absorption spectra of aqueous solutions of sodium sulfamate

Chemical kinetics for the reactions of sulfamate ions (NH 2 SO 3 − ) with the primary products of water and nitric acid radiolysis were measured in aqueous solutions at ambient temperature. Using time-resolved electron pulse radiolysis techniques with a custom multichannel detection system, we examined the reactivity of NH 2 SO 3 − with the hydroxyl radical ( • OH), hydrogen atom (H • ), and nitrate radical (NO 3 • ). The sulfamate ion was found to react with • OH and H • with second-order rate coefficients of k • OH = (5.60 ± 0.04) × 10 6 M −1 s −1 and k H • = (7.96 ± 0.10) × 10 6 M −1 s −1 , respectively, and with NO 3 • with a rate coefficient of k NO 3 • = (1.67 ± 0.06) × 10 7 M −1 s −1 . The reactions of NH 2 SO 3 − with • OH and H • resulted in the formation of two transient radical species, one with maximum absorbance at 300 nm and a second with maxima at both 300 nm and 600 nm. These spectra are tentatively assigned to • NH 2 SO 3 and • NHSO 3 − , respectively. By measuring the absorbance of these radicals as a function of pH, the radical pK a was determined to be 9.5 ± 0.1. Overall, this work has implications for the longevity and performance of ferrous sulfamate, Fe(NH 2 SO 3 ) 2 , as a plutonium reductant in the reprocessing of used nuclear fuel.

Conrad, Jacy K. [Idaho National Laboratory (INL), ↗
DOE OSTI2021

Investigation of post-injection strategies for diesel engine Catalyst Heating Operation using a vapor-liquid-equilibrium-based spray model

Most multidimensional engine simulations spend much time solving for non-equilibrium spray dynamics (atomization, collision, vaporization). However, their accuracy is limited by significant grid dependency, and the need for extensive calibration. This is critical for modeling cold-start diesel fuel post injections, which occur at low temperatures and pressures, far from typical model validation ranges. At the same time, resolving micron-scale spray phenomena would render full Eulerian multiphase calculations prohibitive. In this study, an improved phase equilibrium based approach was implemented and assessed for simulating diesel catalyst heating operation strategies. A phase equilibrium solver based on the model by Yue and Reitz [1] was implemented: a fully multiphase CFD solver is employed with an engineering-size engine grid, and fuel injection is modeled using the standard Lagrangian parcels approach. Mass and energy from the liquid parcels are released to the Eulerian multiphase mixture according to an equilibrium-based liquid jet model. An improved phase equilibrium solver was developed to handle large real-gas mixtures such as those from accurate chemical kinetics mechanisms. The liquid-jet model was improved such that momentum transfer to the Eulerian solver better reproduces the physical spray jet structure. Validation of liquid/vapor penetration predictions showed that the model yields accurate results with very limited tuning and low sensitivity to the few calibration constants. In-cylinder simulations of diesel catalyst heating operation strategies showed that capturing spray structure is paramount when short, transient injection pulses and low temperatures are present. Furthermore, the EP model provides improved predictions of post-injection spray structure and ignitability, while conventional spray modeling does not capture the increase of liquid penetration during the expansion stroke. Finally, the only important EP model calibration constant, Cliq, does not affect momentum transfer, but it changes the local charge cooling distribution through the local energy transfer, which makes it candidate to additional research. The results confirm that non-equilibrium spray processes do not need to be resolved in engineering simulations of high-pressure diesel sprays.

42 ENGINEERING↗
DOE OSTI2020

Combining Mass Spectrometry of Picoliter Samples with a Multicompartment Electrodynamic Trap for Probing the Chemistry of Droplet Arrays

Single droplet levitation provides contactless access to the microphysical and chemical properties of micrometer-sized samples. Most applications of droplet levitation to chemical and biological systems use nondestructive optical techniques to probe droplet properties. To provide improved chemical specificity, we coupled a multicompartment quadrupole electrodynamic trap (QET) with single droplet mass spectrometry. Our QET continuously traps a monodisperse droplet population (tens to hundreds of droplets) and allows for the simultaneous sizing of a single droplet using its Mie scattering pattern. Single droplets are subsequently ejected into the ionization region of an ambient pressure inlet mass spectrometer. We optimized two complementary soft ionization techniques for picoliter aqueous droplets: (1) paper spray (PS) ionization and (2) thermal desorption glow discharge (TDGD) ionization. Both techniques detect oxygenated organic acids in single droplets, with signal-to-noise ratios >100 and detection limits on the order of 10 pg. Sensitivity and reproducibility across single droplets are driven by the droplet deposition location and spray stability in PS-MS and the ionization region humidity and analyte evaporation rate in TDGD-MS. Importantly, the analyte evaporation rate can control the TDGD-MS quantitative capability because high evaporation rates result in significant ion suppression. Additionally, this effect is mitigated by optimizing the vaporization temperature, droplet size range, and analyte volatility. We demonstrate quantitative and reproducible measurements with a droplet internal standard (<10% RSD) and compare the sensitivity of PS-MS and TDGD-MS. Finally, we demonstrate the application of QET-MS to the study of heterogeneous chemical kinetics with the reaction of gas phase O 3 and aqueous maleic acid droplets.

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

Octane Requirements of Lean Mixed-Mode Combustion in a Direct-Injection Spark-Ignition Engine

Here, this study investigates the octane requirements of a hybrid flame propagation and controlled autoignition mode referred to as mixed-mode combustion (MMC), which allows for strong control over combustion parameters via a spark-initiated deflagration phase. Due to the throughput limitations associated with both experiments and 3-D computational fluid dynamics calculations, a hybrid 0-D and 1-D modeling methodology was developed, supported by experimental validation data. This modeling approach relied on 1-D, two-zone engine simulations to predict bulk in-cylinder thermodynamic conditions over a range of engine speeds, compression ratios, intake pressures, trapped residual levels, fueling rates, and spark timings. Those predictions were then transferred to a 0-D chemical kinetic model, which was used to evaluate the autoignition behavior of fuels when subjected to temperature–pressure trajectories of interest. Finally, the predicted autoignition phasings were screened relative to the progress of the modeled deflagration-based combustion in order to determine if an operating condition was feasible or infeasible due to knock or stability limits. The combined modeling and experimental results reveal that MMC has an octane requirement similar to modern stoichiometric spark-ignition engines in that fuels with high research octane number (RON) and high octane sensitivity (S) enable higher loads. Experimental trends with varying RON and S were well predicted by the model for 1000 and 1400 rpm, confirming its utility in identifying the compatibility of a fuel’s autoignition behavior with an engine configuration and operating strategy. However, the model was not effective in predicting (nor designed to predict) operability limits due to cycle-to-cycle variations, which experimentally inhibited operation of some fuels at 2000 rpm. Putting the operable limits and efficiency from MMC in the context of a state-of-the-art engine, the MMC showed superior efficiencies over the range investigated, demonstrating the potential to further improve fuel economy.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Rate constants and third-body collision efficiencies for recombination of Na with OH and O2: Implications for flame inhibition by alkali salts

Flame inhibition by alkali metals has implications for solid fuel combustion, fire safety, and a number of industrial processes. While the mechanism of inhibition is fairly well understood, details related to thermodynamic properties and rate constants are still in question. In the present work, the recombination of Na with OH (R5) and O2 (R7), respectively, was characterized theoretically, and the implications for modeling laminar premixed hydrogen flames doped with sodium species were examined. Third-body collision efficiencies and low-pressure-limit rate constants were obtained using newly fitted ab initio-based potential energy surfaces, classical trajectories, and one-dimensional master equation calculations. The results are consistent with available experimental results and aid in the present modeling study by providing rate information for bath gases and conditions that remain unexplored experimentally. Chemical kinetic modeling of relative Na and absolute H and OH profiles in H2-fueled laminar, premixed flames doped with a sodium salt shows that the most important radical removal cycle is the sequence Na + OH (+M) → NaOH (+M) (R5), NaOH + H → Na + H2O (R12), even under oxidizing conditions. A secondary cycle, Na + O2 (+M) ⇄ NaO2 (+M) (R7), NaO2 + OH → NaOH + O2 (R14), is less important due to the low thermal stability of NaO2. In most flames, reaction R7 is partially equilibrated, and reaction R14 becomes rate-limiting for the second cycle. The flame analysis supports a lower value of k14 than indicated by recent work on KO2 + OH, but more work is required to confirm this.

Jasper, Ahren W.↗
DOE OSTI2023

Radiolytic Evaluation of 3,4,3-LI(1,2-HOPO) in Aqueous Solutions

We report the octadentate hydroxypyridinone ligand 3,4,3-LI(1,2-HOPO) (abbreviated as HOPO) has been identified as a promising candidate for both chelation and f-element separation technologies, two applications that require optimal performance in radiation environments. However, the radiation robustness of HOPO is currently unknown. Here, we employ a combination of time-resolved (electron pulse) and steady-state (alpha self-radiolysis) irradiation techniques to elucidate the basic chemistry of HOPO and its f-element complexes in aqueous radiation environments. Chemical kinetics were measured for the reaction of HOPO and its Nd(III) ion complex ([Nd III (HOPO)] - ) with key aqueous radiation-induced radical transients (eaq - , H · atom, and · OH and NO 3 · radicals). The reaction of HOPO with eaq - is believed to proceed via reduction of the hydroxypyridinone moiety, while transient adduct spectra indicate that reactions with the H · atom and · OH and NO 3 · radicals proceeded by addition to HOPO's hydroxypyridinone rings, potentially allowing for the generation of an extensive suite of addition products. Complementary steady-state 241 Am(III)-HOPO complex ([ 241 Am III (HOPO)] - ) irradiations showed the gradual release of 241 Am(III) ions with increasing alpha dose up to 100 kGy, although complete ligand destruction was not observed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Consequences of 131 I Transmutation in Gas Phase Radioiodine Molecules and Adsorbed on Graphite Surface

In the case of nuclear power plant accidents, several radionuclides are being dispersed into the environment. In particular, the release of radioiodine 131 I represents a major health concern because of its absorption by the thyroid gland. To reduce the release of 131 I into the environment, capture materials, such as activated carbon, provide a mean to efficiently retain radioiodine-bound molecular species. Here, we use ab initio molecular dynamic calculations to investigate the effects of 131 I radioactive decay on the molecular stability of I 2 , acid iodine (HI, HIO, HIO 2 , and HIO 3 ), and methyl-iodine species in both gas phase and adsorbed on graphite. In gas phase, we found that I 2 is the least resistant to 131 I transmutation, followed by the HIO x family of molecules, then HI and CH 3 I which are the most resistant. However, the graphite surface is found to greatly favorize the fragmentation of adsorbed molecular species. While the identification of gas phase fragments can offer new chemical route to consider in chemical kinetics models, the formation reactive fragment products at the graphite surface can favorize their reaction with surface carbon atoms which overtime could reduce the number of adsorption sites of capture materials available, and hence, impacting their capture capability.

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

Point Defects Control Guest Molecule Diffusion in the 1D Pores of Zn(tbip)

Molecular diffusion plays a critical role in metal-organic frameworks (MOFs) within the application of kinetic chemical separations. We carefully study in this work the unexpected role of point defects for short-chain alkanes diffusing in Zn(tbip), an MOF with rigid one-dimensional (1D) channels. Inside a defect-free Zn(tbip), guest molecules are expected to follow single-file diffusion along 1D channels. It has been found previously by Heinke et al. that these parallel 1D channels are connected for molecular diffusion. Our density functional theory (DFT) calculations suggest that linker vacancy defects could arise under experimentally relevant conditions by removing a pair of linkers. Further climbing-image nudged elastic band (cNEB) DFT calculations indicate that hopping of short-chain alkanes between adjacent 1D channels over defect windows can occur at moderate temperatures. In addition to providing insights into connected adjacent 1D channels in Zn(tbip), Heinke et al. also inferred that most 1D pores are blocked from a microkinetic model. Additionally, we explored the influence of hydrolyzed linker created by the formation of linker removal inside 1D pores. Our DFT calculations show that the linkers can effectively block the pores and the linker diffusion in 1D channels is slow. Our results, for the first time, offer a mechanistic explanation of the unexpected molecular diffusion behavior in this MOF with 1D channels.

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

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗
DOE OSTI2024

Source of Processable Vitrimer Viscosities: Swap Frequencies and Steric Factors

Vitrimers exhibit high, processable viscosities, where other polymers do not, and are among the most promising polymers for closed-loop material circularity. We sought to investigate the underlying chemical kinetic factors that result in high viscosities for vitrimers, which are crucial to designing vitrimers with tunable viscosity. To interrogate these factors, we achieved the first simulated predictions of real vitrimer viscosities, using a novel kinetic Monte Carlo molecular dynamics method, overcoming the time and length scale gaps to predict experimental bulk viscosities. The vitrimer architecture investigated is based on poly(dimethylsiloxane) chains and vinylogous urethane bond swaps. We probed the effects of the extent of free swapping groups, %F, the activation energy, E A , and the steric factor, ρ. The steric factor is related to the intrinsic reaction probability for molecules with sufficient energy. All three factors were found to be significant, but the role of ρ was found to be the biggest and also the most underappreciated. The results show that the inclusion of accurate ρ is of critical importance for viscosity predictions, with the evidence suggesting that the typical assumption of ρ = 1 is not valid for vitrimers and that, indeed, very low steric factors are present in bond-swap vitrimers such that values of ρ < 10 –10 may be typical. This greatly influences the bond exchange rates and, ultimately, the viscosities. Recognition of this result is necessary for the prediction of vitrimer viscosities from molecular simulations and to make vitrimers by design from molecular dynamics. We also investigated the effects that E A , ρ, and the number of free swapping groups have upon vitreous range temperatures, TV, with respect to achieving a specific viscosity (η V = 1 × 10 8 Pa·s), as well as for a commonly reported higher viscosity extrapolation (η V = 1 × 10 12 Pa·s). The evidence suggests that vitrimers may follow universal curves for E A vs T V , as a function of ρ. Finally, this study achieves the first of these comparisons of molecular simulations to experiments and reveals critical insights toward creating vitrimers by design, while providing a route for the prediction of T V from kinetic Monte Carlo molecular dynamics simulations.

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

Emergent Behavior at the Calcite–Water Interface during Reactive Transport in a Simple Microfluidic Channel

Geochemical reactive transport processes in natural mineral-fluid systems may produce a wide array of emergent phenomena that are difficult to predict from basic principles and to reproduce in model systems. In this study, we present experimental results obtained from a simple microfluidic system with which we explored the consequences of reacting the calcite (104) cleavage surface with an acidic Pb-bearing solution (pH = 3.5, [Pb] total = 5 mM) as a function of flow rate. This system is relevant to passive remediation systems for Pb-rich acid mine drainage. We observed periodic banding in the amounts of Pb sorption at flow velocities ≥ 926 μm s -1 , where the band spacing was spatially correlated with the amount of calcite dissolution and the development of micropyramidal topography on the calcite (104) surface. The equivalent coverage of Pb deposited in these Pb-rich bands was at least several monolayers per unit cell, yet there was no evidence for precipitation of any secondary Pb phase implying incorporation of Pb within the near-surface calcite lattice. We also observed spatial variations in nucleation and growth of euhedral secondary Pb-carbonate minerals hydrocerusite and cerussite at flow rates ≤ 278 μm s -1 ,. These findings demonstrate potential for exploiting the rich phenomenology afforded by the interplay among transport phenomena and chemical kinetics in experimental systems designed to yield deeper insights into geochemical self-organization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗