Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relative permeability”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Impacts of Permeability Heterogeneity and Background Flow on Supercritical CO 2 Dissolution in the Deep Subsurface

Motivated by CO 2 capture and sequestration (CCS) design considerations, we consider the coupled effects of permeability heterogeneity and background flow on the dissolution of a supercritical CO 2 lens into an underlying deep, confined aquifer. We present the results of a large-scale Monte Carlo simulation study examining the interaction of background flow rate and three parameters describing multi-Gaussian log-permeability fields: mean, variance, and correlation length. Hundreds of high-resolution simulations were performed using the PFLOTRAN finite volume software to model CO 2 dissolution in a kilometer-scale aquifer over 1,000 years. Predictive dimensionless scaling relationships relating CO 2 dissolution rate to heterogeneity statistics, Rayleigh (Ra) and Péclet (Pe) numbers were developed for both gravitationally dominated free convection to background flow-dominated forced convection regimes. An empirical criterion, Pe = Ra 3/4 , was discovered for regime transition. All simulations converged quickly to a quasi-steady, approximately linear dissolution rate. However, this rate displayed profound variability between permeability field realizations sharing the same heterogeneity statistics, even under mild permeability heterogeneity. In general, increased heterogeneity was associated with a lower mean and higher variance of dissolution rate, undesirable from a CCS design perspective. The relationship between dissolution rate and background flow was found to be complex and nonlinear. Dimensionless scaling relationships were uncovered for a number of special cases. Results call into question the validity of the Boussinesq approximation in the context of modest-to-high background flow rates and the general applicability of numerical simulations without background flow.

42 ENGINEERING↗
DOE OSTI2026

Processing-Dependent Structure and Poroelasticity of Nafion in Liquid Water

Ionomers act as the solid polymer electrolyte membrane in many modern electrochemical devices, yet the role of their nanostructure in modulating the poroelastic response remains poorly understood, especially in liquid water, where few techniques can measure simultaneous transport-mechanical properties. Poroelastic Relaxation Indentation (PRI) is uniquely suited for measuring time-dependent transport-mechanical properties of porous solids, specifically hydraulic diffusivity, elastic modulus, Poisson’s ratio, and intrinsic permeability, for porous solids. While ionomers such as Nafion are not porous in the typical sense, Nafion has a nanophase-segregated structure that, when fully swollen in liquid water, behaves as a poroelastic solid with a coupled mechanical-transport response. Using a poroelastic framework, we investigate how casting and pretreatment of Nafion membranes alter their poroelastic response in liquid environments. We characterize both extruded and dispersion-cast Nafion membranes pretreated in water at 25 or 100 °C and relate the mechanical-transport properties to the ionomer structure via hydrophilic and intercrystalline domain spacing measured using Small-Angle X-ray Scattering (SAXS). At 25 °C, dispersion-cast membranes exhibit distinctly lower hydraulic diffusivity and intrinsic permeability than extruded membranes but with comparable mechanical properties. Pretreatment at 100 °C increases hydrophilic domain spacing, improving transport but at the expense of mechanical stiffness. Dispersion-cast membranes respond more strongly to pretreatment than extruded membranes. Using the Carman-Kozeny pore network model and the hydrophilic domain-spacing, we estimate the pore radius and tortuosity to show how pretreatment reduces structure-related tortuosity differences between dispersion-cast and extruded membranes. Here, in this work, we show that nanophase-segregated materials such as Nafion can be rigorously characterized using poroelasticity, resulting in physically meaningful transport-mechanical properties. Coupling PRI with SAXS provides insights into the nanostructural features that govern the coupled mechanical-transport response. By establishing PRI for a nanophase-segregated material, this approach opens avenues for this technique’s application in other hydrated polymeric materials not typically considered to be poroelastic.

Shen, Margaret [University of California, Berkeley↗
DOE OSTI2020

Characterization and Sampling of Ultralow Permeability Geomaterials using Electrokinetics (LDRD Final Report 209234)

This final report on Laboratory Directed Research and Development (LDRD) project 209234 presents background material for electrokinetics at the pore and porous media scales. We present some theoretical developments related to uncoupling electrokinetic flow solutions, from a manuscript recently accepted into Mathematical Geosciences for publication. We present a summary of two pore-scale modeling efforts undertaken as part of the academic alliance with University of Illinois, resulting in one already submitted journal publication to Transport in Porous Media and another in preparation for submission to a journal. We finally show the laboratory apparatus built in Laboratory B59 in Building 823 and discuss some of the issues that occurred with it.

58 GEOSCIENCES↗
DOE OSTI2024

DECOVALEX-2023: Task D Final Report

Task D of DECOVALEX-2023 is focused on the simulation of the coupled thermal hydraulic-mechanical (THM) behaviour in the full-scale engineered barrier system (EBS). The Horonobe EBS experiment is the demonstration of the full-scale EBS in the underground research laboratory (URL) (performed by JAEA in the Horonobe URL in Japan). Task D consisted of the three steps, a preliminary step (Step 0), simulation of the laboratory tests (Step 1) and simulation of the in-situ full-scale EBS experiment (Step 2). Since the Horonobe EBS experiment demonstrates the vertical emplacement option of the EBS, the experiment gallery is also backfilled with the backfill material. Therefore, interaction between the EBS and the backfill material can also be demonstrated, such as deformation (change of density) of the buffer material. The underground water in the Horonobe URL is saline. This fact adds chemical processes to THM behaviour. For example, mechanical properties (such as swelling pressure of the buffer material and backfill material) and hydraulic properties (such as permeability of the buffer material and backfill material) change depending on the water chemistry. Task D was therefore a challenging Task focused on not only the relatively simple THM behaviour but also complex THM behaviour including chemical processes. Six research teams (BGR, CAS, JAEA, KAERI, SNL and Taipower) participated the Task D. BGR, CAS, JAEA, KAERI and Taipower research teams selected a THM approach, while the SNL research team selected a TH approach. Step 1 involved the simulation of laboratory test results and was important to check the numerical codes developed by the research teams. Step 1 was divided into four sub steps. The simulation results through the Step 1 identified the parameters for simulation of the Step 2. Basic parameters of the materials (buffer material, backfill material, rock mass, concrete, sand) were provided by JAEA. Special parameters which research team needed were identified by back analysis of Step 1. Most notably the mechanical behaviour of swelling and displacement depended on the applied model (elastic model or elastoplastic model). Parameters such as Young’s modulus were found to need smaller values than characterised in the fundamental laboratory test results (Step 1-1, 1-2) for the elastic model. Although laboratory experiments are usually simple, test results contained some error. For example, if the saturation level is 100 % or higher, it should be considered an error. This situation was presented in the Step 1-3. A possible reason is that the buffer material is a mixture of bentonite and silica sand. When a specimen is cut to measure volume or weight, sand grains will affect the measurement data. In Step 2, boundary conditions such as temperature on the surface of the simulated overpack, heater power of the electrical heaters installed in the simulated overpack, injection pressure and inflow rate of the test water, were applied. The outer boundary conditions can be selected using measured data (injection pressure and inflow rate of the test water that is controlled by the injection systems installed in the sand layer around the buffer material and in the boundary between backfill material and concrete support). Since such measured data has some noise, research teams developed their own simplified boundary conditions. Inner boundary conditions can be selected using measured data as heater power and temperature on the surface of the simulated overpack. These data also contain some noise, so research teams developed their own simplified developed boundary conditions. Task D validated various approaches thorough the simulation of the in-situ full scale EBS system including backfill of the gallery: variations in the coupling processes (THM or THC), analysis codes, and boundary conditions. Temperature distribution in the buffer material was simulated well by all research teams. This means thermal behaviour is not sensitive to the simulation approaches. Although the water content distribution on the outside of the buffer material was well simulated by all research teams, the simulation results differ from the measured values inside the buffer material (at the centre and inside, near the simulated overpack). The buffer material is made from tap water, but in the in-situ experiment, saline groundwater infiltrates the buffer material. Therefore, the selection of the hydraulic parameters of the buffer material greatly affects the simulation results of the re saturation behaviour of the buffer material. In the Horonobe EBS experiment, measured values suitable for validating the simulation results were not obtained near the simulated overpack. When simulating the pressure and deformation of the buffer material, the measurement data is easily affected by the installation conditions of the measurement sensors, so verifying the measurement data itself remains an issue. Mechanical simulation results differ depending on whether they are considered as elastic or elastoplastic phenomena. The accuracy of measured in-situ data can be assessed by detailed analysis comparing sampling specimen analysis and measured data. The Horonobe EBS experiment is scheduled to be dismantled in the future (FY2026 and 2027). This detailed dismantling investigation will finally confirm the measured data.

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

Convergence of micro-geochemistry and micro-geomechanics towards understanding proppant shale rock interaction: A Caney shale case study in southern Oklahoma, USA

As a direct outcome of economic development coupled with an increase in population, global energy demand will continue to rise in the coming decades. Although renewable energy sources are increasingly investigated for optimal production, the immediate needs require focus on energy sources that are currently available and reliable, with a minimal environmental impact; the efficient exploration and production of unconventional hydrocarbon resources is bridging the energy needs and energy aspirations, during the current energy transition period. The main challenges are related to the accurate quantification of the critical rock properties that influence production, their heterogeneity and the multiscale driven physico-chemical nature of rock–fluid interactions. A key feature of shale reservoirs is their low permeability due to dominating nanoporosity of the clay-rich matrix. As a means of producing these reservoirs in a cost-effective manner, a prerequisite is creation of hydraulic fracture networks capable of the highest level of continued conductivity. Fracturing fluid chemical design, formation brine geochemical composition, and rock mineralogy all contribute to swelling-induced conductivity damage. The Caney Shale is an organic-rich, often calcareous mudrock. Many studies have examined the impact that clay has on different kinds of shale productivity but there is currently no data reported on the Caney Shale in relation to horizontal drilling; all reported data on the Caney Shale is on vertical wells which are shallow, compared to an emerging play that is at double the depth. Here, in this work we develop geochemical–geomechanical integration of rock properties at micro-and nanoscales that can provide insights into the potential proppant embedment and its mitigation. The novel methodology amalgamates the following: computed X-ray tomography, scanning electron microscopy, energy dispersive spectroscopy, micro-indentation, and Raman spectroscopy techniques. Here, our results show that due to the multiscale heterogeneity in the Caney Shale, these geochemical and structural properties translate into a variation in mechanical properties that will impact interaction between the proppant and the host shale rock.

03 NATURAL GAS↗
DOE OSTI2025

Self-Assembled Oligomers Facilitate Amino Acid-Driven CO 2 Capture at the Air–Aqueous Interface

Direct air capture of CO 2 using amino acid absorbents, such as glycine or sarcosine, is constrained by the relatively slow mass transfer of CO 2 through the air–aqueous interface. Our recent study showed a marked improvement in CO 2 capture by introducing CO 2 -permeable oligo-dimethylsiloxane (ODMS-MIM + ) oligomers with cationic (imidazolium, MIM + ) headgroups. Here, in this work, we have employed all-atom molecular dynamics simulations in combination with subensemble analysis using network theory to provide a detailed molecular picture of the behavior of CO 2 and the glycinate anions (Gly – ) at the ODMS-MIM + decorated air–aqueous interfaces. We show that the cationic head groups of the surfactants enhance the concentration and lifetime of Gly – in the interfacial region, while ODMS tails promote the physisorption of CO 2 in the interfacial region. Together, these two factors increase the effective region of contact and the probability of interactions between CO 2 and Gly – compared to that of the pure air–aqueous interface. The fundamental insights gained in this work establish essential foundations for developing hybrid systems with oligomer-decorated interfaces to maximize the overall CO 2 capture rates.

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

Cellulose acetate membranes exhibit exceptional monovalent to divalent cation selectivities

Salt transport properties of cellulose acetate membranes are reported for a series of chloride salts with mono- valent and divalent cations (LiCl, NaCl, MgCl 2 , CaCl 2 ). Measurements include salt permeability and sorption, with diffusivity values calculated from the permeability and sorption results. We report an exceptionally high LiCl/MgCl 2 selectivity of 750:1. Salts with similar valence (LiCl and NaCl; MgCl 2 and CaCl 2 ) have similar transport properties. The high monovalent/divalent selectivity arises from differences in both sorption and diffusion, with a LiCl/MgCl 2 solubility selectivity of about 11 and a diffusivity selectivity of about 70. Atomistic molecular dynamics simulations show that ions tend to reside in isolated clusters of water. Increasing ion charge strengthens ion–water interactions relative to ion–polymer interactions, explaining the reduced sorption of divalent ions. Diffusion of ions through the membrane occurs via hop-like motion between water clusters. Lithium diffuses faster than magnesium due to weaker ion–water coordination for lithium, which allows for greater mobility within water clusters and more frequent hopping. Altogether, our atomistic simulations suggest that the high LiCl/MgCl 2 selectivity is linked to cellulose acetate’s high water/salt selectivity and is a consequence of low water content and relatively uniform water distribution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Formation of Composite SiC/SiC Joints by Embedded Wire Chemical Vapor Deposition

The joining of ceramic monoliths or composites to date has primarily been limited to the formation of brittle monolithic joints using heterogeneous (dissimilar) materials, similar to brazing in metals. The development of a damage-tolerant joint layer by SiC fiber reinforcements is demonstrated here. Tube workpieces made of SiC fiber-SiC matrix composite are joined using a nonwoven SiC fiber mat densified by embedded wire chemical vapor deposition (EWCVD), creating a fiber-reinforced weld-like joint by homogeneous joining. EWCVD uses a localized heating method to target deposition and growth to the joint region specifically, while minimizing thermal damage to the surrounding composite tube material. X-ray computed tomography (XCT) is used to nondestructively characterize as-made joints for relative density, adhesion, and composition. In situ XCT analysis during mechanical testing revealed crack deflections in the bonding layer, which indicates a toughening mechanism typical of ceramic matrix composite phase. Gas permeation testing of these proof-of-concept composite joints identified relatively high leak rates in comparison to fully coated SiC/SiC composite tube workpieces. In conclusion, the novelty of the composite joining method and current technology challenges, including gas permeability, are discussed in comparison with traditional ceramic joints and materials.

SiC↗
DOE OSTI2020

Co-ion specific effect on sodium halides sorption and transport in a cross-linked poly(p-styrene sulfonate-co-divinylbenzene) for membrane applications

Here, sodium halide sorption and transport in a cation exchange membrane based on cross-linked poly(p-styrene sulfonate-co-divinylbenzene) (CR61) were investigated experimentally and theoretically. At fixed co-ion valence, sodium halide partitioning between aqueous electrolyte solutions and ion exchange polymers is only minimally sensitive to co-ion chemical identity. NaCl, NaF and NaI exhibit similar solubility coefficients in CR61, which agrees reasonably well with predictions of the Donnan-Manning model, especially at high external salt concentrations. Ion activity coefficients in the membrane were evaluated and ion sorption thermodynamics were investigated. Ion diffusion coefficients in CR61 were predicted using a recently developed framework based on the Manning diffusion model, which accounts for electrostatic effects on ion diffusion, and the Mackie-Meares model, which accounts for tortuosity/area effects on ion diffusion in swollen polymers. Model predictions agree reasonably well with the experimental data, especially at high external salt concentrations. The observed decrease in salt diffusion coefficients in CR61 relative to values in aqueous solutions is due primarily to tortuosity/area effects. Salt diffusion coefficients correlate with co-ion hydrated radius and membrane water content. Finally, salt permeability coefficients predicted using the solution-diffusion model exhibit good agreement with experimental data, especially at external salt concentration greater than 0.1 mol/L.

42 ENGINEERING↗
DOE OSTI2021

Can we safely go to 200 °C? An integrated approach to assessing impacts to the engineered barrier system in a high-temperature repository

Abstract. This presentation gives on overview of the complex thermo-hydro-mechanical-chemical (THMC) processes occurring during the disposal of heat-producing high-level radioactive waste in geologic repositories. A specific focus is on the role of compacted bentonite, which is commonly used as an engineered backfill material for emplacement tunnels because of its low permeability, high swelling pressure, and radionuclide retention capacity. Laboratory and field tests integrated with THMC modeling have provided an effective way to deepen our understanding of temperature-related perturbations in the engineered barrier system; however, most of this work has been conducted for maximum temperatures around 100 ∘C. In contrast, some international disposal programs have recently started investigations to understand whether local temperatures in the bentonite of up to 200 ∘C could be tolerated with no significant changes to safety relevant properties. Raising the maximum temperature is attractive for economical and safety reasons but faces the challenge of exposing the bentonite to significant temperature increases. Strong thermal gradients may induce complex moisture transport processes while geochemical processes, such as cementation and perhaps also illitization effects may occur, all of which could strongly affect the bentonite and near-field rock properties. Here, we present initial investigations of repository behavior exposed to strongly elevated temperatures. We will start discussing our current knowledge base for temperature effects in repositories exposed to a maximum temperature of 100 ∘C, based on data and related modeling analysis from a large heater experiment conducted for over 18 years in the Grimsel Test Site in Switzerland. We then show results from coupled THMC simulations of a nuclear waste repository in a clay formation exposed to a maximum temperature of 200 ∘C. We also explore preliminary data from a bench-scale laboratory mock-up experiment, which was designed to represent the strong THMC gradients occurring in a “hot” repository, and we finally touch on a full-scale field heater test to be conducted soon in the Grimsel Test Site underground research laboratory in Switzerland (referred to as HotBENT).

Birkholzer, Jens T. (ORCID:0000000279891912)↗
DOE OSTI2023

GeoThermalCloud for EGS – An Open-source, User-friendly, Scalable AI Workflow for Modeling Enhanced Geothermal Systems

Enhanced Geothermal Systems (EGS) offer a vast potential to expand the use of geothermal energy. Heat is extracted from this engineered system by injecting relatively cold water into subsurface fractures, which are in contact with hot dry rock, and brought back to surface through production wells. Creating EGS requires improving the natural permeability of hot crystalline rocks. In this short conference paper, we present a reproducible workflow for modeling EGS. Our workflow called the GeoThermalCloud (GTC) for EGS, leverages recent advances in machine learning, deep learning, and high-performance computing. This GTC framework is currently being made open-source, user-friendly, and reproducible through python scripts as well as Google Colab/Jupyter Notebooks. This GTC for EGS modeling scripts are made available at https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS and will constantly be updated to cater for geothermal community. Current GTC framework provides scripts to train deep learning (DL) models for techno-economics and data worth analysis. The Geothermal Design Tool (https://github.com/GeoDesignTool/GeoDT.git), a fast and simplified multi-physics solver, is used to develop a database for training DL models. This short paper provides details on the scripts to curate, process, and train DL models. The scripts can easily be modified to train on databases generated by other popular open-source simulators such as PFLOTRAN, STOMP, TOUGH, and GEOSX or commercial software such as ResFrac and COMSOL.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Analysis of Orthogonal Efflux and Permeation Properties of Compounds Leads to the Discovery of New Efflux Pump Inhibitors

Optimization of compound permeation into Gram-negative bacteria is one of the most challenging tasks in the development of antibacterial agents. Two permeability barriers–the passive diffusion barrier of the outer membrane (OM) and active drug efflux–act synergistically to protect cells from the antibacterial action of compounds. In Escherichia coli (E. coli) and relatives, these two barriers sieve compounds based on different physicochemical properties that are defined by their interactions with OM porins and efflux pumps, respectively. In this study, we critically tested the hypothesis that the best substrates and inhibitors of efflux pumps are compounds that can effectively permeate the OM and are available at relatively high concentrations in the periplasm. For this purpose, we filtered a large subset of the ZINC15 database of commercially available compounds for compounds containing a primary amine, a chemical feature known to facilitate the uptake through E. coli general porins. The assembled library was screened by ensemble docking to AcrA, the periplasmic component of the AcrAB-TolC efflux pump, followed by experimental testing of the top predicted binders for antibacterial activities, efflux recognition, and inhibition. In this work, we found that the filtered primary amine library is a rich source of compounds with efflux-inhibiting activities and identified efflux pump inhibitors with novel chemical scaffolds effective against E. coli AcrAB-TolC and efflux pumps of multidrug-resistant clinical isolates of Acinetobacter baumannii. However, primary amines are not required for the recognition of compounds by efflux pumps and their efflux-inhibitory activities.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Coreflood Tests to Evaluate Enhanced Oil Recovery Potential of Wettability-Altering Surfactants for Oil-Wet Heterogeneous Carbonate Reservoirs

Summary Oil-wetness and heterogeneity are two main factors that result in low oil recovery (OR) by waterflood in carbonate reservoirs. The injected water is likely to flow through high-permeability regions and bypass the oil in the low-permeability matrix. In this study, systematic coreflood tests were carried out in both “homogeneous” cores and “heterogeneous” cores with a wettability-altering surfactant. The homogeneous coreflood tests were conducted to evaluate surfactant retention, as well as to compare tertiary surfactant flooding with secondary surfactant flooding. The heterogeneous coreflood test was proposed to model bypassing in low-permeability matrix during waterfloods, and dynamic imbibition of surfactant into the low-permeability matrix. Surfactant retention results suggest that retention increases as initial oil saturation decreases. The retention of selected surfactant in the target reservoir cores was measured to be within a range of 0.07–0.12 mg/g-rock, which is economically acceptable. The results of homogeneous coreflood tests showed that both secondary waterflood and secondary surfactant flood can achieve high OR (>50%) from relatively homogeneous oil-wet cores. A shut-in phase after the surfactant injection resulted in a surge in oil production, which suggests that enough time should be given for wettability alteration by surfactants. The results of heterogeneous coreflood tests showed that more oil is bypassed in the tighter matrix by waterflood if the permeability is higher in the flooded layer and this bypassed oil is the target for the wettability-altering surfactant floods. Slow wettability-altering surfactant injection leads to imbibition into bypassed regions. When the oil-wet carbonate reservoirs have large unswept regions after waterflood, wettability-altering surfactants can significantly improve OR if enough time is given for imbibition.

Engineering↗
DOE OSTI2022

Binder jet additive manufacturing of ceramic heat exchangers for concentrating solar power applications with thermal energy storage in molten chlorides

Triply periodic minimal surface (TPMS) geometries can only be fabricated by additive manufacturing methods and are of interest for heat exchangers. Ceramic TPMS heat exchangers can operate at higher temperatures and pressures with superior performance and increased operating efficiencies compared to metal heat exchangers. The properties of ultra-high temperature ceramic (UHTC) materials are also favorable for heat exchangers and potentially suited to concentrated solar power (CSP) systems, such as those based on a molten chloride salt thermal energy storage (TES) medium used to heat CO 2 in a closed-loop Brayton power cycle. We intended to demonstrate binder jet additive manufacturing feasibility of a UHTC-TPMS structure by printing and sintering a null candidate. We aimed to achieve parts with a relative density ≥ 92 % of theoretical and to provide a TPMS part demonstration. The target density indicates the transition from intermediate to final stage sintering, a requirement to inhibit gas permeability and for sintering complex near net shapes to full density with sinter-HIP technology. The goal of the TPMS part demonstration was to determine if printing and sintering parameters developed from test coupons apply to the complex geometries that will eventually be used in a heat exchanger design. Our objective was to print cubic TPMS parts with a 9 cm 3 volume and sinter it without distorting and cracking. We report we were able to sinter ZrB 2 -MoSi 2 composite parts based on the Schwarz-D TPMS and achieve isotropic shrinkage up to 60 % by volume, resulting in densities ranging from 92 % to 96 % of theoretical.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Oxidation and Associated Pore Structure Modification During Experimental Alteration of Granite

Weathering plays a crucial role in a number of environmental processes, and the microstructure and evolution of multi-scale pore space is a critically important factor in the weathering of geological formations. In igneous rocks the infiltration of meteoric water into initially relatively dry material can cause the rock to begin to disaggregate, increasing porosity and surface area, and allowing further disaggregation and weathering to occur. These processes, in turn, allow biota to colonize the rock, further enhancing the weathering rate. In some igneous rocks this may be driven by oxidation of primary minerals. One such mineral, biotite, has been repeatedly mentioned as a cause of cracking during oxidation. However, the scale-dependence of the processes by which this occurs are poorly understood. We cannot, therefore, accurately extrapolate laboratory reaction rates to the field in predictive numerical models. In order to better understand the effects of oxidation and test the hypothesis that fracture and disaggregation are initiated by swelling of oxidizing biotites, we reacted granite cores in a selenic acid-rich aqueous solution at 200°C for up to 438 days. Elevated temperatures and selenic acid were used to provide relatively fast reaction rates and highly oxidizing conditions in sealed reaction vessels. These experiments were analyzed using a combination of imaging, X-ray diffraction, Mössbauer spectroscopy, and small- and ultra-small angle neutron scattering to interrogate porosity and microfracture formation. The experimental results show little observable biotite swelling, but significant transport and growth of iron oxides and/or clays along grain boundaries throughout the sample. Significant increases in porosity were also observed at the sample rim, likely associated with feldspar alteration. Fractures and transport were observed throughout the core, suggesting that stresses due to crystallization pressures caused by the growing iron phases may be the initiating factors in granite weathering, possibly followed by biotite swelling after sufficient permeability is achieved.

Anovitz, Lawrence M.↗
DOE OSTI2024

The vertical-velocity skewness in the atmospheric boundary layer without buoyancy and Coriolis effects

One of the main features of near-neutral atmospheric boundary layer (ABL) turbulence is the positive vertical velocity skewness $Sk_w$ above the roughness sublayer or the buffer region in smooth-walls. The $Sk_w$ variations are receiving renewed interest in many climate-related parameterizations of the ABL given their significance to cloud formation and to testing sub-grid schemes for Large Eddy Simulations (LES). The vertical variations of $Sk_w$ are explored here using wind tunnel and flume experiments collected above smooth, rough, and permeable-walls in the absence of buoyancy and Coriolis effects. These laboratory experiments form a necessary starting point to probe the canonical structure of $Sk_w$ as they deal with a key limiting case (i.e., near-neutral conditions). Diagnostic models based on cumulant expansions, realizability constraints, and constant mass flux approach routinely employed in the convective boundary layer as well as prognostic models based on third-order budgets are used to explain variations in $Sk_w$ for the idealized laboratory conditions. The failure of flux-gradient relations to model $Sk_w$ from the gradients of the vertical velocity variance σ$_w^2$ are explained and corrections based on models of energy transport offered. Novel links between the diagnostic and prognostic models are also featured, especially for the inertial term in the third-order budget of the vertical velocity fluctuation. The co-spectral properties of w′/σ w vs w′ 2 /σ$_w^2$ are also presented for the first time to assess the dominant scales governing $Sk_w$ in the inner and outer layers, where w′ is the fluctuating vertical velocity and σ w is the vertical velocity standard deviation.>

Boundary layer flow↗
DOE OSTI2023

The influence of permeability anisotropy in the upper ocean crust on advective heat transport by a ridge-flank hydrothermal system

Here, in this study, we highlight the importance of permeability anisotropy on the hydrogeological regime of a ridge-flank hydrothermal system. Our study site, North Pond, is a marine sediment pond on ~8 Ma seafloor in the North Atlantic, and represents a low-temperature, end-member ridge-flank hydrothermal system. Previous simulations of North Pond elucidated long-standing hypotheses concerning hydrothermal fluid and heat transport in the upper volcanic crust but failed to fully explain observed patterns of seafloor heat flux in this area. Here we use variography, a geostatistical method, to quantify relations between seafloor heat-flux measurements, and coupled numerical simulations of fluid and heat flow to simulate the hydrogeologic regime. Directional variography shows that heat-flux observations are correlated along-strike of the regional crustal fabric. Three-dimensional simulations that include permeability anisotropy are able to replicate seafloor heat-flux patterns across North Pond. The simulations that result in the best match to thermal data incorporate permeability anisotropy in the horizontal plane. We find that the feedback between permeability anisotropy and the asymmetric geometry of North Pond combine to promote advective removal of heat and mass within the crustal aquifer. These findings suggest that permeability anisotropy in the oceanic crust may influence ridge-flank hydrothermal circulation more broadly.

58 GEOSCIENCES↗
DOE OSTI2022

Effects of time-dependent deformation of shale on the integrity of a geological nuclear waste repository

Safety assessment of geological nuclear waste repositories is essential for the permanent disposal of spent nuclear fuels and high-level radioactive waste. The long-term integrity of the host rock as well as the engineered barrier system (e.g., bentonite buffer) surrounding nuclear waste canisters is of particular importance as decay heat from nuclear waste canisters, which will last over thousands of years, significantly disturbs the thermo-hydromechanical (THM) state of the repository. In this study, THM coupled simulations were carried out to investigate the effect of time-dependent deformation (i.e., creep) of shale on the long-term integrity of a generic subsurface nuclear waste repository. The Norton-Bailey creep model, which is also known as the Lemaitre-Menzel-Schreiner model, was employed to simulate the power-law type creep that is observed in shales. The TOUGH-FLAC simulator was employed for the THM coupled modeling of the repository. The objective of this study is to assess the effect of creep in different shales (i.e., high creep shale vs. low creep shale) on long-term stress and permeability changes in the repository. Results show potential advantages of constructing repositories in the high creep shale, as deviatoric stress levels in the formation decreased to zero in 100 years since the emplacement of nuclear waste canisters and the permeability also decreased to the undamaged, intrinsic levels in 10,000 years. Also, mean effective stress levels in the bentonite buffer increased by 100% in the high creep shale case relative to the low creep shale case at 10,000 years due to creep-induced contraction of the nuclear waste disposal tunnel. However, in earlier periods (e.g., 1000 years), the stress levels in the bentonite buffer were twice smaller in the high creep shale case than in the low creep shale case, which shows a tradeoff between the intermediate- (~1000 years) and long-term (>10,000 years) compaction levels in the bentonite buffer depending on the creep characteristics of the host shale.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗