Search NASA⌕ Search

SEARCH · Search NASA

Results for “Energy sciences”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Realizing the scientific program with polarized ion beams at the future BNL Electron Ion Collider

Polarized ion beams at the Electron Ion Collider (EIC) are essential to address some of the most important open questions at the twenty-first century frontiers of understanding of the fundamental structure of matter. Here, in this work, we summarize the science case and identify polarized 2 H, 3 He, 6 Li, and 7 Li ion beams as critical technology that will enable experiments which address the most important science. Furthermore, we discuss the required ion polarimetry and spin manipulation at the EIC. The current EIC accelerator design is presented. We identify a significant research and development effort across national and international laboratories and universities that is required over about a decade to realize the polarized ion beams and estimate (based on previous experience) that it will require about 20 full-time equivalent (FTE) over 10 years (or a total of about 200 FTE-years) of personnel, including graduate students, postdoctoral researchers, technicians, and engineers. Attracting, educating, and training a new generation of physicists in experimental spin techniques will be essential for the successful realization. Artificial intelligence and machine learning are seen as having significant potential for both acceleration of research and development and amplification of discovery in the optimal realization of this unique quantum technology on a cutting-edge collider. The research and development effort is synergistic with research in atomic physics and fusion energy science.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Support of US Burning Plasma Organization (Final Technical Report)

The objective of the project described here was to support and increase the effectiveness of the US Burning Plasma Organization. This a national organization of scientists and engineers in magnetic fusion energy, established in 2005 at the request of the DOE Office of Fusion Energy Sciences, which serves to “Advance the scientific understanding of burning plasmas and ensure the greatest benefit from a burning plasma experiment by coordinating relevant U.S. fusion research with broad community participation.” It has grown to over 420 members from over 70 institutions, primarily across the US but with some associate members from other countries. The project was divided into two tasks. The first supported leadership activities in the USBPO by the Deputy Director, Dr. Amanda Hubbard. Responsibilities included chairing regular meetings of the USBPO Research Committee, communicating with the USBPO membership via emails and web seminars, forming and overseeing activities of Task Groups, and ensuring good communication about activities of the International Tokamak Physics Activity. A task group report on ‘Recommendations for ITER Experimental Operation, U.S. Team Formation and Participation’ was completed. Hubbard organized a series of highly popular web seminars for the burning plasma community, on a range of topics attracting a large number of participants from diverse institutions. Hubbard ended her term as Deputy Director in May, 2020. She subsequently received the Secretary of Energy’s Appreciation award for her service to USBPO. The second task supported the communications functions of the USBPO, led by Communications Coordinator, Mark London. Communications are critical for the effective operation of the USBPO, given it consists of members at many widely scattered institutions, and conducts its business almost entirely using remote tools. MIT has hosted and maintained the USBPO website, mailing lists, and discussion forums and provided member support, and videoconferencing support. MIT has continued this support following this grant, with funding incorporated into a broader Cooperative Agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Galaxy Clustering with LSST: Effects of Number Count Bias from Blending

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will survey the southern sky to create the largest galaxy catalog to date, and its statistical power demands an improved understanding of systematic effects such as source overlaps, also known as blending. In this work we study how blending introduces a bias in the number counts of galaxies (instead of the flux and colors), and how it propagates into galaxy clustering statistics. We use the 300 deg 2 DC2 image simulation and its resulting galaxy catalog (LSST Dark Energy Science Collaboration et al. 2021) to carry out this study. We find that, for a LSST Year 1 (Y1)-like cosmological analyses, the number count bias due to blending leads to small but statistically significant differences in mean redshift measurements when comparing an observed sample to an unblended calibration sample. In the two-point correlation function, blending causes differences greater than 3σ on scales below approximately 10', but large scales are unaffected. We fit Ω m and linear galaxy bias in a Bayesian cosmological analysis and find that the recovered parameters from this limited area sample, with the LSST Y1 scale cuts, are largely unaffected by blending. Our main results hold when considering photometric redshift and a LSST Year 5 (Y5)-like sample.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

Solidification Using a Baffle in Sealed Ampoules (SUBSA)

Solidification Using a Baffle in Sealed Ampoules (SUBSA) will be the first materials science experiment conducted in the Microgravity Science Glovebox (MSG) Facility at the International Space Station (ISS) Alpha. The launch is schedule for May 31, 2002. Using the specially developed furnace, 10 Te and Zn-doped single crystals of InSb will be directionally solidified in microgravity. A key goal of the SUBSA investigation is to (i) clarify the origin of the melt motion in space laboratories and (ii) to reduce the magnitude of the melt motion to the point that it does not interfere with the transport phenomena. These goals will be accomplished through a special ampoule and furnace design. A disk-shaped baffle, positioned close to the freezing front, is used to reduce melt motion. Furthermore, the solidification will be visualized by using a transparent furnace, with a video camera, continuously sending images to the earth. This allows detection of bubbles and melt de-wetting that could cause surface tension driven convection. In preparation for the space experiments, 30 ground-based experiments were conducted. The results of ground based tests and numerical modeling will be presented. Based on numerical modeling, 12 mm 1D silica ampoules were selected. The small diameter ampoule favors closer placement of the baffle to the interface, without excessive radial segregation caused by forced convection while providing more damping of natural convection. The parts in the silica ampoule include 2 carbon springs made by Energy Science Laboratories, Inc., a pyrocarbon-coated graphite cylinder, pyrocarbon-coated graphite a baffle with the shaft and the InSb charge with the seed crystal grown by W.A. Bonner of Crystallod Inc.

Marin, C.↗
DOE OSTI2026

Synthetic data-driven deep learning for label-free autonomous atomic force microscopy

Atomic force microscopy (AFM) is a widely used tool for nanoscale characterization across materials science, energy research, and biology. However, its adoption in high-throughput materials discovery and statistically driven studies remains limited by a strong dependence on expert operator input and by the scarcity of annotated experimental AFM datasets needed to enable data-driven automation. Here, we introduce SimuScan, a synthetic-data–driven framework that enables reliable AFM feature identification, segmentation, and targeted imaging without requiring large manually labeled experimental datasets. SimuScan generates tunable, high-fidelity synthetic AFM images of defined morphologies while incorporating realistic experimental artifacts, including tip–sample convolution, noise, flattening distortions, and surface debris. These datasets are shown to support scalable, label-free training of modern deep learning models for AFM analysis. When integrated into data-driven AFM workflows, SimuScan-trained models can locate and analyze nanoscale structures across large datasets and guide targeted follow-up imaging. We validate this approach on nanostructured surfaces, DNA assemblies, and bacterial cells, demonstrating robust generalization across diverse sample types with minimal operator intervention. More broadly, this work establishes a general strategy for generating explicitly conditioned, task-relevant synthetic data to improve the reliability of downstream models in autonomous microscopy.

Millan-Solsona, Ruben [Oak Ridge National Laborato↗
DOE OSTI2023

Fundamental Studies of Bifunctional Catalysts for Tandem Reactions

During the past few years of our funding from DOE Basic Energy Sciences, we have focused on the catalytic conversion of biomass-derived 5-hydroxymethyl furfural (HMF) as a platform molecule for sustainable synthesis of value-added chemicals. The combination of the diene and aldehyde functionalities in HMF enables catalytic production of acetalized HMF derivatives with diol or epoxy reactants to allow reversible synthesis of norcantharimide derivatives upon Diels-Alder reaction with maleimides. Whereas the electron-withdrawing nature of aldehyde group inhibits the reactivity of diene group in HMF to undergo Diels-Alder reactions with dienophiles, by converting the aldehyde group in HMF to an electron-donating hydroxyl group, HMF can append maleimide-based chemicals through Diels-Alder coupling. Acetalization of HMF not only produces an electron-rich diene for Diels-Alder reaction, but it also alters the reversibility of acetal formation and thereby allows for the controlled release of chemicals by hydrolysis. Therefore, we synthesized various norcantharimide derivatives from biomass-derived HMF by acetalization over an acid catalyst, Amberlyst-15, followed by Diels-Alder reaction with maleimide. The norcantharimides release the starting materials by retro Diels-Alder reaction that is triggered by acetal hydrolysis under acidic (≤pH 3) conditions. In subsequent work, we explored the synthesis and properties of functional polyurethanes and polyesters with tunable properties from biomass-derived (HMF)-Acetone-HMF (HAH) monomers. HAH can be selectively hydrogenated over Cu and Ru catalysts to produce partially-hydrogenated HAH (PHAH) and fully-hydrogenated HAH (FHAH). HAH contains functional groups, including hydroxyl, furan, enone, and ketone functionalities, that can be exploited to further tune the polymer properties. The π-electron conjugation between the enone and furan groups in HAH monomer stabilizes the furan group to prevent the uncontrolled degradation of furans during reactions, such as etherification. Therefore, HAH is a renewable monomer that can be used to synthesize polymers, having high molecular weight, symmetric functionalities, an inexpensive production price, and being derived from renewable resources without competition with food resources. We demonstrated the synthesis of functional polyurethanes and a renewable polyester from HAH-derived monomers, and 4,4’-methylenebis(phenyl isocyanate) (MDI). We showed that it is essential to be able to selective hydrogenate C=C double bonds in HAH to be able to achieve effective Diels Alder reactions of the furan ring. Therefore, we studied the hydrogenation at temperatures from 313 – 393 K of HAH over Pd, Ru, and Cu based catalysts. HAH was selectively hydrogenated to produce partially-hydrogenated monomers (PHAH) over Cu and Ru catalysts and to fully-hydrogenated HAH monomers (FHAH) over the Ru catalyst. Pd based catalysts yielded a mixture of partially and fully hydrogenated monomers. Reaction kinetics models were employed to quantify the kinetic behavior for hydrogenation over Ru, Cu, and Pd catalysts. A 5-step pathway exhibited over Pd and Ru catalysts consists of both series and parallel reaction steps, where HAH is both converted to fully hydrogenated products sequentially via series reactions of partially hydrogenated intermediates, as well as converted directly in parallel reactions to form the fully hydrogenated products. In contrast, a 3-step pathway over the Cu catalyst consists only of the consecutive reaction steps, where the final product was formed via series reactions of intermediate products. Additionally, reaction over the Cu catalyst did not hydrogenate the furan rings of the HAH molecule and yielded a different final product than those hydrogenation over Pd and Ru catalysts. Using the results from our detailed reaction kinetics studies, we developed universal conditions for the maximum production of the various hydrogenated products for both batch and plug flow systems. In subsequent work, we employed reaction schemes to optimize the yields of products from selective hydrogenations of HAH in isopropanol solvent and in the presence of liquid water. Reaction schemes consisting of 7, 9, and 11 steps were examined to describe the rates of formation of the observed products and reaction intermediates for hydrogenation of HAH over Ru and Pd catalysts, and a 3-step scheme was studied over Cu catalysts. Rate constants and activation energies were calculated using these reaction schemes, and we then apply these reaction schemes to explore the effects of water addition on the hydrogenation pathways. The effects of water addition to isopropanol (IPA) solvents on the hydrogenation of HAH were markedly different over Pd, Ru, and Cu catalysts. Over the Pd catalyst, the addition of water to IPA increased hydrogenation rates and promoted hydrogenation of furan rings. The addition of water to IPA yielded significant carbon losses over the Ru catalyst, and slowed hydrogenation steps over Cu, while significantly inhibiting hydrogenation of the ketone group. This behavior opened routes toward increased production rates of PHAH=O (a partially hydrogenated form of HAH containing a C=O bond), a product in which the diene groups of the furan rings were not hydrogenated. Addition of water also allowed increased feed concentrations of HAH that were previously not possible in pure IPA solvents.

09 BIOMASS FUELS↗
DOE OSTI2024

Aerogel-derived nickel-iron oxide catalysts for oxygen evolution reaction in alkaline media

Anion exchange membrane water electrolyzers (AEMWEs) can generate hydrogen with a pure water feed using noble metal-free catalysts. Here, the development of highly active and stable catalysts for oxygen evolution reaction (OER) is required for improving performance of AEMWEs systems. Ni-Fe (oxy)hydroxides show high OER catalytic activity in alkaline media, but typically have low surface area. In this work, we investigate a series of Ni-Fe oxides with high surface area and disordered morphology, obtained using an aerogel synthesis method. We evaluate the impact of different synthesis variables on the OER activity and demonstrate that heat treatment at high temperatures generates more ordered structure, resulting in a decrease in OER activity. Advanced characterization reveals that maintaining highly disordered and porous structure of the aerogel is essential to achieving high OER activity, as it enables the formation of highly OER-active lamellar structures of the catalyst.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Three influential factors on colloidal nanoparticle deposition for heat conduction enhancement in 3D chip stacks

Thermal management is one of the major challenges facing the development of three-dimensional (3D) chip stacks. Recently, experimental studies have shown that neck-based thermal structure (NTS) between chip layers formed by drying of colloidal suspension in cavity filled with micro-size particles can improve the vertical heat conduction threefold. However, a deep understanding of the mechanisms of neck formation and its influence on heat conduction is still lacking. In this paper, we numerically study the effects of three parameters, i.e., initial nanoparticle concentration, drying temperature and chip surface wettability on neck formation between filler particles and on the resulting heat conduction of the NTS. With increasing nanoparticle concentration, the size and number of necks increase, resulting in an increased effective thermal conductivity (ETC) of NTS. The drying temperature is found to have only little influence on the ETC of resultant NTS, while the neck size and spatial distribution become more uniform at higher drying temperature. When reducing the wettability of the top and bottom surfaces of the cavity, the necks shrink in size until completing evacuating at the top and bottom layers, while the size of the necks between filler particles in the middle height of the cavity expands slowly. In consequence, the ETC of NTS drops at an increasing rate. Being able to reveal the underlying multiple mechanisms of two-phase flow, phase change and heat transport, the current numerical study suggests optimal values for the deposition process, with initial nanoparticle concentration over 0.8%, a drying temperature of 60°C and a uniform contact angle of 30° for practical production of NTS.

3D chip stacks↗
DOE OSTI2021

Physics-Informed Machine Learning Models for Predicting the Progress of Reactive-Mixing

This paper presents a physics-informed machine learning (ML) framework to construct reduced-order models (ROMs) for reactive-transport quantities of interest (QoIs) based on high-fidelity numerical simu-lations. QoIs include species decay, product yield, and degree of mixing. The ROMs for QoIs are applied to quantify and understand how the chemical species evolve over time. First, high-resolution datasets for constructing ROMs are generated by solving anisotropic reaction-di?usion equations using a non-negative finite element formulation for di?erent input parameters. The reactive-mixing model input parameters are: time-scale associated with flipping of velocity, spatial-scale controlling small/large vortex structures of velocity, perturbation parameter of the vortex-based velocity, anisotropic dispersion strength/contrast, and molecular diffusion. Second, random forests, F-test, and mutual information criterion are used to evaluate the importance of model inputs/features with respect to QoIs. We observed that anisotropic dispersion strength/contrast is the most important feature and time-scale associated with flipping of velocity is the least important feature. Third, Support Vector Machines (SVM) and Support Vector Regression (SVR) are used to construct ROMs based on the model inputs. The constructed SVR-ROMs are then used to predict scaling of QoIs. We also present estimates and inequalities on the QoIs, which inform that the species decay, mix, and produce in an exponential fashion. These inequalities also inform that a radial basis function is the most suitable kernel for the SVM/SVR models for QoIs. It is observed that R2-score for SVR-ROMs on unseen data is greater than 0.9, implying that the SVR-ROMs are able to predict the reaction-diffusion system state reasonably well. Finally, in terms of the computational cost, the proposed SVM-ROMs are O(107) times faster than running a high-fidelity finite element simulation for evaluating QoIs. This makes the proposed ML-based ROMs attractive for reactive-transport sensing and real-time monitoring applications as they are significantly faster yet reasonably accurate.

Mudunuru, Maruti K.↗
DOE OSTI2020

Modelling uranium dioxide corrosion under repository conditions: A pore-scale study of the chemical and thermal processes

This paper investigates the corrosion behaviors of spent uranium dioxide (UO 2 ) fuel when placed in geological repositories. We performed pore-scale reactive transport simulations of the UO 2 corrosion process in a defective fuel rod with different orientations of breach or fracture on its clad. It is found that the corrosion rate has strong dependency on the breach orientation. The highest corrosion rate of UO 2 is calculated when the angle between the groundwater flow direction and the clad fracture reaches 180°. UO 2 corrosion simulations with different flow rates are compared. The results show the higher flow rate accelerates the corrosion of UO 2 fuel. The effect of pH on the corrosion process is also determined. It is indicated that higher pH can limit the corrosion of UO 2 by reducing the reaction rate. The dependence of reaction rate and time on reactive surface area is explored. Spent fuel with lower surface area demonstrates longer lifetime under corrosion conditions. A numerical model coupling thermal conduction and chemical reactions is developed to assess the impact of temperature on the process of UO 2 corrosion. The results show that higher temperature leads to larger corrosion rates for UO 2 . Here, the predicted reaction rates are higher in comparison with the isothermal results. Hence, our results can help improve the fundamental understanding of UO 2 corrosion in geological repositories for long-term storage of spent nuclear fuels and provide guidance for the safe operations and selection of appropriate repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Predicting lake surface water phosphorus dynamics using process-guided machine learning

Phosphorus (P) loading to lakes is degrading the quality and usability of water globally. Accurate predictions of lake P dynamics are needed to understand whole-ecosystem P budgets, as well as the consequences of changing lake P concentrations for water quality. However, complex biophysical processes within lakes, along with limited observational data, challenge our capacity to reproduce short-term lake dynamics needed for water quality predictions, as well as long-term dynamics needed to understand broad scale controls over lake P. Here we use an emerging paradigm in modeling, process-guided machine learning (PGML), to produce a phosphorus budget for Lake Mendota (Wisconsin, USA) and to accurately predict epilimnetic phosphorus over a time range of days to decades. In our implementation of PGML, which we term a Process-Guided Recurrent Neural Network (PGRNN), we combine a process-based model for lake P with a recurrent neural network, and then constrain the predictions with ecological principles. We test independently the process-based model, the recurrent neural network, and the PGRNN to evaluate the overall approach. The process-based model accounted for most of the observed pattern in lake P; however it missed the long-term trend in lake P and had the worst performance in predicting winter and summer P in surface waters. The root mean square error (RMSE) for the process-based model, the recurrent neural network, and the PGRNN was 33.0 μg P L -1 , 22.7 μg P L -1 , and 20.7 μg P L -1 , respectively. All models performed better during summer, with RMSE values for the three models (same order) equal to 14.3 μg P L -1 , 10.9 μg P L -1 , and 10.7 μg P L -1 . Although the PGRNN had only marginally better RMSE during summer, it had lower bias and reproduced long-term decreases in lake P missed by the other two models. For all seasons and all years, the recurrent neural network had better predictions than process alone, with root mean square error (RMSE) of 23.8 μg P L -1 and 28.0 μg P L -1 , respectively. The output of PGRNN indicated that new processes related to water temperature, thermal stratification, and long term changes in external loads are needed to improve the process model. By using ecological knowledge, as well as the information content of complex data, PGML shows promise as a technique for accurate prediction in messy, real-world ecological dynamics, while providing valuable information that can improve our understanding of process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A dynamic solvent chamber propagation estimation framework using RNN for warm solvent injection in heterogeneous reservoirs

Warm solvent injection (WSI), injecting low-temperature solvent into formations to reduce the viscosity of heavy oil, is a clean technology for heavy oil production through reducing greenhouse gas emissions and water usage. The success of WSI operation depends on the uniform development and propagation of solvent chambers in reservoirs. However, reservoir heterogeneity stemming from shale barriers plays a detrimental role in the conformance of solvent chamber development and oil production rate. In this work, we developed a novel recurrent neural network (RNN)-based framework with the capability of efficiently tracking and estimating the solvent chamber positions in heterogeneous reservoirs based on only production time-series data. The developed estimation model utilizes the “sequence-to-sequence" mapping methodology to correlate observed production time-series sequence and solvent chamber edge sequence via a long short-term memory (LSTM) algorithm. The trained RNN models exhibit high accuracy, evidenced by the predicted dynamic solvent chamber locations match the corresponding true locations from numerical simulation, with a high coefficient of determination (R 2 ) and a low mean squared error. Specifically, the achieved R 2 values exceed 0.98 on both the training and testing data. The developed RNN-based workflow was tested via several cases from both regularly- and irregularly-shaped shale barriers, and the results were promising. The predicted solvent chambers showed strong agreement with those obtained from numerical simulations. The major benefits of this workflow include reducing computational time and saving overall monitoring and tracking costs for conventional techniques. In conclusion, the present work would provide a good demonstration of the capability of practical integration of machine learning methods in solving engineering problems.

58 GEOSCIENCES↗
DOE OSTI2021

Machine learning to discover mineral trapping signatures due to CO 2 injection

Mineral trapping is pursued as a geological CO 2 sequestration (GCS) mechanism because it permanently stores CO 2 in solid phases or minerals. However, CO 2 mineral-trapping mechanisms are poorly understood due to (1) lack of sufficient field and laboratory data characterizing these complex processes, and (2) challenges to develop site-specific reactive-transport models coupling fluid flow and geochemical reactions occurring at various temporal (from milliseconds to years) and spatial (from pore (millimeters) to field (kilometers)) scales. Reactive transport with additional complexities such as heterogeneity can make the simulation outputs even more difficult to interpret because of complex nonlinearity and multi-scale interdependencies. Furthermore, the values of model outputs such as concentrations can vary by several orders of magnitude, making it harder to correlate and characterize the impact of the variables via traditional data interpretation techniques such as exploratory data analyses. Recently, machine learning (ML) has shown promise in feature discovery and in highlighting hidden mechanisms that cannot be obtained by existing data-analytics and statistical methods. In this study, we applied an unsupervised ML approach, non-negative matrix factorization with custom -means clustering (NMF) to the data generated by reactive-transport simulations of GCS. The reactive-transport data consisted of 19 attributes, including four physio-chemical variables (pH, porosity, aqueous CO 2 , and sequestered CO 2 ), six chemical species (K + , Na + , HCO, Ca 2+ , Mg 2+ , Fe 2+ ), and four carbonate minerals (calcite, dolomite, siderite, and ankerite), a feldspar mineral (albite), and four clay minerals (illite, clinochlore, kaolinite, and smectite) over a period of 200 years of simulation time. Furthermore, the simulation data used was for Morrow B sandstone at the Farnsworth hydrocarbon unit in Texas. Data are sampled at two locations within the model domain: (1) at the injection well and (2) 200 m west of the injection well. The injection was performed for a period of 10 years. Using NMF, we estimated the temporal interdependencies among the 19 attributes over a span of 200 years. We found that NMF was able to identify four reaction stages and their dominant attributes; these cannot be directly discerned through traditional visualization (e.g., line plots, Pareto analysis, Glyph-based visualization methods) or exploratory data analysis tools of the simulation data. The four stages were: reactions in the injection phase followed by short-, mid-, and long-term reactions. The NMF analysis also revealed that 10 among the 19 attributes are dominant. These dominant attributes for mineral trapping include calcite, dolomite at injection well, siderite at 200 m away from the injection well, clinochlore, kaolinite, Na + , K + , Ca 2+ , Mg 2+ , pH, and aqeuous CO 2 . Finally, at late times (65–200 years), our results showed that calcite plays a major role in mineral trapping with insignificant contribution from siderite, ankerite, and clay minerals. These findings make the proposed unsupervised ML-model attractive for reactive-transport sensing towards real-time GCS monitoring.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Prediction and uncertainty quantification of shale well performance using multifidelity Monte Carlo

Uncertainty quantification is an integral component of reservoir management, especially considering the inherent uncertainty in subsurface systems. While a standard practice to estimate the uncertainty, Monte Carlo (MC) simulation is computationally intense when the sampling population comprises high-fidelity simulations. Alternatively, the Multi-fidelity Monte Carlo (MFMC) simulation overcomes this computational intensity by integrating low- and high-fidelity simulations. Our goal is to minimize the number of expensive high-fidelity simulations while maintaining accuracy and using numerous fast and cheap low-fidelity simulations to efficiently sample to input parameter space of interest. We selected gas production from unconventional wells to demonstrate the potential speedups and accuracy of the MFMC approach. The model fidelity usually determines the trade-off between accuracy and efficiency. While the high-fidelity model is more accurate, the low-fidelity model is more efficient. Our high-fidelity simulation includes reservoir simulations of a hydraulically fractured well. On the other hand, our low-fidelity model comprises the parallel-plate flow model. We used differential programming to efficiently solve the 1D flow model, where automatic differentiation is used to efficiently compute the gradients. We matched the production profile of high-fidelity simulations with our low-fidelity simulations. Then, we used a support vector regression to map the high- and low-fidelity input parameters. The mapping function is essential to tune the low-dimensional parameter space of the low-fidelity model to the high-dimensional parameter space of the high-fidelity model. We found that we can use a combination of 9 high fidelity and 10,000 low fidelity simulations to efficiently and accurately simulate pressure management. This method is at least two orders of magnitude faster than only using high-fidelity simulations. Finally, from a broader perspective, MFMC could efficiently estimate the uncertainty of various systems and models, integrating low- and high-fidelity models.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

A robust deep learning workflow to predict multiphase flow behavior during geological C O 2 sequestration injection and Post-Injection periods

Simulation of multiphase flow in porous media is essential to manage the geologic CO 2 sequestration (GCS) process, and physics-based simulation approaches usually take prohibitively high computational cost due to the nonlinearity of the coupled physics. This paper contributes to the development and evaluation of a deep learning workflow that accurately and efficiently predicts the temporal-spatial evolution of pressure and CO 2 plumes during injection and post-injection periods of GCS operations. Based on a Fourier Neural Operator, the deep learning workflow takes input variables or features including rock properties, well operational controls and time steps, and predicts the state variables of pressure and CO 2 saturation. To further improve the predictive fidelity, separate deep learning models are trained for CO 2 injection and post-injection periods due to the difference in primary driving force of fluid flow and transport during these two phases. We also explore different combinations of features to predict the state variables. We use a realistic example of CO 2 injection and storage in a 3D heterogeneous saline aquifer, and apply the deep learning workflow that is trained from physics-based simulation data and emulate the physics process. Through this numerical experiment, we demonstrate that using two separate deep learning models to distinguish post-injection from injection period generates the most accurate prediction of pressure, and a single deep learning model of the whole GCS process including the cumulative injection volume of CO 2 as a deep learning feature, leads to the most accurate prediction of CO 2 saturation. For the post-injection period, it is key to use cumulative CO 2 injection volume to inform the deep learning models about the total carbon storage when predicting either pressure or saturation. The deep learning workflow not only provides high predictive fidelity across temporal and spatial scales, but also offers a speedup of 250 times compared to full physics reservoir simulation, and thus will be a significant predictive tool for engineers to manage the long-term process of GCS.

58 GEOSCIENCES↗
DOE OSTI2025

rcsb-api : Python Toolkit for Streamlining Access to RCSB Protein Data Bank APIs

The Protein Data Bank (PDB) was founded in 1971 as the first open-access digital data resource in biology to serve as the single global archive for three-dimensional (3D) macromolecular structure data. Current PDB holdings exceed 230,000 experimentally determined structures of proteins, nucleic acids, viruses, and macromolecular machines. The RCSB Protein Data Bank RCSB.org research-focused web portal facilitates search, analyses, and visualization of every PDB structure along with more than one million Computed Structure Models from AlphaFold DB and the ModelArchive. It is powered by a set of publicly available Application Programming Interfaces (APIs) that both support RCSB.org users and provide programmatic access to PDB data. Given the breadth and levels of granularity encompassed in this rich data collection, efficiently accessing the information programmatically may be challenging for new users. RCSB PDB has developed a Python software package, rcsb-api , that facilitates easy and efficient use of RCSB PDB APIs within a Python environment. This software tool is designed to streamline access to the extensive corpus of data housed within the PDB, enabling researchers to search, retrieve, and analyze 3D biostructure data seamlessly. Its use will accelerate research in structural biology, molecular biology and biochemistry, drug discovery, and bioinformatics by providing more efficient tools for data integration and analysis. The new toolkit is available on GitHub (github.com/rcsb/py-rcsb-api) and published to the public Python package repository (PyPI) to foster wider usage and support basic and applied research in fundamental biology, biomedicine, and the energy sciences.

FAIR principles↗
DOE OSTI2021

Machine-learning predictions of the shale wells’ performance

The ultra-low permeability nature of shale reservoirs leads to an extended linear flow and necessitates horizontal wells with multi-stage engineered fractures to efficiently extract hydrocarbons resources. These artificially-generated and naturally-occurring fractures form complex networks that create complex flow regimes which control oil production. These fractures are neither identical nor equally-spaced, which leads to a production profile with a masked onset of the boundary-dominated flow. The combination of the extended linear flow with the indeterminate onset of the boundary-dominated flow challenges the current deterministic analytic approaches to forecast the estimated ultimate recovery (EUR). In this work, we propose a novel machine-learning approach which overcomes these challenges and provides reliable EUR estimates based on field-wide analyses. We implement a novel unsupervised machine learning (ML) methodology, which allows for automatic identification of the optimal number of features (signals) present in the data based on non-negative matrix/tensor factorization coupled with k-means clustering incorporating regularization and physics constraints. In the presented analyses, the input data to the ML algorithm is the available (public) production history from the field collected at existing unconventional reservoirs. We validate our approach through hindcasting of the production data, where we achieved an excellent agreement. In addition, our approach is able to identify the poorly-performing wells, which could benefit from early refracing. Our approach provides fast and accurate estimations of the well performance without presumptions about the state of the well or the flow regime.

03 NATURAL GAS↗