SEARCH · Engineering Papers

Results for “flow”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2020

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Effects of Intrusive Confining Units on Groundwater Flow at Pahute Mesa, Nevada National Security Site

Geological, geophysical, hydrogeological, and hydrological evidence indicates that feeder dikes to rhyolitic lava flows may play a role as confining units in the groundwater flow system underlying the Pahute Mesa (PM) area of underground nuclear testing at the Nevada National Security Site (NNSS) (Figure 1). Of geological significance, the PM testing area is largely situated within the Silent Canyon caldera complex (SCCC), which includes both tuffs and lavas that were fed by (and are connected to) intrusive rock units such as feeder dikes prior to and during volcanic eruptions (Orkild et al., 1968; Byers et al., 1976a). Groundwater modeling results presented in this report support extension of “pool-dam” and “breach scenario” concepts presented in Jackson and Fenelon (2018) and Jackson et al. (2021) to low-permeability intrusive confining units (ICUs) such as feeder dikes that are expected to dissect and ring silicic calderas. As laterally extensive, steeply-dipping, and sub-parallel planar structures, ICUs may further compartmentalize groundwater flow in some areas of the volcanic aquifer system at PM. ICUs may explain unusually high horizontal hydraulic gradient southwest of Cheshire. Breaching of an ICU within which the Cheshire test was emplaced may explain changes in water levels before and after the Cheshire test was detonated in the U-20n emplacement hole on 2/14/76 (USDOE, 2015b). This analysis of the effects of ICUs on groundwater flow at PM addresses a need to evaluate uncertainty in the conceptual model or structural aspects of hydrologic source term (HST) or corrective action unit (CAU) models for groundwater flow and contaminant transport (FFACO, 1996, as amended). At Cheshire and possibly other test locations at PM, downgradient extent of radionuclide contamination may be overpredicted by not including the confining effects of ICUs. On the other hand, contaminant extent may be underpredicted by not considering transient flow effects and/or the re-equilibrated flow field caused by breaching of an ICU by an underground nuclear test. Near Cheshire, the effects of ICUs may be most apparent and impactful to the groundwater flow system at PM because of an unusual combination of thickness, extent, and frequency of dikes transecting the saturated CHLFA4 aquifer and feeding into the overlying unsaturated rhyolite lava of Windy Wash. Chapters 2, 3, and 4 present geological, geophysical, hydrogeological, and hydrological evidence for ICUs at PM. Chapters 5 and 6 examine effects of ICUs on groundwater flow including analysis by groundwater modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Hydrogen, Methane, Brine Flow Behavior, and Saturation in Sandstone Cores During H 2 and CH 4 Injection and Displacement

Large-scale underground hydrogen storage (UHS) is a critical component in the emerging hydrogen economy. Knowledge of multiphase flow behavior involving hydrogen in storage reservoir formations is crucial to characterizing hydrogen transport properties and essential for the deliverability and storage operations of UHS. There are still many gaps in fully understanding hydrogen–methane–brine multiphase phase flow that require further investigation. In this work, H 2 and CH 4 were injected through brine-saturated sandstone cores using a tri-axial core holder system fitted with flow rate meters and pressure transducers, while the effluent gas concentrations were analyzed using an online micro gas chromatograph. Brine displacement, permeability, and gas breakthrough curves were measured. We studied the flow behavior of hydrogen and methane in sandstone cores through testing brine displacement by gas injection and comparing the hydrogen displacement of methane with the methane displacement of hydrogen. We also tested the differences between horizontal and vertical flow in brine displacement. The results showed that brine displacement was more efficient in a core with higher permeability and porosity, resulting in a higher initial gas saturation. A higher gas injection rate brought about faster gas breakthrough measured by pore volume and sharper concentration curves. Hydrogen did not exhibit abnormal flow in the sandstone when the flow was horizontal and downward vertical. Gas overriding was observed in brine displacements when the flow was horizontal, with hydrogen showing this behavior more profoundly compared to methane. Downward vertical gas injection induced higher efficiency brine displacement compared to horizontal displacement and resulted in a higher initial gas saturation in the sandstone cores. These findings address critical knowledge gaps regarding gas flow patterns and displacement behaviors during hydrogen injection and recovery phases in UHS facilities using methane as the cushion gas. The insights from this research offer valuable guidance for optimizing UHS systems, ensuring operational efficiency, and advancing sustainable energy solutions in alignment with decarbonization goals.

08 HYDROGEN↗

DOE OSTI2021

X-ray Imaging-Based Void Fraction Measurement in Saturated Flow Boiling Experiments with Seawater Coolant

Recent studies have shown that the presence of dissolved salts in water can exhibit peculiar flow boiling and two-phase flow regimes. Two-phase flow and convective flow boiling are typically characterized with the help of void fraction measurements. To quantitatively improve our understanding of two-phase flow and boiling phenomenon with seawater coolant, void fraction data are needed, which can not be obtained from optical imaging. In this paper, we present experimental void fraction measurements of saturated flow boiling of tap water and seawater using X-ray radiography. X-rays with a maximum energy level of 40 KeV were used for imaging the exit region of the heated test section. At lower heat flux levels, the two phase flow in seawater was bubbly and homogeneous in nature, resulting in higher void fractions as compared to tap water. With an increase in heat flux, the flow regime was similar to slug flow, and void fraction measurements approached similarity with tap water. The predicted pressure drop using the measured void faction shows good agreement with the measured total pressure drop across the test section, demonstrating the validity of the measurement process.

Eckels, Seth↗

DOE OSTI2020

Investigation of the Elliptic Flow Fluctuations of the Identified Particles Using the a Multi-Phase Transport Model

A Multi-Phase Transport (AMPT) model is used to study the elliptic flow fluctuations of identified particles using participant and spectator event planes. The elliptic flow measured using the first order spectator event plane is expected to give the elliptic flow relative to the true reaction plane which suppresses the flow fluctuations. However, the elliptic flow measured using the second-order participant plane is expected to capture the elliptic flow fluctuations. Our study shows that the first order spectator event plane could be used to study the elliptic flow fluctuations of the identified particles in the AMPT model. The elliptic flow fluctuations magnitude shows weak particle species dependence and transverse momentum dependence. Such observation will have important implications for understanding the source of the elliptic flow fluctuations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Anisotropic Turbulent Flows Observed in Above-the-loop-top Regions during Solar Flares

Abstract Solar flare above-the-loop-top (ALT) regions are vital for understanding solar eruptions and fundamental processes in plasma physics. Recent advances in three-dimensional (3D) magnetohydrodynamic (MHD) simulations have revealed unprecedented details on turbulent flows and MHD instabilities in flare ALT regions. Here, for the first time, we examine the observable anisotropic properties of turbulent flows in ALT by applying a flow-tracking algorithm on narrow-band extreme-ultraviolet images that are observed from the face-on viewing perspective. First, the results quantitatively confirm the previous observation that vertical motions dominate and that the anisotropic flows are widely distributed in the entire ALT region with the contribution from both upflows and downflows. Second, the anisotropy shows height-dependent features, with the most substantial anisotropy appearing at a certain middle height in ALT, which agrees well with the MHD modeling results where turbulent flows are caused by Rayleigh–Taylor-type instabilities in the ALT region. Finally, our finding suggests that supra-arcade downflows (SADs), the most prominently visible dynamical structures in ALT regions, are only one aspect of turbulent flows. Among these turbulent flows, we also report the antisunward-moving underdense flows that might develop due to MHD instabilities, as suggested by previous 3D flare models. Our results indicate that the entire flare fan displays group behavior of turbulent flows where the observational bright spikes and relatively dark SADs exhibit similar anisotropic characteristics.

Xie, Xiaoyan (ORCID:0009000705827807)↗

DOE OSTI2020

Catalytic converters having non-linear flow channels

Disclosed is a honeycomb catalyst substrate core having geometrically non-linear flow channels. In an embodiment, the honeycomb catalyst substrate core includes helical flow channels. In another embodiment, the honeycomb catalyst substrate core includes sinusoidal flow channels. In yet another embodiment, the honeycomb catalyst substrate core includes helical plus sinusoidal flow channels. The honeycomb catalyst substrate core comprises a plurality of parallel non-linear flow channels formed along a longitudinal axis of symmetry of the catalyst substrate core, each non-linear flow channel configured such that eddies occurs during engine exhaust gas flow. Also disclosed is a method for manufacturing a ceramic honeycomb having non-linear flow channels, comprising the steps extrusion soft ceramic material through a die whilst the die moves through six degrees of freedom along its axis of symmetry. Disclosure includes a method for manufacturing a ceramic honeycomb having non-linear flow channels using three-dimensional printing.

Masoudi, Mansour↗

DOE OSTI2023

And manufacture of generalized flow profile-producing devices

The present invention provides a process for making a flow conditioning device that transforms an input flow into a desired output flow. The process includes the steps of inputting into a computer program a set of design constraints representative of the input flow and the output flow. The computer program generates a design representative of a flow-conditioning device that transforms the input flow into the output flow. The process then provides the output design to an additive manufacturing or other suitable production system adapted to form a solid representation of the flow-conditioning device.

Ferrar, Anthony M.↗

DOE OSTI2024

Review: Pre-Darcy flows in low-permeability porous media

The widely used Darcy’s law specifies a linear relation between the Darcy velocity of fluid flow and the pressure gradient that drives the flow. However, studies have shown that Darcy velocity can exhibit a nonlinear dependence on the pressure gradient in low-permeability porous media such as clay and shale when the pressure gradient is adequately low. This phenomenon is referred to as low-velocity non-Darcian flow or pre-Darcy flow. This paper provides a comprehensive review of the theories, experimental data, and modeling methods for pre-Darcy flow in low-permeability porous media. The review begins by outlining the fundamental mechanisms underlying pre-Darcy flow that regulate the unique characteristics such as nonlinear dependence of the Darcy velocity on the pressure gradient and its relevance to fluid–rock interactions. The review then proceeds to present a thorough compilation of experimental investigations performed in various low-permeability geomaterials including tight sandstones, shales, and clays. Next, empirical and theoretical models and simulation methods that have been developed to fit and interpret experimental data are reviewed. Finally, the review underscores the challenges encountered in conducting and interpreting pre-Darcy flow experiments and suggests future research directions. By analyzing previous experimental investigations, this review aims to offer a valuable resource for researchers and practitioners seeking to enhance their understanding of fluid dynamics in low-permeability geomaterials. This provides insights into the application of pre-Darcy flow in numerous natural and engineered processes such as shale oil and gas recovery, contaminant transport in low-permeability aquifers, and geological disposal of nuclear waste.

58 GEOSCIENCES↗

DOE OSTI2020

Experimental investigation of turbulent wake flows in a helically wrapped rod bundle in presence of localized blockages

In nuclear sodium fast reactors, bundles of rods are tightly packed into a triangular lattice, enclosed in a hexagonal duct, and each pin is spirally wrapped with a thin wire. Flow blockages can potentially impact the local flow characteristics and heat transfer mechanisms in the bundle due to its small subchannel size. The effects of the blockage on the flow structures and heat transfer mechanisms are important aspects that require an accurate investigation. In this study, the flow-field characteristics in the vicinity of a blockage located in the exterior subchannel of rod bundles with helically wrapped wires were experimentally investigated. The velocity fields in the exterior subchannel were acquired by applying matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) techniques for Reynolds numbers of Re 1 = 4000 and Re 2 = 17 000, i.e., equivalent to Rew 1 = 19 600 and Rew 2 = 83 200, respectively, based on the blockage width. The results from the TR-PIV measurements revealed an arch-shaped vortex with a large flow recirculation and a pair of counter-rotating vortices in the wake region downstream of the blockage, which is commonly observed in the wake flow of bluff bodies. The relative lateral distance and angle between the two vortices decreased when the Reynolds numbers increased. Profiles of maximum turbulence intensity along the shear layers illustrated the transition process including the growth, peak, and decay along the flow direction. From the spectral analysis of the turbulent velocities extracted at points along the shear layer, the Strouhal numbers ( St ) representing the vortex shedding frequency were found to be St = 0.25 and St = 0.56 for the left and right shear layers, respectively. Characteristics of shear layers generated by the blockage in the exterior subchannel were investigated via the two-point cross correlation of fluctuating velocities. The patiotemporal cross correlations of turbulent velocities, computed at points in the region where the left shear layer exhibited rolling effects and vortex breakdowns, were considerably wider and longer. The convection velocity U c was estimated to be ~0.82 U m to 0.93 U m . Proper orthogonal decomposition (POD) analysis was applied to the instantaneous velocity fields to extract the statistically dominant flow structures. It was found that POD modes 2–3 and 4–5 formed the pair modes when the corresponding POD temporal coefficients depicted sinusoidal shapes and exhibited nearly circular orbits in the phase space. Spectral analysis of the POD temporal coefficients confirmed the vortex shedding frequencies detected in the analysis of turbulent velocities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

FreeMHD: Validation and verification of the open-source, multi-domain, multi-phase solver for electrically conductive flows

The extreme heat fluxes in the divertor region of tokamaks may require an alternative to solid plasma-facing components, for the extraction of heat and the protection of the surrounding walls. Flowing liquid metals are proposed as an alternative, but raise additional challenges that require investigation and numerical simulations. Free surface designs are desirable for plasma-facing components, but steady flow profiles and surface stability must be ensured to limit undesirable interactions with the plasma. Previous studies have mainly used steady-state, 2D, or simplified models for internal flows and have not been able to adequately model free-surface liquid metal (LM) experiments. Therefore, FreeMHD has been recently developed as an open-source magnetohydrodynamics (MHD) solver for free-surface electrically conductive flows subject to a strong external magnetic field. The FreeMHD solver computes incompressible free-surface flows with multi-region coupling for the investigation of MHD phenomena involving fluid and solid domains. The model utilizes the finite-volume OpenFOAM framework under the low magnetic Reynolds number approximation. FreeMHD is validated using analytical solutions for the velocity profiles of closed channel flows with various Hartmann numbers and wall conductance ratios. Next, experimental measurements are then used to verify FreeMHD, through a series of cases involving dam breaking, 3D magnetic fields, and free-surface LM flows. These results demonstrate that FreeMHD is a reliable tool for the design of LM systems under free surface conditions at the reactor scale. Furthermore, it is flexible, computationally inexpensive, and can be used to solve fully 3D transient MHD flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

General Correlations for Laminar Flow Friction Loss and Heat Transfer in Plain Rectangular Plate-Fin Cores

In this paper, new generalized correlations for predicting the average fanning friction factor f and average Nusselt number Nu for laminar flow in plain plate-fin compact cores of rectangular cross section are presented. These are based on extended experimental data, as well as three-dimensional computational simulations, obtained for a broad range of fin density and geometrical attributes. The results indicate that while the fully developed forced convection scales only with the interfin channel cross-sectional ratio α (fin spacing by fin height), the entrance region hydrodynamic and thermal performance is additionally a function of the fin-core length L, flow Reynolds number Re, and fluid Prandtl number Pr. The developing flow and convection is further shown to scale as: (fRe)~(L/d h Re) 1/2 , and Nu ~(L/d h Re) 1/2 Pr 1/3 Φ(α), where f, Re, and Nu are all based on the hydraulic diameter dh of the interfin flow channel. Generalized correlations for both (fRe) and Nu are developed by the corresponding scaling of the forced convection behavior and asymptotic matching of the entrance or developing flow (short fin-core flow length) and the fully developed flow (large fin-core flow length) region performance. Finally, the predictions from these correlations are found to be within ±15% of all available experimental data for air, water, and glycol (0.71 ≤ Pr ≤ 10), and fin cores with 0 < α ≤ 1.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Low flow characteristics for regulated and unregulated streams in North Carolina and prediction using climate signals

In this work, low flow statistics of regulated and unregulated streams in the state of North Carolina were updated with streamflow data through 2019. About 22% of the streams considered show a significant downward trend, but considerable low-frequency variability confounding trends. The relationship between regional index time series of groundwater-depth low flow and the Atlantic Multidecadal Oscillation (AMO) was examined to determine whether low-frequency climate modes can account for the long-term pattern in low flows. Consequently, a significant correlation was found between AMO and groundwater-depth low flow, such that positive AMO is associated with lower groundwater-depth low flow and vice versa, particularly for the Piedmont region. Predictive equations for annual low flows at the ecoregion level shows that springtime average streamflow and AMO were selected as the primary predictors of low flow for coastal and Piedmont regions, whereas springtime average streamflow and the November–December–January average Oceanic Niño Index were used as the primary predictors for the mountain region. The relative root mean square error (RMSE) of the disaggregated predictions to US Geologic Survey gauge locations was <23% at 79% of the stations, between 24% and 43% at 10% of the stations, and greater than 44% at 1% of the stations. The remaining 10% of stations showed large RMSEs. This latter percentage is characterized by smaller drainage basins and intermittent flows, suggesting the prediction models are not applicable to drainage basins smaller than roughly 20 square km and intermittent streams.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗

DOE OSTI2024

Continuous Flow Photoelectrochemical Reactor with Gas Permeable Photocathode: Enhanced Photocurrent and Partial Current Density for CO 2 Reduction

Photoelectrochemical (PEC) CO 2 reduction using a photocathode is an attractive method for making valuable chemical products due to its simplicity and lower overpotential requirements. However, previous PEC processes have often been diffusion-limited leading to low production rates of the CO 2 reduction reaction, due to inefficient gas diffusion through the liquid electrolyte to the catalyst surface, particularly at high current densities. In this study, a gas-permeable photocathode in a continuous flow PEC reactor is incorporated, which facilitates the direct supply of CO 2 gas to the photocathode-electrolyte interface, unlike dark reaction-based flow reactors. This concept is demonstrated using Ag-TiO 2 on carbon paper, illuminated through a quartz window and flowing liquid electrolyte. CO 2 supply is managed via pressure and flow control on the non-illuminated side of the carbon paper. The photocurrent density is significantly influenced by the flow rates and pressure of CO 2 gas, and the electrolyte flow rates. Compared to the traditional H-cell, the continuous PEC flow reactor achieves ≈10-fold increase in CO faradaic efficiency, 30-fold increase in production rate and 16-fold increase in stability without catalyst modifications. This work provides essential insights into the design and application of continuous gas-liquid flow PEC reactor systems, highlighting their potential for other PEC reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗

DOE OSTI2024

Impacts of irregularly-distributed acidified brine flow on geo-chemo-mechanical alteration in an artificial shale fracture under differential stress

The efficacy of geological carbon sequestration is reliant on the integrity of the caprock and its resistance to physical and chemical alteration. Caprocks with high abundance of reactive carbonates like calcite are susceptible to acid-promoted dissolution and can result in structural weakening. This work investigates the effect of acidified brine flow through an artificially fractured, high-carbonate (30 % by XRD) shale under differential compressive stress. Cylindrical samples were cut in half vertically and milled to create an artificial fracture with interlocking asperities and open channels. Samples were sheared with a single applied stress in a custom flow cell housed within an industrial CT scanner. Further, either acidic (pH 4) or reservoir-simulated (pH 9.5) brine was flowed through the artificial fracture for 7–8 days under reservoir pressure and room temperature. Model simulations indicate flow mainly occurred in open channels, with limited flow between overlapping asperities. Analysis of fracture surfaces by optical and scanning electron microscopy show increased surface alteration and roughness after exposure to pH 4 versus pH 9.5 brine indicating mineral dissolution/loss, and this effect is greater in areas that receive the highest brine flows. Similarly, surface analysis by scratch testing shows fracture toughness decreases more after exposure to acidic versus reservoir-simulated brine, with the greatest alteration in areas of highest acidic brine flows. Despite weakening, no shear slip occurred. Overall, the results indicate that acidified brine can result in significant physical and geomechanical alteration of irregular fracture surfaces in shale caprock, with greatest effects in preferential flow regions.

58 GEOSCIENCES↗

DOE OSTI2022

Dynamic coevolution of baseflow and multiscale groundwater flow system during prolonged droughts

Field and numerical studies suggest that baseflow is composed of waters from a spectrum of groundwater flow paths termed the Groundwater Flow System (GWFS) – from shallow hillslope contributions to watershed-scale deep circulation originating in headwaters and discharging into lowland rivers. Here, we explore the evolution of the GWFS under prolonged droughts to understand its dynamics and multiscale nature, and to elucidate its role in baseflow generation and recession at the watershed scale. In this work, we consider three drought scenarios of varying severity and simulate groundwater flow in a 2-D cross-section of an idealized watershed with deep permeable bedrock, tracking the evolution of flow paths, baseflow, and residence times during the recession process. We find that baseflow generation at different drainage stages, and within different subwatersheds, is influenced distinctly by flow paths of different scales, depending on the relative strength of the flow paths and the position of the subwatersheds relative to the recharge/discharge zones of the deeper watershed-scale groundwater circulation. Despite having the same local relief, geology, and climate, baseflow from each subwatershed has a distinct recession behavior and time-dependent residence time distribution. Also, the hydraulic and transport characteristics of baseflow generation co-evolve and are strongly affected by the connection state of the water table to subwatersheds. These findings suggest that asynchrony and dissimilarity of baseflow generation from hillslopes under the impact of the watershed-scale groundwater flow, and interactions with local-scale and intermediate-scale groundwater flow, must be taken into account when interpreting baseflow recession data and building conceptual baseflow models at the watershed scale.

54 ENVIRONMENTAL SCIENCES↗