Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reservoir Simulation”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗
NASA NTRS2004

Radicals and Reservoirs in the GMI Chemistry and Transport Model: Comparison to Measurements

The most important use of atmospheric chemistry and transport models is to predict the future composition of the atmosphere. The amounts of gases like chlorofluorcarbons, methyl bromide, nitrous oxide and methane are changing and the stratospheric ozone layer will change because these gases are changing. Methyl bromide, nitrous oxide and methane all have natural sources, and also change because of human activity. Chlorofluorcarbons are man-made gases; these are known to decrease stratospheric ozone and future production is banned. They are long-lived gases, and many decades will pass before they are insignificant in the atmosphere. The models are used to predict changes in ozone and other gases; this is a straightforward application. The models must be also tested using observations for the present day atmosphere. This is a challenging task, because the model contains more than 50 species and more than 150 chemical reactions. Data from satellites, ground stations, aircraft and balloons are used to evaluate the model. Different models that are used in international assessments produce different results; in the most recent assessment some predict that ozone will return to 1980 levels by 2025 and others predict that this will not happen until 2050. Since all the parts of the models are conceptually the same, there must be differences in implementation that produce these differences, This work takes a single model, two different sets of winds and temperatures, and repeats the same prediction for the future. Here we compare the results for these two simulations with many observations. The purpose is to identify differences in the model results for the present atmosphere that will lead to different predictions. This sort of controlled comparison will reduce uncertainty in the predictions for stratospheric ozone.

Douglas, Anne R.↗
DOE OSTI2024

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗
NASA NTRS2017

Synthetic Scenarios from CMIP5 Model Simulations for Climate Change Impact Assessments in Managed Ecosystems and Water Resources: Case Study in South Asian Countries

Increasing population, urbanization, and associated demand for food production compounded by climate change and variability have important implications for the managed ecosystems and water resources of a region. This is particularly true for south Asia, which supports one quarter of the global population, half of whom live below the poverty line. This region is largely dependent on monsoon precipitation for water. Given the limited resources of the developing countries in this region, the objective of our study was to empirically explore climate change in south Asia up to the year 2099 using monthly simulations from 35 global climate models (GCMs) participating in the fifth phase of the Climate Model Inter-comparison Project (CMIP5) for two future emission scenarios (representative concentration pathways RCP4.5 and RCP8.5) and provide a wide range of potential climate change outcomes. This was carried out using a three-step procedure: calculating the mean annual, monsoon, and non-monsoon precipitation and temperatures; estimating the percent change from historical conditions; and developing scenario funnels and synthetic scenarios. This methodology was applied for the entire south Asia region; however, the percent change information generated at 1.5deg grid scale can be used to generate scenarios at finer spatial scales. Our results showed a high variability in the future change in precipitation (-23% to 52%, maximum in the non-monsoon season) and temperature (0.8% to 2.1%) in the region. Temperatures in the region consistently increased, especially in the Himalayan region, which could have impacts including a faster retreat of glaciers and increased floods. It could also change rivers from perennial to seasonal, leading to significant challenges in water management. Increasing temperatures could further stress groundwater reservoirs, leading to withdrawal rates that become even more unsustainable. The high precipitation variability (with higher propensity for localized intense rainfall events) observed in the region can be a key factor for managed ecosystems and water management and could also lead to more incidence of severe urban flooding. The results could be used to assess both mitigation and adaptation alternatives to reduce vulnerabilities in managed ecosystems (agricultural and urban) and water resources.

agriculture↗
DOE OSTI2020

Pore size effect on selective gas transport in shale nanopores

In shale gas production, gas composition may vary over time. To understand this phenomenon, we use molecular dynamics simulations to study the permeation of CH 4 , C 2 H 6 and their mixture from a source container through a pyrophyllite nanopore driven by a pressure gradient. For a pure gas, the flow rate of CH 4 is always higher than that of C 2 H 6 , regardless of pore size. For a 1:1 C 2 H 6 : CH 4 mixture, however, C 2 H 6 :CH 4 flow rate ratio is higher than the compositional ratio in the container (i.e., 1:1) when the pore size is smaller than ~1.8 nm. The selective transport is caused by the competitive adsorption of C 2 H 6 over CH 4 in the nanopore. The selectivity is also determined by the interplay between the surface diffusion of the adsorbed molecules and the viscous flow in the center of the pore, and it diminishes as the viscous flow becomes to dominate the surface diffusion when the pore size becomes larger than 1.8 nm. Our work shows that compositional differentiation of shale gas in production is a consequence of nanopore confinement and therefore a key characteristic of an unconventional reservoir. Finally, the related compositional information can potentially be used for monitoring the status of a production well such as its recovery rate.

03 NATURAL GAS↗
DOE OSTI2022

Experimental study of fracture slip due to stress perturbation in fractured geo-resources

Anthropogenic activities such as underground fluid injection/extraction, well stimulation, and dam impoundment will perturb the state-of-stress in the subsurface, which could activate critically stressed pre-existing fractures and faults. Fault and fracture activation can induce seismicity, change flow behavior, and damage asperities on the fracture surface to alter reservoir permeability. Here, in this study, triaxial shear tests were conducted on saw cut and rough fracture Barre granite specimens to better characterize how a decrease in normal stress and increase in shear stress will influence crystalline rock fracture hydromechanical properties. In a triaxial set up, confining pressure (CP) and differential stress (DS) were used to induce normal and shear stress on a 45° oriented fracture. Fracture slip and the corresponding stress relaxation was observed when the confining stress decreased. Fracture permeability increased when the fracture slipped. Moreover, new micro-cracks observed by x-ray images indicated the creation of new flow paths for fluid movement. The experimental outcomes of this study indicate that decrease in normal stress due to stress perturbations along the faults/fractures in the extensional stress regime can induce the slip along the fracture. This could have a significant influence on slope stability in quarries, stimulation and injection strategy in deep geo-resources, and potentially minimize the risk of induced seismicity. The observed stress dependency of hydraulic aperture in rough fracture during pre- and post-slip stages could provide key inputs for field-scale simulation of geo-resources.

58 GEOSCIENCES↗
DOE OSTI2025

Echo state network for coarsening dynamics of charge density waves

An echo state network (ESN) is a type of reservoir computer that uses a recurrent neural network with a sparsely connected hidden layer. Compared with other recurrent neural networks, one great advantage of ESN is the simplicity of its training process. Yet, despite the seemingly restricted learnable parameters, ESN has been shown to successfully capture the spatial-temporal dynamics of complex patterns. Here we build an ESN to model the coarsening dynamics of charge-density waves (CDWs) in a semiclassical Holstein model, which exhibits a checkerboard electron density modulation at half-filling stabilized by a commensurate lattice distortion. The inputs to the ESN are local CDW order parameters in a finite neighborhood centered around a given site, while the output is the predicted CDW order of the center site at the next time step. Special care is taken in the design of couplings between hidden layer and input nodes to ensure lattice symmetries are properly incorporated into the ESN model. Since the model predictions depend only on CDW configurations of a finite domain, the ESN is scalable and transferrable in the sense that a model trained on dataset from a small system can be directly applied to dynamical simulations on larger lattices. Furthermore, our work opens avenues for efficient dynamical modeling of pattern formations in functional electron materials.

2-dimensional systems↗
DOE OSTI2026

Thermal Reservoir Networks for Modularly Expandable Thermal Microgrids

The Department of Defense (DoD) faces the substantial challenge of cost-effectively retrofitting one to two installations per month, each comprising approximately 1,000 buildings, to improve resilience, reduce energy consumption, and enhance energy supply security. Achieving these objectives requires optimal system selection and effective risk mitigation during system integration. To address this need, we introduce Platform-Based Design (PBD), a structured, hierarchical methodology adapted from other industrial sectors to the domain of energy system retrofits. We demonstrate the effectiveness of PBD through a techno-economic feasibility study comparing geothermal-coupled thermal energy networks (TENs) with conventional energy systems for heating, cooling, and powering 17 buildings at Joint Base Andrews (JBA) in Maryland. Our analysis illustrates that the PBD approach enables rigorous, data-driven, sequential decision making, resulting in a family of Pareto-optimal systems, among which the TEN emerged as the most promising solution. The selected TEN design integrates geothermal borefields, heat recovery heat pumps, photovoltaic (PV) arrays, and battery storage. Compared to the baseline system – gas heating combined with air-source chillers – the proposed TEN reduces annual imported energy by 74% and peak electricity demand by 45%, achieves a levelized cost of energy of $\$0.210$/kWh, and substantially enhances resilience. Life-cycle costs increase by approximately 6%, and initial investment costs are about 2.5 times higher than the baseline. However, if central plant infrastructure, district loops, and utility-scale PV and battery systems are privately funded and operated, the initial investment would fall below the baseline system cost. Critical to achieving these significant performance improvements were detailed nonlinear dynamic simulations coupling geothermal heat transfer, energy system operation, and realistic feedback control logic. These simulations identified essential design modifications and control strategy refinements that substantially reduced energy use, peak demand, and compressor shortcycling, thereby improving durability and reliability—issues that would have been significantly more expensive to resolve during operation. Additionally, the verification step highlighted sensitivities to key design parameters that could reduce initial investment by approximately $\$2$ million and reduce annual life-cycle costs more than $\$300,000$. We recommend adopting the PBD methodology for future feasibility studies and TEN pilot projects to gain valuable operational experience. Furthermore, we recommend that DoD invest in transferring and scaling the PBD methodology to other installations. This entails developing standardized computational frameworks and component libraries as well as training industry in conducting PBD. Such investments would enable rapid, robust, reliable, and cost-effective retrofits, supporting DoD’s ambitious energy system modernization goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Ab initio quantum scattering calculations and a new potential energy surface for the HCl( X 1Σ+)–O2(X3Σg−) system: Collision-induced line shape parameters for O2-perturbed R(0) 0–0 line in H35Cl

The remote sensing of abundance and properties of HCl—the main atmospheric reservoir of Cl atoms that directly participate in ozone depletion—is important for monitoring the partitioning of chlorine between “ozone-depleting” and “reservoir” species. Such remote studies require knowledge of the shapes of molecular resonances of HCl, which are perturbed by collisions with the molecules of the surrounding air. In this work, we report the first fully quantum calculations of collisional perturbations of the shape of a pure rotational line in H35Cl perturbed by an air-relevant molecule [as the first model system we choose the R(0) line in HCl perturbed by O2]. The calculations are performed on our new highly accurate HCl(X1Σ+)–O2(X3Σg−) potential energy surface. In addition to pressure broadening and shift, we also determine their speed dependencies and the complex Dicke parameter. This gives important input to the community discussion on the physical meaning of the complex Dicke parameter and its relevance for atmospheric spectra (previously, the complex Dicke parameter for such systems was mainly determined from phenomenological fits to experimental spectra and the physical meaning of its value in that context is questionable). We also calculate the temperature dependence of the line shape parameters and obtain agreement with the available experimental data. We estimate the total combined uncertainties of our calculations at 2% relative root-mean-square error in the simulated line shape at 296 K. This result constitutes an important step toward computational population of spectroscopic databases with accurate ab initio line shape parameters for molecular systems of terrestrial atmospheric importance.

Chemistry↗
DOE OSTI2021

Development of a Framework for Data Integration, Assimilation, and Learning for Geological Carbon Sequestration (DIAL-GCS) (Final Report)

This project aimed to develop and demonstrate a Data Integration, Assimilation, and Learning framework for geologic carbon sequestration projects (DIAL-GCS). DIAL-GCS is an intelligence monitoring system (IMS) for automating GCS closed-loop management by leveraging recent developments in machine learning technologies, complex event processing (CEP), and reduced-order modeling. The safe and efficient operation of GCS repositories requires integrated monitoring to track the injected CO¬2 as it moves within a storage reservoir. GCS projects are data intensive, as a result of proliferation of digital instrumentation and smart-sensing technologies. GCS projects are also resource intensive, often requiring multidisciplinary teams performing different monitoring, verification, accounting (MVA) tasks throughout the lifecycle of a project to ensure secure containment of injected CO2. The success of GCS thus depends in a large part on our ability to access, assimilate, and analyze heterogeneous data and information sources in a timely manner. This project included a number of meaningful and necessary tasks to transform the human domain knowledge into machine-interpretable rules for automating knowledge extraction and discovery in GCS. The specific technical objectives of the proposed DIAL-GCS project were to develop an ontology-driven GCS data management module for storing, querying, and exchanging GCS data (both historic and live sensor data) from multiple sources and in heterogeneous formats. Incorporate a CEP engine for detecting abnormal situations by seamlessly combining expert knowledge, rule-based reasoning, and machine learning. Enable uncertainty quantification and predictive analytics using a combination of coupled-process modeling, AI/ML methods, and reduced-order modeling, and integrate and demonstrate the system’s capabilities with both real and simulated data. As far as we know, this is one of the first projects aimed to develop intelligent monitoring systems (IMS) targeting the GCS. Under this project, the team had developed a large number of web applications and scientific algorithms that contribute the main theme of intelligent monitoring. The team has published more than a dozen peer reviewed papers and disseminated the research results at multiple technical meetings.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Financial Analysis of the Smallmouth Bass Flows implemented at the Glen Canyon Dam during Water Year 2024

The Glen Canyon Dam (GCD) is a Colorado River Storage Project (CRSP) power resource that is a component of the Salt Lake City Area Integrated Projects (SLCA/IP). The 2016 record of decision (ROD) for the GCD long-term experimental and management plan (LTEMP) final Environmental Impact Statement (EIS) specifies criteria for GCD monthly water releases, daily and hourly operating limits, and experimental releases. This report presents a financial analysis of the Smallmouth bass (Micropterus dolomieu) (SMB) flows implemented at GCD during Water Year (WY) 2024. These bypass flows were introduced by the U.S. Bureau of Reclamation (USBR) as an emergency response to the growing threat posed by invasive SMB in the Colorado River ecosystem downstream of the dam. SMB are a non-native predatory species that pose a significant threat to native fish populations, including the endangered humpback chub (Gila cypha). The thermal regime below GCD, typically cold due to hypolimnetic releases from Lake Powell, has historically served as a thermal barrier limiting SMB establishment. However, persistently low reservoir levels in recent years have reduced stratification in Lake Powell, allowing warmer water to be released downstream. This has enabled SMB to spawn successfully below the dam, prompting urgent ecological concerns. To mitigate the risk of SMB proliferation, the USBR implemented a series of bypass flows in WY 2024. Drawn from a lower elevation than the penstocks, the bypass structures released cooler water downstream. These short-duration bypass flows aimed to keep temperatures cool enough to prevent SMB from spawning, thereby reducing the ecological threat posed by this invasive species. Although motivated by ecological objectives, these bypass flows came with financial tradeoffs. Releasing water through the bypass structures instead of the turbines at GCD reduced hydropower generation, resulting in a significantly lower financial position for Western Area Power Administration (WAPA), which is responsible for marketing the electricity produced by the GCD Powerplant. This report analyzes the financial impact of the SMB flows implemented from July to November 2024. These experimental releases led to an estimated financial cost of approximately $18.9 million, primarily driven by the substantial volume of water diverted through the bypass structures. This study applies an integrated set of tools to estimate WAPA financial impacts by simulating GCD under two types of cases; namely, (1) a “With Experiment” case that mimics the water operations that actually occurred, including the SMB bypass flows, and (2) a “Without Experiment” case that simulates operations under the assumption that the SMB flows did not occur. Both cases comply with LTEMP hourly and daily operating criteria, and the monthly water release volumes are assumed to be identical under both cases. The Colorado River Storage Project Python-based model (CRiSPPy) model was the main modeling tool used to simulate the dispatch of the GCD hydropower plant and associated water releases from Lake Powell. In the modeling process, the research team used extensive data sets and historical information on SLCA/IP power plant characteristics, hydrologic conditions, and WAPA’s power purchases and sales prices.

13 HYDRO ENERGY↗
NASA NTRS2012

Crystallization of Yamato 980459 at 0.5 GPA: Are Residual Liquids Like QUE 94201?

The Martian basaltic meteorites Y980459 and QUE94201 (henceforth referred to as Y98 and QUE respectively) are thought to represent magmatic liquid compositions, rather than being products of protracted crystallization and accumulation like the majority of other martian meteorites. Both meteorite compositions have been experimentally crystallized at 1 bar, and liquidus phases were found to match corresponding mineral core compositions in the meteorites, consistent with the notion that these meteorites represent bona fide melts. They also represent the most primitive and most evolved basaltic martian samples, respectively. Y98 has Mg# (molar Mg/Mg+Fe) approximates 65, and lacks plagioclase; whereas QUE has Mg# approximates 40, and lacks olivine. However they share important geochemical characteristics (e.g. superchondritic CaO/Al2O3, very high epsilon(sub Nd) and low Sr-87/Sr-87) that suggest they sample a similar highly depleted mantle reservoir. As such, they represent likely endmembers of martian magmatic liquid compositions, and it is natural to seek petrogenetic linkages between the two. We make no claim that the actual meteorites themselves share a genetic link (the respective ages rule that out); we are exploring only in general whether primitive martian liquids like Y98 could evolve to liquids resembling QUE. Both experimental and computational efforts have been made to determine if there is indeed such a link. Recent petrological models at 1 bar generated using MELTS suggest that a QUE-like melt can be derived from a parental melt with a Y98 composition. However, experimental studies at 1 bar have been less successful at replicating this progression. Previous experimental crystallization studies of Y98 by our group at 0.5 GPa have produced melt compositions approaching that of QUE, although these results were complicated by the presence of small, variable amounts of H2O in some of the runs owing to the use of talc/pyrex experimental assemblies. Therefore we have repeated the four experiments, augmented with additional runs, all using BaCO3 cell assemblies, which are devoid of water, and these new experiments supersede those reported earlier. Here we report results of experiments simulating equilibrium crystallization; fractional crystallization experiments are currently underway.

Rapp, J. F.↗
DOE OSTI2020

Gaining Perspective on Unconventional Well Design Choices through Play-level Application of Machine Learning Modeling

The recent development of unconventional oil and gas (O&G) reservoirs has led to an abundant hydrocarbon supply, both domestically and globally. However, there is a continued push to develop new and innovative approaches to improve exploration and extraction efficiencies and overall well productivity moving forward. Substantial improvements in unconventional O&G development are expected through optimized well completion and stimulation strategies aimed at maximizing well productivity. Optimizing well designs will require tailoring to the distinctive geologic conditions present for any newly placed well. To better evaluate the impact of well design attributes and their associated interactions on productivity in a major unconventional play, multivariate machine learning-based models that use empirical datasets were developed. A gradient boosted regression tree (GBRT) algorithm was applied. GBRT has been narrowly investigated for O&G applications but enables straightforward parametric importance and influence evaluation, as well as assessment of parameter interaction effects. Models were trained on well design and locational parameters that serve as a proxy for variable geologic conditions to estimate two types of productivity indicator response variables strongly correlated to estimated ultimate recovery (EUR). The dataset utilized consists of over 7,000 well observations that cover the majority of the productive region of the Marcellus Shale. Model performance was evaluated and algorithm parameters tuned by analyzing the goodness-of-fit for simulated results against observed data in a cross-validation approach. Models were found capable of 73–79 percent prediction accuracy on held out testing data of gas equivalent production and can be used to inform future well design and placement decisions for increasing EUR per well and improving overall field-level recovery. Study results indicate that Marcellus well performance improves most with upscaling perforated interval lengths and water and proppant volumes per foot; but relative productivity improvements are spatially dependent across the play. Finally, optimal combinations of water and proppant on well performance were found to vary depending on well location, emphasizing the utility of data-driven models capable of broad application across a play of interest for informing tailored well design approaches prior to their field deployment.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2023

Review and Summary of Oxide Thickness Data for Aluminum-Clad Spent Nuclear Fuel

The U.S. Department of Energy (DOE) owns a large inventory of aluminum-clad spent nuclear fuel (ASNF) in interim storage pending ultimate disposition, with more being generated by currently operating research reactors. Dry storage in sealed DOE Standard Canisters is being investigated as an approach for long-term interim storage and/or disposition in a repository for ASNF. A primary challenge for ASNF storage is the presence of aluminum (oxy)hydroxide layers formed on the cladding surfaces during water exposure in the reactor and in wet storage, which forms a reservoir of chemisorbed water not readily removed at low (<100°C) drying temperatures. Free, physisorbed, and chemisorbed water are all susceptible to radiolytic breakdown under irradiation and can release hydrogen gas. Identifying the likely range of (oxy)hydroxide loadings on ASNF that may be placed in dry storage will help to ensure that the impact of the (oxy)hydroxide is adequately accounted for while avoiding over conservatism and enable mitigation strategies to be implemented where needed. The current report summarizes information on (oxy)hydroxide thicknesses and characteristics from the literature and from recent measurements taken in the present research campaign. In general, corrosion studies have indicated that the corrosion and oxide buildup on aluminum are affected by numerous conditions, including pH, temperature and heat flux, coolant flowrates, irradiation (in-reactor vs. unirradiated tests), duration of water exposure, and the amount of oxide already on the surface. Some of these factors are interrelated (e.g., the local temperature, heat flux, and coolant flowrate). Observations of the impact of irradiation (in-reactor versus ex-reactor measurements) on corrosion and oxide thickness are mixed, but multiple studies indicate that there is a significant difference between in-reactor and ex-reactor corrosion kinetics, even when other operating conditions such as pH, flowrate, heat flux, etc., simulate those in a reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Final Report DE-FE0031785 Mohsen Ahmadian, Ph.D. Demonstration of Proof of Concept of a Multiphysics Approach for Real-Time Remote Monitoring of Dynamic Changes in Pressure and Salinity in Hydraulically Fractured Networks

Hydraulic fracturing has evolved into a multistep process with varying flow rates, carrier fluids (e.g., gel or slickwater), proppant loadings, and proppant grain sizes. As a result, primary recovery from a hydraulically fractured tight-oil reservoir is often a tiny fraction of the original oil in place, ranging between 5 and 10%. As stated in the FOA1990, “part of this problem is due to the inability of current well completion processes to effectively stimulate the entire reservoir volume in contact with the wellbore. Innovative technologies are needed that can help improve the effectiveness of reservoir completion methods, maximize stimulated reservoir volumes, and optimize recovery over the entire producing life span of a well”. We first need to enhance the current fracture diagnostic techniques to improve a well-completion design. However, detecting and delineating a subsurface hydraulic fracture is extremely difficult because the induced fracture network is only fractionally propped, and these propped fractures are generally very thin. Microseismic and tiltmeter monitoring techniques can provide information on the fracture extent but provide little or no information on the movement and final distribution of proppant or production fluids. On the other hand, electromagnetic (EM) imaging has shown the capability to monitor proppant distribution throughout the fracture area, especially in the presence of Electrically Active Proppants (EAPs). A previous EM survey of hydraulic fracturing at the Devine Fracture Pilot Site (DFPS) and subsequent EM code developments demonstrated this survey as a robust technique to remotely interrogate the extent of the EAP-filled hydraulic fracture during its propagation. The objectives of the project were threefold: (1) to capitalize on the material properties of an EAP to demonstrate remote monitoring of relative changes in pressure, pressure, and flow that are commonly encountered during production from a hydraulically fractured reservoir; (2) to evaluate EM imaging tools, to achieve Objective 1 in near real-time; and (3) to develop a multi-physics joint inversion approach to precisely predict flow patterns and physiochemical changes within an EAP-filled fracture network. This research project was built upon our previous work at the Devine Test Site managed by the Bureau of Economic Geology (BEG) at The University of Texas at Austin (UT-Austin). It also leveraged a significant investment from the Advanced Energy Consortium (AEC) to address the DOE's interest in subsurface flow, containment, and characterization by multiple signals. This three-year and three-month project succeeded in demonstrating the feasibility of a real-time dynamic fluid flow mapping technique at Technology Readiness Level 5 (TRL5) by utilizing a commercially available surface-based Controlled-Source Electromagnetic (CSEM) method (Objectives 1, 2). We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, the observed electric field in our study is influenced by bottomhole pressure, flow rate, and salinity, which is demonstrated by comparing these parameters with the electrical field potentials. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce the whole measured electric field magnitudes. Preliminary estimates showed that including Streaming Potential (SP) in our geophysical model is likely needed to reduce the simulation misfit.

02 PETROLEUM↗
NASA NTRS2018

Modeling the Distribution and Type of High-Latitude Natural Wetlands for Methane Studies

High latitude (>50N) natural wetlands emit a substantial amount of methane to the atmosphere, and are located in a region of amplified warming. Northern hemisphere high latitudes are characterized by cold climates, extensive permafrost, poor drainage, short growing seasons, and slow decay rates. Under these conditions, organic carbon accumulates in the soil, sequestering CO2 from the atmosphere. Methanogens produce methane from this carbon reservoir, converting stored carbon into a powerful greenhouse gas. Methane emission from wetland ecosystems depends on vegetation type, climate characteristics (e.g, precipitation amount and seasonality, temperature, snow cover, etc.), and geophysical variables (e.g., permafrost, soil type, and landscape slope). To understand how wetland methane dynamics in this critical region will respond to climate change, we have to first understand how wetlands themselves will change and therefore, what the primary controllers of wetland distribution and type are. Understanding these relationships permits data-anchored, physically-based modeling of wetland distribution and type in other climate scenarios, such as paleoclimates or future climates, a necessary first step toward modeling wetland methane emissions in these scenarios. We investigate techniques and datasets for predicting the distribution and type of high latitude (>50N) natural wetlands from a suite of geophysical and climate predictors. Hierarchical clustering is used to derive an empirical methane-centric wetland model. The model is applied in a multistep process first to predict the distribution of wetlands from relevant geophysical parameters, and then, given the predicted wetland distribution, to classify the wetlands into methane-relevant types using an expanded suite of climate and biogeophysical variables. As the optimum set of predictor variables is not known a priori, the model is applied iteratively, and each simulation is evaluated with respect to observed high-latitude wetlands.

Romanski, Joy↗
NASA NTRS1999

Hyperspectral Remote Sensing and Ecological Modeling Research and Education at Mid America Remote Sensing Center (MARC): Field and Laboratory Enhancement

The purpose of this project was to establish a new hyperspectral remote sensing laboratory at the Mid-America Remote sensing Center (MARC), dedicated to in situ and laboratory measurements of environmental samples and to the manipulation, analysis, and storage of remotely sensed data for environmental monitoring and research in ecological modeling using hyperspectral remote sensing at MARC, one of three research facilities of the Center of Reservoir Research at Murray State University (MSU), a Kentucky Commonwealth Center of Excellence. The equipment purchased, a FieldSpec FR portable spectroradiometer and peripherals, and ENVI hyperspectral data processing software, allowed MARC to provide hands-on experience, education, and training for the students of the Department of Geosciences in quantitative remote sensing using hyperspectral data, Geographic Information System (GIS), digital image processing (DIP), computer, geological and geophysical mapping; to provide field support to the researchers and students collecting in situ and laboratory measurements of environmental data; to create a spectral library of the cover types and to establish a World Wide Web server to provide the spectral library to other academic, state and Federal institutions. Much of the research will soon be published in scientific journals. A World Wide Web page has been created at the web site of MARC. Results of this project are grouped in two categories, education and research accomplishments. The Principal Investigator (PI) modified remote sensing and DIP courses to introduce students to ii situ field spectra and laboratory remote sensing studies for environmental monitoring in the region by using the new equipment in the courses. The PI collected in situ measurements using the spectroradiometer for the ER-2 mission to Puerto Rico project for the Moderate Resolution Imaging Spectrometer (MODIS) Airborne Simulator (MAS). Currently MARC is mapping water quality in Kentucky Lake and vegetation in the Land-Between-the Lakes (LBL) using Landsat-TM data. A Landsat-TM scene of the same day was obtained to relate ground measurements to the satellite data. A spectral library has been created for overstory species in LBL. Some of the methods, such as NPDF and IDFD techniques for spectral unmixing and reduction of effects of shadows in classifications- comparison of hyperspectral classification techniques, and spectral nonlinear and linear unmixing techniques, are being tested using the laboratory.

Cetin, Haluk↗
DOE OSTI2020

From pools to flow: The PROMISE framework for new insights on soil carbon cycling in a changing world

Soils represent the largest terrestrial reservoir of organic carbon, and the balance between soil organic carbon (SOC) formation and loss will drive powerful carbon-climate feedbacks over the coming century. To date, efforts to predict SOC dynamics have rested on pool-based models, which assume classes of SOC with internally homogenous physicochemical properties. However, emerging evidence suggests that soil carbon turnover is not dominantly controlled by the chemistry of carbon inputs, but rather by restrictions on microbial access to organic matter in the spatially heterogeneous soil environment. The dynamic processes that control the physicochemical protection of carbon translate poorly to pool-based SOC models; as a result, we are challenged to mechanistically predict how environmental change will impact movement of carbon between soils and the atmosphere. In this study, we propose a novel conceptual framework to explore controls on belowground carbon cycling: Probabilistic Representation of Organic Matter Interactions within the Soil Environment (PROMISE). In contrast to traditional model frameworks, PROMISE does not attempt to define carbon pools united by common thermodynamic or functional attributes. Rather, the PROMISE concept considers how SOC cycling rates are governed by the stochastic processes that influence the proximity between microbial decomposers and organic matter, with emphasis on their physical location in the soil matrix. We illustrate the applications of this framework with a new biogeochemical simulation model that traces the fate of individual carbon atoms as they interact with their environment, undergoing biochemical transformations and moving through the soil pore space. We also discuss how the PROMISE framework reshapes dialogue around issues related to SOC management in a changing world. We intend the PROMISE framework to spur the development of new hypotheses, analytical tools, and model structures across disciplines that will illuminate mechanistic controls on the flow of carbon between plant, soil, and atmospheric pools.

54 ENVIRONMENTAL SCIENCES↗