Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactivity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

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

An Evaluation of Actinide Reactivity with CO 2 , O 2 , and O 2 /He Gases using Inductively Coupled Plasma Tandem Mass Spectrometry: Application to Simultaneous Measurement of 241 Am/ 241 Pu Ratios in Unseparated Complex Matrices

Accurate actinide measurements are critical within the field of nuclear science. Traditional methods for actinide quantification require time-consuming sample processing prior to analysis. There is a need for rapid analytical techniques that still maintain a high degree of accuracy. In this work, actinide reactivity was assessed for multiple oxygen-containing reaction gases using quadrupole inductively coupled plasma tandem mass spectrometry (Q-ICP-MS/MS) to evaluate actinide analysis in complex sample matrices without analyte-matrix separation. A novel method was developed to measure 241 Am/ 241 Pu in complex sample matrices using O 2 /He reaction gas with no matrix removal or analyte pre-concentration. This inline method reduces matrix-derived polyatomic interferences that complicate traditional ICP-MS analyses by mass-shifting to 241 Am 16 O + and 241 Pu 16 O 2 + , allowing Am and Pu to be mass separated for simultaneous analysis. While mass shifting is efficient, a small portion of Am + (<1.3%) and Pu + (<1.4%) react to from AmO 2 + and PuO + , respectively. Therefore, a mass balance approach was used, in combination with reactivity determined from 242 Pu and 243 Am standard solutions, to correct for residual 241 PuO + and 241 AmO 2 + . The method was validated by measuring 241 Am/ 241 Pu in Pu isotope standards CRM-136 and CRM-137 (separated in March/April 1970 and February 2022, respectively) in both neat solutions and complex matrices containing diluted soil (NIST SRM 2711a, >1000 µg·g -1 ). Method detection limits of 15.9 and 9.6 fg·g -1 were determined for 241 Am and 241 Pu, respectively, and 241 Am/ 241 Pu ratios were measured with accuracies within <3.5%. In conclusion, this work presents the first direct analysis of 241 Am/ 241 Pu in unseparated complex matrices, advancing capabilities for rapid actinide measurements.

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

Butanol-Based Organosolv Lignin and Reactive Modification of Poly(ethylene-glycidyl methacrylate)

Biomass processing industries and modern society are both interested in environmentally friendly plant-based polymers, such as lignin, to replace petrochemical derivatives. In this study, an n-butanol-based solvent system that can be produced through fermentation of cellulosic biomass, was used to fractionate hybrid poplar and isolate high purity lignin. The resulting unique glassy lignin fraction exhibited good melt-flow characteristics but poor mechanical behavior. Therefore, to produce a reusable and tough thermoplastic material, the lignin was reactively blended with glycidyl methacrylate modified polyethylene plastic using a high shear melt-mixing technique at 170 °C. Mechanical and rheological analyses revealed that the reactively processed blends with 30 wt % lignin have good compatibility between both polymers. Scanning electron microscopy, dynamic mechanical thermal analysis, and infrared spectroscopy data indicated improved interaction between lignin and the functional polyethylene matrix. The blends exhibited a thermal relaxation in the temperature range of 10 to 20 °C, because of the branched glycidyl methacrylate moiety. Finally, the optimal composition, consisting of 30% lignin, exhibited ~10 MPa tensile strength with 50% strain to failure; such a lignin-modified polymer may find use as packaging, as textile coatings, and in various horticultural applications.

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

Vibrational Quantum-State-Controlled Reactivity in the O 2 + + C 3 H 4 Reaction

Quantum-state-controlled reactivity is a long-standing goal in the field of physical chemistry. In this work, we explore the vibrational-state-dependent behavior of the ion–molecule reaction between O 2 + in distinct vibrational states and two isomers of C 3 H 4 , allene (H 2 C 3 H 2 ) and propyne (H 3 C 3 H). While most products are formed regardless of the vibrational state of O 2 + , the branching ratios are influenced by vibrational excitation, and a new product, C 2 O + , appears exclusively in the excited-state reactions. This selective formation of C 2 O + demonstrates that vibrational excitation can effectively activate a reaction pathway, providing direct evidence of quantum-state control in the reactivity. These results represent an important step toward the goal of quantum-state-controlled chemistry in molecular systems.

Chemical reactions↗
DOE OSTI2020

Controlling P–C/C–H Bond Cleavage in Nickel Bis(diphosphine) Complexes: Reactivity Scope, Mechanism, and Computations

The synthesis of heteroleptic [Ni(P 2 N 2 )(diphosphine)][BF 4 ] 2 complexes and the cleavage of P–C and C–H bonds of the P 2 N 2 ligand in those complexes are reported here. The products are five-coordinate complexes in which Ni–C and P–H bonds have formed to give a cyclic moiety containing Ni–CH$=$NR 2 . The reactivity of [Ni(P 2 N 2 )(diphosphine)][BF 4 ] 2 complexes is influenced by the rigidity of the diphosphine, the steric effect of the substituents, and length of the carbon linker of the diphosphine ligands. Diphosphine ligands bearing a rigid backbone (e.g., dmpbz, 1,2-bis(dimethylphosphino)benzene) or aromatic substituents (e.g., dppe, 1,2-bis(diphenylphosphino)ethane) react with [Ni(P t Bu 2 N Bn 2 )(CH 3 CN) 2 ][BF 4 ] 2 to give P–C/C–H bond cleavage products. Both [Ni(P t Bu 2 N Bn 2 )(dmpe)(MeCN)][BF 4 ] 2 and [Ni(P t Bu 2 N Bn 2 )(dmpm)(MeCN)][BF 4 ] 2 (dmpm = 1,2-bis(dimethylylphosphino)methane) were prepared by the reaction of [Ni(P t Bu 2 N Bn 2 )(CH 3 CN) 2 ][BF 4 ] 2 with the corresponding diphosphine ligands. [Ni(P t Bu 2 N Bn 2 )(dmpe)(MeCN)][BF 4 ] 2 readily undergoes P–C/C–H bond cleavage in nitromethane. In sharp contrast, [Ni(P t Bu 2 N Bn 2 )(dmpm)][BF 4 ] 2 is stabilized by dmpm, a diphosphine with small bite angle, and does not show P–C/C–H bond cleavage reactivity. Computational results show that for complexes bearing less bulky diphosphine ligands, such as dmpm, the barriers for the rate-determining transition states are in some examples higher than 30 kcal/mol with the M06 functional, higher than those for complexes bearing more rigid or more bulky ligands, consistent with experimental studies. The calculated barriers for the first transition state correlated with increased values of the dihedral angle formed by the two NiP 2 planes.

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

Ultrafast Reactive Laser Sintering of Highly Conductive Garnet-Type LLZTO Solid Electrolytes

Rapid and scalable fabrication of garnet-type solid electrolytes remains a major challenge for the practical deployment of lithium metal batteries. Here, we report reactive laser sintering (RLS) as an ultrafast and potentially scalable strategy for fabricating garnet-type Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 (LLZTO) solid electrolytes. RLS of LLZTO enables simultaneous reaction and densification, achieving ∼95% relative density while minimizing lithium loss and suppressing secondary phase formation. Compared to conventional furnace sintering, RLS promotes enhanced grain growth and improved densification, leading to improved ionic conductivity (0.36 ± 0.08 mS cm −1 ) while maintaining comparable activation energies for Li + transport. Structural characterization by X-ray diffraction (XRD), Raman spectroscopy, and solid-state 6 Li/ 7 Li NMR confirms the formation of cubic garnet LLZTO with homogeneous microscale elemental distribution. In addition, nanoindentation measurements demonstrate that RLS preserves the mechanical properties of the garnet framework despite ultrafast localized thermal processing. By integrating simultaneous reaction and densification with tunable microstructural control, reactive laser sintering provides a promising manufacturing pathway for high-performance garnet solid electrolytes toward next-generation solid-state batteries.

CO2 laser↗
DOE OSTI2025

Chemical Recycling of Polycaprolactones via Reactive Melt Processing

Chemical recycling is a promising technology for the deconstruction of waste plastics into monomers or other chemical intermediates, which can be converted into other value-added products or repolymerized into polymers. In this work, polycaprolactone and poly(4-propylcaprolactone) were chemically recycled by exploiting ring-closing depolymerization (RCD). To demonstrate the utility of this approach, reactive distillation experiments were initially conducted in a heated round-bottom flask with a distillation head to explore the effect of temperature and catalyst loading on RCD. Inspired by industrial devolatilization equipment, analogous experiments were conducted using a twin-screw melt compounder that was modified to simultaneously heat and mix polymer in the presence of a catalyst while continuously removing the evolved monomer vapors with a vacuum system and collecting them in a cold trap. Yield and selectivity of the recovered monomer were characterized by 1 H nuclear magnetic resonance spectroscopy and gas chromatography-mass spectrometry. Importantly, the more scalable twin-screw melt compounder approach produced values of monomer selectivity (∼90%) and yield (∼80%) similar to those of the round-bottom flask experiments. The recovered monomers were also successfully repolymerized into the same neat polymer with comparable molar masses without purification steps or adding initiator. Here, the results of this study promote the use of reactive devolatilization extrusion as a scalable approach to polyester depolymerization.

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

CH 4 Activation over Perovskite Catalysts: True Density and Reactivity of Active Sites

The high thermal stability of perovskites has drawn attention toward their applications for catalytic CH 4 activation and conversion, typically occurring at high temperatures. The reaction rates of perovskite catalysts for CH 4 combustion, however, trail behind those of noble metal catalysts. Ways to optimize the performance of perovskite catalysts are destined to trial-and-error approaches unless their complex reconstructed surfaces are correlated with fundamental kinetic parameters. Discerning the intrinsic activity of surface catalytic sites and the density of those sites is crucial to rationally envision complex metal oxides with enhanced catalytic performance. Here, the present work presents a detailed kinetic analysis of catalytic CH 4 combustion over a set of seven perovskites (SrTiO 3 , SrZrO 3 , SrFeO 3 , LaFeO 3 , LaInO 3 , LaCoO 3 , LaMnO 3 ) with various surface terminations. Steady-state isotopic transient kinetic analysis was employed to measure turnover frequency (TOF) and density of surface intermediates (N) under operando conditions. Top surface characterization elucidated performance-structure relationships between near-monolayer surface composition and intrinsic reactivity of the catalysts. By using a chemical etching procedure to expose Fe-sites at the top surface of LaFeO 3 (LaFeO 3 , HNO 3 ), its TOF was increased 4-fold, compared with the unmodified sample, although N on the surface of LaFeO 3 , HNO 3 decreased. Density functional theory simulations corroborated that surface Fe-termination and La-Fe termination offer lower energetic barriers for CH 4 activation when compared with La-termination. In general, surface reconstruction is shown as a tool to tune TOF and N to improve reaction rates. This work fills a gap in current kinetic studies of perovskites through a careful assessment and discussion of the density and intrinsic reactivity of active sites for methane combustion over well-characterized reconstructed perovskite surfaces.

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

A Thioether-Ligated Cupric Superoxide Model with Hydrogen Atom Abstraction Reactivity

The central role of cupric superoxide intermediates proposed in hormone and neurotransmitter biosynthesis by noncoupled binuclear copper monooxygenases like dopamine-β-monooxygenase has drawn significant attention to the unusual methionine ligation of the Cu M (“Cu B ”) active site characteristic of this class of enzymes. The copper–sulfur interaction has proven critical for turnover, raising still-unresolved questions concerning Nature’s selection of an oxidizable Met residue to facilitate C–H oxygenation. We describe herein a model for Cu M , [( TMG N 3 S)Cu I ] + ([1] + ), and its O 2 -bound analog [( TMG N 3 S)Cu II (O 2 •– )] + ([1·O 2 ] + ). The latter is the first reported cupric superoxide with an experimentally proven Cu–S bond which also possesses demonstrated hydrogen atom abstraction (HAA) reactivity. Introduction of O2 to a precooled solution of the cuprous precursor [1]B(C 6 F 5 ) 4 (-135 °C, 2-methyltetrahydrofuran (2-MeTHF)) reversibly forms [1·O 2 ]B(C 6 F 5 ) 4 (UV/vis spectroscopy: λ max 442, 642, 742 nm). Resonance Raman studies (413 nm) using 16 O 2 [ 18 O 2 ] corroborated the identity of [1·O 2 ] + by revealing Cu–O (446 [425] cm –1 ) and O–O (1105 [1042] cm –1 ) stretches, and extended X-ray absorption fine structure (EXAFS) spectroscopy showed a Cu–S interatomic distance of 2.55 Å. HAA reactivity between [1·O 2 ] + and TEMPO–H proceeds rapidly (1.28 × 10 –1 M –1 s –1 , -135 °C, 2-MeTHF) with a primary kinetic isotope effect of kH/kD = 5.4. Comparisons of the O 2 -binding behavior and redox activity of [1] + vs [2] + , the latter a close analog of [1]+ but with all N atom ligation (i.e., N 3 S vs N 4 ), are presented.

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

Enhanced CO 2 Reactive Capture and Conversion Using Aminothiolate Ligand–Metal Interface

Metallic catalyst modification by organic ligands is an emerging catalyst design in enhancing the activity and selectivity of electrocatalytic carbon dioxide (CO 2 ) reactive capture and reduction to value-added fuels. However, a lack of fundamental science on how these ligand–metal interfaces interact with CO 2 and key intermediates under working conditions has resulted in a trial-and-error approach for experimental designs. With the aid of density functional theory calculations, we provided a comprehensive mechanism study of CO 2 reduction to multicarbon products over aminothiolate-coated copper (Cu) catalysts. Our results indicate that the CO 2 reduction performance was closely related to the alkyl chain length, ligand coverage, ligand configuration, and Cu facet. The aminothiolate ligand–Cu interface significantly promoted initial CO 2 activation and lowered the activation barrier of carbon–carbon coupling through the organic (nitrogen (N)) and inorganic (Cu) interfacial active sites. Experimentally, the selectivity and partial current density of the multicarbon products over aminothiolate-coated Cu increased by 1.5-fold and 2-fold, respectively, as compared to the pristine Cu at –1.16 V RHE , consistent with our theoretical findings. Furthermore, this work highlights the promising strategy of designing the ligand–metal interface for CO 2 reactive capture and conversion to multicarbon products.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Unraveling Reactivity Differences: Room-Temperature Ring-Opening Metathesis Polymerization (ROMP) versus Frontal ROMP

Here, in this study, we explore the distinct reactivity patterns between frontal ring-opening metathesis polymerization (FROMP) and room-temperature solventless ring-opening metathesis polymerization (ROMP). Despite their shared mechanism, we find that FROMP is less sensitive to inhibitor concentration than room-temperature ROMP. By increasing the initiator-to-monomer ratio for a fixed inhibitor/initiator quantity, we find reduction in the ROMP background reactivity at room temperature (i.e., increased resin pot life). At elevated temperatures where inhibitor dissociation prevails, accelerated frontal polymerization rates are observed because of the concentrated presence of the initiator. Surprisingly, the strategy of employing higher initiator loading enhances both pot life and front speeds, which leads to FROMP rates exceeding prior reported values by over 5 times. This counterintuitive behavior is attributed to an increase in the proximity of the inhibitor to the initiator within the bulk resin and to whether the temperature favors coordination or dissociation of the inhibitor. A rapid method was developed for assessing resin pot life, and a straightforward model for active initiator behavior was established. Modified resin systems enabled direct ink writing of robust thermoset structures at rates much faster than previously possible.

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

Reactive Transport Modeling of Mineral Precipitation and Carbon Trapping in Discrete Fracture Networks

In this study we use numerical experiments to analyze reactive flow and transport behavior in discrete intersecting fracture networks, focusing on (a) how reaction-induced changes in physical and chemical properties affect flow connectivity and (b) how fracture networks developed in the Earth's critical zone contribute to carbon sequestration via mineral weathering reactions. In the first part of the study, we used two-dimensional reactive flow and transport simulations to analyze the impacts of mixing in a natural discrete fracture network. We concluded that reaction-induced changes can substantially alter the flow connectivity, especially at fracture intersections. The second set of simulations considered the problem of natural weathering of fractured mafic and ultramafic rocks in the partially saturated Earth's critical zone as a function of infiltration rates, fracture permeability, and partially saturated flow parameters. As a model system, we considered an incongruent reaction network with dissolution of forsterite and precipitation of magnesite. The behavior is complex in terms of the rate-controlling processes because of the multicomponent nature of the system as shown by the grid Peclet number: the CO 2 behavior is gas diffusion-controlled in the partially saturated zone, while the rate of water flow via the Damkӧhler number controls Mg 2+ transport through the fracture network. The amounts of carbon that can be trapped are modest, but the naturally fractured domain considered here provides a useful “base case” against which various engineered solutions can be compared.

58 GEOSCIENCES↗
DOE OSTI2026

Identifying mechanistic differences between co-fed CO 2 hydrogenation and reactive CO 2 capture using Ru and Pd dual function materials

Dual function materials (DFMs) enable reactive carbon capture (RCkeC), an intensified approach to carbon dioxide capture and utilization for cost and energy input reductions. Yet, there is a fundamental lack of understanding of mechanisms around CO 2 adsorption and subsequent conversion on these materials, hindering further development. Herein, we investigated several supported alkaline metal oxides for their CO 2 adsorption characteristics to find that Na/Al 2 O 3 had the highest CO 2 adsorption capacity, accompanied by a variety of CO 2 adsorption geometries as identified by in situ DRIFTS and computational modeling. The addition of catalytic metals (Ru, Pd) increased the adsorption capacity of Na/Al 2 O 3 without altering binding modes. In the subsequent reactive desorption step, acetate and formate intermediates were observed. Notably, this mechanistic investigation identified that the formation of acetate species was unique to RCC on a DFM, as these species were not observed in co-fed hydrogenation over the DFM or RCC over a Na-free catalyst.

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

Active learning for robust, high-complexity reactive atomistic simulations

Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. In this paper, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn–Sham density functional theory in terms of structure, dynamics, and speciation.

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

Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li 2 + Na reaction

A first-principles based quantum dynamics study of the Li + LiNa(ν = 0, j = 0) → Li 2 (ν', j') + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li 2 Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Furthermore, of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.

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

A physics-informed deep learning description of Knudsen layer reactivity reduction

A physics-informed neural network (PINN) is used to evaluate the fast ion distribution in the hot spot of an inertial confinement fusion target. The use of tailored input and output layers to the neural network is shown to enable a PINN to learn the parametric solution to the Vlasov–Fokker–Planck equation in the absence of any synthetic or experimental data. As an explicit demonstration of the approach, the specific problem of Knudsen layer fusion yield reduction is treated. Here, the predictions from the Vlasov–Fokker–Planck PINN are used to provide a non-perturbative solution of the fast ion tail in the vicinity of the hot spot, thus allowing the spatial profile of the fusion reactivity to be evaluated for a range of collisionalities and hot spot conditions. Excellent agreement is found between the predictions of the Vlasov–Fokker–Planck PINN and the results from traditional numerical solvers with respect to both the energy and spatial distribution of fast ions and the fusion reactivity profile, demonstrating that the Vlasov–Fokker–Planck PINN provides an accurate and efficient means of determining the impact of Knudsen layer yield reduction across a broad range of plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electric fields imbue enzyme reactivity by aligning active site fragment orbitals

It is broadly recognized that intramolecular electric fields, produced by the protein scaffold and acting on the active site, facilitate enzymatic catalysis. This field effect can be described by several theoretical models, each of which is intuitive to varying degrees. In this contribution, we show that a fundamental effect of electric fields is to generate electrostatic potentials that facilitate the energetic alignment of reactant frontier orbitals. We apply this model to demystify the impact of electric fields on high-valent iron–oxo heme proteins: catalases, peroxidases, and peroxygenases/monooxygenases. Specifically, we show that this model easily accounts for the observed field-induced changes to the spin distribution within peroxidase active sites and explains the transition between epoxidation and hydroxylation pathways seen in Cytochrome P450 active site models. Thus, for the intuitive interpretation of the chemical effect of the field, the strategy involves analyzing the response of the orbitals of active site fragments, and their energetic alignment. We note that the energy difference between fragment orbitals involved in charge redistribution acts as a measure for the chemical hardness/softness of the reactive complex. This measure, and its sensitivity to electric fields, offers a single parameter model from which to quantitatively assess the effects of electric fields on reactivity and selectivity. Thus, the model provides an additional perspective to describe electrostatic preorganization and offers ways for its manipulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Thermal-chemical instability of weakly ionized plasma in a reactive flow

The control of plasma instability of weakly ionized plasma in a reactive flow is of great importance in plasma-assisted combustion, catalysis, fuel reforming and material synthesis. In this work, we propose a new concept of plasma chemical instability and analyze the mechanisms and impact of plasma chemical instability on the transition from a uniform discharge to a contracted state of a self-sustained glow discharge in a reactive H 2 –O 2 –N 2 mixture. Here, a one-dimensional numerical model for plasma chemical instability was developed which accounted for convective heat loss, Joule heating of plasma, and major non-equilibrium plasma-assisted combustion kinetic pathways including electron-impact ionization, vibrational energy transfer, electron attachment, combustion reactions, and heat release. The results showed that plasma chemical instability significantly modified the onset plasma current of plasma thermal instability. Specifically, the critical condition of the instability transition was strongly influenced by electron-impact reactant ionization, electron attachment to oxygen, endothermic/exothermic chemical reactions, and the formation of reaction products.

reactive flow↗
DOE OSTI2021

Global Sensitivity Analysis of a Reactive Transport Model for Mineral Scale Formation During Hydraulic Fracturing

Injection of water-based hydraulic fracturing fluid (HFF) into tight shale gas/oil formations can increase formation permeability and enhance production rates, but this process frequently causes mineral scale formation that can occlude pore space and hinder flow. To identify the most important factors that control the formation of mineral scales, we applied a novel global sensitivity analysis method—distance-based generalized sensitivity analysis (DGSA)—to a reactive transport model (RTM) that was previously built and calibrated to simulate precipitation of barite [BaSO4] and iron (hydr)oxide [Fe(OH) 3 ] in shale matrices and on fracture surfaces. Reactive transport simulations were run with model parameters randomly sampled based on assigned uncertainties. Modeling results for barite and Fe(OH)3 formation were clustered using machine-learning algorithms. A list of ranked critical input parameters was obtained after statistical quantification of cumulative distribution functions of input parameters. We found that barite formation is most sensitive to the rate of sulfate ion generation, which is determined by the pyrite dissolution rate coefficient and oxidant availability. In addition, barite formation is sensitive to the initial amounts of barite in HFF and shale, followed by barite thermodynamics/kinetics. For Fe(OH) 3 formation, the ranked factors are Fe(OH)3 precipitation rate coefficients, initial HFF pH, initial Fe(OH) 3 amount in HFF, and oxidant availability. Overall, our results provide insights into managing mineral scale formation during hydraulic fracturing to enhance production. Meanwhile, this study serves as an example of global sensitivity analysis of RTMs using the efficient, straightforward, and open-source DGSA method.

58 GEOSCIENCES↗