Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “underground”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling of thermal pressurization in tight claystone using sequential THM coupling: Benchmarking and validation against in-situ heating experiments in COx claystone

We apply thermoporoelasticity and a sequentially coupling technique for modeling thermally-driven coupled Thermo-Hydro-Mechanical (THM) processes in tight claystone. A THM benchmark case with a corresponding analytic solution for thermoporoelasticity under a constant heat loading verifies the model. Thereafter, two in situ heating experiments are simulated for model validation: a smaller-scale heating experiment (TED experiment) and a larger-scale experiment (ALC experiment) in Callovo-Oxfordian (COx) claystone at the Meuse/Haute-Marne underground research laboratory in France. The model exhibits good performance to match the observed temperature and pore pressure evolution for the smaller-scale TED experiment. For the larger-scale ALC experiment, general trends of thermal-pressurization are captured in the modeling, but pressure is underestimated at some monitoring points during cool-down. This indicates that the THM response in the field may be affected by the variability of rock's properties or irreversible or time-dependent mechanical processes that are not included in the current thermoporoelastic model. The main contributions of this work are as follows: (1) we verify and validate the numerical simulator, TOUGH-FLAC, to be a valuable coupled THM modeling tool; (2) prove that the laboratory determined material parameters can be used as reference values for upscaling experiments. However, to better identify and quantify THM processes with modeling of in situ tests, more emphasize should be dedicated to obtaining high-quality mechanical deformation data.

58 GEOSCIENCES↗
DOE OSTI2024

Modeling nuclear waste disposal in crystalline rocks at the Forsmark and Olkiluoto repository sites – Evaluation of potential thermal–mechanical damage to repository excavations

We conduct coupled thermo-hydro-mechanical modeling of a KBS-3V repository design in crystalline rocks, using data and conditions from the Forsmark in Olkiluoto repository sites in Sweden and Finland. The study focuses on repository performance related to the impact of thermal and hydraulic evolution on the potential for thermal–mechanical damage to underground repository excavations. For the designs and conditions considered at the Forsmark and Olkiluoto repository sites, the simulations show a peak temperature well under the adopted performance target of a 100°C maximum temperature, whereas there is still a high potential for thermal–mechanical damage to the KBS-3V waste deposition holes. The thermal–mechanical damage is much more likely if rock permeability is so low that it delays saturation and swelling of bentonite-clay-based backfill beyond the time for the thermal–mechanical peak, which occurs 50 to 100 years after nuclear waste deposition. We also found that sidewalls of the KBS-3V emplacement tunnels are vulnerable to tensile fracturing due to the combined effect of thermal stressing and backfill swelling. The study highlights a strong interaction between bentonite-based backfill and host rock through capillary suction along with induced rock desaturation. A careful design and selection of the bentonite-clay-based backfill materials for KBS-3V tunnels and deposition holes can facilitate a timely saturation and backfill swelling that in turn can minimize thermal–mechanical damage.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Nonenergy Biomass Carbon Removal and Storage (BiCRS): Assessing Durability of Nongaseous Carbon Products Across Terrestrial Storage Fates

Biomass Carbon Removal and Storage, or BiCRS, pathways use plants or algae that remove carbon dioxide from the atmosphere through photosynthesis and store it underground or in long-lived products. While some BiCRS approaches generate an energy product, all BiCRS approaches generate a carbon product. A new subset of BiCRS approaches focus on the storage of these raw or converted carbon products for generation of carbon credits. However, the durability of these approaches is highly variable as carbon products vary widely in their “form” and the conditions of their “fate.” We organize our thinking about carbon products and their durability around these two primary axes. The durability of carbon product “forms” is mediated by chemical recalcitrance and ranges substantially across agricultural residues, municipal solid waste, woody biomass, and nongaseous products of thermochemical conversion (e.g., biochars and bio-oils). Meanwhile, terrestrial storage “fates” vary in the mechanism employed to stall decay, including surface storage, dry storage, shallow anoxic storage, and deep or geologic anoxic storage (or injection). Each mechanism has different implications for suitability with different feedstock forms as well as long-term risks. We present a framework for assessing durability of solid or liquid raw and conversion carbon products under terrestrial storage fates, highlighting knowns, unknowns, and research priorities moving forward.

09 BIOMASS FUELS↗
DOE OSTI2017

Adsorption Characteristics of Coals Pyrolyzed at Slow Heating Rates

Underground coal pyrolysis is a process similar to in situ oil shale production wherein heat is applied to deep coal formations to produce light hydrocarbons. Like enhanced coal bed methane, the injection of CO 2 post thermal treatment can improve hydrocarbon recovery and serves as a means of carbon sequestration. Little information exists pertaining to coals pyrolyzed to temperatures expected with this process. This work examines the development of meso- and micropores and their influence on methane and carbon dioxide adsorption and permeability on a thermally treated Utah bituminous coal. Residual tars were found to affect the pore size distributions in pyrolyzed coals. Generally, with increasing treatment temperature, there are more meso- and micropores. A correlation was established between the prevalence of mesopores in pyrolyzed coals and adsorption and, to a lesser extent, permeability. Finally, the treatment temperature of this particular coal is directly related to the amount of CO 2 said coal can store and how the plume of injected CO 2 moves through the formation.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2018

Computational Evaluation of Mg–Salen Compounds as Subsurface Fluid Tracers: Molecular Dynamics Simulations in Toluene–Water Mixtures and Clay Mineral Nanopores

Molecular tracers that can be selectively placed underground and uniquely identified at the surface using simple on-site spectroscopic methods would significantly enhance subsurface fluid monitoring capabilities. To ensure their widespread utility, the solubility of these tracers must be easily tuned to oil- or water-wet conditions as well as reducing or eliminating their propensity to adsorb onto subsurface rock and/or mineral phases. In this work, molecular dynamics simulations were used to investigate the relative solubilities and mineral surface adsorption properties of three candidate tracer compounds comprising Mg–salen derivatives of varying degrees of hydrophilic character. Simulations in water–toluene liquid mixtures indicate that the partitioning of each Mg–salen compound relative to the interface is strongly influenced by the degree of hydrophobicity of the compound. Simulations of these complexes in fluid-filled mineral nanopores containing neutral (kaolinite) and negatively charged (montmorillonite) mineral surfaces reveal that adsorption tendencies depend upon a variety of parameters, including tracer chemical properties, mineral surface type, and solvent type (water or toluene). Simulation snapshots and averaged density profiles reveal insight into the solvation and adsorption mechanisms that control the partitioning of these complexes in mixed liquid phases and nanopore environments. As a result, this work demonstrates the utility of molecular simulation in the design and screening of molecular tracers for use in subsurface applications.

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

Acid Erosion of Carbonate Fractures and Accessibility of Arsenic-Bearing Minerals: In Operando Synchrotron-Based Microfluidic Experiment

Underground flows of acidic fluids through fractured rock can create new porosity and increase accessibility to hazardous trace elements such as arsenic. In this study, we developed a custom microfluidic cell for an in operando synchrotron experiment using X-ray attenuation. The experiment mimics reactive fracture flow by passing an acidic fluid over a surface of mineralogically heterogeneous rock from the Eagle Ford shale. Over 48 h, calcite was preferentially dissolved, forming an altered layer 200–500 μm thick with a porosity of 63–68% and surface area >10$\times$ higher than that in the unreacted shale as shown by xCT analyses. Calcite dissolution rate quantified from the attenuation data was 3 $\times$ 10 –4 mol/m 2 s and decreased to 3 $\times$ 10 –5 mol/m 2 s after 24 h because of increasing diffusion limitations. Erosion of the fracture surface increased access to iron-rich minerals, thereby increasing access to toxic metals such as arsenic. Quantification using XRF and XANES microspectroscopy indicated up to 0.5 wt % of As(-I) in arsenopyrite and 1.2 wt % of As(V) associated with ferrihydrite. This study provides valuable contributions for understanding and predicting fracture alteration and changes to the mobilization potential of hazardous metals and metalloids.

58 GEOSCIENCES↗
DOE OSTI2016

Determination of Organic Partitioning Coefficients in Water-Supercritical CO 2 Systems by Simultaneous in Situ UV and Near-Infrared Spectroscopies

CO2 injected into depleted oil or gas reservoirs for long-term storage has the potential to mobilize organic compounds and distribute them between sediments and reservoir brines. Understanding this process is important when considering health and environmental risks, but little quantitative data currently exists on the partitioning of organics between supercritical CO2 and water. In this work, a high-pressure, in situ measurement capability was developed to assess the distribution of organics between CO2 and water at conditions relevant to deep underground storage of CO2. The apparatus consists of a titanium reactor with quartz windows, near-infrared and UV spectroscopic detectors, and switching valves that facilitate quantitative injection of organic reagents into the pressurized reactor. To demonstrate the utility of the system, partitioning coefficients were determined for benzene in water/supercritical CO2 over the range 35-65 °C and approximately 25-150 bar. Density changes in the CO2 phase with increasing pressure were shown to have dramatic impacts on benzene's partitioning behavior. Our partitioning coefficients were approximately 5-15 times lower than values previously determined by ex situ techniques that are prone to sampling losses. The in situ methodology reported here could be applied to quantify the distribution behavior of a wide range of organic compounds that may be present in geologic CO2 storage scenarios.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2023

Lanthanum-Doped Graphene for Electrocatalytic Reduction of Nitrogen Monoxide

Electrocatalytic conversion of nitrogen oxides is a promising approach to address nitrogen pollution in underground water. Nitrogen oxides, such as nitrate and nitrite, can be reduced to N 2 and NH 3 over an electrocatalyst, with NO as the key intermediate, subsequent reaction of which controls the overall activity and selectivity. Here, we report density functional theory calculations of NO electrocatalytic reduction (NOER) over lanthanum embedded in graphene, through which we show that the localized states in La drive the reaction favorably due to more pronounced molecular adsorption and charge transfer than transition metals such as Co. The free energy profiles are compared between the La-based single-atom catalyst (SAC) and Co-based SAC, for producing N 2 and NH 3 . We find that the La-SAC intensifies the electron transfer to the adsorbed NO, promoting the first protonation step of NO, which is the potential-limiting step on the Co-SAC. Also, an intriguing effect of the water solvent is revealed on the La-SAC. In addition to stabilizing the intermediate species, water molecules that are coordinated with La participate directly in the protonation steps, enhancing the catalyst activity. Furthermore, this study reveals the unique mechanism of NOER over rare-earth-based catalysts, highlighting the potential application of atomically dispersed f-block elements for electrocatalytic reactions.

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

Initial Laboratory Measurements Probing Hydrogen Interactions with Eagle Ford Shale and Pyrite: Potential Implications for Subsurface Hydrogen Storage

Hydrogen (H 2 ) has the potential to be a transformative technology as an enabler to a low-carbon future and promoter of renewable energy. When H 2 is injected and stored in the subsurface, it has the potential to interact with the caprock (usually shale) which overlies and seals the storage reservoir. Here, this study examines geochemical reactions or changes in surface morphology to Eagle Ford Shale, a proxy for caprock, upon exposure to H 2 at 50°C and 10.3 MPa. Reactions were also performed with N 2 to provide an experimental control. Scanning Electron Microscopy with Energy Dispersive Spectroscopy (SEM-EDS), Atomic Force Microscopy (AFM), and Optical Photothermal-Infrared (O-PTIR) Spectroscopy were applied to quantify changes on the microscale and nanoscale level. Fluid chemistry changes were monitored with Ion Chromatography (IC) and Inductively Coupled Plasma Mass Spectrometry (ICP-MS). Exposure of Eagle Ford Shale to H 2 gas alone did not result in any alterations to the shale chemically or any changes in the surface morphology. Exposure of Eagle Ford Shale to both H 2 and water as well as N 2 and water resulted in changes to the surface morphology because of gypsum dissolution and reprecipitation, thus indicating that H 2 is not necessary to promote changes. Pure pyrite was the most reactive with H 2 possibly resulting in a reduction to pyrrhotite. These initial studies suggest that the extent of reactions activated by hydrogen with caprock are minor under the temperature and pressure conditions that would represent underground hydrogen storage.

08 HYDROGEN↗
DOE OSTI2022

Porous Colloidal Nanoparticles as Injectable Multimodal Contrast Agents for Enhanced Geophysical Sensing

Injecting fluids into underground geologic structures is crucial for the development of long-term strategies for managing captured carbon and facilitating sustainable energy extraction operations. Here, we have previously reported that the injection of metal–organic frameworks (MOFs) into the subsurface can enhance seismic monitoring tools to track fluids and map complex structures, reduce risk, and verify containment in carbon storage reservoirs because of their absorption capacity of low-frequency seismic waves. Here, we demonstrate that water-based Cr/Zn/Zr MOF colloidal suspensions (nanofluids) are multimodal geophysical contrast agents that enhance near-wellbore logging tools. Based on experimental fluid-only measurements, MIL-101(Cr), ZIF-8, and UiO-66 nanofluids have distinct complex conductivity and/or low-field nuclear magnetic resonance (NMR) signatures that are relevant to field-deployed technologies, implying the potential to enhance near-wellbore monitoring of CO 2 injection and associated processes with downhole logging tools. Small- and wide-angle X-ray scattering characterization of ~0.5 wt % MIL-101(Cr) suspensions confirmed phase stability and provided insight into the fractal nature of colloidal nanoparticles. Finally, low-field (2 MHz) NMR measurements of MIL-101(Cr) nanofluid injection into a prototypical Berea sandstone demonstrate how paramagnetic high-surface area MOFs may dominate the relaxation times of hydrogen-bearing fluids in porous geologic matrices, enhancing the mapping of near-surface and near-wellbore transport pathways and advancing sustainable subsurface energy technologies.

NMR↗
DOE OSTI2019

Computationally Guided Discovery and Experimental Validation of Indole-3-acetic Acid Synthesis Pathways

Elucidating the interaction networks associated with secondary metabolite production in microorganisms is an ongoing challenge made all the more daunting by the rate at which DNA sequencing technology reveals new genes and potential pathways. Developing the culturing methods, expression conditions, and genetic systems needed for validating pathways in newly discovered microorganisms is often not possible. Therefore, new tools and techniques are needed for defining complex metabolic pathways. Here, we describe an in vitro computationally assisted pathway description approach that employs bioinformatic searches of genome databases, protein structural modeling, and protein–ligand-docking simulations to predict the gene products most likely to be involved in a particular secondary metabolite production pathway. This information is then used to direct in vitro reconstructions of the pathway and subsequent confirmation of pathway activity using crude enzyme preparations. As a test system, we elucidated the pathway for biosynthesis of indole-3-acetic acid (IAA) in the plant-associated microbe Pantoea sp. YR343. This organism is capable of metabolizing tryptophan into the plant phytohormone IAA. BLAST analyses identified a likely three-step pathway involving an amino transferase, an indole pyruvate decarboxylase, and a dehydrogenase. However, multiple candidate enzymes were identified at each step, resulting in a large number of potential pathway reconstructions (32 different enzyme combinations). Our approach shows the effectiveness of crude extracts to rapidly elucidate enzymes leading to functional pathways. Furthermore, results are compared to affinity purified enzymes for select combinations and found to yield similar relative activities. Further, in vitro testing of the pathway reconstructions revealed the “underground” nature of IAA metabolism in Pantoea sp. YR343 and the various mechanisms used to produce IAA. Importantly, our experiments illustrate the scalable integration of computational tools and cell-free enzymatic reactions to identify and validate metabolic pathways in a broadly applicable manner.

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

Modeling the Ionic Strength Effect on Diffusion in Clay. The DR-A Experiment at Mont Terri

Solute diffusion in compacted clays depends on ionic strength through its control on the thickness of the electrical double layer (EDL) on the charged clay surfaces. In the DR-A field experiment (Mont Terri, Switzerland), synthetic porewater was circulated through a borehole for 189 days, leading to the out-diffusion of a variety of tracers into the Opalinus Clay. The borehole solution was then replaced with a higher-salinity solution for an additional 540 days, leading to the diffusion of Cs + , Ca 2+ , Mg 2+ , and Sr 2+ back into the borehole and to an increase in the out-diffusion of anions (I - , Br - ) and 3 H. The experimental results were interpreted using the CrunchClay code, which includes a mean electrostatic potential model for the EDL. The EDL corresponds to a second continuum in addition to bulk electrically neutral porewater. Species-specific diffusion (Nernst-Planck equation) occurs through both domains. A 1D radial model considered a single pore diffusion coefficient (Dp = 10-9 m2/s) for cations and 3H in the bulk porosity, and a smaller D p (3 × 10 -10 m 2 /s) for anions. Dp values in the EDL were smaller (10 -11 m 2 /s), except for Cs + and K + (5 × 10 -10 and 2 × 10 -10 m 2 /s, respectively). The model reproduced well the experimental results and showed the capability to consider temporal changes in geochemical conditions affecting the transport and retention of potentially important radionuclide contaminants (e.g., 137 Cs + , 90 Sr 2+ , 129 I - ) in underground geological nuclear waste repositories. Finally, coupled multicomponent diffusion together with the electrostatic properties of the charged surfaces are essential in the development of predictive models for ion transport in clays.

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

Molecular Insights into Geochemical Reactions of Iron-Bearing Minerals: Implications for Hydrogen Geo-Storage

Here, this study investigates the reaction of hydrogen (H 2 ) with pyrite (FeS 2 ), focusing on how temperature and the presence of water influence the reaction pathways and kinetics. Utilizing computational molecular simulations and kinetic analyses, we explore the impact of these factors on the formation of hydrogen sulfide (H 2 S) and related species. First, grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations reveal that physical H 2 adsorption occurs in distinct layers on the pyrite surface. In addition, increased temperatures reduce the absolute adsorption capacities. Reactive MD simulations demonstrate that H 2 interacts differently with pyrite under varying conditions. At 298 K, H 2 reacts with pyrite to form HS – , leading to the formation of HS – through covalent bonding with sulfur of pyrite. However, no H 2 S is produced at this temperature, suggesting that a kinetic barrier (i.e., activation energy) may prevent this reaction. At higher temperatures, H 2 S production significantly increases. The presence of water introduces additional complexity to the reaction mechanism. Unlike dry conditions, water enhances H 2 S generation, even at low temperatures. Water also facilitates the formation of additional products, such as SOH, indicating a more intricate chemical environment on the pyrite surface. Our findings identify the association of HS – ions to form H 2 S as the rate-limiting step, with temperature influencing this process. This finding suggests that while the presence of water can create a more dynamic reaction environment, the overall mechanisms leading to H 2 S formation remain consistent. These outcomes suggest the need for developing targeted strategies to manage and control H 2 S emissions within the context of underground hydrogen storage.

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

In Situ Study of CO2 and H2O Partitioning Between Na-Montmorillonite and Variably Wet Supercritical Carbon Dioxide

Shale formations play fundamental roles in large-scale geologic carbon sequestration (GCS) aimed primarily to mitigate climate change, and in smaller-scale GCS targeted mainly for CO2-enhanced gas recovery operations. In both technologies, CO2 is injected underground as a supercritical fluid (scCO2), where interactions with shale minerals could influence successful GCS implementation. Reactive components of shales include expandable clays, such as montmorillonites and mixed-layer illite/smectite clays. In this work, we used in situ X-ray diffraction (XRD) and in situ infrared (IR) spectroscopy to investigate the swelling/shrinkage and water/CO2 sorption of a pure montmorillonite, Na-SWy-2, when the clay is exposed to variably hydrated scCO2 at 50 °C and 90 bar. Measured interlayer spacings and sorbed water concentrations at varying levels of scCO2 hydration are similar to previously reported values measured in air at ambient pressure over a range of relative humidities. IR spectra show evidence of both water and CO2 intercalation, and variations in peak shapes and positions suggest multiple sorbed types with distinct chemical environments. Based on the intensity of the asymmetric CO stretching band of the CO2 associated with the Na-SWy-2, we observed a significant increase in sorbed CO2 as the clay expands from a 0W to a 1W state, suggesting that water props open the interlayer so that CO2 can enter. However, as the clay transitions from a 1W to a 2W state, CO2 desorbs sharply. These observations were placed in the context of two conceptual models concerning hydration mechanisms for expandable clays and were also discussed in light of recent theoretical studies on CO2-H2O-clay interactions. The swelling/shrinkage of expandable clays could affect solid volume, porosity, and permeability of shales. Consequently, the results from this work could aid predictions of shale caprock integrity in large-scale GCS, as well as methane transmissivity in enhanced gas recovery operations.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2018

Analytical Model to Detect Fault Permeability Alteration Induced by Fault Reactivation in Compartmentalized Reservoirs

Underground fluid injection induces changes in in situ stress condition of the target formation and local faults that can potentially lead to fault reactivation, which may result in the leakage of injected and/or native fluids into neighboring formations. In this paper, we introduce an analytical method to detect fault reactivation caused by fluid injection into deep faulted aquifers considering across-fault leakage. The fundamental assumption made in our model is that fault permeability will be altered upon fault slip. Therefore, we model fault reactivation as a sudden change in fault permeability at the onset of fault slip. The fault is modeled as a linear interface between two permeable formations with equal rock and fluid properties. The governing equations are coupled through the fault interface and are solved using the Laplace-Fourier integral transform technique. Based on the analytical solution, we find the characteristic bottomhole pressure and pressure derivative responses that enable detecting fault reactivation using diagnostic plots. We observe that pressure derivative undergoes a rapid change at the onset of fault slip followed by a late-time trend to attain a new equilibrium governed by the altered fault permeability. Furthermore, we discuss the evolution of the across-fault leakage rate upon and after fault slippage. Here, the results from this study are presented in the form diagnostic plots and type curves that may be used for reservoir and fault characterization purposes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Predicting Effective Diffusion Coefficients in Mudrocks Using a Fractal Model and Small-Angle Neutron Scattering Measurements

The determination of effective diffusion coefficients of gases or solutes in the water–saturated pore space of mudrocks is time consuming and technically challenging. Yet, reliable values of effective diffusion coefficients are important to predict migration of hydrocarbon gases in unconventional reservoirs, dissipation of (explosive) gases through clay barriers in radioactive waste repositories, mineral alteration of seals to geological CO 2 storage reservoirs and contaminant migration through aquitards. In this study, small angle and very small angle neutron scattering techniques have been utilized to determine a range of transport properties in mudrocks, including porosity, pore size distributions and surface and volume fractal dimensions of pores and grains, from which diffusive transport parameters can be estimated. Using a fractal model derived from Archie's Law, we calculate effective diffusion coefficients from these parameters and compare them to laboratory–derived effective diffusion coefficients for CO 2 , H 2 , CH 4 and HTO on either the same or related mudrock samples. The samples include Opalinus Shale from the underground laboratory in Mont Terri, Switzerland; Boom Clay from a core drilled in Mol, Belgium and a marine claystone cored in Utah, USA. The predicted values were compared to laboratory diffusion measurements. The measured and modeled diffusion coefficients show good agreement, differing generally by less than factor 5. Furthermore, neutron or X–ray scattering analysis is therefore proposed as a novel method for fast, accurate estimation of effective diffusion coefficients in mudrocks, together with simultaneous measurement of multiple transport parameters including porosity, pore size distributions and surface areas, important for (reactive) transport modelling.

58 GEOSCIENCES↗
DOE OSTI2019

Revisiting the Analytical Solutions of Heat Transport in Fractured Reservoirs Using a Generalized Multirate Memory Function

Numerous analytical solutions have been developed for modeling thermal perturbations to underground formations caused by deep-well injection of fluids. Each solution has been derived for a specific boundary value problem and a simplified flow network with one set of parallel fractures. In this paper, new generalized solutions G*( x, s ) are developed using (existing) global transfer functions G$^*_0(x,s)$ and a new memory function g*( s ), where x and s are the space and Laplace variable. The memory function represents the solutions of conductive heat exchange between fractures and matrix blocks and between fractured aquifers and unfractured aquitards. The memory function is developed to account for multirate exchange induced by different shapes, sizes, properties, and volumetric fractions of matrix blocks bounded by multiple sets of orthogonal fractures with different spacing. The global transfer functions represent the fundamental solutions to convective, convective-conductive, and convective-dispersive heat transport in fractures (or aquifers) without exchange and are available for various (1-D linear, 1-D radial, 2-D dipole, and single-well injection-withdrawal) flow fields. The new solutions with exchange are developed using $G^*(x,s) = B^*(s)G^*_0{x, s[1 + Θg^*(s)]}$, thereby greatly simplifying solution development in a novel way, where Θ and B*( s ) are a fracture-matrix scaling factor and the boundary condition function. As a result, the new solutions are applied to several example problems, showing that heat transport in fractured aquifers is significantly impacted by (1) thermal dispersion in fractures that is rarely considered, (2) multirate heat exchange with a wide range of size and anisotropy of rectangular matrix blocks, and (3) heat exchange between aquifers and aquitards.

58 GEOSCIENCES↗
DOE OSTI2019

Volcanoes Erupt Stressed Quartz Crystals

Volcanic eruptions are energetic events driven by the imbalance of magmatic forces. The magnitudes of these forces remain poorly constrained because they operate in regions that are inaccessible, either underground or dangerous to approach. New techniques are needed to quantify the processes that drive eruptions and to probe magma storage conditions. Here we present X-ray microdiffraction measurements of volcanic stress imparted as lattice distortions to the crystal cargo of magma from Yellowstone and Long Valley eruptions. Elevated residual stresses between 100 and 300 MPa are preserved in erupted quartz. Multiple volcanic forces could be culpable for the deformation so we analyzed crystals from pyroclastic falls, pyroclastic density currents, and effusive lavas. In conclusion, stresses are preserved in all quartz but cannot be attributed to differences in eruption style. Instead, lattice deformation likely preserves an in situ measurement of the deviatoric stresses required for the brittle failure of viscous, crystal-bearing glass during ascent.

58 GEOSCIENCES↗