Search NASA⌕ Search

SEARCH · Search NASA

Results for “Normalizing flows”

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

Kinematic Flow and the Emergence of Time

Perhaps the most basic question we can ask about cosmological correlations is how their strength changes as we smoothly vary kinematic parameters. The answer is encoded in differential equations that govern this evolution in kinematic space. In this Letter, we introduce a new perspective on these differential equations. We show that, in the simplified setting of conformally coupled scalars in power-law Friedmann-Robertson-Walker spacetimes, the equations for arbitrary tree-level processes can be obtained from a small number of simple combinatorial rules. While this “kinematic flow” is defined purely in terms of boundary data, it reflects the physics of bulk time evolution. The unexpected regularity of the equations suggests the existence of an autonomously defined mathematical structure from which cosmological correlations and the time evolution of the associated spacetime emerge.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On Alfvénic turbulence of solar wind streams observed by Solar Orbiter during March 2022 perihelion and their source regions

It has been recently accepted that the standard classification of the solar wind solely according to flow speed is outdated, and particular interest has been devoted to the study of the origin and evolution of so-called Alfvénic slow solar wind streams and to what extent such streams resemble or differ from fast wind. In March 2022, Solar Orbiter completed its first nominal phase perihelion passage. During this interval, it observed several Alfvénic streams, allowing for characterization of fluctuations in three slow wind intervals (AS1-AS3) and comparison with a fast wind stream (F) at almost the same heliocentric distance. This work makes use of Solar Orbiter plasma parameters from the Solar Wind Analyzer (SWA) and magnetic field measurements from the magnetometer (MAG). The magnetic connectivity to the solar sources of selected solar wind intervals was reconstructed using a ballistic extrapolation based on measured solar wind speed down to the (spherical) source surface at 2.5 R s below which a potential field extrapolation was used to map back to the Sun. The source regions were identified using SDO/AIA observations. A spectral analysis of in situ measured magnetic field and velocity fluctuations was performed to characterize correlations, Alfvénicity, normalized cross-helicity, and residual energy in the frequency domain as well as intermittency of the fluctuations and spectral energy transfer rate estimated via mixed third-order moments. A machine learning technique was used to separate proton core, proton beam, and alpha particles and to study v − b correlations for the different ion populations in order to evaluate the role played by each population in determining the Alfvénic content of solar wind fluctuations. The comparison between fast wind and Alfvénic slow wind intervals highlights the differences between the two solar wind regimes: The fast wind is characterized by larger amplitude fluctuations, and magnetic and velocity fluctuations are closer to equipartition of energy. In fact the Alfvénic slow wind streams appear to be on a spectrum of wind types, with AS1, originating from open field lines neighboring active regions and displaying similarities with the fast wind in terms of fluctuation amplitude and turbulence characteristics, but not with respect to the alpha particles and proton beams. The other two slow streams differed both in their sources as well as plasma characteristics, with AS2 coming from the expansion of a narrow coronal hole corridor and AS3 from a region straddling a pseudostreamer. The latter displayed the coldest and highest density but the slowest stream with the smallest fluctuation amplitude and greatest magnetic energy excess. It also showed the largest scatter in proton beam speeds and the greatest difference in speed between proton beam and alpha particles. This study shows how the old fast–slow solar wind dichotomy, already called into question by the observations of slower Alfvénic solar wind streams, should further be refined, as the Alfvénic slow wind, originating in different solar wind regions, show significant differences in density, temperature, and proton and alpha-particle properties in the inner heliosphere. The observations presented here provide the starting point for a better understanding of the origin and evolution of different solar wind streams as well as the evolving turbulence contained within.

magnetohydrodynamics (MHD)↗
DOE OSTI2026

Operational Parameter Database for Molten Salt Thermal Energy Storage Tank Modeling

The second generation of concentrated solar power (CSP) plants is characterized by the use of a central receiver (either cavity or external), two molten nitrate salt tanks (60 wt.% NaNO 3 and 40 wt.% KNO 3 ), and a steam Rankine power-generation cycle connected through a primary heat exchanger. Molten salt thermal energy storage (TES) tanks have been widely deployed in commercial CSP plants worldwide and have been essential for increasing plant dispatchability and capacity factor, while also reducing the levelized cost of electricity (LCOE). These systems enable energy storage at the gigawatt-hour scale, typically providing 6 to 17 hours of storage duration. Despite being a commercial technology, the multiple failures observed after only a few months or years of operation in plants around the world demonstrate the technology's relative infancy and highlight the need for further research to improve its reliability. The National Laboratory of the Rockies (NLR), in collaboration with industry partners and academic and research institutions, has been leading multiple projects funded by the U.S. Department of Energy (DOE). These projects focus on addressing molten salt tank failures by improving tank design and welding fabrication practices, evaluating new alloys and weld fillers, and providing guidelines for tank commissioning and safe operation. In particular, this report presents modeling results on the effect of key tank operation parameters during 60 minutes of operation, including the mass flow rate and temperature of the salt inflow, tank salt inventory temperature, and inventory level for a representative molten salt tank design. These results form a database of tank operation behaviors that captures the effects of each specific parameter during charging, charging/discharging, and discharging processes.

14 SOLAR ENERGY↗
DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗
DOE OSTI2025

Operation of helium sub-atmospheric multistage cryogenic centrifugal compressor trains: Part 2 – Transient modeling and pump-down path selection

Low-pressure conditions required for operation of helium cryogenic systems below the normal boiling point (i.e. 4.2 K) are established through a transient process, commonly referred to as ‘pump-down’. This process is defined as the transition from pressures above atmospheric conditions to the saturation pressure which corresponds to a specified operational temperature. The FRIB 2 K system consists of five cryogenic centrifugal compressors which are operated in series. Historically, the pump-down process path has been established through empirical methods and system operator experience. Investigation into the pump-down process at FRIB aimed to develop a pump-down methodology which relies on theoretical model predictions rather than empirically developed process paths. Ensuring stable operation during the pump-down process involved application of a centrifugal compressor performance prediction model, which is described in Part 1 of this paper. Compressor performance maps can be directly used to evaluate the stability of a selected pump-down path and anticipate the overall reliability of the selected path. In conjunction with the compressor performance maps, a system pressure model was developed to estimate the transient pressure response during the pump-down process. Lastly, an explicit equation was developed to establish a mass flow rate profile for the pump-down process. Implementation of the presented methodology (including the developed models) allows for the system operator to determine a continuous pump-down path which maintains compressor stability while conforming to overall system capabilities. Altogether, the methodology presented has resulted in simplification of transient pump-down operations and increased the reliability, stability and efficiency of the pump-down process.

Compressor train control↗
DOE OSTI2024

A Q-band frequency tunable Doppler backscattering (DBS) system for pedestal and scrape-off layer density fluctuation and flow measurements in the DIII-D tokamak

We present the design and laboratory tests for a new Q-band frequency tunable Doppler backscattering (DBS) system suitable for probing poloidal wavenumber kñ = 6–8 cm−1 density fluctuations and their flow velocities in the pedestal and scape-off layer (SOL) of the DIII-D tokamak. This system will provide new measurements in the increasingly important and under-diagnosed far pedestal and SOL plasma regions. These results are important for experimental transport studies and necessary for the validation of transport models, both of which are important to fusion energy research. The use of a single tunable frequency reduces the complexity and potential failure points as compared to a multichannel system. This new system utilizes a 33–50 GHz tunable source and will be integrated into the current V-band DBS in DIII-D using a broadband Q- and V-band multiplexer. A full-scale mockup of the quasi-optical system was used to test and optimize the performance. These tests include beam profile measurements at different distances (and angles) from a paraboloidal focusing and steering mirror. The measurements cover the full frequency range 33–75 GHz of the integrated/combined Q–V band DBS system and target a large radial coverage of the low-field side of the plasma from ρ = 1.1 to ρ = 0.5, where ρ is the normalized flux surface radial coordinate.

Instruments & Instrumentation↗
DOE OSTI2025

Differential equations for cosmological correlators

Cosmological fluctuations retain a memory of the physics that generated them in their spatial correlations. The strength of correlations varies smoothly as a function of external kinematics, which is encoded in differential equations satisfied by cosmological correlation functions. In this work, we provide a broader perspective on the origin and structure of these differential equations. As a concrete example, we study conformally coupled scalar fields in a power-law cosmology. The wavefunction coefficients in this model have integral representations, with the integrands being the product of the corresponding flat-space results and “twist factors” that depend on the cosmological evolution. Similar twisted integrals arise for loop amplitudes in dimensional regularization, and their recent study has led to the discovery of rich mathematical structures and powerful new tools for computing multi-loop Feynman integrals in quantum field theory. The integrals of interest in cosmology are also part of a finite-dimensional basis of master integrals, which satisfy a system of first-order differential equations. We develop a formalism to derive these differential equations for arbitrary tree graphs. The results can be represented in graphical form by associating the singularities of the differential equations with a set of graph tubings. Upon differentiation, these tubings grow in a local and predictive fashion. In fact, a few remarkably simple rules allow us to predict — by hand — the equations for all tree graphs. While the rules of this “kinematic flow” are defined purely in terms of data on the boundary of the spacetime, they reflect the physics of bulk time evolution. We also study the analogous structures in tr ϕ 3 theory, and see some glimpses of hidden structure in the sum over planar graphs. This suggests that there is an autonomous combinatorial or geometric construction from which cosmological correlations, and the associated spacetime, emerge.

Cosmological models↗
DOE OSTI2025

WHONDRS River Corridor Sediment and Water Geochemistry and In Situ Sensor Data from 7 Perennial and 7 Intermittent Streams across San Antonio, Texas (v3)

This dataset supports a broader study examining the effects of intermittency on sediment respiration. The dataset provides sediment and surface water geochemistry and in situ sensor data from 7 perennial and 7 intermittent streams in San Antonio, Texas. Each stream/site was visited both in summer during base flow (July-September 2023) and winter during peak flow (January-February 2024). Related data were collected and will be published separately in collaboration with A. Veach. The data package was originally published in April 2025. It was updated in June 2025 (v2; modified and new files) and September 2025 (v3; modified files). See the change history section in the readme for more details. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of two folders of field photos and videos, one folder of raw Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) data and one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) readme; (5) international generic sample number (IGSN) mapping file; (6) field protocol; (7) a subfolder with sample data; and (8) a subfolder with sensor data. The sample data subfolder contains (1) surface water and sediment dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data and averages; (2) surface water and sediment total nitrogen data and averages; (3) sediment grain size data; (4) sediment iron (II) data and averages; (5) wet sediment mass, dry sediment mass, water mass, and wet sediment volume in incubation and sediment ICR vials; (7) sediment incubation respiration rate data and averages; (8) normalized respiration rate data and averages; (9) methods codes; (10) sediment percent carbon and nitrogen; (11) sediment X-ray diffraction (XRD) data; (12) gravimetric moisture and averages; (13) a subfolder with sediment incubation respiration data, scripts, and plots; (14) surface water and sediment FTICR methods; and (15) a subfolder of 9.4 Tesla (9.4T) FTICR-MS data. This folder contains five subfolders, one containing the sediment .xml data files, one containing the water .xml files, one containing the sediment CoreMS output files, one containing the water CoreMS output files, and the other containing instructions and scripts for processing the files in CoreMS (https://github.com/EMSL-Computing/CoreMS). The sensor data subfolder contains (1) a subfolder with miniDOT dissolved oxygen and temperature data and plots; (2) miniDOT dissolved oxygen and temperature summary data; and (3) miniDOT installation methods. All files are .csv, .pdf, .R, .xml, .d, .html, .Rmd, .py, .cal, .json, .jpg, .jpeg, .png, .mov, or .mp4. CORRECTION: The data processing methods for FTICR described in “v3_WHONDRS_AV1_Methods_Codes.csv” mistakenly indicate that users should process the data in Formultitude. The corrected description should read: “Both unprocessed and processed data are provided to allow users flexibility in data processing. Instructions and scripts for processing the data using CoreMS are included.” CORRECTION: Carbon and nitrogen content are reported as percentages. The current column headers "01395_C_percent_per_mg" and "01397_N_percent_per_mg" are incorrect. These should read "01395_C_percent" and "01397_N_percent" and will be corrected in the next version of this data package.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of an x-ray multilayer mirror for radiography at sub-1-keV photon energies

An x-ray multilayer mirror on a spherical substrate designed for near-normal incidence with a photon energy of ~738 eV (F He α ) was procured and tested. This device is intended to be used for in-flight radiography of the shell in inertial confinement fusion experiments with cryogenic targets on the OMEGA laser at the Laboratory for Laser Energetics. Experiments in self-emission on a small (~10 J) laser system showed that the reflectivity of the mirror is high enough to record an image at laser energies as low as 0.1 J. Furthermore, a second set of tests in backlighting geometry on a larger (kJ)-scale, short-pulse laser yielded usable radiographs with laser energies as low as 40 J with a spatial resolution of ~10 μm.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Is There a Scalar Atmospheric Surface Layer Within a Convective Boundary Layer? Implications for Flux Measurements

Top‐down entrainment shapes the vertical gradients of sensible heat, latent heat, and CO 2 fluxes, influencing the interpretation of eddy covariance (EC) measurements in the unstable atmospheric surface layer (ASL). Using large eddy simulations for convective boundary layer flows, we demonstrate that decreased temperature gradients across the entrainment zone increase entrainment fluxes by enhancing the entrainment velocity, amplifying the asymmetry between top‐down and bottom‐up flux contributions. These changes alter scalar flux profiles, causing flux divergence or convergence and leading to the breakdown of the constant flux layer assumption (CFLA) in the ASL. As a result, EC‐measured fluxes either underestimate or overestimate “true” surface fluxes during divergence or convergence phases, contributing to energy balance non‐closure. The varying degrees of the CFLA breakdown are a fundamental cause for the non‐closure issue. These findings highlight the underappreciated role of entrainment in interpreting EC fluxes, addressing non‐closure, and understanding site‐to‐site variability in flux measurements.

eddy covariance fluxes of scalars↗
DOE OSTI2026

Investigation of Flux Spreading in a Light-Trapping, Planar-Cavity Receiver for Enclosed Solar Particle Heating

Concentrating solar thermal power (CSP) technology development has recently focused on increasing the operating temperatures to accommodate high efficiency power cycles and thermochemical processes. Inert solid particles as heat transfer media enable solar receivers to operate above 700 degrees Celsius resulting in increased system thermal efficiency compared to the conventional molten salt based CSP system. An open-cavity falling-particle solar receiver that can efficiently heat particles by direct heating from concentrated solar radiation faces challenges with large particle losses from wind and unable to support thermochemical reactions. A light-trapping, planar cavity reiver (LTPCR) where particles are indirectly heated can significantly minimize the particle losses during the operation, support thermochemical reactions, and offer scalability potential. The LTPCR features an array of vertical planar receiver/absorber panels arranged within a cavity configuration. Concentrated solar radiation from heliostats is focused onto the receiver walls, where heat is indirectly transferred to solid particles flowing inside the receiver channels. Heat transfer occurs through direct contact between the receiver panel walls and particles, and can be enhanced by fluidizing particles with air. This fluidization increases particle-wall contact and extends particle residence time, maximizing heat transfer efficiency. The unique vertical planar receiver structure originated from a near-blackbody tubular light absorber, effectively distributing the incoming solar beam spread across the panel walls and trapping light. This flux spreading effect, driven by cosine projection, converts high incident solar flux into a lower, more uniform heat flux on the panel walls. This redistribution enhances heat transfer efficiency between particle-wall or reaction gases-wall, while preventing localized overheating of the receiver panel. Indirect planar cavity solar receivers completely separate solid particles from the ambient environment that can greatly reduce the thermal losses in heated particles resulting in high efficiency at high temperatures above 700 degrees Celsius. This design ensures no particle losses to the environment during the operation while open-cavity designs can experience significant particle losses from wind. An experimental investigation was conducted to observe flux spreading on the receiver panel wall. A lab-scale prototype planar receiver, fabricated using Haynes 230 alloy, was tested under direct concentrated solar radiation using the high-flux solar furnace (HFSF) facility at NREL. The experiment was performed under normal peak radiative heat fluxes ranging from 800 to 1900 kW/m2. A temperature distribution on the panel wall was measured using a thermal imaging camera (FLIR A 6600). To prevent overheating at the receiver front tip, prism-shaped heat shields (Zircar UNIFROM C1) were placed in front of the receiver, and their influence on flux spreading was also studied. Absorbed flux distribution on the panel wall was modeled using SolTrace. The total solar power and flux distributions delivered from HFSF were determined based on the heliostat mirror optical properties, direct normal irradiance (DNI) on the on-sun testing days, peak flux measurement during the on-sun testing, and shutter/attenuator settings Due to the large incident angles of the solar beam on the panel wall, the angular optical properties of Haynes 230 alloy and Zircar heat shields were incorporated into the model. This flux distribution model was then integrated into a computational fluid dynamics (CFD) simulation to predict the receiver panel wall temperature, which was compared with the experimental measurements. Both prediction and measurements identified a temperature hotspot at the backside of the panel, indicating that the incident solar beam can fully reach to the rear of the receiver. The heat shields positioned at the front of the receiver effectively reduced the excessive temperature rise at the receiver front tip. Overall, the temperature was well distributed over the panel wall, with a minor hotspot at the back of the receiver. The model slightly overpredicted the temperature, possibly due to discrepancies in optical properties of the panel and an underprediction of thermal loss in the receiver. The advancement of the particle LTPCR offers a viable alternative to open-cavity receivers by addressing particle loss issues. Additionally, it presents a pathway for enabling solar thermochemical processes, extending CSP technology beyond power generation to fuel and chemical production.

14 SOLAR ENERGY↗
DOE OSTI2024

Feasibility of an experiment on clumping induced by the Crow instability along a shocked cylinder

The growth of three-dimensional perturbations subject to the Crow instability along a vortex dipole resulting from the passage of a shock wave through a heavy gaseous cylinder is examined numerically. A linear stability analysis is performed based on geometric parameters extracted from two-dimensional simulations to determine the range of unstable wavenumbers, which is found to extend from 0.0 to 1.3 when normalized by the core separation distance. The analysis is then verified by comparison to three-dimensional simulations, which clearly show the development of the instability and the pinch-off of the vortex dipole into isolated vortex rings, which manifest as clumps of the original cylinder material. A scaling law is developed to determine the relevant spatiotemporal scales of the instability development, which is then used to assess the feasibility of a high-energy-density experiment visualizing clump formation. Specifically, a shocked cylinder with an initial diameter of 100 μm consisting of a perturbation of approximate wavelength and amplitude of 600 and 10 μm, respectively, is expected to form clumps resulting from the Crow instability approximately 40 ns after it is shocked, with dynamics which can be readily visualized on the Omega EP laser facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Higher-order symmetry plane correlations in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The correlations between event-by-event fluctuations of symmetry planes are measured in Pb-Pb collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV recorded by the ALICE detector at the Large Hadron Collider. This analysis is conducted using the Gaussian estimator technique, which is insensitive to biases from correlations between different flow amplitudes. The study presents, for the first time, the centrality dependence of correlations involving up to five different symmetry planes. The correlation strength varies depending on the harmonic order of the symmetry plane and the collision centrality. Comparisons with measurements from lower energies indicate no significant differences within uncertainties. Additionally, the results are compared with hydrodynamic model calculations. Although the model predictions provide a qualitative explanation of the experimental results, they overestimate the data for some observables. This is particularly true for correlators that are sensitive to the nonlinear response of the medium to initial-state anisotropies in the collision system. As these new correlators provide unique information—independent of flow amplitudes—their usage in future model developments can further constrain the properties of the strongly interacting matter created in ultrarelativistic heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Anisotropic jet broadening and jet shape

In this paper, we explore the use of jet substructure as a way of probing phenomena which break the isotropic behavior of jets, such as jet propagation through an anisotropically flowing quark-gluon plasma or spin correlations. We introduce two novel observables for this purpose: the azimuthal-dependent jet broadening and the azimuthal-dependent jet shape, which generalize the traditional isotropic substructure studies. Using Soft-Collinear Effective Theory, we explicitly calculate the jet functions associated with these observables with a standard jet axis and with a Winner-Take-All jet axis in both the resummed and fixed order limits. While our analysis first and foremost establishes the formalism for the azimuthal-dependent jet substructure, it also brings to light new results for jet substructure in the azimuthally integrated case, such as the semi-inclusive jet function and the exclusive jet shape for the Winner-Take-All axis, and the jet broadening in the fixed order region. As an illustrative example for the new formalism we demonstrate that the azimuthal-dependent jet broadening can be used as a direct probe of the transversity parton distribution function in deep inelastic scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

High-Voltage DC MULTI-terminal SIMulation (HVDC MULTISIM): Technical Program Summary

Expansion of the power grid in the USA is essential to fulfilling the rising energy demand of the country. In particular, the transmission grid, a network of electrical energy corridors that enable the flow of large amounts of energy from the point of generation to the point of consumption, needs significant expansion to cater to this growth of electrification. The current infrastructure is also dated, and any new transmission corridor should be based on a vision of building a modern infrastructure that is futureproof. High voltage DC (HVDC) transmission technology falls in this category: it is the most economical way to build large transmission lines and relies on sophisticated electronics and flexible controls, rather than just passive components like switchgear and transformers. This project is aimed at establishing the technical and economic feasibility of a network of HVDC lines that are interconnected to form a multi-terminal DC network. HVDC transmission is not a new technology and United States has a number of such lines; however, these are point to point transmissions and do not form a DC grid. Modeling and analyzing DC grids is challenging since there are no established modeling tools and the traditional ways of modeling AC grids fall short of providing the fidelity required the fast dynamics of the DC grid. This project aims to build a software-in-the-loop simulator for a multi-terminal HVDC grid, called MULTISIM and, establish new control and protection algorithms to operate such a connected DC system. NLR, a key partner, established the core framework of the MULTISIM simulator and validated its functionality. The GE Vernova team, created a full-scale model of a four-terminal HVDC network using PSCAD (Power Systems Computer Aided Design) software and established the baseline performance of the system during normal and fault operation. The project was started on 10/12024 with a kick-off meeting held on 12/11/2024 and ran through two quarters till termination. All the tasks, milestones and deliverables during this period were met. Several interim reports detailing the various tasks were submitted. The following sections provide details on the completed tasks till the project pause and termination after Q2.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Origin of the ring ellipticity in the black hole images of M87*

We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT)'s 2018 array, including improved (u,v)-coverage from the Greenland Telescope, we measured the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of τ = 0.08−0.02+0.03 with a position angle of ξ = 50.1−7.6+6.2 degrees. To interpret this measurement, we compared it to general relativistic magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including the thermal or nonthermal electron distribution function, spins, and ion-to-electron temperature ratios in both low- and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in nonthermal models and models with higher jet emission. These results indicate that the ellipticity measured from the M87* emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects.Key words: accretion, accretion disks / black hole physics / gravitation / galaxies: active** These authors have contributed equally to this work.*** NASA Hubble Fellowship Program, Einstein Fellow.† Deceased.

Dahale, Rohan [Granada U., Theor. Phys. Astrophys.↗
DOE OSTI2025

Global River Topology (GRIT): A Bifurcating River Hydrography

Existing global river networks underpin a wide range of hydrological applications but do not represent channels with divergent river flows (bifurcations, multi‐threaded channels, canals), as these features defy the convergent flow assumption that elevation‐derived networks (e.g., HydroSHEDS, MERIT Hydro) are based on. Yet, bifurcations are important features of the global river drainage system, especially on large floodplains and river deltas, and are also often found in densely populated regions. Here we developed the first raster and vector‐based Global RIver Topology that not only represents the tributaries of the global drainage network but also the distributaries, including multi‐threaded rivers, canals and deltas. We achieve this by merging a 30 m Landsat‐based river mask with elevation‐generated streams to ensure a homogeneous drainage density outside of the river mask for rivers narrower than approximately 30 m. Crucially, we employ the new 30 m digital terrain model, FABDEM, based on TanDEM‐X, which shows greater accuracy over the traditionally used SRTM derivatives. After vectorization and pruning, directionality is assigned by a series of elevation, flow angle and continuity approaches. The new global network and its attributes are validated using gauging stations, comparison with existing networks, and randomized manual checks. The new network represents 19.6 million km of streams and rivers with drainage areas greater than 50 km 2 and includes 67,495 bifurcations. With the advent of hyper‐resolution modeling and artificial intelligence, GRIT is expected to greatly improve the accuracy of many river‐based applications such as flood forecasting, water availability and quality simulations, or riverine habitat mapping.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗