Search NASA⌕ Search

SEARCH · Search NASA

Results for “interactive data visualization”

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

UNH TDP - Concurrent Measurements of Inflow, Power Performance, and Loads for a Grid-Synchronized Vertical Axis Cross-Flow Turbine Operating in a Tidal Estuary

This data was collected between October 12 and December 15 of 2021 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). This data set includes over 29 days of grid connected turbine operation during this 65 day time frame. The priority for this measurement campaign was to collect data while the turbine was electrically connected to the grid by means of a rectifier and inverter. The Fall_2021_UNH_Measurement_Timeline.png highlights when each instrument was functioning and the Fall_2021_UNH_Test_Log.jpg indicates the four main regions for analysis available from this measurement campaign. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal forcing with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Kaelin Chancey (Assessment Of The Localized Flow And Tidal Energy Conversion System At An Estuarine Bridge - UNH MS Thesis 2019). The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. This project was pursued to quantify device performance for cross flow turbines operating in a marine environment. Accurate physical models, to characterize cross flow turbine performance, require real operational data sets due to the complexity of blade fluid interactions. This data can help support model development which will help predict turbine performance when analyzing perspective project locations in the future. Instrumentation was deployed to measure; water speed/direction, electrical power output, turbine shaft speed, turbine thrust force, and platform motion. Concurrent measurements of these parameters allow for correlations (cause and affect) to be inferred, allowing for characterization of device performance over a range of operating conditions. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) directly upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier and the Shark power meter. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. This submission includes zipped folders for each instrument containing quality controlled (QC'd) data in daily .csv files for the relevant duration specific to each instrument, along with separate .csv file that contains the units for each variable. Some instrument daily files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. Note the degree of QC varied between each instrument due to time constraints. Particular time and attention was given to perform quality control tests on the acoustic based instruments that are particularly suscep...

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Quantification of the Crack Evolution Process by Extracting Relevant Signal Components from Wave Propagation and Diffusive Transport Front Measurements

Wave propagation and diffusive transport phenomena in a geological rock sample undergoing crack evolution process are expected to interact with the mechanical discontinuities in the medium. The measurements of the signals associated with these phenomena can be used to assess and monitor the crack-driven micromechanical alterations in the rock. Different wave/diffusion phenomena, such as sonic propagation, pressure diffusion, and acoustic emission (AE), are sensitive to different elements of the mechanical discontinuities generated during the evolution of the crack clusters from initiation to coalescence. Sonic propagation, AE, and pressure diffusion monitoring have the potential to map the crack evolution because the transmitter-receiver arrays can be designed, arranged and tuned to (1) achieve maximum recovery of the scattered waveforms and travel times, (2) capture the later arrivals and multiple reflections, and (3) illuminate large rock volume. However, the structural/topological complexities of the mechanical discontinuities, complex distribution of the stress fields, complex mechanical alterations in media, and fluid redistribution in the crack system pose serious challenges for the detection and modeling of the crack evolution process (from here on, we will use the term ‘crack evolution process’ to mean that the crack evolution occurred under shallow crustal conditions). For purposes of accurately accounting such complexities and heterogeneities in the absence of reliable physical laws, simulation methods, and signal processing techniques, my early-career research proposal will develop and apply novel data-driven machine learning methods to: (1) extract signal components relevant to the various phases of crack evolution and (2) generate a 2D visual map of the crack evolution process.

58 GEOSCIENCES↗
DOE OSTI2026

From 2D to 4D: a containerized workflow and browser to explore dynamic chromatin architecture

Background Characterizing the physical organization of the genome is essential for understanding long-range gene regulation, chromatin compartmentalization, and epigenetic accessibility. Hi-C experiments generate two-dimensional (2D) genome-wide contact maps of chromatin interactions by capturing the spatial proximity between genomic loci, which reveal interaction frequencies but lack the spatial resolution needed to interpret the three-dimensional (3D) genome structure(s). Emerging evidence suggests that epigenetic regulation is closely linked to 3D genome architecture, and that structural changes over time (4D) drive key biological processes in development, disease, and environmental response. Thus, integrating 3D structure with functional data is critical for a more complete understanding of genome regulation. Previous work, most notably the 4DHiC chromosome modeling framework, has shown that physical multi-dimensional modeling approaches rooted in polymer physics and molecular dynamics can resolve these structures at biologically meaningful resolutions by integrating temporal Hi-C data with physical constraints to uncover dynamic chromosome reorganization. Thus, molecular dynamics simulations, constrained by Hi-C contact matrices, can resolve fine-scale structural changes and reveal functionally significant transitions in chromatin conformation. Results Herein, we present the 4D Genome Browser Workflow (4DGBWorkflow) and the 4D Genome Browser (4DGB). The algorithm is based on the 4DHiC method, and the containerized tool is an end-to-end workflow that can transform, filter, and view 4D epigenomics and chromatin datasets, allowing non-specialists to apply three-dimensional modeling principles to diverse datasets and experimental conditions. The software executes on a laptop running macOS, Linux or Windows. From input Hi-C files (.hic), the 4DGBWorkflow produces 3D reconstructions of chromosomes, integrates the reconstruction with track data (e.g., epigenetic marks, transcriptome profiles), and provides comparative visualization of the results in a single workflow. Conclusions The 4DGBWorkflow and 4D Genome Browser are open-source tools for comparative analysis and visualization of 4D chromosome datasets, including chromatin architecture and epigenomic signals. Automatic integration of Hi-C data with molecular dynamics democratizes the construction of time resolved 3D genome structures, simplifying complex simulations and data integration schemes.

3D Genome Browser↗
DOE OSTI2023

Post-starburst properties of post-merger galaxies

Post-starburst galaxies (PSBs) are transition galaxies showing evidence of recent rapid star formation quenching. To understand the role of galaxy mergers in triggering quenching, we investigate the incidence of PSBs and resolved PSB properties in post-merger galaxies using both SDSS single-fibre spectra and MaNGA resolved IFU spectra. We find post-mergers have a PSB excess of 10–20 times that relative to their control galaxies using single-fibre PSB diagnostics. A similar excess of ∼ 19 times is also found in the fraction of central (C)PSBs and ring-like (R)PSBs in post-mergers using the resolved PSB diagnostic. However, 60 per cent of the CPSBs + RPSBs in both post-mergers and control galaxies are missed by the single-fibre data. By visually inspecting the resolved PSB distribution, we find that the fraction of outside-in quenching is seven times higher than inside-out quenching in PSBs in post-mergers while PSBs in control galaxies do not show large differences in these quenching directions. In addition, we find a marginal deficit of HI gas in PSBs relative to non-PSBs in post-mergers using the MaNGA-HI data. The excesses of PSBs in post-mergers suggest that mergers play an important role in triggering quenching. Resolved IFU spectra are important to recover the PSBs missed by single-fibre spectra. The excess of outside-in quenching relative to inside-out quenching in post-mergers suggests that AGNs are not the dominant quenching mechanism in these galaxies, but that processes from the disc (gas inflows/consumption and stellar feedback) play a more important role.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Build and Execute Environment

At exascale, the challenge becomes to develop applications that run at scale and use exascale platforms reliably, efficiently, and flexibly. Workflows become much more complex because they must seamlessly integrate simulation and data analytics. They must include down-sampling, post-processing, feature extraction, and visualization. Power and data transfer limitations require these analysis tasks to be run in-situ or in-transit. We expect successful workflows will comprise multiple linked simulations along with tens of analysis routines. Users will have limited development time at scale and, therefore, must have rich tools to develop, debug, test, and deploy applications. At this scale, successful workflows will compose linked computations from an assortment of reliable, well-defined computation elements, ones that can come and go as required, based on the needs of the workflow over time. We propose a novel framework that utilizes both virtual machines (VMs) and software containers to create a workflow system that establishes a uniform build and execution environment (BEE) beyond the capabilities of current systems. In this environment, applications will run reliably and repeatably across heterogeneous hardware and software. Containers, both commercial (Docker and Rocket) and open-source (LXC and LXD), define a runtime that isolates all software dependencies from the machine operating system. Workflows may contain multiple containers that run different operating systems, different software, and even different versions of the same software. We will run containers in open-source virtual machines (KVM) and emulators (QEMU) so that workflows run on any machine entirely in user-space. On this platform of containers and virtual machines, we will deliver workflow software that provides services, including repeatable execution, provenance, checkpointing, and future proofing. We will capture provenance about how containers were launched and how they interact to annotate workflows for repeatable and partial re-execution. We will coordinate the physical snapshots of virtual machines with parallel programming constructs, such as barriers, to automate checkpoint and restart. We will also integrate with HPC-specific container runtimes to gain access to accelerators and other specialized hardware to preserve native performance. Containers will link development to continuous integration. When application developers check code in, it will automatically be tested on a suite of different software and hardware architectures.

Guan, Qiang↗
DOE OSTI2017

Build and Execute Environment

At exascale, the challenge becomes to develop applications that run at scale and use exascale platforms reliably, efficiently, and flexibly. Workflows become much more complex because they must seamlessly integrate simulation and data analytics. They must include down-sampling, post-processing, feature extraction, and visualization. Power and data transfer limitations require these analysis tasks to be run in-situ or in-transit. We expect successful workflows will comprise multiple linked simulations along with tens of analysis routines. Users will have limited development time at scale and, therefore, must have rich tools to develop, debug, test, and deploy applications. At this scale, successful workflows will compose linked computations from an assortment of reliable, well-defined computation elements, ones that can come and go as required, based on the needs of the workflow over time. We propose a novel framework that utilizes both virtual machines (VMs) and software containers to create a workflow system that establishes a uniform build and execution environment (BEE) beyond the capabilities of current systems. In this environment, applications will run reliably and repeatably across heterogeneous hardware and software. Containers, both commercial (Docker and Rocket) and open-source (LXC and LXD), define a runtime that isolates all software dependencies from the machine operating system. Workflows may contain multiple containers that run different operating systems, different software, and even different versions of the same software. We will run containers in open-source virtual machines (KVM) and emulators (QEMU) so that workflows run on any machine entirely in user-space. On this platform of containers and virtual machines, we will deliver workflow software that provides services, including repeatable execution, provenance, checkpointing, and future proofing. We will capture provenance about how containers were launched and how they interact to annotate workflows for repeatable and partial re-execution. We will coordinate the physical snapshots of virtual machines with parallel programming constructs, such as barriers, to automate checkpoint and restart. We will also integrate with HPC-specific container runtimes to gain access to accelerators and other specialized hardware to preserve native performance. Containers will link development to continuous integration. When application developers check code in, it will automatically be tested on a suite of different software and hardware architectures.

Guan, Qiang↗
DOE OSTI2025

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗
DOE OSTI2022

Flow structures of the cross-flow over a five-layer helically coiled steam generator geometry

The helically-coiled heat exchanger (HCSG) offers advantages over straight tubes such as compactness in geometry, increased heat transfer coefficients, and capability to absorb thermal expansion. In this study, the time-resolved velocity field data of shell-side crossflow over a five-layer helically coiled steam generator geometry was obtained using Particle Image Velocimetry (PIV) between adjacent rods in four regions with up to 5000 frames per second (5 kHz) at approximately Re u = 3600. The slant rod bundle created different flow patterns in different regions. Proper Orthogonal Decomposition (POD) revealed high-energy-mode flow structures. Due to the confined flow channel geometry, the large flow structures were not advected resulting in no strong POD mode paring. Continuous Wavelet Transform (CWT) visualized the characteristics of flow fluctuation in time and frequency domain simultaneously, which ranged in Strouhal number, S u , between 0.03 and 1.17. In the wake regions of the rods, flows changed in both pattern and magnitude over time. Once a pattern formed, it continued for relatively shorter than straight tube arrays, before changing the pattern again. The histograms of the POD time coefficients presented this multimodal flow characteristics. The current study focuses more on whole flow field analyses in the selected regions rather than local point-wise analyses. Finally, this study revealed that various Strouhal numbers can exist in the helically coiled heat exchanger geometry, which was not observed in straight tube arrays. In particular, Strouhal numbers smaller than 0.1 do not seem to be caused by vortex shedding but seem to be the result of the combination of complex geometry and possible multimodal trends. Based on the current observation, further works would involve intensive instantaneous point-wise analyses on the interactions between important locations such as separation points, stagnation points, shear layers, etc.

42 ENGINEERING↗
DOE OSTI2025

In situ Detection of Plasma Induced Surface Interaction based on Deep Learning based Visual Diagnostics (Technical Report)

It is characteristic for many plasma devices to undergo plasma-material interaction leading to surface erosion. These processes, often not easily detectable, lead to changes in device performance and lifespan. State-of-the-art lifetime tests and wear experiments require over 1000s hours. A self-consistent model for accurately predicting the erosion's effects is not available. In situ detection of these processes is not a trivial task since the surface variations at the early stages have a micron scale. Such limitations not only restrict testing and prediction capabilities but also slow the development of new thrusters and limit mission duration. To address these challenges, an in-situ diagnostic for real-time erosion assessment has been developed, aiming to expedite lifetime testing and broaden experimental campaigns. Several works were dedicated to real-time and in situ monitoring of material erosion during plasma exposure using laser holography, microscopy, and with telemicroscopes. However, the applicability of these approaches is limited due to complexity, cost and less flexibility as they often require placing diagnostic equipment inside the vacuum chamber. In collaboration with Princeton Collaborative Research Facility (PCRF), Princeton Plasma Physics Laboratory (PPPL), a new diagnostic approach is developed, where geometry modifications to the ceramic channel walls were introduced that would result in accelerated channel erosion. We employed Long-distance microscope (LDM) imagery, combined with Deep-Learning based Shape from focus or depth from focus (DFF or SFF) approach, that provides an accessible and cost-effective solution. LDM employs focus variation techniques to continuously capture multiple images of the target object at distinct focal planes. DFF, an optical focus variation method, generates a 3D topographical surface depth map from a sequence of variably focused images. Combined with the developed diagnostic, this approach offers a controllable means to study erosion under accelerated conditions. In this work, we develop Neural Network-based DFF algorithm applicable for LDM data to quantitatively evaluate plasma induced surface modification from LDM data. Next, we develop Deep Learning-based super-resolution depth map image reconstruction technique to increase the resolution of depth maps obtained from DFF algorithm to improve the accuracy of erosion measurements. Thirdly, we develop several image processing techniques to remove noise and improve the quality of depth map image. Here we report the results of initial tests for this approach. An experimental setup designed and built in PPPL was employed that consists of a 3-cm gridded ion source that produces a neutralized argon beam with energies up to 600 eV. A hexagonal boron nitride (h-BN) ceramic target, designed based on computational predictions, was used. Tests were conducted to reconstruct the complex geometry of the target under the lighting conditions of the operated ion source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Initial Assessment of Ruthenium Removal Systems for Tritium Pretreatment Off-Gas

The processing of used nuclear fuel releases volatile and semivolatile radionuclides into the off-gas streams associated with the various processing steps. Ongoing efforts to develop methods to control the release of gaseous radionuclides to the environment have identified ruthenium (specifically 106 Ru) as one of the semivolatile nuclides requiring a high degree of abatement. The tritium pretreatment process will release a fraction of the ruthenium as volatile ruthenium tetroxide (RuO 4 ). Scoping tests were conducted to characterize the adsorption/deposition of RuO 4 from dry gas streams by two different substrates: steel wool and silica gel. A simple test system was assembled to produce RuO 4 through oxidation of powdered ruthenium metal, followed by a temperature-controlled sorption column. Tests were conducted with no sorbent present in the column to obtain baseline data on the evolution of RuO 4 and deposition of ruthenium on the internal surfaces of the column. Silica gel and steel wool were tested as adsorbent media. Gaseous effluent from the column was scrubbed with caustic solution to remove any remaining ruthenium. Ultraviolet/visible spectroscopy (UV/Vis) and inductively coupled plasma–mass spectroscopy analyses were used to measure the amount of ruthenium in the scrub solutions. The UV/Vis analysis proved to be problematic because of precipitation of ruthenates, ruthenium oxides, and hydrated ruthenium oxides that occurred over time. An effective stabilization agent that does not adversely affect the desired spectra needs to be identified if this method is to be used in the future. There were indications that oxidation of powdered ruthenium metal did not generate a stream with steady concentration of RuO 4 , so improved techniques are needed for follow-on experiments designed to infer deposition rate. Both steel wool and silica gel are able to decrease the effluent Ru concentration by a factor of 10–100 during the 5-hour tests conducted. These tests were conducted using a dilute RuO 4 -bearing dry gas stream and the short test duration did not result in sorbent saturation or in an estimation of Ru breakthrough time. RuO 4 desorbs from silica gel at room temperature. As interactions between RuO 4 and silica gel are not irreversible, it cannot function as a disposal media without pretreatment to reduce RuO 4 to a less mobile species. Silica gel may be useful as a RuO 4 sorbent or delay media in some applications. Based on visual observations, steel wool was an effective deposition substrate with no evidence of RuO 4 release. The deposition form should be identified and methods to assay the ruthenium content need to be developed.

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

High repetition rate mapping of the interaction between a laser plasma and magnetized background plasma via laser induced fluorescence

In this work, the laminar coupling of energy between a laser-produced plasma and a background magnetized plasma was investigated via planar laser induced fluorescence diagnostic and magnetic flux probes. Experiments performed on the Large Plasma Device at the University of California, Los Angeles, mapped out the two-dimensional spatiotemporal evolution of the laser-plasma (debris) ion velocity distribution function (VDF) to assess debris-background coupling in a sub-Alfvénic regime. The acquisition of these data necessitates high repetition rate (1 Hz) as each dataset is the accumulation of thousands of laser shots, which would not be feasible in single-shot experiments. Fully kinetic, three-dimensional particle-in-cell simulations are compared to the measured VDFs to provide a framework in which we can understand the coupling of a sub-Alfvénic plasma flow through a preformed, magnetized plasma. The simulations display the same departure from the expected gyromotion of the debris plasma as observed in the experimental data, and in conjunction with the measured magnetic field traces, have led to the direct observation of the collisionless coupling via laminar fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Supporting Information - How Can Crop Production Adapt to growing groundwater restrictions in the West?

Abstract This data plan outlines the structure and content of datasets generated and utilized in the research presented in the manuscript abstract. Groundwater overdraft has led to serious water supply issues in the US West. Most western states manage groundwater use through a permitting system, but California has only recently begun to restrict groundwater use statewide with the 2014 passage of groundwater restrictions which target the elimination of groundwater overdraft practices by 2042. With groundwater extraction curtailed, crop production in the US West (a $95 billion industry annually) will be affected, and appropriate response strategies will be needed to ensure minimal disruption to food production and the regional economy. In this paper, we explore the adoption of alternative adaptive responses: (a) deficit irrigation; (b) switching to less water-intensive crops; (c) changing the extent of irrigated land (including fallowing); and (d) geographically shifting crop production. Employing an integrated modeling approach, we explicitly capture the interactions and feedbacks between local hydrology, changes in crop yields, crop and land use decision-making, changes in crop prices, and regional shifts in crop production. We find that the optimal adaptive response is spatially heterogeneous and comprises a portfolio of strategies. Southwestern states and California will be the most impacted by groundwater restrictions. The optimal responses in these states are to both adopt deficit irrigation strategies and reduce a portion of their irrigated croplands, resulting in a shift in crop production to northwestern states with a larger supply of water. The datasets detailed in this repository represent the output from these integrated models, specifically designed to support the analysis and visualization presented in the manuscript and supplemental information. These datasets, used in conjunction with the scripts available at our associated GitHub repository (https://github.com/pches/Femeena_etal_How_can_crop_production_adapt), enable the reproduction of figures and facilitate a comprehensive understanding of the optimal adaptive strategies in agriculture for mitigating water stress under varying groundwater extraction scenarios, as outlined in the abstract.

DNDCe↗
DOE OSTI2025

Supporting Information - How Can Crop Production Adapt to growing groundwater restrictions in the West?

Abstract This data plan outlines the structure and content of datasets generated and utilized in the research presented in the manuscript abstract. Groundwater overdraft has led to serious water supply issues in the US West. Most western states manage groundwater use through a permitting system, but California has only recently begun to restrict groundwater use statewide with the 2014 passage of groundwater restrictions which target the elimination of groundwater overdraft practices by 2042. With groundwater extraction curtailed, crop production in the US West (a $95 billion industry annually) will be affected, and appropriate response strategies will be needed to ensure minimal disruption to food production and the regional economy. In this paper, we explore the adoption of alternative adaptive responses: (a) deficit irrigation; (b) switching to less water-intensive crops; (c) changing the extent of irrigated land (including fallowing); and (d) geographically shifting crop production. Employing an integrated modeling approach, we explicitly capture the interactions and feedbacks between local hydrology, changes in crop yields, crop and land use decision-making, changes in crop prices, and regional shifts in crop production. We find that the optimal adaptive response is spatially heterogeneous and comprises a portfolio of strategies. Southwestern states and California will be the most impacted by groundwater restrictions. The optimal responses in these states are to both adopt deficit irrigation strategies and reduce a portion of their irrigated croplands, resulting in a shift in crop production to northwestern states with a larger supply of water. The datasets detailed in this repository represent the output from these integrated models, specifically designed to support the analysis and visualization presented in the manuscript and supplemental information. These datasets, used in conjunction with the scripts available at our associated GitHub repository (https://github.com/pches/Femeena_etal_How_can_crop_production_adapt), enable the reproduction of figures and facilitate a comprehensive understanding of the optimal adaptive strategies in agriculture for mitigating water stress under varying groundwater extraction scenarios, as outlined in the abstract.

DNDCe↗
DOE OSTI2025

Supporting Information - How Can Crop Production Adapt to growing groundwater restrictions in the West?

Abstract This data plan outlines the structure and content of datasets generated and utilized in the research presented in the manuscript abstract. Groundwater overdraft has led to serious water supply issues in the US West. Most western states manage groundwater use through a permitting system, but California has only recently begun to restrict groundwater use statewide with the 2014 passage of groundwater restrictions which target the elimination of groundwater overdraft practices by 2042. With groundwater extraction curtailed, crop production in the US West (a $95 billion industry annually) will be affected, and appropriate response strategies will be needed to ensure minimal disruption to food production and the regional economy. In this paper, we explore the adoption of alternative adaptive responses: (a) deficit irrigation; (b) switching to less water-intensive crops; (c) changing the extent of irrigated land (including fallowing); and (d) geographically shifting crop production. Employing an integrated modeling approach, we explicitly capture the interactions and feedbacks between local hydrology, changes in crop yields, crop and land use decision-making, changes in crop prices, and regional shifts in crop production. We find that the optimal adaptive response is spatially heterogeneous and comprises a portfolio of strategies. Southwestern states and California will be the most impacted by groundwater restrictions. The optimal responses in these states are to both adopt deficit irrigation strategies and reduce a portion of their irrigated croplands, resulting in a shift in crop production to northwestern states with a larger supply of water. The datasets detailed in this repository represent the output from these integrated models, specifically designed to support the analysis and visualization presented in the manuscript and supplemental information. These datasets, used in conjunction with the scripts available at our associated GitHub repository (https://github.com/pches/Femeena_etal_How_can_crop_production_adapt), enable the reproduction of figures and facilitate a comprehensive understanding of the optimal adaptive strategies in agriculture for mitigating water stress under varying groundwater extraction scenarios, as outlined in the abstract.

DNDCe↗
DOE OSTI2022

Testing the role of AGN on the star formation and metal enrichment of ‘twin galaxies’

Here, we explore the effect of AGN activity on the star formation history of galaxies by analysing the stellar population properties of 10 pairs of nearby twin galaxies – selected as being visually similar except for the presence of an AGN. The selection of such twin samples represents a method to study AGN feedback, as recently proposed by del Moral Castro et al. We use integral field unit (IFU) data from CALIFA, stacked within three fixed apertures. AGN galaxies in a twin pair suggest more evolved stellar populations than their non-AGN counterpart 90 per cent of the time, regardless of aperture size. A comparison with a large sample from SDSS confirms that most twins are representative of the general population, but in each twin the differences between twin members is significant. A set of targeted line strengths reveal the AGN member of a twin pair is older and more metal-rich than the non-AGN galaxy, suggesting AGN galaxies in our sample may either have an earlier formation time or follow a different star formation and chemical enrichment history. These results are discussed within two simple, contrasting hypotheses for the role played by AGN in galaxy evolution, which can be tested in the future at a greater detail with the use of larger data sets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Automated detection of double nuclei galaxies using GOTHIC and the discovery of a large sample of dual AGN

We present a novel algorithm to detect double nuclei galaxies (DNG) called Gothic (Graph-bOosTed iterated HIll Climbing) – that detects whether a given image of a galaxy has two or more closely separated nuclei. Our aim is to test for the presence of dual/multiple active galactic nuclei (AGN) in galaxies that visually represent a DNG. Although galaxy mergers are common, the detection of dual AGN is rare. Their detection is very important as they help us understand the formation of supermassive black hole (SMBH) binaries, SMBH growth and AGN feedback effects in multiple nuclei systems. There is thus a need for an algorithm to do a systematic survey of existing imaging data for the discovery of DNGs and dual AGNs. We have tested Gothic on an established sample of DNGs with a 100 per cent detection rate and subsequently conducted a blind search of 1 million SDSS DR16 galaxies (with spectroscopic data available) lying in the redshift range of z = 0 to 0.75. From the list of candidate DNGs found, we have detected 159 dual AGNs, of which 2 are triple AGN systems. Our results show that dual AGNs are not common, and triple AGN even rarer. The colour (u–r) magnitude plots of the DNGs indicate that star formation is quenched as the nuclei come closer and as the AGN fraction increases. The quenching is especially prominent for dual/triple AGN galaxies that lie at the extreme end of the red sequence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Scalable edge clustering of dynamic graphs via weighted line graphs

Timestamped relational datasets consisting of records (or connections) between pairs of entities are ubiquitous in network science. For applications like peer-to-peer communication, email, various social network interactions, and computer network security, it is useful to organize these records into groups based on how and when they are occurring. Weighted line graphs offer a natural way to model how records are related in such datasets but for large real-world graph topologies, building and utilizing the line graph is prohibitively expensive. Here, we present the framework to cluster the edges of a dynamic graph via the associated line graph that contains two major contributions. The first is a method to work with the line graph implicitly and the second is a distributed scale implementation of an agglomerative hierarchical graph clustering algorithm. We outline a novel hierarchical dynamic graph edge clustering approach that efficiently breaks massive relational datasets into small sets of edges containing events at various timescales. This is in stark contrast to traditional graph clustering algorithms that prioritize highly connected (clique-like) community structures. Our approach relies on constructing a sufficient subgraph of a weighted line graph and applying a hierarchical agglomerative clustering. This approach is related to scalable techniques from spatial clustering, nonlinear-dimension reduction, topological data analysis, and draws particular inspiration from HDBSCAN. As an edge clustering, this method yields an overlapping node clustering. Our algorithm is parallelizable and we demonstrate efficient clustering of a billion-scale, real-world dynamic graph into small edge sets that correlate in topology and time. The entire clustering process for a graph with tens of billions of edges takes just a few minutes of run time on 256 nodes of a distributed compute environment. We argue how the output of the edge clustering is useful for a multitude of data visualization and powerful machine learning tasks, both involving the original massive dynamic graph data and metadata associated with the nodes and edges. Finally, we describe how this approach can be extended to dynamic hypergraphs and dynamic graphs/hypergraphs with unstructured data living on vertices and edges.

Data Analysis↗
DOE OSTI2022

Findings on Subtask 3.1 - Bakken Rich Gas Enhanced Oil Recovery Project

Total in-place oil for the Bakken petroleum system (BPS) (which includes the Bakken and Three Forks Formations) has been estimated to be 600 billion barrels (bbl). However, BPS wells have decline rates as high as 85% over the first 3 years of their lives, and primary recovery factors typically range from 3% to 10% of original oil in place. Given the low initial recovery rates, even small incremental productivity improvements could dramatically increase technically recoverable oil in the BPS. One potential solution is enhanced oil recovery (EOR) using gas injection, such as carbon dioxide (CO2) or hydrocarbon (HC) gases. While commonly used in conventional reservoirs, CO2 EOR in unconventional tight oil reservoirs has been limited to pilot tests. EOR using rich gas (mixture of methane, ethane, and propane) has also been employed in numerous pilots in several unconventional plays and has recently been successfully applied in the Eagle Ford play. If successful, large-scale gas-based EOR in the BPS could dramatically increase oil productivity and recovery factors and extend the life of the play for decades. While CO2 may be a technically suitable working fluid for EOR in the BPS, supplies are limited and costs for using CO2 in EOR pilots are prohibitively high. Meanwhile, produced gas flaring has presented challenges for BPS operators in North Dakota. Analysis conducted by the North Dakota Pipeline Authority indicates that the current gas-gathering infrastructure in North Dakota is insufficient to accommodate all of the associated gas that is produced from the BPS. The geographically isolated location of North Dakota relative to large natural gas markets, combined with suppressed natural gas prices, has made it economically challenging for industry to invest capital in expanding gas-gathering infrastructure in the state. These circumstances led to a research program conducted by the Energy & Environmental Research Center (EERC) in partnership with Liberty Resources Management Company LLC (LR) to examine the potential to use rich gas injection for EOR and mitigate flaring. A rich gas EOR pilot test was designed and executed by LR at its Stomping Horse development area in Williams County, North Dakota. From July 2018 through May 2019, a total of 160 million standard cubic feet (MMscf) of rich produced gas was injected into the BPS using five different wells in a sequential injection strategy. LR’s Leon–Gohrick drill spacing unit (DSU) was used as the test site. Regulatory oversight was provided by the North Dakota Industrial Commission (NDIC). Technical support was provided by the EERC through a series of laboratory, modeling, and field-based activities, and additional post-pilot research activities incorporated learnings from the test, developed new laboratory data, improved fracture modeling methods, and developed machine learning and big data analytics. The results from the Stomping Horse rich gas EOR pilot activities indicate that developing an effective, economical EOR approach for the BPS will require more field tests. Another key lesson learned from the Stomping Horse tests is that detailed pre- and posttest data on reservoir conditions and fluids production are essential. Robust reservoir characterization provides information that is crucial to creating realistic geomodels and conducting valid dynamic simulations of potential EOR scenarios. A detailed understanding of the completions and production history of offset wells is also necessary for valid test result interpretations. This knowledge is essential to designing the operational parameters of injectivity tests and interpreting the results. A conformance control strategy is also essential to success. Laboratory-based examinations of rich gas interactions with reservoir fluids and rocks were conducted, with an emphasis on determining the ability to mobilize oil in the tight reservoir rocks and shales of the BPS. Injection fluid composition was shown to have a positive impact on reducing reservoir oil minimum miscibility pressure (MMP), reducing interfacial tension (IFT), and altering wettability. IFT and contact angle measurements demonstrated that wettability can be altered in the presence of rich gas, suggesting the potential to improve oil recovery. Iterative modeling of surface infrastructure and reservoir performance using data generated by the various project activities was conducted. A geologic model of the Stomping Horse area was built; history-matched oil, gas, and water production was used in simulations of various EOR scenarios. Early programmatic modeling results were used to support LR’s design and operation of the EOR pilot and to provide insight regarding optimization of future commercial-scale BPS EOR design and operations. Post-pilot modeling focused on alternative methods of understanding complex fracture networks and accelerating simulation time. These led to improved simulation run times and provide excellent history-matching results. Several of these iterative models were used as the bases for developing algorithms into machine learning and big data analytics. History matching in reservoir simulation is time-consuming and computer processing-intensive. Machine learning algorithms were created, and an automated history-matching tool was developed. A large set of synthetic reservoir simulations were created to generate well responses (oil, gas, and water production, well bottomhole pressure [BHP], and tracer or propane breakthrough) for a set of EOR operating parameters that included offset well status (open or closed), injectate (rich gas or propane), injection rate, and injection well BHP. A user interface was developed to provide real-time visualization. Machine learning-based models were developed to provide rapid forecasting of well performance given a set of user-defined EOR operating parameters. These predictive models allow the user to modify the offset well status, injection rate, and injection well BHP and rapidly forecast future production performance. The combination of real-time visualization tools with real-time forecasting tools provides a framework for real-time control—operational changes that the EOR site operator can enact (e.g., changing gas injection rates) to affect the observed performance and potentially improve the EOR outcome. There is great reason to be optimistic about the future of EOR in the Bakken. The results of the laboratory studies suggest significant potential for high rates of oil mobilization using produced field gas injection under the right conditions. The results of the lab studies, combined with rigorous statistical analysis of well production data and associated modeling efforts, confirm the notion that fluid mobility within the reservoir is controlled by fractures. As more knowledge is gained about the nature and distribution of fracture networks in the Bakken, the industry will be in a better position to predict and, ultimately, influence fluid mobility. New field tests are necessary to develop a more complete understanding of those conditions. Thoughtful and creatively engineered field tests within a well-characterized geologic setting will yield the fundamental knowledge needed to take Bakken oil production to the next level. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission’s Oil and Gas Research Program and Computer Modelling Group.

04 OIL SHALES AND TAR SANDS↗