Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydraulics”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Code validation of SAM using natural-circulation experimental data from the compact integral effects test (CIET) facility

The primary objective of this study is to validate the system analysis code, SAM, using experimental data from the Compact Integral Effects Test (CIET) experimental loop. SAM is a modern system analysis code being developed at Argonne National Laboratory for safety analysis of designs for advanced non-light water reactors (non-LWRs), such as sodium-cooled fast reactors, high-temperature gas-cooled reactors, and fluoride salt-cooled high-temperature reactors (FHRs). To support SAM code development for the wide range of non-LWR applications, it is of paramount importance to validate the code against experiments highly relevant to these reactor concepts. Additionally, the CIET facility, which was designed to reproduce the thermal-hydraulics response of FHRs under both forced- and natural-circulation conditions, has been identified and selected as one of the benchmark test facilities for SAM code validation. In this study, two sets of available CIET tests were selected for SAM code validation purposes, namely, power step change transient tests and steady-state natural-circulation tests. For all selected tests, SAM-predicted results show very good agreement with experimental data. The successful validation of SAM against these selected CIET experiments demonstrates that the computer code is well suited for thermal-hydraulics analysis of FHR designs.

42 ENGINEERING↗
DOE OSTI2022

An experimental investigation on the characteristics of heat pipes with annular type composite wick structure

Here, an experimental investigation has been conducted on the hydraulic characteristics of annular type wick structures for heat pipes. An experimental facility which can measure porosity, permeability, and effective pore radius of the wick structures in a vacuum condition is established. Nine different types of multi-layered (6 layers in total) composite screen meshes are characterized. Based on the measurement results, a wick structure composed of one layer of 100×100 mesh, three layers of 400×400 mesh, and two layers of 60×60 mesh is determined to have the highest permeability to effective pore radius ratio $(K/r_{eff})$ and was selected as a targeted sample. The wick-to-wall gap effect on the hydraulic characteristics of annular type wick structure is also investigated by measuring the permeability and effective pore radius of the sample wick structure with varying gap widths. The result shows that the permeability increases as the gap increases from 0 mm to 1.2 mm. After a peak at 1.2 mm, the permeability decreases as the gap increases and converges to the value of the case measured without a wall structure. The effective pore radius becomes smaller as the gap increases, making a peak at a gap distance of 1.2 mm. This result implies that there is an optimal point in gap distance which is determined to be 1.2 mm for the selected composite mesh structure. A capillary limitation correlation with a multiplying factor which explains the enhanced performance of the wick due to the gap is suggested. An annular wick type heat pipe is constructed and tested. The capillary limitation of the heat pipe showed good agreement with the suggested correlation.

42 ENGINEERING↗
DOE OSTI2022

Simulation of natural circulation cartridge loop experiments and application to molten salt reactors

This work uses the TRAC/RELAP Advanced Computational Engine (TRACE) thermal hydraulics (TH) code to model natural circulation cartridge loop experiments previously conducted at Oak Ridge National Laboratory (ORNL) using water and compares the simulated and experimental results. TRACE is also used to characterize natural circulation in the cartridge loop vehicle using FLiNaK as the working fluid. The experimental vehicle is a buoyancy-aided, annular cartridge loop, referred to as a thermosyphon, and is designed to aid in qualifying liquid–fueled and/or liquid–cooled irradiation experiments for the Versatile Test Reactor (VTR), which is currently being designed in the United States. Out-of-pile water experiments have been conducted using the cartridge and the Thermosyphon Test Loop facility at ORNL, and future experiments are anticipated that would use other molten salt surrogates as the working fluid, followed by eventual insertion of a similar cartridge into VTR. Additionally, this work aims to determine how well TRACE can replicate the natural convection conditions that were observed experimentally; this serves as an initial step for validating the modeling tool for design and safety calculations to support future irradiation experiments in VTR. Initial predictions of potential experiments were made using FLiNaK as the natural circulation fluid to demonstrate the relevance of the cartridge design to molten salt reactors (MSRs). Results from this study indicate that TRACE can accurately capture natural convection phenomena in the thermosyphon and that several design changes to the current cartridge vehicle are necessary to achieve hydraulic conditions similar to those expected in MSRs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Coregistered positron emission particle tracking (PEPT) and X-ray computed tomography (CT) for engineering flow measurements

Increasingly, fully 3D experimental measurements of flow in complex engineering geometries are required to validate computational fluid dynamics models that support and inform reactor design and licensing. One barrier to such measurements is the complexity of typical reactor components and subsequent lack of optical access in these systems. To overcome this, the deployment of coregistered positron emission particle tracking (PEPT) and X-ray computed tomography (CT) is explored for flow measurement in reactor thermal hydraulic components and model (scaled) systems. Through this methodology, fully 3D flow information (via PEPT) and detailed internal geometry (via CT) are captured in opaque systems such as pipes, rod bundles, packed beds, etc. The reconstructed flow field and geometry can then be overlain to reveal detailed flow features around internal structures within a given test section. This is enabled through the use of a combined preclinical PET/CT scanner with overlapping PET and CT fields of view. Such measurements are useful for characterizing flow inside such intricate nuclear thermal hydraulic components as core geometries and heat exchangers, among others, and providing valuable 3D validation data for CFD models. In this work, basic tests of this 3D flow/geometry mapping are presented, and the implications of such measurements are discussed. Further, preliminary measurements are made with both point sources and flow in a simple pipe flow geometry to evaluate the capabilities of this technique. PEPT and CT features are coregistered with up to 0.1 mm precision, and pipe flow mean velocity and Reynolds stresses are reconstructed with similar accuracy to previous PEPT demonstrations. The utility of PEPT/CT is shown herein, and suggestions for future measurements are made.

3D flow measurement↗
DOE OSTI2025

A hybrid surrogate modeling framework for the Digital Twin of a Fluoride-salt-cooled High-temperature Reactor (FHR)

While nuclear energy is a non-greenhouse-gas emitting energy source, expensive operational costs due to the high-level of safety requirements decreases their competitiveness in the sustainable energy market. Advanced reactor concepts paired with Digital Twins aim to increase the commercialization gains of nuclear energy by reducing operational costs, increasing reactor reliability and enhancing power generation. To support Digital Twin tasks such as real-time autonomous control, proactive maintenance monitoring or optimizing power demand operations, a fast and accurate virtual representation of the Nuclear Power Plant (NPP) is required. The computational cost of high-fidelity, physics-based models are unsuitable for real-time analysis or scalability. Here, in this work, a hybrid surrogate modeling framework is developed fora Fluoride-salt-cooled High-temperature Reactor (FHR) that leverages physics-inspired models for key reactor components and uses data-driven methods for rapid system state space prediction. The Xenon reactivity feedback model is integrated to inform the surrogate model about the reactor core and the homologous pump theory model is the basis for representing pump degradation. Using a detailed, two dimensional thermal hydraulics model to generate data on the FHR, we train a network of Vectorized Autoregressive Moving-Average with eXogenous input (VARMAX) models to predict the remaining state values. The result is a surrogate model that provides a detailed reactor state representation of 41 system states and a pump degradation analysis. The framework is applied to Load Follows profiles, yielding high accuracy and a speedup that is more than 4000x faster compared to the higher- fidelity thermal hydraulics model, enabling real-time operational intelligence and applications in long horizon predictions. While the surrogate model framework is demonstrated for the particular case of FHR, the hybrid physical/data-driven modeling approach including the network of surrogates and the underlying modularity has the potential to be applied to other physical asset systems.

Digital Twins↗
DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗
DOE OSTI2021

Surfactants are Ineffective for Reducing Imbibition of Water-Based Fracturing Fluids in Deep Gas Reservoirs

Minimizing loss of injected hydraulic fracturing fluids into shale along fracture-matrix boundaries is desired because imbibed water restricts gas production and wastes valuable water resources. This problem has motivated the addition of surfactants into water-based hydraulic fracturing fluids in order to reduce the capillary driving force for imbibition. Here, we show that reduction in interfacial tension and wettability alteration has negligible ability to reduce imbibition in deep gas reservoirs. The effectiveness of altering capillary forces acting at the wetting front also depends on the injection pressure acting at the fracture-matrix boundary. The pressure at the interface between the fracture and the shale matrix is constrained between the reservoir pore pressure and formation pressure (rock fracture pressure, also known as breakdown pressure of the rock) and increases with depth to magnitudes that greatly exceed that of capillary pressures. Furthermore, the analyses presented here show that even maximum alteration of interfacial properties that result in strongly hydrophobic interactions between the fracturing fluid and reservoir rock is incapable of significantly reducing imbibition in deep reservoirs. Instead of using surfactants, this analysis points to decreases in wellbore shut-in pressures and shut-in times as practical options for reducing imbibition losses of water-based fluids.

03 NATURAL GAS↗
DOE OSTI2025

Inorganic Alterations in Unconventional Shale Reservoirs: Importance of Additive and Base Fluid Chemistry

The effective development of unconventional petroleum systems requires the use of significant water resources. In an effort to reduce the consumption of freshwater resources for hydraulic fracturing, highly saline produced waters are increasingly recycled for use as a base fluid. However, there are significant knowledge gaps regarding potential water–rock interactions resulting from the introduction of produced waters and associated additives into shale reservoirs such as formation and deposition of mineral scale, which can negatively affect hydrocarbon production through wellbore restriction and damage to hydraulically generated fractures. To assess the impacts of field stimulation practices in the subsurface, a series of laboratory experiments were completed using (a) three distinct sedimentary rock formations of the Midland Basin (Texas, USA) and (b) additives with two different base fluids: municipal fresh water and clean brine. The experimental approach used relevant injection sequences and mixing ratios in specialized reactors for 3 weeks. Static pressurized experiments and nonpressurized time-resolved experiments were undertaken. The resulting solids and liquids were analyzed by using a variety of laboratory- and synchrotron-based techniques. The use of an acid spearhead (15% HCl) resulted in texturing of both clay-rich and calcareous shales, which can temporarily enhance porosity but subsequently result in mineral scale deposition. The primary matrix scale was Fe(III)-bearing phases, which occurred in all experiments regardless of base fluid chemistry. Additionally, strontium sulfate (SrSO 4 ) precipitated on shale surfaces when clean brines were used. It was concluded that clean brine was the main source of Sr 2+ species, while persulfate breaker degradation and oxidation of pyrite were the sources of SO 4 2– . Sulfate scaling was more pronounced in clay-rich shales, suggesting that Sr sorption is important for promoting celestite formation. This work demonstrates that mineral scale deposition is a complex phenomenon, whereby the type and proportions of various mineral phases are determined from reservoir alteration processes and coprecipitation of constituents from injection fluids. In conclusion, the experimental results shown here should be considered when evaluating different base fluids and additives in order to mitigate mineral precipitation in unconventional shale reservoirs, which could result in reservoir degradation.

Jew, Adam D. [SLAC National Accelerator Laboratory↗
DOE OSTI2023

Structure–Transport Relationships of Water–Organic Solvent Co-transport in Carbon Molecular Sieve (CMS) Membranes

We explore the effects of the carbon molecular sieve (CMS) microstructure on the separation performance and transport mechanism of water–organic mixtures. Specifically, we utilize PIM-1 dense films and integrally skinned asymmetric hollow fiber membranes as polymer precursors for the CMS materials. The PIM-1 membranes were pyrolyzed under several different pyrolysis atmospheres (argon, carbon dioxide, and diluted hydrogen gas) and at multiple pyrolysis temperatures. Detailed gas physisorption measurements reveal that membranes pyrolyzed under 4% H 2 and CO 2 had broadened ultramicropore distributions (pore diameter <7 Å) compared to Ar pyrolysis, and pyrolysis under CO 2 increased ultramicropore volume and broadened micropore distributions at increased pyrolysis temperatures. Gravimetric water and p-xylene sorption and diffusion measurements reveal that the PIM-1-derived CMS materials are more hydrophilic than other CMS materials that have been previously studied, which leads to sorption-diffusion estimations showing water-selective permeation. Water permeation in the vapor phase, pervaporation, and liquid-phase hydraulic permeation reveal that the isobaric permeation modes (vapor permeation and pervaporation) are reasonably well predicted by the sorption-diffusion model, whereas the hydraulic permeation mode is significantly underpredicted (>250×). Conversely, the permeation of p-xylene is well predicted by the sorption-diffusion model in all cases. The collection of pore size analysis, vapor sorption and diffusion, and permeation in different modalities creates a picture of a combined transport mechanism in which water–under high transmembrane pressures–permeates via a Poiseuille-style mechanism, whereas p-xylene solutes in the mixture permeate via sorption-diffusion.

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

Hybrid Sorbents for 129 I Capture from Contaminated Groundwater

Radioiodine ( 129 I) poses a risk to the environment due to its long half-life, toxicity, and mobility. It is found at the U.S. Department of Energy Hanford Site due to legacy releases of nuclear wastes to the subsurface where 129 I is predominantly present as iodate (IO 3 – ). To date, a cost-effective and scalable cleanup technology for 129 I has not been identified, with hydraulic containment implemented as the remedial approach. Here, novel high-performing sorbents for 129 I remediation with the capacity to reduce 129 I concentrations to or below the US Environmental Protection Agency (EPA) drinking water standard and procedures to deploy them in an ex-situ pump and treat (P&T) system are introduced. Furthermore, this includes implementation of hybridized polyacrylonitrile (PAN) beads for ex-situ remediation of IO 3 – -contaminated groundwater for the first time. Iron (Fe) oxyhydroxide and bismuth (Bi) oxyhydroxide sorbents were deployed on silica substrates or encapsulated in porous PAN beads. In addition, Fe–, cerium (Ce)–, and Bi–oxyhydroxides were encapsulated with anion-exchange resins. The PAN–bismuth oxyhydroxide and PAN–ferrihydrite composites along with Fe- and Ce-based hybrid anion-exchange resins performed well in batch sorption experiments with distribution coefficients for IO 3 – of >1000 mL/g and rapid removal kinetics. Of the tested materials, the Ce-based hybrid anion-exchange resin was the most efficient for removal of IO 3 – from Hanford groundwater in a column system, with 50% breakthrough occurring at 324 pore volumes. The functional amine groups on the parent resin and amount of active sorbent in the resin can be customized to improve the iodine loading capacity. These results highlight the potential for IO 3 – remediation by hybrid sorbents and represent a benchmark for the implementation of commercially available materials to meet EPA standards for cleanup of 129 I in a large-scale P&T system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2019

The Community Land Model Version 5: Description of New Features, Benchmarking, and Impact of Forcing Uncertainty

The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.

58 GEOSCIENCES↗
DOE OSTI2020

Changes in Physical Properties of Everglades Peat Soils Induced by Increased Salinity at the Laboratory Scale: Implications for Changes in Biogenic Gas Dynamics

Abstract Saltwater intrusion due to sea level rise is a major concern for the Florida Everglades because it may induce shifts in ecosystem productivity and physical soil properties. However, the effects of saline water intrusion into the current carbon gas dynamics of the Everglades (particularly in terms of biogenic gas production and emissions, i.e., CH 4 and CO 2 ) are still uncertain. In this work, we present a laboratory‐based study to simulate how sea level rise may alter the physical properties (i.e., hydraulic conductivity) of peat soils from the Everglades and consequently affect the accumulation and release of biogenic gases within the peat matrix. Peat monoliths collected from the Everglades were subjected to progressive increases in salinity from a NaCl solution, and changes to the biogenic gas dynamics regime were simultaneously monitored using a combination of time‐lapse ground‐penetrating radar measurements, manometers, time‐lapse photography, and gas traps. Physical changes to the peat matrix at each salinity interval were assessed using constant head permeameter tests. Consistent with previous research, results show that a progressive increase in salinity (from fresh to saltwater) results in (1) a progressive increase in hydraulic conductivity and (2) a progressive decrease in gas content within the peat matrix (i.e., production) and gas releases. This work has implications for better understanding the potential effects of saltwater intrusion into freshwater peatland systems in the Everglades, particularly in terms of carbon gas dynamics.

Sirianni, Matthew J.↗
DOE OSTI2021

Effect of Glacial/Interglacial Recharge Conditions on Flow of Meteoric Water Through Deep Orogenic Faults: Insights Into the Geothermal System at Grimsel Pass, Switzerland

Many meteoric-recharged, fault-hosted geothermal systems in amagmatic orogenic belts have been active through the Pleistocene glacial/interglacial climate fluctuations. The effects of such climate-induced recharge variations on fluid flow patterns and residence times of the thermal waters are complex and may influence how the geothermal and mineralization potential of the systems are evaluated. We report systematic thermal–hydraulic simulations designed to reveal the effects of recharge variations, using a model patterned on the orogenic geothermal system at Grimsel Pass in the Swiss Alps. Fault-bound circulation of meteoric water through the Grimsel fault is driven to depths of ~10 km by the high alpine topography. Simulations suggest that the current single pass flow is typical of interglacial periods, during which 1) meteoric recharge into the fault is high (above tens of centimeters per year), 2) conditions are at or somewhat below the critical Rayleigh number, and 3) the hydraulic connectivity along the fault plane is extensive (an extent of at least 10 km into increasingly higher terrain is required to explain the 10 km penetration depth). The subcritical condition constrains the bulk fault permeability to < 1e-14 m2. In contrast, the limited recharge during the numerous Pleistocene glaciation events likely induced a layered flow system with single pass flow confined to shallow depths while non-Rayleigh convection occurred in the deep fault. The same layering can be observed at low aspect ratios (length/depth) of the fault plane, when the available recharge area limits the flux through the fault.

58 GEOSCIENCES↗
DOE OSTI2021

Regional Subsidence Analysis Through a Multi-Scale Modeling Framework Based on Breakage Mechanics

Although poroelastic models are often used to explain the delay between subsurface fluid depletion and ground subsidence, inelastic compaction involving permanent changes of rock microstructure may exacerbate hydro-mechanical coupling, thus influencing the interpretation of measurements and long-term forecasts. Here, a multi-scale modeling approach is discussed, which accounts for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. A constitutive model built within the framework of breakage mechanics is proposed to link the hydraulic conductivity of granular rocks with inelastic deformations and changes in grading caused by injection-depletion cycles at stress levels far from yielding. The proposed model has been incorporated into large-scale simulation frameworks, thus enabling the spatiotemporal mapping of regional subsidence through a hybrid, semi-analytical approach. Numerical results based on this strategy show that the model allows isolating near-field and far-field effects into the computation of land subsidence and can generate forecasts for different modeling scenarios (e.g., elastic and inelastic compaction, constant permeability, and concurrent change of compressibility and permeability). In particular, examples of simulations for the case of the Groningen gas field are discussed, showing the model capabilities to use both field measurements and laboratory tests for the generation of reasonable subsidence maps, without expensive computational costs. Results indicate that ignoring coupled inelastic effects has major consequences on the predicted timescale of subsidence. Specifically, while all the model scenarios produced similar long-term ground settlements, those ignoring breakage-dependent permeability changes result in a variation of the temporal window of residual subsidence of the order of several decades.

58 GEOSCIENCES↗
DOE OSTI2022

Physics‐Informed Machine Learning Method for Large‐Scale Data Assimilation Problems

Abstract We develop a physics‐informed machine learning approach for large‐scale data assimilation and parameter estimation and apply it for estimating transmissivity and hydraulic head in the two‐dimensional steady‐state subsurface flow model of the Hanford Site given synthetic measurements of said variables. In our approach, we extend the physics‐informed conditional Karhunen‐Loéve expansion (PICKLE) method to modeling subsurface flow with unknown flux (Neumann) and varying head (time‐dependent Dirichlet) boundary conditions. We demonstrate that the PICKLE method is comparable in accuracy with the standard maximum a posteriori (MAP) method, but is significantly faster than MAP for large‐scale problems. Both methods use a mesh to discretize the computational domain. In MAP, the parameters and states are discretized on the mesh; therefore, the size of the MAP parameter estimation problem directly depends on the mesh size. In PICKLE, the mesh is used to evaluate the residuals of the governing equation, while the parameters and states are approximated by the truncated conditional Karhunen‐Loéve expansions with the number of parameters controlled by the smoothness of the parameter and state fields, and not by the mesh size. For a considered example, we demonstrate that the computational cost of PICKLE increases near linearly (as N 1.15 ) with the number of grid nodes N , while that of MAP increases much faster (as N 3.28 ). We also show that once trained for one set of Dirichlet boundary conditions (i.e., one river stage), the PICKLE method provides accurate estimates of the hydraulic head for any value of the Dirichlet boundary conditions (i.e., for any river stage).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Theory of an Automatic Seepage Meter and Ramifications for Applications

A new approach for measuring fluxes across surface water—groundwater interfaces was recently proposed. The Automatic Seepage Meter (ASM) is equipped with a precise water level sensor and digital memory that analyzes water level time series in a vertical tube inserted into a streambed. The ability to infer flux values with high temporal resolution relies on an accurate interpretation of water level dynamics inside the tube. Here, we reduce the three-dimensional hydrodynamic problem that describes the ASM water level in a variety of field conditions to a single ordinary differential equation. This novel general analytical solution for estimating ASM responses is more comprehensive and flexible than previous approaches and is applicable to the entire range of field conditions, including steady or transient stream stages, evaporation, rainfall, and noise. For example, our analysis determines the timing of the nonmonotonic ASM response to a monotonic linear stream stage variation and explains previously used empirical parabolic approximation for estimating fluxes. We present algorithms for simultaneous inference of vertical interface flux and hydraulic conductivity values together with an example code. We quantify how the accuracy of parameter estimation depends on test duration and noise amplitude and propose how our analysis can be used to optimize field test protocols. On this basis, changing the ASM geometry by increasing the radius and decreasing tube insertion depth may enable ASM field test protocols that estimate interface flux and hydraulic conductivity faster while maintaining desired accuracy. Potential applications of joint parameter estimation are suggested.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Quartz Dissolution Effects on Flow Channelization and Transport Behavior in Three‐Dimensional Fracture Networks

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Hyman, Jeffrey D.↗