Search NASA⌕ Search

SEARCH · Search NASA

Results for “high energy density science”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

MX Anti-MXenes from Non-van der Waals Bulks for Electrochemical Applications: The Merit of Metallicity and Active Basal Plane

Two-dimensional transition-metal compounds (2DTMCs) are promising materials for electrochemical applications, but 2DTMCs with metallicity and active basal planes are rare. In this work, we proposed a simple and effective strategy to extract 2DTMCs from non-van der Waals bulk materials and established a material library of 79 2DTMCs, which we named as anti-MXenes since they are composed of one M atomic layer sandwiched by two X atomic layers. By means of density functional theory computations, 24 anti-MXenes were confirmed to be thermodynamically, dynamically, mechanically, and thermally stable. The metallicity and active basal plane endow these anti-MXenes with potential as excellent electrode materials, for example, as electrocatalysts for hydrogen evolution reactions (HER). Among the noble-metal free anti-MXenes with favorable H-binding, CuS can boost HER at the whole range of H coverages, while CoSi, FeB, CoB, and CoP show promise for HER at some specific H coverages. The active sites are the tetra-coordinating nonmetal atoms at the basal planes, thus rendering a very high density of active sites for these materials. CoB is also a promising anode material for lithium-ion batteries, showing low Li diffusion energy barriers, a very high capacity, and a suitable open circuit voltage. Overall, this work promotes the “computational exfoliation” of 2D materials from non-van der Waals bulks and exemplifies the applications of anti-MXenes in various electrochemical processes.

36 MATERIALS SCIENCE↗
NASA NTRS2015

Conformal Ablative Thermal Protection System for Small and Large Scale Missions: Approaching TRL 6 for Planetary and Human Exploration Missions and TRL 9 for Small Probe Missions

In 2011, NASAs Aeronautics Research Mission Directorate (ARMD) funded an effort to develop an ablative thermal protection system (TPS) material that would have improved properties when compared to Phenolic Impregnated Carbon Ablator (PICA) and AVCOAT. Their goal was a conformal material, processed with a flexible reinforcement that would result in similar or better thermal characteristics and higher strain-to-failure characteristics that would allow for easier integration on flight aeroshells than then-current rigid ablative TPS materials. In 2012, NASAs Space Technology Mission Directorate (STMD) began funding the maturation of the best formulation of the game changing conformal ablator, C-PICA. Progress has been reported at IPPW over the past three years, describing C-PICA with a density and recession rates similar to PICA, but with a higher strain-to-failure which allows for direct bonding and no gap fillers, and even more important, with thermal characteristics resulting in half the temperature rise of PICA. Overall, C-PICA should be able to replace PICA with a thinner, lighter weight, less complicated design. These characteristics should be particularly attractive for use as backshell TPS on high energy planetary entry vehicles. At the end of this year, the material should be ready for missions to consider including in their design, in fact, NASAs Science Mission Directorate (SMD) is considering incentivizing the use of C-PICA in the next Discovery Proposal call. This year both scale up of the material to large (1-m) sized pieces and the design and build of small probe heatshields for flight tests will be completed. NASA, with an industry partner, will build a 1-m long manufacturing demonstration unit (MDU) with a shape based on a mid LD lifting body. In addition, in an effort to fly as you test and test as you fly, NASA, with a second industry partner, will build a small probe to test in the Interactive Heating Facility (IHF) arc jet and, using nearly the same design, build the aeroshell and TPS, with instrumentation, for a small probe flight test article, due to fly in 2017. At the end of the year, the C-PICA will be at TRL 5+, and with the flight data in 2017, it will be at TRL 9 for missions needs with C-PICA at a small scale (12 diameter). The scale-up and small probe efforts will be de-scribed in this presentation.

Conformal Ablative↗
DOE OSTI2021

First-principles calculation of lattice distortions in four single phase high entropy alloys with experimental validation

Exceptional properties of high entropy alloys (HEAs) are attributed to the disordered random solid solution of multiple alloying elements. Despite numerous studies, the fundamental understanding at the electronic and atomic levels is still missing. We report a comparative study of four fcc HEAs of NiFeCoCr and NiFeCoCrX (X = Mn, Cu, or Pd) based on ab initio calculations using large supercells with 500 atoms in equal composition. After fully optimizing their structures using the VASP package, their electronic structure, interatomic bonding, partial charge distribution, and mechanical properties are calculated and compared, revealing the intricate interdependence among them. A novel parameter based on the quantum mechanical metric for internal cohesion, the total bond order density (TBOD), is used to interpret the calculated properties. The highest TBOD is found in Cantor alloy NiFeCoCrMn but lower in NiFeCoCrPd. The atomic radii vary depending on their local chemical environment. The resulting lattice distortions is validated experimentally in NiFeCoCrMn and NiFeCoCrPd. Moreover, modeling of Cu and Pd clustering in the supercell shows they have lower total energy in agreement with the observation of Cu-enhanced nano-participates in NiFeCoCrCu and Pd-induced concentration wave in NiFeCoCrPd HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

2D Ni 0.25 Mn 0.75 O 2 : A high-performance cathode for multivalent ion batteries

Although Li-ion batteries have driven portable energy storage in recent decades, there is increasing concern about their safety, cost, and abundance of constituents. Multivalent ion batteries (MVIBs) have the potential to remedy these issues, but they are limited by the currently known MVIB cathodes, which fail to deliver unanimously favorable voltage, energy density, and diffusion kinetics. Here we used density functional theory (DFT) to model the performance of Li, Na, Mg, Ca, and Al ions when paired with 2D Ni 0.25 Mn 0.75 O 2 , a novel cathode that uses increased layer separation to improve on the kinetics of its 3D analog. Our calculations yielded maximum voltages of 3.38 V for Na and 2.7 V for Ca, outperforming 2D Na x MnO 2 and Na x NiO 2 . Diffusion barriers for Li, Na, and Ca are below 300 meV, comparable to existing battery technology and the endpoint 2D cathodes; meanwhile, Mg and Al have prohibitively high diffusion barriers, implying their incompatibility with this cathode. Lastly, density of states calculations and Bader charge analysis show that the cathode becomes conducting following ion adsorption, which is necessary for high-rate performance. 2D Ni 0.25 Mn 0.75 O 2 maintains performance seen with other 2D transition metal oxides while increasing cathode conductivity, indicating that it is a promising candidate for experimental investigation with Li, Na, and Ca ions.

25 ENERGY STORAGE↗
DOE OSTI2021

Progress on continuum modeling of lithium–sulfur batteries

While lithium–sulfur batteries are a promising next-generation chemistry devices due to their high theoretical energy density, commercialization has been slow due to low coulombic efficiency and poor cycle life. This review explores the ways in which continuum modeling contributes to the understanding of lithium–sulfur (LiS) battery mechanisms and cell-level performance through the lens of micro- and macroscale phenomena. We examine different approaches to modeling important physical phenomena such as reaction mechanisms, cathode microstructure, shuttling, nucleation and precipitation, and transport limitations. This paper also emphasizes the significance and challenge of connecting typical modeling parameters and assumptions to systems-level metrics of a standard state-of-art high performing lithium–sulfur cell. Particularly important, the considerations for high energy density cells and the areas where continuum models can facilitate better collaboration are discussed. Here, we also summarize a few selected works to highlight experimentally-driven modeling, use of electroanalytical techniques, and parameter identification approaches to enable model-based design and advanced battery management systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Laser-driven flash x-ray radiography of a shocked metallic foil

Characterizing hydrodynamic instability evolution in millimeter-scale, high-Z foils is crucial for understanding complex phenomena in high-energy-density physics. Here, we demonstrate a proof-of-concept, laser-driven flash x-ray radiography platform tailored for two-dimensional linear density mapping in shocked high-Z foils. Using chromium (Cr) foils with internal shockwaves (∼100 μm width), our platform achieves a spatial resolution of 59.8 ± 1.4 μm by employing a broadband x-ray source extending into the hundreds of keV range. The setup combines a compound parabolic concentrator cone with a tantalum wire target, a magnetic field to deflect residual transmitted electrons, and a copper casing to shield the sides and rear of the image plate pack. By varying the delay of the short-pulse beam driving the flash x-ray source, we resolve shockwave dynamics, specifically the velocity, position, width, and density profile, within the Cr foil. Reported experimental results are consistent with the corresponding hydrodynamics and radiation transport simulations, which accurately reproduce the measured electron and x-ray source terms. These developments enable the conversion of shockwave radiographs into two-dimensional density maps, enhancing interpretability for hydrodynamic instability evolution applications and validating the simulation approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-coercivity ThMn 12 -type monocrystalline Sm–Zr–Fe–Co–Ti particles by high-temperature reduction diffusion

The ThMn 12 -type (Sm,Zr)1(Fe,Co,Ti) 12 compounds have the potential of powerful permanent magnets. Magnetically hard and anisotropic powders of such compound have been prepared by subjecting elemental oxides and Co mixed with Ca and CaO dispersant to a succession of high-energy ball-milling, reduction diffusion at 990–1220 °C and repeated washing. The size of the resulting ThMn 12 -type crystallites, their coercivity and fraction of monocrystalline particles were all found to increase with the reduction-diffusion temperature. Particles synthesized at 1220 °C were highly monocrystalline with a mean size of 0.54 μm and, after a magnetic-field alignment, exhibited a coercivity of 1.26 T and a full-density-projected maximum energy product of at least 209 kJ/m 3 (26.3 MGOe). Here, the strong positive effect of the reduction-diffusion temperature on the coercivity has been attributed to separation of the crystallites and to decrease in the incidence of structural defects.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The effect of microstructure on the dynamic shock response of 1045 steel

Herein, the effect of microstructure on the shock response of 1045 steel is investigated via plate impact experiments and postmortem characterization. Three unique microstructures are explored: ferrite-pearlite, martensite, and ferrite with spheroidal cementite (i.e. spheroidized). Two spall recovery experiments, at approximate peak pressures of 3.2 and 3.5 GPa, are conducted to assess the Hugoniot elastic limit (HEL), spall strength, and damage morphology of the various microstructures. The ferrite-pearlite and martensite microstructures exhibit continuous yielding at both quasi-static and dynamic rates, while the spheroidized condition displays discontinuous yielding. Discontinuous yielding of the spheroidized microstructure is attributed to a combined low initial dislocation density coupled with a low dislocation nucleation rate. The spall strength of ferrite-pearlite is consistently lower than the spheroidized microstructure, attributed to elongated cementite that is more susceptible to cracking than more spherical cementite precipitates. Despite a high density of boundaries, martensite exhibits the highest spall strength. A large percentage of the boundaries within the martensite microstructure are found to be low energy (i.e. Σ3 or low angle), and are thus less susceptible to spall damage. Overall, the high spall strength of martensite is likely linked to traditional strengthening mechanisms that limit dislocation motion.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Microstructure-sensitive modeling of high temperature creep in grade-91 alloy

Predicting the effects of microstructure on high-temperature creep responses of steel components is critical to minimize the risks of failure and maximize economic viability in the energy sector. In this research, a recently developed advanced mechanistic constitutive model is employed to study the effect of microstructure on the creep responses and to rationalize experimentally reported variability in the performances of grade-91 alloy. Hundreds of experimental creep tests from the literature are used to assess the predictability of the model. The model proposed is shown to accurately predict the steady-state creep rates and also the less commonly considered, yet important, primary-to-secondary transients for a wide range of creep conditions. Further, the model is capable of extrapolating the creep response of grade 91 alloy even outside of the calibration regime. Using this model, the roles of initial microstructure described in terms of grain size, dislocation density, and precipitate content, on the creep responses are investigated. The effect of initial microstructure on the steady-state creep rates and creep-strain is found to vary with stress and temperature. Grain size plays a significant role in the low-stress creep, where the diffusional creep is dominant. On the other hand, the effects of dislocation densities and precipitate content are significant in the dislocation-plasticity-dominated high-stress and temperature regimes. Furthermore, by comparing model predictions against a large experimental creep database, it is found that variability in the creep responses can be rationalized on the basis of differences in the initial microstructure. Overall, this work provides a robust pathway for microstructure-aware engineering design, verification, and validation of metallic components for energy applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Non-planar platinum group metal-free fuel cell cathodes for enhanced oxygen transport and water rejection

Proton exchange membrane fuel cells (PEMFC) with cathodes using platinum group metal-free (PGM-free) catalysts could significantly reduce costs, but the lower volumetric oxygen reduction reaction (ORR) activity requires thick electrodes that suffer from liquid water flooding and increased oxygen transport resistance. To address these challenges, we developed a 3D gas diffusion electrode (GDE) architecture to enhance liquid water removal through the diffusion media and reduce cathode saturation. The cathode features a uniform catalyst layer adjacent to the membrane for high ORR activity and then pillars of the hydrophilic catalyst layer that pass through the microporous layer (MPL), providing a low capillary pressure barrier pathway to the carbon fiber paper layer and channel. The non-planar cathode also increases the interfacial area between the catalyst layer and hydrophobic MPL for a greater fraction of the cathode with high O 2 concentration. Our studies included parametric experimental study of the pillar density to identify the optimum pitch between pillars. Our measurements show significant improvements in the mass transport region of the polarization curve with 3D structured electrodes leading to an 8% increase in maximum current density, 19% increase in maximum power density, and 16% increase in current density at 0.67 V with air.

25 ENERGY STORAGE↗
DOE OSTI2023

Detection of hot electrons originating from an upper valley at ∼1.7eV above the Γ valley in wurtzite GaN using electron emission spectroscopy

Using electron emission spectroscopy, measurement and analysis were conducted on the energy distribution of vacuum emitted electrons from electrically driven InGaN/GaN green (peak wavelengths 𝜆 ≈ 515nm) light-emitting diodes (LEDs) with and without a prewell superlattice (SL). Here, we report on the detection of a high-energy upper valley at ∼ 1.7eV above the Γ valley from samples with no prewell SL. We propose that these upper valley electrons originate predominantly from trap-assisted Auger recombination (TAAR) in green LEDs, as the intensity of these peaks is found to have quadratic dependence on the carrier density 𝑛 [see Espenlaub et al., J. Appl. Phys. 126, 184502 (2019)]. The high-energy upper valley peak was not observed in the sample with a prewell SL which is attributed to gettering by the prewell SL of still unidentified impurities that act as TAAR centers.

carrier generation & recombination↗
DOE OSTI2020

3D‐Printed Structure Boosts the Kinetics and Intrinsic Capacitance of Pseudocapacitive Graphene Aerogels

Abstract The performance of pseudocapacitive electrodes at fast charging rates are typically limited by the slow kinetics of Faradaic reactions and sluggish ion diffusion in the bulk structure. This is particularly problematic for thick electrodes and electrodes highly loaded with active materials. Here, a surface‐functionalized 3D‐printed graphene aerogel (SF‐3D GA) is presented that achieves not only a benchmark areal capacitance of 2195 mF cm −2 at a high current density of 100 mA cm −2 but also an ultrahigh intrinsic capacitance of 309.1 µF cm −2 even at a high mass loading of 12.8 mg cm −2 . Importantly, the kinetic analysis reveals that the capacitance of SF‐3D GA electrode is primarily (93.3%) contributed from fast kinetic processes. This is because the 3D‐printed electrode has an open structure that ensures excellent coverage of functional groups on carbon surface and facilitates the ion accessibility of these surface functional groups even at high current densities and large mass loading/electrode thickness. An asymmetric device assembled with SF‐3D GA as anode and 3D‐printed GA decorated with MnO 2 as cathode achieves a remarkable energy density of 0.65 mWh cm −2 at an ultrahigh power density of 164.5 mW cm −2 , outperforming carbon‐based supercapacitors operated at the same power density.

asymmetric supercapacitors↗
DOE OSTI2021

A first-principles study of hydrogen surface coverage on δ-Pu (100), (111), and (110) surfaces

Hydriding corrosion of plutonium leads to surface cracking, pitting, and ultimately structural failure. Laboratory experiments demonstrate that hydriding begins on the surface or near the subsurface of plutonium. However, there has not yet been a systematic evaluation of hydrogen surface coverage on plutonium. In this work, we compute the surface energies of the low facet surfaces of face-centered cubic δ-Pu. The adsorption free energies of expected hydrogen structures at low and high coverage are presented along with the likely progression for filling sites as the H 2 partial pressure increases. Implications for near-equilibrium pressure hydride nucleation and non-equilibrium millibar pressure hydriding are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Synthesis of high-density olivine LiFePO 4 from paleozoic siderite FeCO 3 and its electrochemical performance in lithium batteries

The lithium-ion cathode material olivine LiFePO 4 (LFP) has been synthesized for the first time from natural paleozoic iron carbonate (FeCO 3 ). The ferrous carbonate starting material consists of the mineral siderite at about 92 wt. % purity. Because FeCO 3 has divalent iron, the reaction with lithium dihydrogen phosphate (LiH 2 PO 4 ) provides a unique method to develop iron-(II) containing LFP in an inert atmosphere. Since siderite FeCO 3 is a common mineral that can be directly mined, it may, therefore, provide an inexpensive route for the production of LFP. After carbon-coating, the LFP yields a capacity in the range of 80–110 mAh g -1 LFP (in one chosen specimen sample), which is lower than commercially available LiFePO 4 (150–160 mAh g -1 LFP ). However, the tap density of LFP derived from siderite is noticeably high at 1.65 g cm -3 . The material is likely to be improved with powder purification, nanosized processing, and more complete carbon-coating coverage with increased optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗
DOE OSTI2026

The Maximal Entanglement Limit in Statistical and High-energy Physics

These lectures advocate the idea that quantum entanglement provides a unifying foundation for both statistical physics and high-energy interactions. I argue that, at sufficiently long times or high energies, most quantum systems approach a Maximal Entanglement Limit (MEL) in which phases of quantum states become unobservable, reduced density matrices acquire a thermal form, and probabilistic descriptions emerge without invoking ergodicity or classical randomness. Within this framework, the emergence of probabilistic parton model, thermalization in the break-up of confining strings and in high-energy collisions, and the universal small-x behavior of structure functions arise as direct consequences of entanglement and geometry of high-dimensional Hilbert space.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Plasma electron analysis: Voyager plasma science experiment

The Plasma Science Experiment (PLS) on the Voyager spacecraft provide data on the plasma ions and electrons in the interplanetary medium and the magnetospheres of the giant planets Jupiter and Saturn. A description of the analysis used to obtain electron parameters (density, temperature, etc.) from the plasma science experiment PLS electron measurements which cover the energy range from 10 eV to 5950 eV is presented. The electron sensor (D cup) and its transmission characteristics are described. A derivation of the fundamental analytical expression of the reduced distribution function F(e) is given. The electron distribution function F(e), used in the moment integrations, can be derived from F(e). Positive ions produce a correction current (ion feedthrough) to the measured electron current, which can be important to the measurements of the suprathermal electron component. In the case of Saturn, this correction current, which can either add to or subtract from the measured electron current, is less than 20% of the measured signal at all times. Comments about the corrections introduced by spacecraft charging to the Saturn encounter data, which can be important in regions of high density and shadow when the spacecraft can become negatively charged are introduced.

Sittler, E. C., Jr.↗
DOE OSTI2021

On the feasibility of using physics-informed machine learning for underground reservoir pressure management

In this work, we evaluate the feasibility of using physics-informed machine learning (PIML) for underground energy-related pressure management. To this end, we develop a PIML framework to manage underground reservoir pressures by training neural networks to determine fluid extraction rates for dedicated extraction wells during fluid injection operations given a range of reservoir conditions (e.g., transmissivity and storativity). We implement an automatically-differentiable analytical physics model of fluid flow in porous media within the PIML framework as a proxy for more complicated models. This allows us to execute a sufficient number of training scenarios to fully evaluate the feasibility of using PIML to support pressure management activities. We quantify the number of physics-model parameters required for automatic differentiation to become more efficient than finite-difference gradient calculations. We use a simple scenario with a single injector, extractor, and critical location for our feasibility analysis. We evaluate the effect of the size of the training dataset (i.e., the number of reservoir condition samples) on the accuracy and efficiency of the PIML framework. For an equivalent number of model evaluations, the larger training dataset took less time to train and produced a neural network that was able to more accurately manage reservoir pressures. We also evaluate the effect of the training dataset batch size (i.e., number of reservoir condition samples used to update the neural network coefficients during training; i.e., how the training dataset is partitioned). While training ran faster with larger batch sizes, they produced neural networks that managed pressures less accurately. We demonstrate the approach on a more complex scenario involving 10 injectors, 10 extractors, and 4 critical locations (a relatively high well density of 20 wells/km2). We provide the number of forward and adjoint model evaluations required in each case as an indication of the feasibility of using PIML for pressure management when more complicated physics models with longer execution times are used.

54 ENVIRONMENTAL SCIENCES↗