Search NASA⌕ Search

SEARCH · Search NASA

Results for “REACTOR 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 271 records · Page 15

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗

Reduced-order condensed-phase kinetic models for polyethylene, polypropylene and polystyrene thermochemical recycling

Thermochemical recycling of plastic waste (PW) into chemicals and energy vectors requires coupling particle and reactor-scale simulations to accurate condensed phase pyrolysis mechanisms for each constituent. This work proposes a methodology to derive reduced-order condensed-phase kinetic models from validated semi-detailed kinetic mechanisms. Two types of kinetic models are obtained for polyethylene (PE), polypropylene (PP) and polystyrene (PS): reduced semi-detailed models and multi-step fully lumped ones. These families offer different compromises between accuracy and computational cost. The former employ 50–100 gas + liquid species and describe both the radical degradation and the detailed carbon distribution of the products. Conversely, the latter involves 5–10 species per polymer tracking only the main petroleum cuts. The kinetic mechanisms are complemented by the definition of thermochemical properties of gas, liquid, and solid-phase species, accounting for phase-transitions through pseudo-chemical reactions. Model validations are performed by comparison with experimental data and the original semi-detailed mechanisms in terms of mass loss, heat fluxes and product distribution profiles. The resulting CHEMKIN-like condensed-phase models are attached as Supplementary Material and as a GitHub repository. Extending the proposed approach to other polymers and coupling it with existing subsets in the CRECK kinetic framework (e.g., biomass, PVC, PET) offers a powerful tool to model thermochemical recycling of PW and biomass/PW mixtures.

kinetics↗

Experimental and Kinetic Modeling Study of 3-Methyl-2-butenol (Prenol) Oxidation

Longer chain alcohols with 4–5 carbon atoms are attractive alternative fuels as they can be derived from biological sources and since their combustion leads to lower exhaust gas levels of NO x and soot compared to commercial fossil fuels. The auto-ignition behavior of fuels that contain both a hydroxyl group and a C=C double bond in their molecular structure is not well established in the literature. Understanding the influence of these functional groups on the ignition behavior of fuels is critical in the development of tailor-made fuels for advanced combustion engines. In this study, ignition delay times of an unsaturated alcohol, 3-methyl-2-butenol (prenol), are measured using a high-pressure shock tube and a rapid compression machine at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 in “air” in the temperature range 600–1400 K. Furthermore, a detailed kinetic model is developed and validated using the new experimental data in this study. In addition, speciation data in a jet-stirred reactor, ignition delay times, and laminar burning velocities available in the literature were also used to validate the new kinetic model. Fuel flux and sensitivity analyses are performed using this new model to determine the important fuel consumption pathways and critical reactions that affect the reactivity of prenol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling the Pyrolytic Behavior and Kinetics of Single Polymers and Plastic-Rich Municipal Solid Waste Using Thermal Analysis

Pyrolysis is a highly promising thermochemical recycling technology for converting heterogenous plastic waste into sustainable fuels in a single step. Therefore, understanding the pyrolysis mechanism is essential for enabling rational reactor design and enhancing efficient recycling techniques. In this study, the thermal degradation behaviors and corresponding kinetics of pure polymers (PE, PP, and PET) and plastic-rich MSW were examined using simultaneous thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC). Experiments were carried out in the temperature range of 30-800°C with variable heating rates from 5°C/min to 20°C/min in an ultra-high purity Argon atmosphere. Our results indicated that the plastic pyrolysis was an endothermic process, with varying decomposition temperature ranges depending on their structure and composition. Various iso-conversional model-free methods (Friedman, Flynn-Wall-Ozawa, Starink, and Kissinger-Akahira-Sunose) were utilized to determine the apparent activation energy of the plastic degradation, which increased in the following order: PET (214 kJ/mol), PP (218 kJ/mol), PE (245 kJ/mol), and MSW (249 kJ/mol). Finally, Criado’s master plots were employed to identify the best-fitting reaction model and the pre-exponential factor was subsequently determined.

Bashir, Muhammad Aamir↗

An experimental and kinetic modeling study of cyclopentane and dimethyl ether blends

Cyclopentane is a suitable naphthene, or cycloalkane, in a palette for multi-component gasoline surrogate fuels due to its presence in market fuels and its relevance to alkyl substituted cyclopentanes also present. However, the previous oxidation studies of cyclopentane have primarily focused on neat mixtures. Blending cyclopentane with dimethyl ether in this work therefore serves to inform our understanding of, and improve predictive models for, multi-component mixtures. Here, the auto-ignition of cyclopentane/dimethyl ether blends was studied in a high-pressure shock tube and in a rapid compression machine. A wide range of temperatures (650 – 1350 K) and elevated pressures of 20 and 40 bar were studied at equivalence ratios of 0.5, 1.0 and 2.0 in air for two blending ratios (30/70 and 70/30 mole% cyclopentane/di-methyl ether mixtures). A detailed kinetic model for cyclopentane was revised to capture the measured ignition delay times and apparent heat release rates in this study. Literature ignition delay time, jet-stirred reactor, and laminar burning velocity measurements of neat cyclopentane were used as additional validation. Improvements to the kinetic model were based on recent literature studies related to sub-models including cyclopentene and cyclopentadiene which allowed the removal of previous local rate-constant optimizations. Low temperature reactivity of cyclopentane was found to be controlled by the branching ratio between concerted elimination of HO͘ 2 and the strained formation of Q&#775OOH radicals in agreement with previous studies. In this study, the low branching ratio of Q&#775OOH formation increases the influence of a competing consumption pathway for cyclopentyl-peroxy (CPTO͘ 2 J) radicals. The sensitivity of the simulated ignition delay times to the formation of cyclopentyl hydroperoxide (CPTO 2 H), from CPTO͘ 2 J and HO͘ 2 , is discussed. The current model is used to analyze the influence of dimethyl ether on the reactivity of cyclopentane in the context of previous literature studies of dimethyl ether binary blends with ethanol and toluene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-Phase Turbulent Kinetic Energy Budget Computation in Co-Current Taylor Bubble Flow

Simulation of two-phase flows is relevant for reactor design and safety at normal operation or during accident scenarios. Often, the two-phase flow is in a regime in which slugs are formed or where the flow stratifies. Modeling such situations using standard single-phase Reynolds-averaged Navier-Stokes (RANS) turbulence models fails due to an overestimation of the eddy viscosity at the resolved two-phase interface. To solve this, an ad hoc turbulence damping term has been proposed in the literature that reduces the turbulence production locally at a two-phase interface, analogously to turbulence wall functions. However, this approach must be tailored to the specific setting and does not consider physical contributions such as surface tension or flow topology. Therefore, the problem of two-phase interfacial turbulence must be studied more in-depth. In this work, we consider co-current turbulent Taylor bubble flow using high-fidelity numerical simulation. The Basilisk code is used to simulate a Taylor bubble rising in a vertical pipe. By simulating the bubble in a moving frame of reference, we may study the turbulent kinetic energy (TKE) budgets ahead of the bubble, in its wake, and across the interface. The implementation of the TKE budget computation and the underlying averaging techniques are validated for the single-phase region ahead of the Taylor bubble using reference direct numerical simulation data. The analysis of the TKE budgets in the setting of Taylor bubble flow allows for the study of how turbulence behaves due to the presence of a two-phase interface and, in turn, supports the improvement of two-phase RANS models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Study on ammonia and dimethyl ether oxidation and kinetic interaction up to 100 atm

NH 3 and DME dual fuel oxidation and kinetic coupling are experimentally studied by using a supercritical pressure jet-stirred reactor (SP-JSR) at 20 and 100 atm, over a temperature range of 500–900 K, and at fuel-lean and stoichiometric conditions with NH 3 to DME molar ratios of 4 and 0.62, respectively. An HP-Mech model for high-pressure NH3/DME oxidation is developed based on our previous studies and it shows generally better performance than other models in the literature on high-pressure oxidation experiments. Due to DME’s strong low-temperature reactivity, it dramatically enhances the oxidation of NH3 at low temperature. However, the effect of NH 3 on DME oxidation varies with temperature. At low temperatures, NH 3 inhibits low-temperature reactivity by consuming OH radicals. In addition, the NH 2 and NO x formation from NH 3 further suppresses the low-temperature DME reactivity by reducing alkylperoxyl and O 2 QOOH radicals via RO 2 + NH 2 = RO + H 2 NO, RO 2 + NO = RO + NO 2 and O 2 QOOH + NO = 2CH 2 O + HO 2 + NO 2 . At intermediate temperatures, due to enhanced kinetic coupling of NH 2 /NO x /HO 2 chemistry, DME oxidation is significantly promoted. It is found that there are two major NH 2 /NO x /HO 2 coupling pathways for OH radical production brought by NH 3 : (1) NH 2 + HO 2 = H 2 NO + OH and (2) NH 2 + NO 2 = H 2 NO + NO with NO + HO 2 = OH + NO 2 . Moreover, the resulting H 2 NO will further contribute to OH production via the H 2 NO/NO x /HONO coupling pathways: H 2 NO + NO 2 = HNO + HONO, HNO + NO 2 = HONO + NO, and HONO (+M) = OH + NO (+M). These new NH 2 /NO x /HO 2 and H 2 NO/NO x /HONO pathways play a critical role in promoting DME oxidation at intermediate temperatures. Novelty and Significance Statement: Blending NH 3 with DME can significantly enhance the reactivity of NH 3 and facilitate its practical application in advanced internal combustion engines. NH 3 /DME oxidation is studied in a novel supercritical pressure jet-stirred reactor up to 100 atm, a much higher pressure than previous studies. An interesting non-monotonic effect of NH3 addition on fuel oxidation has been found. Unique kinetic couplings of NH2/HO2/NOx and H2NO/NOx/HONO brough by NH3 and DME blending are identified to play important roles in radical generation and reactivity promotion. The valuable experimental data and key kinetics revealed in the current work can tremendously improve our understanding of NH 3 /DME oxidation via NO x /RO 2 and NH 2 / H 2 NO/HO 2 /NO x kinetic couplings at low to intermediate temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetic modelling of an environmentally friendly carbamazepine synthesis via urea and iminostilbene in batch and continuous processes

Accurate kinetic models for reaction systems allow for improved process understanding and greater quality control, which is particularly beneficial as the pharmaceutical industry shifts from batch to continuous manufacturing (CM). In this work, a first principles kinetic model has been developed for the synthesis of carbamazepine (CBZ) from iminostilbene and urea, starting in a batch reactor and subsequently in a continuous flow reactor. An eco-friendly reaction pathway using urea was selected to avoid the toxic reagents that are typically used for synthesis of CBZ. The kinetic parameters determined from batch reactions were utilized in a MATLAB based kinetic model to simulate the yield for the continuous process. Overall, good agreement between the model prediction and corresponding experimental values was observed for the batch and continuous reaction systems within the full factorial design space. However, the model slightly overpredicted the yield of the continuous reaction system for higher conversion values (>60%) since it did not account for the reverse reaction that can occur at the studied reaction conditions. Here, the use of broken order kinetics was compared with whole number orders, and it was determined that the whole number orders resulted in better agreement at all conversion values for the continuous system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Directly monitoring the shift in corrosion mechanisms of a model FeCrNi alloy driven by electric potential

Abstract Stainless steels are used in a myriad of engineering applications, including construction, automobiles, and nuclear reactors. Developing accurate, predictive mechanistic models for corrosion and electrochemical corrosion kinetics of stainless steels has been a topic of research studies over many decades. Herein, we quantified the aqueous corrosion kinetics of a model austenitic Fe–18Cr–14Ni (wt%) alloy in the presence and absence of applied potential using systematic in situ electrochemical atomic force microscopy (EC-AFM) and transmission electron microscopy (TEM). Without an applied bias, vertical dissolution of corrosion pits is controlled by the surface kinetics/diffusion hybrid mechanism, whereas lateral dissolution is diffusion controlled. When an electric bias is applied, the increase in corrosion rate is dominated by the nucleation of new pits. These insights gained by the in situ EC-AFM will allow applications of this method for a quantitative understanding of corrosion of a wider class of materials.

36 MATERIALS SCIENCE↗

Aqueous-phase effects on ethanol decomposition over Ru-based catalysts

The effects of an aqueous phase on ethanol decomposition for hydrogen production over a Ru(0001) catalyst surface model have been investigated from first principles. Solvent effects on the reaction mechanism and kinetic parameters have been quantified with the help of a microkinetic reactor model, density functional theory, and an implicit solvation scheme (iSMS). Our calculations indicate that in both vapor- and aqueous-phase reaction environments, the ethanol decomposition starts with acetaldehyde formation on the surface, some of which further dehydrogenates to ketenyl species (CHCO), where the C–C bond cleaves to form methylidyne (CH) and CO. In the vapor phase, adsorbed CH gets hydrogenated to methane, and CO desorbs or undergoes methanation reducing the amount of hydrogen produced. In contrast, under aqueous phase reaction conditions, the methanation is inhibited, and the water–gas shift (WGS) reaction is accelerated, leading to complete conversion of CO to CO 2 and H 2 . Furthermore, calculations indicate that the observed reaction behavior under aqueous phase reforming conditions originates primarily from the higher water chemical potential, and implicit solvent models predict only a small solvation effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integrated control of individual plasma scalars with simultaneous neoclassical tearing-mode suppression

A novel integrated-control architecture has been tested in nonlinear, one-dimensional simulations using the control-oriented transport simulator (COTSIM©) and in DIII-D experiments. Integrated architectures that can perform continuous-mission control while also handling off-normal events will be vital in future reactor-grade tokamaks. Continuous-mission controllers for individual magnetic and kinetic scalars (thermal stored-energy (W), volume-average toroidal rotation (Ω Φ ), and safety factor profile (q) at different spatial locations) have been integrated in this work with event-triggered neoclassical tearing-mode (NTM) suppression controllers by combining them into an architecture augmented by a supervisory and exception handling (S&EH) system and an actuator management (AM) system. Here, the AM system, which enables the integration of competing controllers, solves in real time a nonlinear optimization problem that takes into account the high-level control priorities dictated by the S&EH system. The resulting architecture offers a high level of integration and some of the functionalities that will be required to fulfill the advanced-control requirements anticipated for ITER. Initial simulations using COTSIM suggest that the plasma performance and its MHD stability may be improved under integrated feedback control. In addition, the integrated-control architecture has been implemented in the DIII-D plasma control system and tested experimentally for the first time ever in DIII-D in a high-q min scenario, which is a candidate for steady-state operation in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A computational study of the thermodynamic conditions leading to autoignition in nanosecond pulsed discharges

Nanosecond pulsed discharges have attracted the attention of engine manufacturers due to the possibility of attaining distributed ignition sites that accelerate burn rates while resulting in very little electrode erosion. Multidimensional modeling tools currently capture the electrical structure of such discharges accurately, but resolving the chemical structure remains a challenging problem owing to the disparity of time-scales in streamer propagation (nanoseconds) and ignition phenomena (microseconds). The purpose of this study is to extend multidimensional results towards resolving the chemical structure in the wake of streamers (or the afterglow) by using a batch reactor model. This can afford the use of very detailed chemical kinetic information. The full non-equilibrium nature of the electrons is taken into account, along with fast gas heating, shock wave propagation, and thermal diffusion. The results shed light on ignition phenomena brought about by such discharges.

ignition, plasma, kinetics↗

SPERLU Selective Process for Efficient Removal of Lignin and Upgrading (Final Technical Report)

To meet the growing demand for bio-based chemicals and improve the profitability and efficiency of the emerging biorefinery industry, lignin is an abundant and attractive feedstock. Spero Renewables has developed a proprietary and patent-pending technology, the Selective Process for Efficient Removal of Lignin and Upgrading (SPERLU). The technology is a catalytic deconstruction and upgrading of lignin in lignocellulosic biomass or commercially available technical lignin to produce multi-functional phenol (MFP) products. The SPERLU process has been demonstrated and vetted on a lab scale to convert commercially available technical lignin to MFPs in yields of 60-90% (based on lignin). The MFP products can be separated and further upgraded through catalytic or biological means. Spero performed a detailed kinetic study of the SPERLU process and used the resulting kinetic data to design and construct a mini-pilot-scale reactor for routine production of hundreds of grams of MFP products. Spero’s lignin-based MFPs have been demonstrated for use as a drop-in replacement for bisphenol A (BPA) in the synthesis of thermoset polymers. Spero’s lignin-based thermoset polymers exhibit excellent mechanical properties and can be used in many applications. Additionally, a collaboration with the National Renewable Energy Lab (NREL) was used to evaluate the feasibility of upgrading SPERLU products through biological means into valuable chemicals. A comprehensive technoeconomic analysis (TEA) and Life-cycle analysis (LCA) have been completed and support the case for commercialization of the SPERLU technology. Through this project Spero has significantly de-risked the SPERLU technology and is actively planning larger scale pilot projects.

09 BIOMASS FUELS↗

The hard start phenomena in hypergolic engines. Volume 4: The chemistry of hydrazine fuels and nitrogen tetroxide propellant systems

The various chemical reactions that occur and that could possibly occur in the RCS engines utilizing hydrazine-type fuel/nitrogen tetroxide propellant systems, prior to ignition (preignition), during combustion, and after combustion (postcombustion), and endeavors to relate the hard-start phenomenon to some of these reactions are discussed. The discussion is based on studies utilizing a variety of experimental techniques and apparatus as well as current theories of chemical reactions and reaction kinetics. The chemical reactions were studied in low pressure gas flow reactors, low temperature homogeneous- and heterogeneous-phase reactors, simulated two-dimensional (2-D) engines, and scaled and full size engines.

Miron, Y.↗

The kinetics of the reactions of ground-state NH with N 2 O and implications for ammonia combustion

Quantum chemistry and canonical transition state theory have been applied to derive rate constants k for reactions of 3 Σ − imidogen with nitrous oxide. 3 channels were quantified. The fastest leads to NNH+NO and may be summarized as k = 2.8 × 10 7 T 1.68 exp.(−23.1 kcal mol −1 /RT) cm 3 mole −1 s −1 . Two much slower pathways lead to N 2 + 3 HNO and spin-forbidden 1 HNO. Modeling of literature data, from an NH 3 /O 2 flame and jet-stirred reactor measurements on dilute NH 3 /N 2 O mixtures, with the new kinetic information suggests that these channels are too slow to consume N 2 O significantly, and confirms that a widely used proposed rate constant for NH+N 2 O overestimates the reactivity.

Marshall, Paul [University of North Texas, Denton,↗

Radiation damage effects in beryllium for next generation neutrino beam targetry (Final Technical Report)

Current and future high-power accelerators put severe requirements on materials used for target and beam windows and target facilities have been recognized as a critical challenge in development of future particle accelerators. In accelerators, window and target materials are exposed to extreme conditions, which include bombardment with very high energy protons (1- 100 GeV) and thermomechanical shock waves. Radiation can cause direct damage in the material, and it leads to production of transmutation products (especially helium), both phenomena having a potential adverse effect on the stability and durability of the target/window material. At high enough temperatures, He can aggregate to form gas bubbles, which in turn cause significant dimensional changes (swelling), enable easy crack propagation, and eventually cause failure by fracture. On the other hand, if the temperature is too low, radiation damage accumulates in the form of internal defects (e.g., dislocations), leading to hardening and a decreased ductility of the material. In this project, we will focus on beryllium since it is considered to be one of the candidate materials for beam windows and targets in the next-generation proton accelerators, e.g., the Long Baseline Neutrino Facility (LBNF). Radiation effects in Be have been studied in the context of nuclear fusion reactor applications. However, key differences exist between reactor and accelerator conditions, including neutron vs. proton irradiation, continuous vs. pulsed beam flux, much higher energies of bombarding particles in accelerators, and higher operating temperatures for typical reactors. For example, the impact of beam pulsing on the radiation damage and the He bubble kinetics is largely unknown. While results obtained on Be from fusion research might not be directly transferrable to understanding target materials, there is an opportunity to bring state-of-the-art tools from materials research in nuclear reactors to aid design of target and beam window materials in high-power accelerators. To this end, the overarching goal of this project are to develop an experimentally-validated computational framework capable of predicting radiation damage evolution in beryllium relevant to beam window and target conditions, focusing on He bubble formation and growth as a function of irradiation temperature. Our model will be based on the cluster dynamics formalism, where size distribution of defects and He bubbles is simulated as a function of time, temperature, and radiation dose. Parameters for the model will be taken from published experiments and from high-fidelity atomistic simulations proposed in this project. In addition, we will carry out a series of targeted ex-situ and in-situ dual-beam experiments using low-energy protons to provide critical data for validation of the model on the effects of radiation on He clustering, He bubble distribution, and dislocation loop density/size in proton irradiated Be.

36 MATERIALS SCIENCE↗

Gas Hydrate Film Growth in Microfluidic Channels for Carbon Dioxide Capture and Sequestration Applications

Gas or clathrate hydrates are a solid, crystalline compound composed of water and guest molecules that typically form at high pressure and low temperature conditions. Carbon dioxide (CO2) hydrates may be involved in several carbon dioxide capture and sequestration (CCS) applications, including CO2 pipeline transportation and CO2 offshore sequestration. Within these applications, the formation mechanism and kinetics must be well understood to manage the CCS processes, either by preventing or promoting hydrate formation. In this work, a high-pressure glass microfluidic reactor is used in tandem with visual microscopy and in-situ Raman spectroscopy to study both the morphological and kinetic behavior of gas hydrate crystals. Subcooling, pressure, and CO2 flow rate are investigated for their impact on the thickening behavior of pure CO2 hydrates, with flow rate being the only parameter to have a significant effect. Visual and Raman spectroscopy evidence show that both a dense hydrate layer and a porous hydrate layer form, and the latter may provide a path for mass transfer to continue hydrate crystallization. A first principles mass transfer model is developed to describe CO2 hydrate crystal thickening at the interface between gas and water. The impacts of gas impurities and channel wettability are also studied. This method is further applied to investigate the conversion of methane hydrate to CO2 hydrate for combined energy recovery and methane hydrate formation. The authors acknowledge the US Department of Energy Basic Energy Science award # DE-SC0022162.

Wadsworth, Lindsey [Colorado School of Mines, Gold↗