Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analysis and 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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Community Energy Operations and Planning System: Concept, Use cases, Metrics, and Benefits

Community and city leaders are interested in achieving sustainability goals, providing resilient energy infrastructure, and improving economic competitiveness. Community-level data acquisition and analysis can provide energy and associated benefits that are not possible at the single building level. However, there is a lack of organizational structure, common semantic data models, interoperable systems, and methods to support data-driven decision-making for community-scale energy supply and demand systems. We explored the need and opportunity for a Community Energy Operations and Planning System (Community EOPS), a potential data exchange platform. We conducted “customer discovery” interviews, and reviewed literature, public tools, and technology platforms to identify key energy data “users” and use cases in communities. The key users of the Community EOPS could be developers of mixed-use districts, corporate, defense and university campus energy managers, and city managers of cities that own their energy utility. The value could be for community planning and reporting (for energy data-integrated land use planning and community infrastructure investments in microgrids, storage, district heating and cooling), energy efficiency (leveraging optimizations for community scale energy supply and demand), flexible load management (grid-edge load management to offset, shift, and flatten loads for multiple buildings and EV fleets), cost savings and revenue generation (participating in grid services), and social benefits such as energy resilience, equity, and awareness. We developed a conceptual Community EOPS architecture with recommendations for streamlined and prioritized data acquisition, sharing, and integration driven by prioritized use cases, common metrics, and actionable visualizations that can provide value to a community’s users.

Singh, Reshma↗
DOE OSTI2021

Online Monitoring System for Concrete Structures Affected by Alkali-Silica Reaction

This report presents a comprehensive study of developing ultrasonic wave and acoustic emission (AE) techniques for long-term monitoring of alkali-silica reaction (ASR) development in concrete. Small, medium, and full-scale concrete specimens were cast, conditioned, and monitored for various periods, from at least one year to 2.4 years. The reactive concrete specimens contain coarse or fine reactive aggregates to study the effects of different types of reactive aggregates. Confinements were also designed to simulate the 2-dimensional confinement effects of reinforcement in the shield building of nuclear power plants. The concrete specimens were stored in a environmental chamber with high humidity and high temperature to accelerate the ASR development. The ultrasonic monitoring data shows high sensitivity to ASR development and could detect cracking initiation well before visible surface cracks occurred. However, the linear ultrasonic analysis based on wave velocity is strongly affected by temperature variation. Therefore, a nonlinear ultrasonic method was proposed to measure thermally induced nonlinear acoustic responses of concrete (thermal modulation of ultrasonic wave). The measured nonlinear acoustic parameter shows a high correlation with ASR expansion across specimens with different reactive aggregates and confinement conditions. The same conclusion was obtained from nonlinear resonance tests on small concrete prisms. Compared to the linear acoustic methods, the nonlinear tests show high sensitivities to ASR damage from internal microcracking initiation at the early stage to visible cracks at the late stage of ASR. The attributes of the thermal modulation of nonlinear ultrasonic method include high sensitivity, immunity to temperature effects, and strong correlation with ASR expansion, which demonstrate great potentials of the nonlinear ultrasonic method for diagnosis of ASR damage and prediction of concrete deterioration process. Acoustic emission is a passive sensing technique for damage assessment, and access to only one surface is needed even for thick and heavily reinforced elements such as the walls utilized for nuclear shield building. Additionally, relatively few sensors are required to monitor the progression of the damage process. Results indicate that the AE data can be related to the damage rating index, which is a petrography-based means of assessing damage due to ASR in reinforced concrete. Furthermore, AE is capable of detecting ASR damage long before surface cracking is visually noticeable. Boundary conditions play in important role in the progression of ASR damage and differences in boundary conditions are reflected in the AE data. Entropy based data assessment methods provide a means to assess the damage state, and convolutional neural network based data assessment provided a means to assess the damage state in real-time. Results indicate that artificial neural network models may be used as a means to predict volumetric expansion based on AE data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

NEMESIS: Exoplanet Transit Survey of Nearby M-dwarfs in TESS FFIs. I

In this work, we present an analysis of 33,054 M-dwarf stars, located within 100 parsecs, via the Transiting Exoplanet Survey Satellite (TESS) full-frame images (FFIs) of observed sectors 1–5. We present a new pipeline called NEMESIS, developed to extract detrended photometry, and to perform transit searches of single-sector data in TESS FFIs. As many M-dwarfs are faint, and are not observed with a two-minute cadence by TESS, FFI transit surveys can provide an empirical validation of how many planets are missed, using the 30-minute cadence data. In this work, we detect 183 threshold crossing events, and present 29 candidate planets for sectors 1–5, 24 of which are new detections. Our sample contains orbital periods ranging from 1.25 to 6.84 days, and planetary radii from 1.26 to 5.31 R {sub ⊕}. With the addition of our new planet candidate detections, along with detections previously observed in sectors 1–5, we calculate an integrated occurrence rate of 2.49 ± 1.58 planets per star, for the period range ∈ [1, 9] days, and planet radius range ∈ [0.5,11] R {sub ⊕}. We project an estimated yield of 122 ± 11 transit detections of nearby M-dwarfs. Of our new candidates, 23 have signal-to-noise ratios >7, transmission spectroscopy metrics >38, and emission spectroscopy metrics >10. We present all of our data products for our planet candidates via the Filtergraph data visualization service, located at https://filtergraph.com/NEMESIS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Date Release Report for Large Surface Explosion Coupling Experiment (LSECE) Nevada National Security Site

The DTRA sponsored 2020-2022 Large Surface Explosion Coupling Experiment (LSECE) consists of two large ground surface chemical explosions, data collection, analysis and modeling carried out in 2020-2022. The LSECE chemical explosions were carried out at the site of prior NNSA sponsored buried chemical explosions that were part of the Source Physics Experiment (SPE) Phase II in Dry Alluvium Geology (DAG). This report described the data collected under LSECE that is being made publicly available. The prior buried explosion DAG data are described in a separate report (https://www.osti.gov/biblio/1825534) and are also publicly available. The LSECE explosions, data analysis and numerical modeling work sponsored by DTRA were intended to address three objectives: 1) Generate seismo-acoustic data to test and improve numerical models of explosion energy coupling in dry alluvium geology; 2) To improve seismo-acoustic yield estimation techniques as a function of depth and medium properties; 3) To study and improve acoustic propagation modeling under two different atmospheric conditions. The LSECE data that were collected have expanded the prior DAG buried explosion dataset to include surface chemical explosions recorded on a common set of seismo-acoustic stations. They have allowed a more detailed study of above and below surface seismo-acoustic energy coupling with applications to explosion monitoring and assessment. The two approximately 1-ton TNT equivalent yield LSECE explosions were detonated at different times of the day to explore the effects of the different atmospheric conditions. The first chemical explosion “Artemis” was conducted before dawn when temperature inversions were present. The second chemical explosion “Apollo” was conducted on a sunny afternoon when the temperature gradient was more linear. The LSECE chemical explosion data were also collected across a variety of different sensor types to allow evaluation of their effectiveness in recovering useful information. The LSECE instrumentation included fiber optic or Distributed Acoustic Sensing (DAS), a dense array or Large-N array of seismometers, borehole accelerometers, DAS and a velocity meter, airborne acoustic instruments and a variety of visual and remote sensing data.

58 GEOSCIENCES↗
DOE OSTI2024

Solid–liquid separation of lignocellulosic sugars from biomass by rotating ceramic disc filtration

In many biomass conversion processes, the separation of cellulosic sugars from residual, lignin-rich solids is a critical step, and achieving high recovery yields of sugars by conventional tangential crossflow and vacuum filtration is challenged by the presence of biomass solids, which rapidly foul filters, resulting in decreased throughput. Considering the performance limitations of traditional filtration methods, dynamic filtration, which generates high shear at the membrane surface to decrease fouling, is emerging as a viable alternative for demanding solid–liquid separations. For high solids separations, there is little available information regarding the performance, limitations, and energy consumption of dynamic filtration. To that end, here we characterized the performance of a dynamic filtration module, specifically a rotating ceramic disc (RCD) filter, for the aseptic recovery of cellulosic sugars from biomass solids following pretreatment and enzymatic hydrolysis. We show how RCD rotational velocity and percent biomass solids impact the filter throughput. Additionally, we used computational fluid dynamics (CFD) simulations to estimate the shear rate at the membrane surface and to visualize hydrodynamic profiles within the module. With the combined CFD simulations and experimental results, we estimated the energy demand and operating expenses for a viable dynamic filtration system operating with a lignocellulosic feed slurry. Our results indicate that an RCD filter can achieve ≥95% recovery of sugars and produce a retentate slurry containing 12 wt% insoluble solids with low energy consumption (a 2.2-fold improvement over cross-flow filtration) and low operating costs ($\$$0.06 per kg sugars). These results show a viable path towards operationally reliable, energy efficient, and cost-effective separations of sterilized cellulosic sugars from biomass solids and highlight the potential of dynamic filtration systems for challenging solid–liquid separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The edges of galaxies: Tracing the limits of star formation

The outskirts of galaxies have been studied from multiple perspectives for the past few decades. However, it is still unknown if all galaxies have clear-cut edges similar to everyday objects. We address this question by developing physically motivated criteria to define the edges of galaxies. Based on the gas density threshold required for star formation, we define the edge of a galaxy as the outermost radial location associated with a significant drop in either past or ongoing in situ star formation. We explore ~1000 low-inclination galaxies with a wide range in morphology (dwarfs to ellipticals) and stellar mass (10 7 M ⊙ < M * < 10 12 M ⊙ ). The location of the edges of these galaxies (R edge ) were visually identified as the outermost cutoff or truncation in their radial profiles using deep multi-band optical imaging from the IAC Stripe82 Legacy Project. We find this characteristic feature at the following mean stellar mass density, which varies with galaxy morphology: 2.9 ± 0.10 M ⊙ pc -2 for ellipticals, 1.1 ± 0.04 M ⊙ pc -2 for spirals, and 0.6 ± 0.03 M ⊙ pc -2 for present-day star-forming dwarfs. Additionally, we find that R edge depends on its age (colour) where bluer galaxies have larger R edge at a fixed stellar mass. The resulting stellar mass–size plane using R edge as a physically motivated galaxy size measure has a very narrow intrinsic scatter (≲0.06 dex). These results highlight the importance of new deep imaging surveys to explore the growth of galaxies and trace the limits of star formation in their outskirts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Clustering-Based Scenario Generation Framework for Power Market Simulation with Wind Integration

A critical step in stochastic optimization models of power system analysis is to select a set of appropriate scenarios and significant numbers of scenario generation methods exist in the literature. This paper develops a clustering based scenario generation method, which aims to improve the performance of existing scenario generation techniques by grouping a set of correlated wind sites into clusters according to their cross-correlations. Copula based models are utilized to model spatiotemporal correlations and the Gibbs sampling is then used to generate scenarios for day-ahead markets. Our results show that the generated scenarios based on clustered wind sites outperform existing approaches in terms of reliability and sharpness and can reduce the total computational time for scenario generation and reduction significantly. The clustering-based framework can therefore provide a better support for real-world market simulations with high wind penetration.

data visualization↗
DOE OSTI2022

X-ray nano-imaging of defects in thin film catalysts via cluster analysis

Functional properties of transition-metal oxides strongly depend on crystallographic defects; crystallographic lattice deviations can affect ionic diffusion and adsorbate binding energies. Scanning x-ray nanodiffraction enables imaging of local structural distortions across an extended spatial region of thin samples. Yet, localized lattice distortions remain challenging to detect and localize using nanodiffraction, due to their weak diffuse scattering. Here, in this study, we apply an unsupervised machine learning clustering algorithm to isolate the low-intensity diffuse scattering in as-grown and alkaline-treated thin epitaxially strained SrIrO 3 films. We pinpoint the defect locations, find additional strain variation in the morphology of electrochemically cycled SrIrO 3 , and interpret the defect type by analyzing the diffraction profile through clustering. Our findings demonstrate the use of a machine learning clustering algorithm for identifying and characterizing hard-to-find crystallographic defects in thin films of electrocatalysts and highlight the potential to study electrochemical reactions at defect sites in operando experiments.

42 ENGINEERING↗
DOE OSTI2023

Data-driven selection and spectral classification of white dwarf stars

The next generation of spectroscopic surveys is expected to provide spectra for hundreds of thousands of white dwarf (WD) candidates in the upcoming years. Currently, spectroscopic classification of white dwarfs is mostly done by visual inspection, requiring substantial amounts of expert attention. We propose a data-driven pipeline for fast, automatic selection, and spectroscopic classification of WD candidates, trained using spectroscopically confirmed objects with available Gaia astrometry, photometry, and Sloan Digital Sky Survey (SDSS) spectra with signal-to-noise ratios ≥9. The pipeline selects WD candidates with improved accuracy and completeness over existing algorithms, classifies their primary spectroscopic type with ≥ 90 % accuracy, and spectroscopically detects main sequence companions with similar performance. We apply our pipeline to the Gaia Data Release 3 cross-matched with the SDSS Data Release 17 (DR17), identifying 424 096 high-confidence WD candidates and providing the first catalogue of automated and quantifiable classification for 36 523 WD spectra. Both the catalogue and pipeline are made available online. Such a tool will prove particularly useful for the undergoing SDSS-V survey, allowing for rapid classification of thousands of spectra at every data release.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Full Scale HEPA Filter Encapsulation in Ultra-High-Performance Grout - Proof of Concept

The Hanford Tank Waste Treatment and Immobilization Plant (WTP) currently being constructed to treat radioactive waste, includes vitrification facilities for both the high-level waste (HLW) and low activity waste (LAW) fractions. Operation of the WTP will produce contaminated high-efficiency particulate air filters (HEPA), as part of the solid secondary waste (SSW) stream. The HEPA filters receive off-gas from the vessel vent header and primary off-gas treatment system in LAW Facility and remove particulate contaminants including 99 Tc and 129 I salts. The current disposal method for HEPA filters is encapsulation in metal containers using cementitious material (CM) and disposal in the Integrated Disposal Facility (IDF). Results from the 2017 IDF Performance Assessment (PA) WRPS (2018) demonstrated that while compliance is maintained for the 1000 year compliance period mandated by DOE O 435.1 and its accompanying manual, release of constituents from the HEPA filters result in exceedance of the performance objective imposed as the groundwater regulatory limit at later times. For example, at about 1500 years post-closure, solid secondary waste (SSW), including HEPA filters is predicted to become a dominant contributor to 99 Tc release and over the 10,000-year sensitivity analysis period, SSW is the dominant contributor of 129 I release to the groundwater. The estimated release could potentially be reduced if the HEPA filters are not compacted and waste containers could be distributed throughout a large space, thereby diluting the contaminant release. Additionally, a better cementitious material could be used to encapsulate the HEPA filters. One alternative method for disposal of contaminated HEPA filters is encapsulation of the filters in ultra high-performance grout (UHPG). UHPG is a variation of ultra-high-performance concrete (UHPC) is commonly used in the prestressed concrete industry for large structural members. Recent studies of UHPG show it has excellent properties for containing radionuclides such as 99 Tc and 129 I Nichols and Kaplan (2021). This report presents the results of the first attempt to encapsulate a clean, full-size HEPA filter in UHPG and evaluate the effectiveness of the immobilization process and final waste form. A full-scale proof-of-concept simulated waste form was prepared by encapsulating a HEPA filter in a 110-gallon stainless steel (SS) drum using UHPG. A change from Type I/II PC to Type 1L PLC was made after American Rock Products informed the team that they would no longer be using Type I/II by the end of 2024 and the northwest was phasing out Type I/II PC overall. Type I/II PC used in previous studies of UHPG for encapsulation (Nichols and Kaplan 2021). After the UHPG was cured both the scaled mockup and the full-scale simulated waste forms were sectioned for visual examination. UHPG completely encapsulated the filters and bonded to the external surfaces of materials comprising the filters. No cracks were observed in the sectioned waste forms.

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

Comparative study on the overall energy performance between photovoltaic and Low-E insulated glass units

In this report a novel semi-transparent building integrated photovoltaic (BIPV) laminate was developed and introduced in this paper. It was produced by cutting standard mono-crystalline silicon solar cells into small strips and then making electrical connections between each strip before laminating the cells between two layers of glass. The overall energy performance and energy saving potential of the BIPV insulated glass unit (IGU) under real world conditions were identified through a side by side comparative study. Compared to the reference IGU, the BIPV IGU had lower solar heat gain coefficient (SHGC) but much higher U-factor. The average HVAC electricity saving of the BIPV IGU was about 10% relative to the reference IGU. Daylighting measurement and analysis were carried out to evaluate the trade-offs associated with the BIPV IGU between daylight, glare, and lighting energy use. The results indicated that the BIPV IGU is better than the reference IGU in reducing discomfort glare. However, if the most conservative viewpoint near the window is used for the assessment, a lower transmittance BIPV IGU is required to bring the overall discomfort levels below the perceptible level. Lastly, the net energy saving potential associated with the novel BIPV IGU was identified based on the power, thermal and daylighting performance. On average, the BIPV IGU saved 16.8% of the total electricity use of the room. Further studies and improvement on the energy conversion efficiency of solar cells, the optimal transmittance as well as the thermal properties would make this technology more energy-efficient and affordable.

14 SOLAR ENERGY↗
DOE OSTI2022

Induced microearthquakes predict permeability creation in the brittle crust

Predicting the evolution of permeability accurately during stimulation at the reservoir scale and at the resolution of individual fractures is essential to characterize the fluid transport and the reactive/heat-transfer characteristics of reservoirs where stress exerts significant control. Here, we develop a hybrid machine learning (ML) model to visualize in situ permeability evolution for an intermediate-scale (~10 m) hydraulic stimulation experiment. This model includes an ML model that was trained using the well history of flow rate and wellhead pressure and MEQ data from the first three stimulation episodes to predict average permeability from the statistical features of the MEQs alone for later episodes. Moreover, a physics-inspired model is integrated to estimate in situ fracture permeability spatially. This method relates fracture permeability to fracture dilation and scales dilation to the equivalent MEQ magnitude, according to laboratory observations. The seismic data are then applied to define incremental changes in permeability in both space and time. Our results confirm the excellent agreement between the ground truth and model-predicted permeability evolution. The resulting permeability map defines and quantifies flow paths in the reservoir with the averaged permeability comparing favorably with the ground truth of permeability.

induced seismicity↗
DOE OSTI2023

Multi-Scale Characterization of Porosity and Cracks in Silicon Carbide Cladding after Transient Reactor Test Facility Irradiation

Silicon carbide (SiC) ceramic matrix composite (CMC) cladding is currently being pursued as one of the leading candidates for accident-tolerant fuel (ATF) cladding for light water reactor applications. The morphology of fabrication defects, including the size and shape of voids, is one of the key challenges that impacts cladding performance and guarantees reactor safety. Therefore, quantification of defects’ size, location, distribution, and leak paths is critical to determining SiC CMC in-core performance. This research aims to provide quantitative insight into the defect’s distribution under multi-scale characterization at different length scales before and after different Transient Reactor Test Facility (TREAT) irradiation tests. A non-destructive multi-scale evaluation of irradiated SiC will help to assess critical microstructural defects from production and/or experimental testing to better understand and predict overall cladding performance. X-ray computed tomography (XCT), a non-destructive, data-rich characterization technique, is combined with lower length scale electronic microscopic characterization, which provides microscale morphology and structural characterization. This paper discusses a fully automatic workflow to detect and analyze SiC-SiC defects using image processing techniques on 3D X-ray images. Following the XCT data analysis, advanced characterizations from focused ion beam (FIB) and transmission electron microscopy (TEM) were conducted to verify the findings from the XCT data, especially quantitative results from local nano-scale TEM 3D tomography data, which were utilized to complement the 3D XCT results. In this work, three SiC samples (two irradiated and one unirradiated) provided by General Atomics are investigated. The irradiated samples were irradiated in a way that was expected to induce cracking, and indeed, the automated workflow developed in this work was able to successfully identify and characterize the defects formation in the irradiated samples while detecting no observed cracking in the unirradiated sample. These results demonstrate the value of automated XCT tools to better understand the damage and damage propagation in SiC-SiC structures for nuclear applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

The geometry of signal space: a case study of direct mapping between seismic signals and event distribution

SUMMARY Under favourable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal-space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection and scaling, with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning and possible application to studies of repeating events in induced, volcanic and glacial seismicity. Our results place limits on the widely held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 2. Evaluating Controls on Flow Persistence in an Urbanized Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in an urbanized catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, distributed temperature sensing (DTS), continuous self-potential (SP) monitoring, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Field_Application subfolder contains the ATS XML input scripts, data files, output data for the SP site. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. The flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.m can only be used with COMSOL with MATLAB) is executed using the ATS output data to simulate the potential field. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) DTS Contains collated DTS data including raw Stokes and anti-Stokes measurement (provided as .h5 file). It also includes DTS processing.ipynb, a Jupyter notebook for calibrating the DTS data using dts_calibration Python package. cooler_calibration.csv is the DTS calibration CSV used in the calibration sequence. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion. 6) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 7) SP Contains the SP data collected in field at the SP sites (provided as CSV files). 8) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). Note: Code files (.ipynb, .py, .xml) can be opened in any standard code editor, .exo file can be viewed using Paraview, .h5 files can be opened using HDFView software and h5py Python package, and .resipy file can be opened with the open-source ResIPy software.

ATS↗
DOE OSTI2023

Molecular Hypernetworks for Exploration of Multi-Dimensional Metabolomics Data (Chyper)

Orthogonal separations of data from high-resolution mass spectrometry can provide insight into sample composition and help address the challenge of complete annotation of molecules in untargeted metabolomics. “Molecular networks” (MNs), as used, for example, in the Global Natural Products Social Molecular Networking platform, are an increasingly popular computational strategy for exploring and visualizing molecular relationships and improving annotation. MNs use graph representations to show the relationships between measured multidimensional data features. MNs also show promise for using network science algorithms to automatically identify targets for annotation candidates and to dereplicate features associated to a single molecular identity. However, more advanced methods may better represent the complexity present in samples. Our work aims to increase confidence in annotation propagation by extending molecular network methods to include “molecular hypernetworks” (MHNs), able to natively represent multiway relationships among observations supporting both human and analytical processing. In this paper we first introduce MHNs illustrated with simple examples, and demonstrate how to build them from liquid chromatography- and ion mobility spectrometry- separated MS data. We then describe a method to construct MHNs directly from existing MNs as their “clique reconstructions”, demonstrating their utility by comparing examples of previously published graph-based MNs to their respective MHNs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗
DOE OSTI2025

Geometric Interpretation of the Cluster Location Problem Part I: Theory

We present a new framing of the seismic location problem using principles drawn from differential geometry. Our interpretation relies upon the common assumption that travel times observed across a network are continuous, differentiable functions of source location. In consequence, travel‐time functions constitute a differentiable map between the source region and a Riemannian manifold. The manifold is said to be the image of the source region embedded in a generally high‐dimension travel‐time vector space. A cluster of events in the source region has an image of discrete points on the manifold, that, except in the simplest cases, cannot be viewed directly. However, it is possible to project the image of a cluster into a tangent space of the manifold for direct visualization. The projection operator can be computed directly from the data without a velocity model, but produces a distorted rendering of the cluster geometry. With a model we can predict the distortions and correct them to estimate cluster geometry. We develop these points with the simplest possible example, one for which direct visualization of the manifold is possible, using the example as an introduction to the relevant concepts from differential geometry in a familiar setting. The tangent space, a local linearization of the manifold, plays a key role. We develop a metric to estimate the limits of linearization, that is, to determine when the curvature of the manifold invalidates the linear assumption. We also examine the interplay of model error, inadequate network geometry, and pick error. We then generalize our results from the simple case to the general case of 3D source regions observed by general networks. Although we do suggest a new “project and correct” method for location, we do not develop it into a practical algorithm. In conclusion, our intention rather is to highlight new analytical methods grounded in differential geometry.

East Pacific Ocean Islands↗