Results for “Pressure effects” — Page 40

DOE OSTI2024

Simulation of triggering and evolution of ELM by pellet injection in EAST under BOUT++ framework

A BOUT ++ three-field magnetohydrodynamic model is employed to study the triggering and evolution of edge localized mode (ELM) by Li pellets injected along the outer mid-plane in the EAST configuration. The linear simulation shows that compared with a large deposition on the pedestal top (scenario I), a smaller deposition within the steep-gradient pedestal region (scenario II) can stimulate much larger linear growth rates of all-n peeling-ballooning modes (PBMs). The nonlinear simulation shows that there exists a pellet size threshold for ELM triggering for two deposition locations; the threshold for scenario I predicted in the present study matches the EAST observation well. Comparison of the two scenarios reveals that a smaller deposition is sufficient to trigger an ELM in a much shorter time in scenario II, whose ELM size is comparable to that in scenario I. This conclusion confirms previous DIII-D and ASDEX-Upgrade observations, suggesting that the steep-gradient pedestal region is a favorable deposition location for ELM triggering with minimum pellet size. Simulation analyses also find that the positive radial gradient of the hump-like pressure profile in the outer mid-plane induced by the pellet deposition plays a different role in the two scenarios. In scenario I, the force resulting from the gradient hinders the outflow of core plasmas and in return, the perturbation is suppressed from spreading inwards after ELM crashes. In scenario II, with a sizable deposition, the gradient results in another competitive perturbation growth region during the linear phase, thus dispersing the free energy and reducing the efficiency of destabilizing PBMs by pellet injection. The suppressing effect of saturated zonal flow on other modes, the short ELM fast crash phase, and the restricting transport effect of the positive radial pressure gradient work together to constrain the pedestal energy loss, especially when the pellet deposition amount is high.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

High Pressure Infiltration/Expulsion of Aqueous NaCl in Planar Hydrophobic Nanopores

Pressure-driven permeation of water in a poorly wettable material results in a conversion of mechanical work into surface free energy representing a new form of energy storage, or absorption. When water is replaced by a concentrated electrolyte solution, the storage capacity of a nanoporous medium becomes comparable to high-end supercapacitors. The addition of salt can also reduce the hysteresis of the infiltration/expulsion cycle. Our molecular simulations provide a theoretical perspective into the mechanisms involved in the process, and underlying structures and interactions in compressed nanoconfined solutions. Specifically, we consider aqueous NaCl in planar confinements of widths of 1.0 nm and 1.64 nm and pressures of up to 3 kbar. Open ensemble Monte Carlo simulations utilizing fractional exchanges of molecules for efficient additions/removal of ions have been utilized in conjunction with pressure-dependent chemical potentials to model bulk phases under pressure. Confinements open to these pressurized bulk, aqueous electrolyte phases show reversibility at narrow pore sizes and irreversibility in wider ones, consistent with experiment, as well as strong hysteresis at both pore size. The addition of salt results in significant increases in the solid/liquid interfacial tension in narrower pores and associated infiltration and expulsion pressures. Furthermore, these changes are consistent with strong desalination effects at the lower pore size observed irrespective of external pressure and initial concentration.

25 ENERGY STORAGE↗

DOE OSTI2023

Fiber Optic Temperature Sensor System Using Air-Filled Fabry–Pérot Cavity with Variable Pressure

We report a high-resolution fiber optic temperature sensor system based on an air-filled Fabry–Pérot (FP) cavity, whose spectral fringes shift due to a precise pressure variation in the cavity. The absolute temperature can be deduced from the spectral shift and the pressure variation. For fabrication, a fused-silica tube is spliced with a single-mode fiber at one end and a side-hole fiber at the other to form the FP cavity. The pressure in the cavity can be changed by passing air through the side-hole fiber, causing the spectral shift. We analyzed the effect of sensor wavelength resolution and pressure fluctuation on the temperature measurement resolution. A computer-controlled pressure system and sensor interrogation system were developed with miniaturized instruments for the system operation. Experimental results show that the sensor had a high wavelength resolution (<0.2 pm) with minimal pressure fluctuation (~0.015 kPa), resulting in high-resolution (±0.32 ℃) temperature measurement. It shows good stability from the thermal cycle testing with the maximum testing temperature reaching 800 ℃.

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Scaled ILU Smoothers for Navier-Stokes Pressure Projection

Incomplete LU (ILU) smoothers are effective in the algebraic multigrid (AMG) V-cycle for reducing high-frequency components of the error. However, the requisite direct triangular solves are comparatively slow on GPUs. Previous work has demonstrated the advantages of Jacobi iteration as an alternative to direct solution of these systems. Depending on the threshold and fill-level parameters chosen, the factors can be highly nonnormal and Jacobi is unlikely to converge in a low number of iterations. We demonstrate that row scaling can reduce the departure from normality, allowing us to replace the inherently sequential solve with a rapidly converging Richardson iteration. There are several advantages beyond the lower compute time. Scaling is performed locally for a diagonal block of the global matrix because it is applied directly to the factor. Further, an ILUT Schur complement smoother maintains a constant GMRES iteration count as the number of MPI ranks increases, and thus parallel strong-scaling is improved. Our algorithms have been incorporated into hypre, and we demonstrate improved time to solution for linear systems arising in the Nalu-Wind and PeleLM pressure solvers. For large problem sizes, GMRES+AMG executes at least five times faster when using iterative triangular solves compared with direct solves on massively parallel GPUs.

algebraic multigrid↗

DOE OSTI2021

Characterization of Oil Well Cement Performance during Early Hydration under Simulated Borehole Conditions

Summary Experiments on oil well cement (OWC) slurries were performed using the newly developed laboratory-scale wellbore simulation chamber (WSC). The WSC can simulate hydrostatic pressure reduction in the cemented annulus and possible gas migration under representative conditions. Forensic analysis shows that pressurized fluids can result in porous cement and gas channeling during cement slurry gelation. The effects of different factors on slurry pore pressure were also studied, including formation permeability, initial overburden pressure (OBP) representing the depth of interest, wellbore temperature, water/cement (w/c) ratio, cement composition, and the use of a calcium chloride (CaCl2)-based accelerator. By analyzing the temperature history of hydrating cement using degree of hydration, the evolution of cement hydration was characterized for slurry designs cured at different hydration rates. This provides the opportunity to parameterize the slurry designs and other important factors associated with wellbore conditions.

Engineering↗

DOE OSTI2023

Low Cost Glass-Ceramic Matrix Composite Heat Exchanger

As part of ARPA-E’s High Intensity Thermal Exchange through Materials and Manufacturing Processes (HITEMMP) program, this project sought to develop novel heat exchanger (HX) capabilities to enable efficient and power dense power generation cycles. This class of HX comes under the category of ceramic/composite materials with the higher temperature goal in the program of ≥1100 °C inlet temperature operation. The enabling capability of this effort is the use of glass-ceramic matrix composite (GCMC) material which provides the high temperature durability of a ceramic, the flaw tolerance of a composite, a significantly faster and lower cost manufacturing process than conventional matrix CMCs and very low porosity levels < 0.5%. For thin-walled HX structures and the need to minimize leakage, the low porosity differentiator is particularly important. RTRC has prior experience with this material system and in the current project advanced the component design and manufacturing methods into new territory to produce features required for effective heat exchange under high pressures. In this approach, silicon carbide fiber is fabricated into a fiber preform using various textile processes. Graphite tooling is used both during the build-up of the fiber preform (interior tooling) and after the fiber preform has been completed (exterior tooling). This tooling assembly is heated to high temperature in an environment that has been evacuated and backfilled with inert gas. A reservoir of specialty glass is present and once the desired temperature has been reached to achieve the desired glass viscosity, an actuator distributes the glass throughout the fiber preform using passageways which are part of the tooling design in a process known as glass transfer molding. After the tooling has been removed, the composite is heat treated to convert the amorphous glass to a crystalline ceramic, providing improved properties. The project was divided into three phases focusing on the following: 1) 10 kW HX design and coupon-level tube sheet fabrication, 2) 10 kW HX fabrication, 3) 50 kW HX fabrication. During Budget Period 1 (BP1), additional risks were encountered and the need for additional funds was agreed upon by ARPA-E program leadership. Due to a variety of factors, the contract modification required nominally 18 months to execute at which time the HITEMMP program was effectively concluding. Because of this and the time that would be required to perform BP2 tasks, it was decided to conclude the project at the end of BP1. During the design of the 10 kW HX, manufacturing constraints were learned and incorporated, leading to a revised configuration for the fiber preform and HX. Heat exchange and pressure drop predictions also played a role in modifying the original design concept to be a higher aspect ratio shell-and-tube HX, simplifying the manufacturing process and improving the heat exchanger performance. Good gravimetric and volumetric thermal power densities of 11.2 kW/kg and 10,200 kW/m3 for the entire HX were projected that involved thermo-structural Finite Element Analysis to determine the structural mass needed for the high operation pressures of 250 bar cold inlet and 80 bar hot inlet. Fiber preforms using textile processes were produced for multiple headered tube sheets. Additional challenges were encountered during the glass transfer molding step for which solutions were identified, but programmatics did not allow them to be implemented in BP1. While complete HX test articles were not fabricated, the benefits of this GCMC material for a variety of high temperature applications remain.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2026

Experimental and Computational Study of Weld Drop-Through Hydraulic Effects in Multi-Layer Concentric Pipes Servicing Cryogenic Liquid Hydrogen

The Cryogenic Moderator System (CMS) at the Oak Ridge National Laboratory (ORNL)-Spallation Neutron Source (SNS) supplies liquid hydrogen, through vacuum insulated transfer lines to the moderators. The welding process used to make the transfer lines produces weld drop-through, which can constrict the flow of hydrogen and/or cause thermal bridges between adjacent layers. To ensure proper operation of the CMS it is therefore necessary to conduct a flow test to verify the pressure drop caused by the welds. Measurements in as-built piping sections intended to provide liquid hydrogen service to the SNS-CMS system were performed to test the effect of weld penetrations into the fluid flow area focusing on pressure drop effects. The reference test sections were designed using prototypical piping dimensions, materials and expected penetration welds. A numerical study using Computational Fluid Dynamics (CFD) tools in the transfer lines provided a conversion factor from the pressure drop measured in the lines with water to the expected pressure drop using liquid cryogenic hydrogen. The results show evidence of a maximum pressure drop modification in water of 78 KPa in the tested parts, which scales to liquid hydrogen as PHydrogen = 5.2 KPa. The maximum measured equivalent blockage area percentage was found in the as-built TDM transfer line spool 2 with a value of 31%. The measured values obtained in the as-built transfer lines provided accurate values to estimate the expected total pressure drop in the CMS system providing operational limits for the current recirculatory capacity. The effect of the weld drop-through in hydrogen evacuation capacity scenarios was better evaluated as a result of this study.

Dominguez-Ontiveros, Elvis [ORNL] (ORCID:000000018↗

DOE OSTI2021

Effects of Hydrogen on the Phase Relations in Fe-FeS at Pressures of Mars-Sized Bodies

The large radius, and therefore low density, of the Martian core found in the InSight mission data analysis highlights the importance of considering other light elements besides sulfur (S), which has been considered as the main light element for Mars for decades. Hydrogen (H) is abundant in the solar system and becomes siderophile at high pressures. Although Fe-S and Fe-H systems have been studied individually, the Fe-S-H ternary system has only been investigated up to 16 GPa and 1723 K. We have investigated the Fe-S-H system at pressures and temperatures (P-T) relevant to the cores of Mars-sized planets (up to 45 GPa and well above the melting temperature of FeS) in the laser-heated diamond anvil cell combined with in situ synchrotron X-ray diffraction. Here we found that sufficient hydrogen leads to the disappearance of Fe 3 S at high P-T. Instead, separate Fe-H and Fe-S phases appear at 23–35 GPa. At pressures above 35 GPa, we found a new phase appearing while Fe-S phases disappear and Fe-H phases remain. Our analysis indicates that the new phase likely contains both S and H in the crystal structure (tentatively FeS x H y where x ≈ 1 and y ≈ 1). The observed pressure-dependent changes in the phase relation may be important for understanding the structure and dynamics of the Martian core and the cores of Mars-sized exoplanets.

58 GEOSCIENCES↗

DOE OSTI2024

Quantum effects in H-bond symmetrization and in thermodynamic properties of high pressure ice

Not provided.

Materials Science↗

DOE OSTI2021

Tritium Transport Phenomena in Molten-Salt Reactors: Molten Salt Tritium Transport Experiment Design

Tritium is produced from neutron interactions with both lithium and beryllium. Large quantities of tritium are generated in Molten Salt Reactors (MSRs) which use LiF/BeF2 (FLiBe) as the fuel salt. Tritium is unique among the radionuclide hazards as it readily permeates through metal structural materials at high temperatures. All metal surfaces are potential release paths for tritium. For adequate safety analysis and eventual licensing of new reactors, predictive models for tritium transport and release from MSRs must be developed. These models must account for the multiple transport phenomena involved with tritium: fuel salt phase mass transport, dissociation/recombination reactions on metal surfaces, interstitial diffusion through the metal structure, and salt or gas phase mass transport in the downstream fluid. These models also must also be validated with representative experiments. Our previous report outlined tritium transport phenomena involved in MSRs, made suggestions on gaps in the transport dataset, and proposed an experimental test stand to test combined transport effects – tritium transport through pipe walls in a convective salt flow. In this report, we summarize an updated analysis framework for tritium transport in MSRs, report our results on hydrogen and deuterium permeation through Hastelloy N, and describe the final design of the Molten Salt Tritium Transport Experiment (MSTTE, pronounced “misty”). The MSRE provides the only wholistic experimental data set for tritium transport in MSRs and understanding the transport phenomena involved in the MSRE is crucial for future model development. One set of parameters in our analysis framework was unknown for the MSRE—surface reaction rates for tritium on Hastelloy N. This warranted our hydrogen and deuterium permeation campaign to assess the permeability, diffusivity, and solubility of hydrogen isotopes in clean Hastelloy N. Surface reaction rate constants were probed by low pressure measurements, however, no surface effects were observed in the limits of our permeation apparatus. Permeation experiments on oxidized Hastelloy N were not performed for this report but are planned in future work. The experimental test stand, MSTTE, measures combined transport properties of the salt-metal system. MSTTE is a forced convection FLiBe loop with custom designed test section to measure tritium transport through candidate structural materials. We use MSRE relevant dimensionless numbers to design and scale the test section. Hastelloy N is a candidate loop and test section material due to the relevance for the MSRE and related designs, however, other metals are being considered (e.g. 316H SS) which may better align with current vendor concepts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Hydraulic Response to Thermal Stimulation Efforts at Raft River Based on Stepped Rate Injection Testing

The injection well stimulation project at the Raft River geothermal field tests the effect of long-term cold water injection and high pressure injection on well injectivity, improvements to which could reduce operating costs. The primary data for analysis and interpretation of the injection test are step-rate flow tests run before each new phase of the injection. These tests were analyzed using a combination of standard pump-test analytical solution methods and methods developed expressly for the observed conditions. The stepped rate injection tests, combined with long-term flow and pressure response data suggest that the well is located within a fractured formation of low transmissivity but high storativity. These calculated parameters appeared to increase with pressure during the first injection test and the higher values were reproduced during the second stepped rate test. Calculated transmissivity and storativity are on the order of 4E-5 m cm 2 and 1E-4 m Pa -1 , respectively. The apparent pressure dependence of fitted hydraulic parameters may reflect near-well fracture compliance that increased the effective radius of the wellbore during the first test. While the type curve fit analysis also suggests that the reservoir behaves as a uniformly fractured reservoir with a radial flow regime, the hydraulic parameters indicate that condition may exist only a very limited distance (<10 m) from the well. Longer-term pressure response suggests that flow in the system effectively reaches steady state in a period of less than a day, which may reflect pressure stabilization resulting from pressure-dependent permeability or a region of much higher permeability located with a few meters of the well. The transmissivity estimates obtained from this analysis, converted to approximate fracture density and aperture, provide useful constraints on the distance to which the thermal front may migrate from the well during the cold water injection phase of the stimulation project. We estimate that the cooling front will migrate less than a tenth of a kilometer over an approximately one-year injection period. Here, the effects of that cooling, however, may be substantial, because increases in permeability have maximum effect nearest the well.

cold water injection↗

DOE OSTI2026

Solvent selection for a biomass-to-bioproduct pipeline through integrated reductive catalytic fractionation and microbial funneling

The growing significance of lignin-first biorefineries, which focus on upgrading the aromatics resulting from lignin depolymerization, presents opportunities for bioproduct synthesis using microbial strains capable of funneling a diverse array of phenolics into a single commodity chemical. In this study, we evaluated a biomass-to-bioproduct pipeline involving the reductive catalytic fractionation (RCF) of poplar biomass followed by biological funneling with a Novosphingobium aromaticivorans strain that produces 2-pyrone-4,6-dicarboxylic acid (PDC), a potential bioplastic precursor. Considering the impact of solvent on RCF reactor operating pressure, and the potential inhibitory effects of solvent on downstream microbial funneling, we performed an analysis of six pure solvents, namely methanol, ethanol, isopropanol, isobutanol, 1,4-dioxane and ethylene glycol, and different variations of their aqueous mixtures comprising 5 to 50 vol% water. For each pure solvent and solvent/water system, we measured phenolic monomer yields in the RCF process and PDC yields from the phenolic monomers. We then developed correlation models that relate phenolic monomer yields from RCF-derived samples to Hansen solubility parameters to determine solvent descriptors that contribute to high yields. Furthermore, we developed an integrated biorefinery system to estimate the minimum selling price (MSP) of PDC and the associated carbon footprint to identify solvent systems with better costs and sustainability metrics. These analyses resulted in the 50 vol% methanol/water system being identified as optimal because it reduces RCF reactor pressure and is compatible with microbial funneling with N. aromaticivorans. This solvent system produced 63 g PDC per kg biomass (264 g PDC per kg lignin) from 85 g phenolic monomers per kg biomass at a reduced reactor pressure of 48 bar (reduced by 26% compared to our previous poplar-to-PDC pipeline). The MSP for this system is $\$$13.98 per kg of purified PDC (carbon footprint of 1.47 kg CO 2 e per kg), which is about 24% lower than a previously described poplar-to-PDC pipeline and 46% lower than a lignin-to-PDC pipeline that used pure methanol as the solvent. The results from this study illustrate improvements that can be made in lignocellulosic biorefineries that are compatible with the hybrid chemical and biological processes needed to gain value from lignin.

Sripada, Sarada [Great Lakes Bioenergy Research Ce↗

DOE OSTI2025

In situ observation of irradiation-induced enhancement to the desorption pressure of zirconium hydride in a nuclear reactor

We quantify the effect of a nuclear-reactor environment on the hydrogen isotope equilibrium vapor pressure over pure zirconium and zirconium hydride. A vacuum-sealed capsule containing a zirconium foil with 6 atom% deuterium was irradiated at a neutron flux of ~10 14 cm -2 s -1 at the University of Missouri Research Reactor (MURR). The internal stainless-steel (SS) sample holder acted as the heat source via gamma absorption. To measure low desorption pressures in a high-flux environment, we developed a method to transduce pressure from the measured sample temperature during irradiation, calibrating with known deuterium pressures in unirradiated capsules at various heating powers using an internal filament-heated system designed to mimic irradiation-induced heating. Our temperature-pressure transduction method operates similarly to a Pirani or thermocouple pressure gauge. The in-reactor measurements revealed a roughly 4-fold enhancement in desorption pressure after only 6 h of irradiation (~2 × 10 18 cm -2 neutron fluence) compared to thermal desorption in control experiments, indicating a nonthermal contribution from neutron irradiation. The slower temperature/pressure stabilization rate in the reactor suggests that desorption pressure enhancement increases with neutron fluence. Further, this enhancement signifies increased solubility of hydrogen isotopes in zirconium during irradiation. We propose that high-energy neutron collisions with hydrogen isotopes in hydrides lead to their decomposition at lower temperatures, supersaturating the surrounding αZr lattice and resulting in higher desorption pressure, which continues to rise as more hydrides dissolve with increasing neutron fluence.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗

DOE OSTI2023

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗

DOE OSTI2024

Viscosity Measurements at High Pressures: A Critical Appraisal of Corrections to Stokes' Law

Abstract Fluids and melts in planetary interiors significantly influence geodynamic processes from volcanism to global‐scale differentiation. The roles of these geofluids depend on their viscosities ( η ). Constraining geofluid η at relevant pressures and temperatures relies on laboratory‐based measurements and is most widely done using Stokes' Law viscometry with falling spheres. Yet small sample chambers required by high‐pressure experiments introduce significant drag on the spheres. Several correction schemes are available for Stokes' Law but there is no consensus on the best scheme(s) for high‐pressure experiments. We completed high‐pressure experiments to test the effects of (a) the relative size of the sphere diameter to the chamber diameter and (b) the top and bottom of the chamber, that is, the ends, on the sphere velocities. We examined the influence of current correction schemes on the estimated viscosity using Monte Carlo simulations. We also compared previous viscometry work on various geofluids in different experimental setups/geometries. We find the common schemes for Stokes' Law produce statistically distinct values of η . When inertia of the sphere is negligible, the most appropriate scheme may be the Faxén correction for the chamber walls. Correction for drag due to the chamber ends depends on the precision in the sinking distance and may be ineffective with decreasing sphere size. Combining the wall and end corrections may overcorrect η . We also suggest the uncertainty in η is best captured by the correction rather than propagated errors from experimental parameters. We develop an overlying view of Stokes' Law viscometry at high pressures.

Geochemistry & Geophysics↗

DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Exact solution of burst pressure for thick-walled pipes using the flow theory of plasticity

This paper develops an exact solution of burst pressure for defect-free, thick-walled pipes using the flow theory of plasticity in terms of the average shear stress yield criterion (or Zhu-Leis yield criterion). The pipe steel is assumed to obey the power-law strain hardening rule, and large plastic deformation is described by the finite strain theory. On this basis, internal pressure is obtained as a power series function of the effective strains on the inside and outside surfaces of the thick-walled pipe with use of the second Bernoulli numbers. At burst failure, the Zhu-Leis flow solution of burst pressure is determined as a power series solution, and the burst effective strains, the burst effective stresses, and the burst pressure are functions of the diameter ratio (D o / D i ), strain hardening exponent (n), and ultimate tensile strength (UTS). Similarly, the Tresca and von Mises flow solutions of burst pressure are also determined as a power series solution in terms of the Tresca and von Mises yield criteria. A general closed-from exact solution of burst pressure is then proposed for the three yield criteria, and the results showed that the proposed exact solution matches well with the power series solution for each yield criterion. Moreover, the von Mises flow solution is an upper bound prediction, the Tresca flow solution is a lower bound prediction, and the Zhu-Leis flow solution is an intermediate prediction that agrees well with the finite element analysis results of burst pressure for thick-walled pipes. Furthermore, two datasets of full-scale burst tests are then utilized to evaluate and validate the proposed flow solutions of burst pressure for both thin and thick-walled pipes.

42 ENGINEERING↗