Search NASA⌕ Search

SEARCH · Search NASA

Results for “resource partitioning”

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 OSTI2021

Influences on Discharge Partitioning on a Large River Delta: Case Study of the Mississippi-Atchafalaya Diversion, 1926-1950

The modern Mississippi River (M.R.) Delta is plumbed by the Mississippi and Atchafalaya rivers, setting water and sediment dispersal pathways for Earth’s fifth-largest river system. The Atchafalaya River’s (A.R.) partial annexation of discharge from the M.R., particularly between 1926 and 1950, prompted warnings of a rapid river avulsion and the construction of the Old River Control Structure to regulate flow. Natural and anthropogenic causes of this discharge annexation are difficult to disentangle. Here, we develop and validate a hydrodynamic model of flow partitioning through the historic channel network. We then isolate how several key changes to the system affected discharge partitioning and stage at the A.R.-M.R. diversion. Simulations show that erosion of the upper A.R. can account for 73% of the water discharge increase. Dredging in the lower A.R. between 1932 and 1950 can account for 35% of the water discharge increase, and was also an important control on shear stress distribution. The lower M.R. was slightly erosional during this period, and therefore hindered the A.R. discharge increase. Significant lacustrine delta deposition in A.R. had little effect on partitioning. These findings highlight the importance of A.R. enlargement processes on avulsion dynamics at this site. Given the essential nature of this river junction to the society, transportation, and commerce of the United States, improved attribution of discharge increases may lead to future management strategies that are broadly impactful.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Continental Scale Hydrostratigraphy: Basin‐Scale Testing of Alternative Data‐Driven Approaches

Abstract Integrated hydrological modeling is an effective method for understanding interactions between parts of the hydrologic cycle, quantifying water resources, and furthering knowledge of hydrologic processes. However, these models are dependent on robust and accurate datasets that physically represent spatial characteristics as model inputs. This study evaluates multiple data‐driven approaches for estimating hydraulic conductivity and subsurface properties at the continental‐scale, constructed from existing subsurface dataset components. Each subsurface configuration represents upper (unconfined) hydrogeology, lower (confined) hydrogeology, and the presence of a vertical flow barrier. Configurations are tested in two large‐scale U.S. watersheds using an integrated model. Model results are compared to observed streamflow and steady state water table depth (WTD). We provide model results for a range of configurations and show that both WTD and surface water partitioning are important indicators of performance. We also show that geology data source, total subsurface depth, anisotropy, and inclusion of a vertical flow barrier are the most important considerations for subsurface configurations. While a range of configurations proved viable, we provide a recommended Selected National Configuration 1 km resolution subsurface dataset for use in distributed large‐and continental‐scale hydrologic modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluating asynchronous Schwarz solvers on GPUs

With the commencement of the exascale computing era, we realize that the majority of the leadership supercomputers are heterogeneous and massively parallel. Even a single node can contain multiple co-processors such as GPUs and multiple CPU cores. For example, ORNL’s Summit accumulates six NVIDIA Tesla V100 GPUs and 42 IBM Power9 cores on each node. Synchronizing across compute resources of multiple nodes can be prohibitively expensive. Hence, it is necessary to develop and study asynchronous algorithms that circumvent this issue of bulk-synchronous computing. In this study, we examine the asynchronous version of the abstract Restricted Additive Schwarz method as a solver. We do not explicitly synchronize, but allow the communication between the sub-domains to be completely asynchronous, thereby removing the bulk synchronous nature of the algorithm. We accomplish this by using the one-sided Remote Memory Access (RMA) functions of the MPI standard. We study the benefits of using such an asynchronous solver over its synchronous counterpart. We also study the communication patterns governed by the partitioning and the overlap between the sub-domains on the global solver. Finally, we show that this concept can render attractive performance benefits over the synchronous counterparts even for a well-balanced problem.

Nayak, Pratik↗
DOE OSTI2022

Chemical Thermodynamic Modeling of Molten Salts to Support Off-Gas Abatement Systems

The reprocessing of used nuclear fuel by any means will liberate gaseous fission products, such as hydrogen ( 3 H), carbon ( 14 C), noble gases ( 85 Kr), and halogens ( 129 I), from the irradiated fuel. These elements will distribute through chemical processing operations and partition into process off-gas streams, and the specific volatile release fractions will be dictated by the chemical and physical properties of the system. A recent assessment found that there were significant knowledge gaps regarding the release quantities of volatile radionuclides from individual unit operations. These knowledge gaps limited the ability to determine what dedicated off-gas treatment technologies could be required for electrochemical-based reprocessing facilities. Unfortunately, experimental efforts to quantify release fractions are limited by the challenges associated with performing experiments using irradiated fuel. This report documents preliminary thermodynamic predictions of iodine and tritium release from chloride-based molten salts as part of an effort to better direct resources toward those experiments (both simulant and irradiated) that will be of the greatest impact. It was predicted that less than 0.5% of tritium was expected to be released and that nearly all of that amount would be released as H 2 . Thermochemical data for hydrogen (H 2 ) are of high fidelity, and no additional validation is recommended. Less than 0.05% of iodine was predicted to be released, with the primary volatile species being Cs 2 I 2 . Unlike the tritium predictions, the data underpinning the iodine release predictions are of low quality. It is recommended that a limited experimental program be dedicated to expanding the current physical property and thermodynamic property data for iodine in electrochemical processing and molten salt conditions, including vapor pressure measurements for the dimerized species predicted to comprise the majority of iodine release. Validating and improving the thermodynamic properties used in predictive modeling can reduce the need for expensive testing with irradiated fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Benchmarking quantum trial wavefunctions for phaseless auxiliary-field quantum Monte Carlo

The phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) method is a stochastic imaginary-time projection technique for computing ground-state properties of strongly correlated quantum systems, with accuracy that depends critically on the choice of trial wavefunction. Here, we investigate ph-AFQMC with trial states prepared using parameterized quantum circuits. In this work, we present a comprehensive benchmarking study of quantum trial wavefunctions spanning unitary coupled-cluster, Hamiltonian-informed, Jastrow-inspired, and adaptively constructed ansatze. The benchmarking evaluates accuracy, expressibility, and scalability of these ansatze within the QC-AFQMC framework. We test these ansatze on linear hydrogen chains under bond stretching and find that several ansatz families produce chemically accurate ph-AFQMC energies across the dissociation curve. We have performed simulations using the CUDA-Q quantum development platform on the GPU partition of the Perlmutter supercomputer. When comparing ansatze at similar numbers of variational parameters, we find that different ansatz families yield comparable ph-AFQMC results despite exhibiting substantially different variational energies, optimization costs, and circuit depths. Our results indicate that the variational energy of an ansatz is not always a reliable indicator of its quality for ph-AFQMC and reveal instances of over-parameterization. In the strongly correlated regime, trial wavefunctions obtained from adaptive ansatze, exemplified here by ADAPT-VQE with the UCCSD operator pool, can outperform their fixed-ansatz counterparts (UCCSD) in terms of projected energies while using substantially more compact circuits, providing a flexible route to optimize quantum resources within the ph-AFQMC framework.

Rofougaran, Rod [LBNL, Berkeley; Columbia U.; PNL,↗
DOE OSTI2021

Dimensionality reduction of the many-body problem using coupled-cluster subsystem flow equations: classical and quantum computing perspective

We discuss reduced-scaling strategies employing recently introduced sub-system embedding sub-algebras coupled-cluster formalism (SES-CC) to describe many-body systems. These strategies utilize properties of the SES-CC formulations where the equations describing certain classes of sub- systems can be integrated into a computational flows composed coupled eigenvalue problems of reduced dimensionality. Additionally, these flows can be defined at the level of the CC Ansatz defined by selected classes of cluster amplitudes, which define the wave function ”memory” of possible partitionings of the many-body system into constituent sub-systems. One of the possible ways of solving these coupled problems is through implementing procedures, where the information is passed between the sub-systems in a self-consistent manner. As a special case, we consider local flow formulations where the so-called local character of correlation effects can be closely related to properties of sub-system embedding sub-algebras employing localized molecular basis. We also generalize flow equations to the time domain and to downfolding methods utilizing double exponential unitary CC Ansatz (DUCC), where reduced dimensionality of constituent sub-problems offer a possibility of efficient utilization of limited quantum resources in modeling realistic systems.

Electron correlation, quantum chemistry, quantum c↗
DOE OSTI2024

HERO WEC 2024 Hydraulic Configuration Deployment Data

The following submission includes raw and processed data from the in water deployment of NREL's Hydraulic and Electric Reverse Osmosis Wave Energy Converter (HERO WEC), in the form of parquet files, TDMS files, CSV files, bag files and MATLAB workspaces. This dataset was collected in March 2024 at the Jennette's pier test site in North Carolina. This submission includes the following: - Data description document (HERO WEC FY24 Hydraulic Deployment Data Descriptions.doc) - This document includes detailed descriptions of the type of data and how it was processed and/or calculated. - Processed MATLAB workspace - The processed data is provided in the form of a single MATLAB workspace containing data from the full deployment. This workspace contains data from all sensors down sampled to 10 Hz along with all array Value Added Products (VAPs). - MATLAB visualization scripts - The MATLAB workspaces can be visualized using the file "HERO_WEC_2024_Hydraulic_Config_Data_Viewer.m/mlx". The user simply needs to download the processed MATLAB workspaces, specify the desired start and end times and run this file. Both the .m and .mlx file format has been provided depending on the user's preference. - Summary Data - The fully processed data was used to create a summary data set with averages and important calculations performed on 30-minute intervals to align with the intervals of wave resource data reported from nearby CDIP ocean observing buoys located 20km East of Jennette's pier and 40km Northeast of Jennette's pier. The wave resource data provided in this data set is to be used for reference only due the difference in water depth and proximity to shore between the Jennette's pier test site and the locations of the ocean observing buoys. This data is provided in the Summary Data zip folder, which includes this data set in the form of a MATLAB workspace, parquet file, and excel spreadsheet. - Processed Parquet File - The processed data is provided in the form of a single parquet file containing data from all HERO WEC sensors collected during the full deployment. Data in these files has been down sampled to 10 Hz and all array VAPs are included. - Interim Filtered Data - Raw data from each sensor group partitioned into 30-minute parquet files. These files are outputs from an intermediate stage of data processing and contain the raw data with no Quality Control (QC) or calculations performed in a format that is easier to use than the raw data. - Raw Data - Raw, unprocessed data from this deployment can be found in the Raw Data zip folder. This data is provided in the form of TDMS, CSV, and bag files in the original format output by the MODAQ system. - Python Data Processing Script - This links to an NREL public github repository containing the python script used to go from raw data to fully processed parquet files. Additional documentation on how to use this script is included in the github repository. This data set has been developed by the National Renewable Energy Laboratory, operated by Alliance for Sustainable Energy, LLC, for the U.S. Department of Energy (DOE) under Contract No. DE-AC36-08GO28308. Funding provided by the U.S. Department of Energy Office of Energy Efficiency and Renewable Energy Water Power Technologies Office.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

A Task Based Approach for Co-Scheduling Ensemble Workloads on Heterogeneous Nodes

Scientific workflows consist of multiple, connected applications, with data and results flowing from one to another in a pipeline. Traditionally, such workflows are executed in sequential order, storing intermediate data in storage disks. Co-scheduling application workflows concurrently on the same compute nodes would greatly reduce the cost of moving data to/from storage and allow real-time analysis of intermediate results. Nevertheless, most parallel programming runtimes do not allow seamless integration of various applications in a scientific workflow, in part due to the complexity of managing data and resources. The situation is even more complicated for heterogeneous systems. In this work we extend the Minos Computing Library (MCL) runtime to accelerate pipe-lined and parallel workloads where multiple applications are running in the same system. MCL’s asynchronous task library and runtime dynamically manages resources to allow co-scheduling of multiple processes sharing heterogeneous resources. In addition, we design a custom ex- tension of the Open Compute Language (OpenCL) to enable multiple processes to share device memory. We enable MCL to coordinate these shared buffers to allow for easy, fast data sharing between applications. Using malleable micro-benchmarks and two application workflows that combine scientific simulation and AI-based analysis, we show that our method outperforms traditional approaches.

Index Terms—Parallel systems, Scheduling and Task ↗
NASA NTRS1993

Volatility in the lunar crust: Trace element analyses of lunar minerals by PIXE proton microprobe

In situ determination of mineral compositions using microbeam techniques can characterize magma compositions through mineral-melt partitioning, and be used to investigate fine-grained or rare phases which cannot be extracted for analysis. Abundances of Fe, Mn, Sr, Ga, Zr, Y, Nb, Zn, Cu, Ni, Se, and Sb were determined for various mineral phases in a small number of lunar highlands rocks using the PIXE proton microprobe. Sr/Ga ratios of plagioclase and Mn/Zn ratios of mafic silicates show that the ferroan anorthosites and Mg-suite cumulates are depleted in volatile lithophile elements to about the same degree compared with chondrites and the Earth. This links the entire lunar crust to common processes or source compositions. In contrast, secondary sulfides in Descartes breccia clasts are enriched in chalcophile elements such as Cu, Zn, Ni, Se, and Sb, and represent a potential resource in the lunar highlands.

Norman, M. D.↗
DOE OSTI2021

Fast and Scalable Sparse Triangular Solver for Multi-GPU Based HPC Architectures

Designing efficient and scalable sparse linear algebra kernels on modern multi-GPU based HPC systems is a daunting task due to significant irregular memory references and workload imbalance across the GPUs. This is particularly the case for \textit{Sparse Triangular Solver (SpTRSV)} which introduces additional two-dimensional computation dependencies among subsequent computation steps. Dependency information is exchanged and shared among GPUs, thus warrant for efficient memory allocation, data partitioning, and workload distribution as well as fine-grained communication and synchronization support. In this work, we demonstrate that directly adopting unified memory can adversely affect the performance of SpTRSV on multi-GPU architectures, despite linking via fast interconnect like NVLinks and NVSwitches. Alternatively, we employ the latest NVSHMEM technology based on Partitioned Global Address Space programming model to enable efficient fine-grained communication and drastic synchronization overhead reduction. Furthermore, to handle workload imbalance, we propose a malleable task-pool execution model which can further enhance the utilization of GPUs. By applying these techniques, our experiments on the NVIDIA multi-GPU supernode V100-DGX-1 and DGX-2 systems demonstrate that our design can achieve on average 3.53x (up to 9.86x) speedup on a DGX-1 system and 3.66x (up to 9.64x) speedup on a DGX-2 system with 4-GPUs over the Unified-Memory design. The comprehensive sensitivity and scalability studies also show that the proposed zero-copy SpTRSV is able to fully utilize the computing and communication resources of the multi-GPU system.

Xie, Chenhao↗
DOE OSTI2023

Bark beetle impacts on forest evapotranspiration and its partitioning

Insect outbreaks affect forest structure and function and represent a major category of forest disturbance globally. However, the resulting impacts on evapotranspiration (ET), and especially hydrological partitioning between the abiotic (evaporation) and biotic (transpiration) components of total ET, are not well constrained. As a result, we combined remote sensing, eddy covariance, and hydrological modeling approaches to determine the effects of bark beetle outbreak on ET and its partitioning at multiple scales throughout the Southern Rocky Mountain Ecoregion (SRME), USA. At the eddy covariance measurement scale, 85 % of the forest was affected by beetles, and water year ET as a fraction of precipitation (P) decreased by 30 % relative to a control site, with 31 % greater reductions in growing season transpiration relative to total ET. At the ecoregion scale, satellite remote sensing masked to areas of >80 % tree mortality showed corresponding ET/P reductions of 9–15 % that occurred 6–8 years post-disturbance, and indicated that the majority of the total reduction occurred during the growing season; the Variable Infiltration Capacity hydrological model showed an associated 9–18 % increase in the ecoregion runoff ratio. Long-term (16–18 year) ET and vegetation mortality datasets extend the length of previously published analyses and allowed for clear characterization of the forest recovery period. During that time, transpiration recovery outpaced total ET recovery, which was lagged in part due to persistently reduced winter sublimation, and there was associated evidence of increasing late summer vegetation moisture stress. Overall, comparison of three independent methods and two partitioning approaches demonstrated a net negative impact of bark beetles on ET, and a relatively greater negative impact on transpiration, following bark beetle outbreak in the SRME.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Proteome reallocation enables the selective de novo biosynthesis of non-linear, branched-chain acetate esters

The one-carbon recursive ketoacid elongation pathway is responsible for making various branched-chain amino acids, aldehydes, alcohols, and acetate esters in living cells. Controlling selective microbial biosynthesis of these target molecules at high efficiency is challenging due to enzyme promiscuity, regulation, and metabolic burden. In this study, we present a systematic modular design approach to control proteome reallocation for selective microbial biosynthesis of branched-chain acetate esters. Through pathway modularization, we partitioned the branched-chain ester pathways into four submodules including keto-isovalerate submodule for converting pyruvate to keto-isovalerate, ketoacid elongation submodule for producing longer carbon-chain keto-acids, ketoacid decarboxylase submodule for converting ketoacids to alcohols, and alcohol acyltransferase submodule for producing branched-chain acetate esters by condensing alcohols and acetyl-CoA. By systematic manipulation of pathway gene replication and transcription, enzyme specificity of the first committed steps of these submodules, and downstream competing pathways, we demonstrated selective microbial production of isoamyl acetate over isobutyl acetate. In this work, we found that the optimized isoamyl acetate pathway globally redistributed the amino acid fractions in the proteomes and required up to 23-31% proteome reallocation at the expense of other cellular resources, such as those required to generate precursor metabolites and energy for growth and amino acid biosynthesis. From glucose fed-batch fermentation, the engineered strains produced isoamyl acetate up to a titer of 8.8 g/L (> 0.25 g/L toxicity limit), a yield of 0.22 g/g (61% of maximal theoretical value), and 86% selectivity, achieving the highest titers, yields and selectivity of isoamyl acetate reported to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1996

Design and Implementation of Replicated Object Layer

One of the widely used techniques for construction of fault tolerant applications is the replication of resources so that if one copy fails sufficient copies may still remain operational to allow the application to continue to function. This thesis involves the design and implementation of an object oriented framework for replicating data on multiple sites and across different platforms. Our approach, called the Replicated Object Layer (ROL) provides a mechanism for consistent replication of data over dynamic networks. ROL uses the Reliable Multicast Protocol (RMP) as a communication protocol that provides for reliable delivery, serialization and fault tolerance. Besides providing type registration, this layer facilitates distributed atomic transactions on replicated data. A novel algorithm called the RMP Commit Protocol, which commits transactions efficiently in reliable multicast environment is presented. ROL provides recovery procedures to ensure that site and communication failures do not corrupt persistent data, and male the system fault tolerant to network partitions. ROL will facilitate building distributed fault tolerant applications by performing the burdensome details of replica consistency operations, and making it completely transparent to the application.Replicated databases are a major class of applications which could be built on top of ROL.

Koka, Sudhir↗
DOE OSTI2025

Water Stable Isotopes in Precipitation, Rivers, and Groundwater Across an Elevation Gradient in the Sierra Nevada Mountains (USA) Reflect Source Elevation

Understanding watershed processes is critical to predict the impacts of climate change and forest management on water resources. However, collecting hydrological data in mountainous terrain is challenging. Precipitation, river water, and groundwater H and O stable isotope data can provide insights into processes occurring at the mountain range scale. Water δ 2 H and δ 18 O values in precipitation vary with terrain elevation; thus, the resulting isotopic lapse rates of precipitation, groundwater, and river water have the potential to elucidate watershed processes and source elevations of major rivers. We analysed H and O stable isotope data of precipitation, groundwater, and river water over the course of one Water Year (Oct 2016—Oct 2017) in the Sierra Nevada mountains of California, USA. We calculated elevation-dependent isotopic lapse rates of these waters to estimate the source elevation of major rivers draining the west flank of the Sierra Nevada mountains. We also investigated the Cosumnes River's watershed in more detail to determine how river flow may be more fully partitioned. Here, we found that H and O stable isotopes in precipitation are temporally variable, but isotopic lapse rates are generally consistent with prior studies. However, groundwater samples across an elevation gradient provide a more consistent and accessible isotopic lapse rate to predict river water source elevations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

TrustDER: Trusted, Private and Scalable Coordination of Distributed Energy Resources

In this project, the Stanford and SLAC Teams have developed a Trusted, Private and Scalable platform for coordinating Coordination of Distributed Energy Resources (TrustDER). This is a layered system that ensures private, trusted and scalable coordination and monitoring of DERs. It accommodates a variety of resources, such as solar generation, gensets and loads, with a particular focus on battery systems-based resources, as they are a transformational technology experiencing fast growth in adoption by large critical facilities. The platform can be used as standalone or added to existing aggregation systems to enable trust, privacy and resilience. TrustDER consists of layers that address each of the shortcomings of the existing state of the art. Each layer in the platform can operate independently but provides information to the layers above it to enable a novel form of overall coordination architecture. The project consists of several tasks, with each task dedicated to the design of each layer. Task 2 Resource Virtualization defined a software abstraction layer for distributed energy resources (DERs). The goal of this abstraction was to simplify the implementation of algorithms utilizing cooperation of DERs resources in a variety of use cases. Task 3 is on Secure ID for Asset Authentication. Identity Management Systems (IDMS) are a foundational infrastructure for interactions between entities (organizations, users, devices, and services). Secure ID is blockchain-based a distributed identity management system allowing (1) identity provisioning, (2) authentication, (3) authorization, and (4) identity data sharing for IoT-enabled assets on the electricity grid. In this project, the SLAC team focused on designing and testing Keymaker, a protocol for authenticating device identity managed by Secure ID. Task 5 Private and Safe Integration is focused on the design and evaluation of a DER cooperation scheme which allows for the aggregation of DERs without impacting network reliability. The approach is designed based on realistic assumptions regarding data availability, communication infrastructure limitations, and privacy. Task 6 Scalable Distributed Privacy for Information explored how virtualized batteries could be managed privately. Specifically, it examined the case in which a principal provides a partitioned battery to multiple clients. Task 7 Use Cases was to ensure that this technology was applied in relevant situations and scenarios. Primarily, this means that virtualization needed to be employed in a manner that either improved flexibility, bolstered security or privacy, or decreased costs.

25 ENERGY STORAGE↗
DOE OSTI2021

Local indistinguishability and incompleteness of entangled orthogonal bases: Method to generate two-element locally indistinguishable ensembles

Highlights: • Local indistinguishability of orthogonal quantum states. • Bipartite and multipartite quantum ensembles. • Unextendible entangled bases and uncompletable entangled bases. • Two-element locally indistinguishable quantum ensembles. • Multiparty unextendible entangled bases with unextendibility in all partitions. We relate the phenomenon of local indistinguishability of orthogonal states with the properties of unextendibility and uncompletability of entangled bases for bipartite and multipartite quantum systems. We prove that all two-qubit unextendible entangled bases are of size three and they cannot be perfectly distinguished by separable measurements. We identify a method of constructing two-element orthogonal ensembles, based on the concept of unextendible entangled bases, that can potentially lead to information sharing applications. Two-element ensembles form the fundamental unit of ensembles, and yet does not offer locally indistinguishable ensembles for pure state elements. Going over to mixed states does open this possibility, but can be difficult to identify. The method provided using unextendible entangled bases can be used for their systematic generation. In multipartite systems, we find a class of unextendible entangled bases for which the unextendibility property remains conserved across all bipartitions. We also identify nonlocal operations, local implementation of which require entangled resource states from a higher-dimensional quantum system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A Linear-Complexity Tensor Butterfly Algorithm for Compressing High-Dimensional Oscillatory Integral Operators

This paper presents a multilevel tensor compression algorithm called tensor butterfly algorithm for efficiently representing large-scale and high-dimensional oscillatory integral operators, including Green's functions for wave equations and integral transforms such as Radon transforms and Fourier transforms. The proposed algorithm leverages a tensor extension of the so-called complementary low-rank property of existing matrix butterfly algorithms. The algorithm partitions the discretized integral operator tensor into subtensors of multiple levels and factorizes each subtensor at the middle level as a Tucker-type interpolative decomposition, whose factor matrices are formed in a multilevel fashion. For a d-dimensional (d > 1) integral operator discretized into a 2d-mode tensor with n2d entries, the overall CPU time and memory requirement scale as O(nd), in stark contrast to the O(nd log n) complexity of existing matrix algorithms such as matrix butterfly algorithms and fast Fourier transforms (FFTs), where n is the number of points per direction. When comparing with other tensor algorithms such as quantized tensor train (QTT), the proposed algorithm also shows superior CPU and memory performance for tensor contraction. Remarkably, the tensor butterfly algorithm can efficiently model high-frequency Green's function interactions between two unit cubes, each spanning 512 wavelengths per direction, which represents problems of scale over 512× larger than that existing butterfly algorithms can handle, with the same amount of computation resources. On the other hand, for a problem representing 64 wavelengths per direction, which is the largest size existing algebraic matrix algorithms can handle, our tensor butterfly algorithm exhibits 200x speedups and 30× memory reduction compared with existing ones. Moreover, the tensor butterfly algorithm also permits O(nd)-complexity FFTs and Radon transforms up to d = 6 dimensions.

Kielstra, P Michael↗
DOE OSTI2024

Formation of a Large Cold Groundwater Mantle Helium Anomaly and High Temperature Geothermal Resources in Response to Bimodal Magmatism Near Roosevelt Hot Springs and Utah FORGE, Milford Valley, Southwest Utah

Abstract A large mantle helium anomaly and separate domains of high heat flow are the predominant manifestations of bimodal magmatic activity in the Milford valley. The mantle helium anomaly (1.9–2.6 R/Ra) covers 270 km 2 and is subdivided into two separated domains: a cold shallow groundwater regime and high temperature hydrothermal activity. The zone of anomalous heat flow covers >100 km 2 and is also subdivided into two adjacent domains, comprising hydrothermal activity at Roosevelt Hot Springs (RHS) (3–7 W/m 2 ) and conductive heat flow (100–180 mW/m 2 ). While the transfer of heat and mantle helium at RHS are coupled, heat and helium transfer are decoupled in the adjacent cold groundwater regime to the west. Both the mantle helium and geothermal anomalies are attributed to recent mafic‐felsic magmatic intrusions of >400 km 3 , however, the absence of volcanic eruptions <500,000 years indicates magmas stall before rising to shallow crustal level <10 km depth. Deep level magmatism produces a felsic composition melt, which is inferred to be responsible for the widespread and near uniform range of diluted mantle helium values. A thick and impermeable mass of crystalline granitic basement rock at the mid‐crustal level divides the ascent of mantle helium into separate flow paths. It may also impede the rise of buoyant magma trapping thermal energy that facilitates partial melting, slow cooling, and development of a thick thermal aureole. Partitioning of convective and conductive thermal regimes and independent flow paths supplying deeply derived helium characterize the development of a large long‐lived magma‐related geothermal system.

Geochemistry & Geophysics↗