Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydraulic characterization”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning

We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.

Geomechanics↗
DOE OSTI2023

Physics-embedded inverse analysis with algorithmic differentiation for the earth’s subsurface

Abstract Inverse analysis has been utilized to understand unknown underground geological properties by matching the observational data with simulators. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented with embedded physics and a technique known as algorithmic differentiation. We use a physics-embedded generative model, which takes statistically simple parameters as input and outputs subsurface properties (e.g., permeability or P-wave velocity), that embeds physical knowledge of the subsurface properties into inverse analysis and improves its performance. We tested the application of this approach on four geologic problems: two heterogeneous hydraulic conductivity fields, a hydraulic fracture network, and a seismic inversion for P-wave velocity. This physics-embedded inverse analysis approach consistently characterizes these geological problems accurately. Furthermore, the excellent performance in matching the observational data demonstrates the reliability of the proposed method. Moreover, the application of algorithmic differentiation makes this an easy and fast approach to inverse analysis when dealing with complicated geological structures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Discrete Measurements of the Least Horizontal Principal Stress from Core Data: An Application of Viscoelastic Stress Relaxation

Summary The emerging Paradox Oil Play in southeastern Utah is among the most significant unconventional plays in the western USA. The mean total undiscovered oil resources within just the Pennsylvanian Cane Creek interval of the Paradox Basin are believed to exceed 215 million barrels. However, to date, less than 5% (~9 million barrels) of the total Cane Creek resource has been produced from fewer than 40 wells, and only approximately one-half of those are horizontal wells. More than 95% of production is from the central Cane Creek Unit (CCU). Natural fractures are a key feature of many production wells, but stimulation by induced hydraulic fractures is not consistently successful. We hypothesize that more effective production in this play will rely on better fundamental characterization, especially on better quantification of the state of stress. Approximately 110 ft of core, well logs, and a diagnostic fracture injection test (DFIT) were acquired from the State 16-2 well within the CCU. With these data, we applied two methods to constrain and clarify the state of stress. The first technique, the Simpson’s coefficient method, provides lower bounds on the two horizontal principal stresses and relies on only limited data. Alternatively, the viscoelastic stress relaxation (VSR) method is used to estimate the least horizontal principal stress, building on observations that principal stresses become more isotropic as the viscous behavior of a rock is more pronounced. Results of these two methods support the hypothesis that the state of stress in the CCU of the Paradox Basin is nearly lithostatic and isotropic. Other factors consistent with this hypothesis include high formation pore pressure, which tends to reduce the possible stress states by changing the frictional failure equilibrium; lack of induced fractures in the core, which should be present in the case of stress anisotropy; and interbedded halite layers, which given their high degree of ductility, probably lead to greater VSR for the entire sedimentary package.

Energy & Fuels↗
DOE OSTI2021

Analysis of convection heat transfer on multiscale rough superhydrophobic and liquid infused surfaces

Multiscale rough superhydrophobic or slippery liquid infused porous surfaces have gained much interest in recent years for their improved transport phenomena properties. While there have been several studies on drag reduction and condensation on non-wetting surfaces, convection heat transfer that is important in many thermal and thermochemical applications has not been addressed systematically. Here, this article utilizes a fractal description of rough surface topographies to develop analytical models for the Nusselt number and the thermal hydraulic factor for fluid flow and heat transfer inside a cylinder with non-wetting surfaces. For air-infused superhydrophobic surfaces, the model considers the dynamic stability of the air/fluid interface in the asperities. Using the analytical formulations and the stability criteria, systematic studies are presented on the effects of the fractal surface parameters, cylinder radius and Reynolds number on the convective heat transfer characteristics, from which surface texture design maps are developed for maximizing the convection heat transfer. It is shown that multiscale non-wetting surfaces are most effective in the range of lower Reynolds number and small cylinder radius for achieving the best convective heat transfer and thermal hydraulic performance. Applying the models to actual non-wetting surface topographies fabricated using electrodeposition and chemical etching, it is shown that contrary to prevailing notion, superhydrophobicity, characterized by the highest contact angles, does not always lead to the maximum convective heat transfer performance, and that under certain fluid flow conditions, hydrophobic surfaces may offer a greater thermal performance.

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

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2024

Impacts of Topography‐Driven Water Redistribution on Terrestrial Water Storage Change in California Through Ecosystem Responses

Abstract Lateral subsurface flow plays an essential role in sustaining the terrestrial ecosystem, but it is not explicitly represented in most Earth System Models. In this study, we implemented an explicit lateral saturated flow model into the E3SM land model (ELM). The model explicitly describes lateral flow in the saturated zone by representing, for each model grid, an idealized hillslope consisting of five hydrologically connected soil columns. We conducted three model experiments driven by 0.125° atmospheric forcing data during 1980–2015 over California using models of the default ELM, a modified version of ELM to enhance infiltration, and the model with the lateral saturated flow model. The simulated runoff, evapotranspiration, and terrestrial water storage anomaly (TWSA) from the three simulations were evaluated against available observations, and the model explicitly representing lateral flow performs best. The new model produces greater gridcell‐averaged evapotranspiration especially over the mountainous regions with moderate relief and seasonally dry climates. Most importantly, it improves the modeled seasonal variations, interannual variabilities, and the recent decadal decline of TWSA. Many of these improvements can be attributed to the enhanced ecosystem resilience to droughts as demonstrated by transpiration increases caused by lateral flow. Model sensitivity experiments suggest that subsurface runoff is most sensitive to the ratio between horizontal and vertical saturated hydraulic conductivity, followed by hillslope planforms (convergent, divergent, and uniform), number of columns, and lower boundary conditions. Future work should effectively characterize hillslopes in global models and explore the long‐term influences of lateral water movement on modeled biogeochemical cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2022

Estimating biofuel contaminant concentration from 4D ERT with mixing models

Here, we present the results of a lab-scaled feasibility study to assess the performance of electrical resistivity tomography for detection, characterization, and monitoring of fuel grade ethanol releases to the subsurface. Further, we attempt to determine the concentration distribution of the ethanol from the electrical resistivity tomography data using mixing-models. Ethanol is a renewable fuel source as well as an oxygenate fuel additive currently used to replace the known carcinogen methyl tert-butyl ether; however, ethanol is preferentially biodegraded and a cosolvent. When introduced to areas previously impacted by nonethanol-based fuels, it will facilitate the persistence of carcinogenic fuel compounds like benzene and ethylbenzene, as well as remobilize them to the ground water. These compounds would otherwise be retained in the soil column undergoing active or passive remediation processes such as soil vapor extraction or natural attenuation. Here, we introduce ethanol to a saturated Ottawa sand in a tank instrumented for four-dimensional geoelectrical measurements. Forward model results suggest pure phase ethanol released into a water saturated silica sand should present a detectable target for electrical resistivity tomography relative to a saturated silica sand only. We observe the introduction of ethanol to the closed hydraulic system and subsequent migration over the duration of the experiment. One-dimensional and three–dimensional temporal data are assessed for the detection, characterization, and monitoring of the ethanol release. Results suggest one-dimensional geoelectrical measurements may be useful for monitoring a release, while three-dimensional geoelectrical field imaging would be useful to characterize, monitor, and design effective remediation approaches for an ethanol release, assuming field conditions do not preclude the application of geoelectrical methods. We then attempt to use predictive mixing models to calculate the distribution of ethanol concentration within the measurement domain. For this study we examine four different models: a nested parallel mixing model, a nested cubic mixing model, the complex refractive index model (CRIM), and the Lichtenecker-Rother (L-R) model. The L-R model, modified to include an electrical formation factor geometry term, provided the best agreement with expected EtOH concentrations.

09 BIOMASS FUELS↗
DOE OSTI2021

Coupled Thermal-Hydraulic Analysis and Species Mass Transport in a Versatile Experimental Salt Irradiation Loop (VESIL) Using CTF

With the resurgence of interest in molten salt reactors, there is a need for new experiments and modeling capabilities to characterize the unique phenomena present in this fluid fuel system. A Versatile Experimental Salt Irradiation Loop (VESIL) is currently under investigation at Idaho National Laboratory to be placed in the Advanced Test Reactor (ATR). One of the key phenomena this proposed experiment plans to elucidate is fission product speciation in the fuel-salt and the subsequent effects this has on the fuel-salt properties, source term generation, and corrosion control. Specifically, noble gases (Xe & Kr) will bubble out to a plenum or off-gas system, and noble metals (Mo, Tc, Te, etc.) will precipitate and deposit in specific zones in the loop. This work extends the mass transfer and species interaction models in CTF (Coolant-Boiling in Rod Arrays—Two Fluids) and applies these models to give a preliminary estimation of fission product behavior in the proposed VESIL design. A noble metal–helium bubble mass transfer model is coupled with the thermal-hydraulic results from CTF to determine the effectiveness of this insoluble fission product (IFP) extraction method for VESIL. Amounts of IFP species extracted to the off-gas system and species distributions in VESIL after a 60-day ATR cycle are reported.

Walker, Samuel A.↗
DOE OSTI2025

In-situ neutron radiography of pore-water movement in portland cement mortar under subzero temperature exposure: Insights from D 2 O-H 2 O systems

Neutron radiography was employed in this study to monitor real-time pore water movement in portland cement mortar exposed to subzero temperatures, providing novel insights into freeze-thaw (F-T) damage mechanisms. Deuterated water (D 2 O) was used to hydrate the cement in the mortar specimens, and the pore structure was saturated with water (H 2 O), allowing for a contrast that enabled detailed tracking of water movement predominantly in the pore structure. Furthermore, the mortar sample was placed on a chiller plate set to subzero temperatures, thereby providing a one-dimensional temperature gradient vertically in the sample. Neutron radiography revealed that unfrozen water migrated upward as the freezing front advanced, which can create hydraulic pressure within the mortar microstructure and may contribute to the initiation and propagation of microcracks during F-T cycles. Additional characterization of the concrete made with D 2 O showed a marked contrast to the samples made with H 2 O, including delayed hydration, diminished compressive strength, and greater porosity. Therefore, while mortar made with D 2 O allowed for greater neutron contrast to the pores filled with H 2 O, further study is required to produce mortars with similar properties to those made with H 2 O. In addition, further refinement is needed to the radiography experiment to allow for greater control of temperatures and potential quantification of the water concentration.

Air-void system↗
DOE OSTI2020

Permeability Hysteresis From Microchannels Opening During Dissolution/Reprecipitation Cycle

Abstract Permeability is a critical parameter for geological resources characterization. Its evolution with respect to porosity is particularly interesting and many research initiatives focus on deriving such relationships, to understand some hydraulic impacts of microstructure alteration. Permeability evolution from chemical reactions, for instance, can become complex as flow channels may open during rock dissolution. In this contribution, we show that this phenomenon can lead to irregular porosity‐permeability curves and permeability hysteresis after reprecipitation. Current approaches describing permeability as a simple function of porosity can therefore not capture this behavior, and we advocate instead the use of dynamic modeling for such scenarios. We demonstrate our approach by modeling a dissolution/precipitation cycle for a unit cell pore channel and quantify the process at larger scale on three different rock samples, whose microstructures are reconstructed from segmented micro‐Computerised Tomography scans.

Lesueur, Martin↗
DOE OSTI2021

Water balance response of permafrost-affected watersheds to changes in air temperatures

Abstract Observations show increases in river discharge to the Arctic Ocean especially in winter over the last decades but the physical mechanisms driving these changes are not yet fully understood. We hypothesize that even in the absence of a precipitation increase, permafrost degradation alone can lead to increased annual river runoff. To test this hypothesis we perform 12 millennium-long simulations over an idealized hypothetical watershed (1 km 2 ) using a distributed, physically based water balance model (Water flow and Balance Simulation Model, WaSiM). The model is forced by both a hypothetical warming defined by an air temperature increase of 7.5 ∘ C over 100 years, and a corresponding cooling scenario. To assess model sensitivity we vary soil saturated hydraulic conductivity and lateral subsurface flow configuration. Under the warming scenario, changes in subsurface water transport due to ground temperature changes result in a 7%–14% increase in annual runoff accompanied by a 6%–20% decrease in evapotranspiration. The increase in runoff is most pronounced in winter. Hence, the simulations demonstrate that changes in permafrost characteristics due to climate warming and associated changes in evapotranspiration provide a plausible mechanism for the observed runoff increases in Arctic watersheds. In addition, our experiments show that when lateral subsurface moisture transport is not included, as commonly done in global-scale Earth System Models, the equilibrium water balance in response to the warming or cooling is similar to the water balance in simulations where lateral subsurface transport is included. However, the transient changes in water balance components prior to reaching equilibrium differ greatly between the two. For example, for high saturated hydraulic conductivity only when lateral subsurface transport is considered, a period of decreased runoff occurs immediately after the warming. This period is characterized by a positive change in soil moisture storage caused by the soil moisture deficit developed during prior cooling.

Debolskiy, Matvey V. (ORCID:0000000296343627)↗
DOE OSTI2023

Signal Decomposition for Intrusion Detection in Reliability Assessment in Cyber Resilience (Summary Report)

The complexity of assuring cyber resilience for physical process interactions in connected systems such as energy grids increases dramatically as the coupling between processes becomes more direct and responsive. An example of this growing complexity is provided by Integrated Energy Systems (IES), in which various processes such as nuclear heat generation and commodity production are being directly coupled for increased responsiveness to highly variable signals such as market pricing or electricity demand. As such, the potential attack surface of the coupled processes is larger than the two processes independently. Securing these complex systems requires two-fold monitoring: cybersecure monitoring for potential malicious incursion, and physics monitoring for system tampering. Physics monitoring includes analyzing the behavior of the signals within the system for anomalous behavior. This analysis has been shown to be insufficient if approached by only data-driven machine learning and artificial intelligence (MLAI) techniques or only low-level model comparison. Previous efforts at Purdue University suggested combining high-fidelity models with MLAI algorithms as a basis for a software tool for detecting anomalies in physical processes. This work built on that suggestion, developing an advanced library for signal decomposition and analysis using both MLAI and high-fidelity physics algorithms for greatly improved anomaly detection, especially false data injection. This software can be used as part of a secure imbedded intelligence (SEI) system designed under Consequence-driven Cyber-informed Engineering (CCE) for complex coupled systems. This library established a foundation for online and posteriori analysis of digital signals for the purpose of detecting potential malicious tampering in digital signals representing physical processes. Demonstrations carried out throughout the development highlight the effective use of characterization algorithms to detect signal perturbations, particularly triangle attack-style perturbations, in three wide-ranging applications: seismic monitoring, nuclear thermal hydraulics system simulation, and custom manufacturing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Characterization of an electrically conductive proppant for fracture diagnostics

Fracture diagnosis with electromagnetic (EM) and electrical tools requires proppants with high electrical conductivity and mechanical strength. Lab measurements of the electrical and hydraulic conductivity of proppants are critical for selecting the best candidates. Such measurements greatly benefit simulations, field tests, and the ultimate application of such proppants in the field. To that end, a new lab protocol is developed for measuring the electrical and hydraulic conductivity of proppants. The lab setup, which mainly includes a resistivity core holder and a Hassler sleeve core holder, allows for simulation of realistic pressure and temperature conditions when making measurements. Petroleum coke (PC) is proposed as a candidate proppant because of its widespread availability and low cost. Lab measurements show that the effective electrical conductivity of pure PC in a model fracture is approximately 5000 S/m, under a closure stress greater than [Formula: see text] (4000 psi). When PC is mixed with sand, the effective electrical conductivity of the mixture decreases with an increasing weight percentage of sand. Although sand degrades the contact between PC particles, the electrical conductivity stays reasonably high (approximately 1700 S/m) when 50% sand is added. Hydraulic conductivity measurements show that when a fracture is propped with pure PC, the measured fracture conductivity is greater than [Formula: see text] ([Formula: see text]) (dimensionless fracture conductivity greater than 100 for a shale with [Formula: see text] or 100 nD permeability) under a confining pressure of [Formula: see text] (6000 psi). This means that a fracture propped with PC is infinitely conductive in a typical shale formation. When sand is added, the fracture’s hydraulic conductivity becomes even higher, which clearly shows PC’s ability of sustaining high stresses. The proposed protocol provides a robust and effective method that can be generalized for lab testing for other candidate proppants. The data presented clearly show that PC has the potential for field-scale applications in EM hydraulic fracture diagnostics.

Geochemistry & Geophysics↗
DOE OSTI2023

Wave Energy Power Take-Off Validation with a Hydraulicly Actuated Rotary Dynamometer and a Bidirectional High-Power DC Supply: Preprint

There are many organizations working towards the commercialization of wave energy converter technologies and are advancing their designs through the technology readiness levels (TRLs). A critical step before the field deployment of prototype wave energy converters is the validation of the subsystems and components that are contained in the wave energy converter through laboratory testing and performance characterization. In 2021, the National Renewable Energy Laboratory (NREL) developed and demonstrated a system for testing power takeoffs (PTO) with a low-speed, high-torque dynamometer and a grid-tied high-power DC power source and sink before field deployment. The hydraulic dynamometer allows for the simulation of PTO actuation from wave motion and is capable of a wide range of wave periods and heights which are represented as various speeds and torques from the dynamometer. The high-power bidirectional power supply allows for hardware in the loop and controller in the loop testing to be conducted on WEC power electronics. This presentation was made to describe the methods used by NREL research staff to test all components and sub-systems in the PTO of a novel wave energy converter before field deployment.

dynamometer↗
DOE OSTI2023

Wave Energy Power Take-off Validation with a Hydraulicly Actuated Rotary Dynamometer and a Bi-Directional High-Power DC Supply

There are many organizations working towards the commercialization of wave energy converter technologies and are advancing their designs through the technology readiness levels (TRLs). A critical step before the field deployment of prototype wave energy converters is the validation of the subsystems and components that are contained in the wave energy converter through laboratory testing and performance characterization. In 2021, the National Renewable Energy Laboratory (NREL) developed and demonstrated a system for testing power takeoffs (PTO) with a low-speed, high-torque dynamometer and a grid-tied high-power DC power source and sink before field deployment. The hydraulic dynamometer allows for the simulation of PTO actuation from wave motion and is capable of a wide range of wave periods and heights which are represented as various speeds and torques from the dynamometer. The high-power bidirectional power supply allows for hardware in the loop and controller in the loop testing to be conducted on WEC power electronics. This presentation was made to describe the methods used by NREL research staff to test all components and sub-systems in the PTO of a novel wave energy converter before field deployment.

dynamometer↗
DOE OSTI2020

Machine learning for geophysical characterization of brittleness: Tuscaloosa Marine Shale case study

Brittleness is one of the most important reservoir properties for unconventional reservoir exploration and production. Better knowledge about the brittleness distribution can help to optimize the hydraulic fracturing operation and lower costs. However, there are very few reliable and effective physical models to predict the spatial distribution of brittleness. We have developed a machine learning-based method to predict subsurface brittleness by using multidiscipline data sets, such as seismic attributes, rock physics, and petrophysics information, which allows us to implement the prediction without using a physical model. The method is applied on a data set from Tuscaloosa Marine Shale, and the predicted rock physics template is close to the calculated value from conventional inverted elastic parameters. Therefore, the proposed method helps determine areas of the reservoir that have optimal geomechanical properties for successful hydraulic fracturing.

Geochemistry & Geophysics↗