Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed Processing”

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 397 records · Page 22

Distributed Multi-GPU Community Detection on Exascale Computing Platforms

Community detection is a fundamental operation in graph mining, and by uncovering hidden structures and patterns within complex systems it helps solve fundamental problems pertaining to social networks, such as information diffusion, epidemics, and recommender systems. Scaling graph algorithms for massive networks becomes challenging on modern distributed-memory multi-GPU (Graphics Processing Unit) systems due to limitations such as irregular memory access patterns, load imbalances, higher communication-computation ratios, and cross-platform support. We present a novel algorithm HiPDPL-GPU (Distributed Parallel Louvain) to address these challenges. We conduct experiments involving different partitioning techniques to achieve an optimized performance of HiPDPL-GPU on the two largest supercomputers: Frontier and Summit. Remarkably, HiPDPL-GPU processes a graph with 4.2 billion edges in less than 3 minutes using 1024 GPUs. Qualitatively, the performance of HiPDPL-GPU is similar or better compared to other state-of-the-art CPU- and GPU-based implementations. While prior GPU implementations have predominantly employed CUDA, our first-of-its-kind implementation for community detection is cross-platform, accommodating both AMD and NVIDIA GPUs.

Sattar, Naw Safrin↗

Deep kernel methods learn better: from cards to process optimization

Abstract The ability of deep learning methods to perform classification and regression tasks relies heavily on their capacity to uncover manifolds in high-dimensional data spaces and project them into low-dimensional representation spaces. In this study, we investigate the structure and character of the manifolds generated by classical variational autoencoder (VAE) approaches and deep kernel learning (DKL). In the former case, the structure of the latent space is determined by the properties of the input data alone, while in the latter, the latent manifold forms as a result of an active learning process that balances the data distribution and target functionalities. We show that DKL with active learning can produce a more compact and smooth latent space which is more conducive to optimization compared to previously reported methods, such as the VAE. We demonstrate this behavior using a simple cards dataset and extend it to the optimization of domain-generated trajectories in physical systems. Our findings suggest that latent manifolds constructed through active learning have a more beneficial structure for optimization problems, especially in feature-rich target-poor scenarios that are common in domain sciences, such as materials synthesis, energy storage, and molecular discovery. The Jupyter Notebooks that encapsulate the complete analysis accompany the article.

97 MATHEMATICS AND COMPUTING↗

The Utility of Macroecological Rules for Microbial Biogeography

Macroecological rules have been developed for plants and animals that describe large-scale distributional patterns and attempt to explain the underlying physiological and ecological processes behind them. Similarly, microorganisms exhibit patterns in relative abundance, distribution, diversity, and traits across space and time, yet it remains unclear the extent to which microorganisms follow macroecological rules initially developed for macroorganisms. Additionally, the usefulness of these rules as a null hypothesis when surveying microorganisms has yet to be fully evaluated. With rapid advancements in sequencing technology, we have seen a recent increase in microbial studies that utilize macroecological frameworks. Here, we review and synthesize these macroecological microbial studies with two main objectives: (1) to determine to what extent macroecological rules explain the distribution of host-associated and free-living microorganisms, and (2) to understand which environmental factors and stochastic processes may explain these patterns among microbial clades (archaea, bacteria, fungi, and protists) and habitats (host-associated and free living; terrestrial and aquatic). Overall, 78% of microbial macroecology studies focused on free living, aquatic organisms. In addition, most studies examined macroecological rules at the community level with only 35% of studies surveying organismal patterns across space. At the community level microorganisms often tracked patterns of macroorganisms for island biogeography (74% confirm) but rarely followed Latitudinal Diversity Gradients (LDGs) of macroorganisms (only 32% confirm). However, when microorganisms and macroorganisms shared the same macroecological patterns, underlying environmental drivers (e.g., temperature) were the same. Because we found a lack of studies for many microbial groups and habitats, we conclude our review by outlining several outstanding questions and creating recommendations for future studies in microbial ecology.

59 BASIC BIOLOGICAL SCIENCES↗

Opening Up Transactive Systems: Introducing TESS and Specification in a Field Deployment

Transactive energy systems (TS) use automated device bidding to access (residential) demand flexibility and coordinate supply and demand on the distribution system level through market processes. In this work, we present TESS, a modularized platform for the implementation of TS, which enables the deployment of adjusted market mechanisms, economic bidding, and the potential entry of third parties. TESS thereby opens up current integrated closed-system TS, allows for the better adaptation of TS to power systems with high shares of renewable energies, and lays the foundations for a smart grid with a variety of stakeholders. Furthermore, despite positive experiences in various pilot projects, one hurdle in introducing TS is their integration with existing tariff structures and (legal) requirements. In this paper, we therefore describe TESS as we have modified it for a field implementation within the service territory of Holy Cross Energy in Colorado. Importantly, our specification addresses challenges of implementing TS in existing electric retail systems, for instance, the design of bidding strategies when a (non-transactive) tariff system is already in place. We conclude with a general discussion of the challenges associated with “brownfield” implementation of TS, such as incentive problems of baseline approaches or long-term efficiency.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Predictions of the glass transition temperature and viscosity of organic aerosols from volatility distributions

Abstract. Volatility and viscosity are important properties of organic aerosols (OA), affecting aerosol processes such as formation, evolution, and partitioning of OA. Volatility distributions of ambient OA particles have often been measured, while viscosity measurements are scarce. We have previously developed a method to estimate the glass transition temperature (Tg) of an organic compound containing carbon, hydrogen, and oxygen. Based on analysis of over 2400 organic compounds including oxygenated organic compounds, as well as nitrogen- and sulfur-containing organic compounds, we extend this method to include nitrogen- and sulfur-containing compounds based on elemental composition. In addition, parameterizations are developed to predict Tg as a function of volatility and the atomic oxygen-to-carbon ratio based on a negative correlation between Tg and volatility. This prediction method of Tg is applied to ambient observations of volatility distributions at 11 field sites. The predicted Tg values of OA under dry conditions vary mainly from 290 to 339 K and the predicted viscosities are consistent with the results of ambient particle-phase-state measurements in the southeastern US and the Amazonian rain forest. Reducing the uncertainties in measured volatility distributions would improve predictions of viscosity, especially at low relative humidity. We also predict the Tg of OA components identified via positive matrix factorization of aerosol mass spectrometer (AMS) data. The predicted viscosity of oxidized OA is consistent with previously reported viscosity of secondary organic aerosols (SOA) derived from α-pinene, toluene, isoprene epoxydiol (IEPOX), and diesel fuel. Comparison of the predicted viscosity based on the observed volatility distributions with the viscosity simulated by a chemical transport model implies that missing low volatility compounds in a global model can lead to underestimation of OA viscosity at some sites. The relation between volatility and viscosity can be applied in the molecular corridor or volatility basis set approaches to improve OA simulations in chemical transport models by consideration of effects of particle viscosity in OA formation and evolution.

54 ENVIRONMENTAL SCIENCES↗

Production of b b ¯ at forward rapidity in p + p collisions at s = 510 GeV

The cross section of bottom quark-antiquark ( b b ¯ ) production in p + p collisions at s = 510 GeV is measured with the PHENIX detector at the Relativistic Heavy Ion Collider. The results are based on the yield of high mass, like-sign muon pairs measured within the PHENIX muon arm acceptance ( 1.2 < | y | < 2.2 ). The b b ¯ signal is extracted from like-sign dimuons by utilizing the unique properties of neutral B meson oscillation. We report a differential cross section of d σ b b ¯ → μ ± μ ± / d y = 0.16 ± 0.01 ( stat ) ± 0.02 ( syst ) ± 0.02 ( global ) nb for like-sign muons in the rapidity and p T ranges 1.2 < | y | < 2.2 and p T > 1 GeV / c , and dimuon mass of 5 – 10 GeV / c 2 . The extrapolated total cross section at this energy for b b ¯ production is 13.1 ± 0.6 ( stat ) ± 1.5 ( syst ) ± 2.7 ( global ) μ b . The total cross section is compared to a perturbative quantum chromodynamics calculation and is consistent within uncertainties. The azimuthal opening angle between muon pairs from b b ¯ decays and their p T distributions are compared to distributions generated using ps pythia 6, which includes next-to-leading order processes. The azimuthal correlations and pair p T distribution are not very well described by pythia calculations, but are still consistent within uncertainties. Flavor creation and flavor excitation subprocesses are favored over gluon splitting.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluation of UKESM aerosol size and composition using ATom measurements indicates missing marine aerosol formation mechanisms

Atmospheric aerosols influence climate through their interactions with radiation and clouds, yet large uncertainties remain in their simulation by global models. This study evaluates the United Kingdom Earth System Model version 1.1 (UKESM1.1) using global-scale aircraft observations from the Atmospheric Tomography (ATom) mission, focusing on aerosol lifecycle processes in the remote marine atmosphere. We assess model performance in simulating aerosol precursor vapours, number size distributions, chemical composition, and environmental conditions. Several process improvements are tested, including sulfuric acid-ammonia nucleation, ammonium nitrate scheme, methanesulfonic acid condensation, and low-temperature isoprene-derived secondary organic aerosol formation. Model biases differ significantly between the upper troposphere (UT) and the marine boundary layer (MBL). In the UT, UKESM1.1 overestimates nucleation and Aitken mode particles while underestimating accumulation mode, indicating insufficient growth. In the MBL, the model overestimates primary aerosols (e.g. seasalt) and precursor gases but underestimates nucleation and Aitken mode particles, even after incorporating updated nucleation and ammonium nitrate scheme. The persistence of low aerosol number concentrations, despite overestimated precursors, suggests missing formation pathways likely involving other species such as iodine, amines, and organic vapours. These limitations result in an unbalanced cloud condensation nuclei budget that over-relies on primary emissions. Sensitivity tests reveal that model outputs are strongly influenced by dimethyl sulfide emissions and vapour condensation schemes. Our results highlight the need for future model development to prioritise mechanistic representation of currently missing aerosol sources, rather than relying on empirical tuning, to improve aerosol-climate interaction estimates.

He, Xu-Cheng [Univ. of Cambridge (United Kingdom);↗

The role of supernova convection for the lower mass gap in the isolated binary formation of gravitational wave sources

ABSTRACT Understanding the astrophysical phenomena involving compact objects requires an insight about the engine behind the core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of the growing populations of the detected black hole–black hole, black hole–neutron star, and neutron star–neutron star mergers. The time-scale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. We adopt new formulas for the relation between the pre-SN star properties and their remnants and check how they impact the population of double compact object (DCO) mergers formed via the isolated binary evolution. The new formulas give one the ability to test a wide spectrum of assumptions on the convection growth time. In particular, the different variants allow for a smooth transition between having a deep mass gap and a remnant mass distribution filled by massive neutron stars and low-mass black holes. We present the distribution of masses, mass ratios, and the local merger rate densities of DCO for the different variants of new formulas and test them together with different approaches to other highly uncertain processes. We find that the mass distribution of DCO mergers is sensitive to the adopted assumption on the SN convection growth time-scale up to $m_1+m_2 \lesssim 35 \, \mathrm{M}_{\odot }$. Between the two extreme tested variants the probability of compact object formation within the mass gap may differ by up to approximately two orders of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum Sensing for Energy Applications: Review and Perspective

On its revolutionary threshold, quantum sensing is creating potentially transformative opportunities to exploit intricate quantum mechanical phenomena in new ways to make ultrasensitive measurements of multiple parameters. Concurrently, growing interest in quantum sensing has created opportunities for its deployment to improve processes pertaining to energy production, distribution, and consumption. Safe and secure utilization of energy is dependent upon addressing challenges related to material stability and function, secure monitoring of infrastructure, and accuracy in detection and measurement. A summary of well-established and emerging quantum sensing materials and techniques, as well as the corresponding sensing platforms that have been developed for their deployment is provided here. Specifically, the enhancement of existing advanced sensing technologies with quantum methods and materials is focused on, enabling the realization of an unprecedented level of sensitivity, placing an emphasis on relevance to the energy industry. The review concludes with a discussion on high-value applications of quantum sensing to the energy sector, as well as remaining barriers to sensor deployment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗

Record acceleration of the two-dimensional Ising model using a high-performance wafer-scale engine

The versatility and wide-ranging applicability of the Ising model, originally introduced to study phase transitions in magnetic materials, have made it a cornerstone in statistical physics and a valuable tool for evaluating the performance of emerging computer hardware. Here, we present a novel implementation of the two-dimensional Ising model on Cerebras Wafer-Scale Engine (WSE) – a revolutionary processor that is opening new frontiers in computing. In our deployment of the checkerboard algorithm, we optimized the Ising model to take advantage of the unique WSE architecture. Specifically, we employed a compressed bit representation storing 16 spins on each int16 word, and efficiently distributed the spins over the processing units enabling seamless weak scaling and limiting communications to only immediate neighboring units. Our implementation can handle up to 754 simulations in parallel, achieving an aggregate of over 61.8 trillion flip attempts per second for Ising models with up to 200 million spins. This represents a gain of up to 148 times over previously reported single-devices with a highly optimized implementation on NVIDIA V100 and up to 88 times in productivity compared to NVIDIA H100. Our findings highlight the significant potential of the WSE in scientific computing, particularly in the field of materials modeling.

Ising model↗

An immersed interface method for microstructure-scale electrochemical battery models: numerical formulation and performance portable implementation

We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.

25 ENERGY STORAGE↗

Life Cycle Analysis of Natural Gas Supply Chain and End Use Applications in the United States

Natural gas (NG) plays a crucial role in current and future energy systems in the United States due to its abundance and affordability. In this study a life cycle analysis of the NG supply chain in the United States was conducted using Argonne's R&D GREET model, examining stages from recovery to distribution using reported field data processed and documented by National Energy Technology Laboratory. Supply chain emissions were evaluated across multiple spatial scales, including national average, overall regional production, region-to-region, and basin-to-region scenarios. The GHG intensity of the U.S. average NG supply chain was estimated at 10.3 kg CO 2 e/MMBtu (lower heating value), with a range across regions from 7.8 kg CO 2 e/MMBtu (Northeast) to 15.1 kg CO 2 e/MMBtu (Pacific). The analysis further assessed how upstream NG emissions influence the life cycle GHG emissions of key end-use applications, including electricity generation (0.044–0.086 kg CO 2 e/kWh from upstream NG in combined cycle facilities), hydrogen production (1.04–2.20 kg CO 2 e/kg H 2 for steam methane reforming [SMR] and 1.06–2.23 kg CO 2 e/kg for autothermal reforming [ATR]), and transit bus operation utilizing compressed natural gas fuel (0.19–0.37 kg CO 2 e/mile) and hydrogen fuel (0.12–0.25 kg CO 2 e/mile for hydrogen produced in SMR and ATR).

compression↗

Optimization Modeling for Advanced Syngas to Olefin Reactive Systems

Reactor designs with mixed catalysts play an important role in transforming a multiple reactor system to single-shot reactors. In addition to savings in capital and ease of implementation, single-shot reactors are useful to break equilibrium limitations, thereby increasing the yield and selectivity of desired product as shown in previous studies. However, the nonlinear and highly exothermic nature of mixed-catalyst systems makes it difficult for commercial process simulation and optimization tools to optimize these systems. This study describes the development and application of optimization strategies for mixed-catalyst, single-shot reactors for syngas to olefin (STO) processes. Finding the optimal catalyst distribution is challenging and requires advanced solution strategies for singular optimal control problems, which are poorly conditioned and often lead to flat response surfaces. The graded bed and partial-moving finite-element approaches are used to find the optimal catalyst distribution that maximizes the olefins yield. A 1.3% increase in the yield is observed from one zone to three zones. The yield further improves from three zones to the exact infinite dimensional solution by 0.2%. This improvement can be realized in practice by changing only the catalyst distribution, without any extra investment. Lastly, the results suggest that a mixed-catalyst single shot reactor bed can be applied to other reaction mechanisms to increase reactor performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon Coating Influence on the Formation of Percolating Electrode Networks for Silicon Anodes

Previous studies have demonstrated that chemical vapor deposition carbon coating on silicon (Si@C) can enhance the electrochemical performance of lithium-ion batteries with Si-based anodes. However, the underlying mechanisms contributing to this improvement have not been fully explored. Here, we address this knowledge gap by applying a suite of characterization methods to evaluate Si@C anodes prepared by reducing acetylene on ball-milled Si particles. Raman mapping measurements show that the C coating (<5 nm thick) enables a homogeneous Si and carbon distribution during the slurry casting process, thereby promoting Si utilization during cycling. The coating’s microstructure and morphology were evaluated using X-ray photoelectron spectroscopy (XPS), scanning transmission electron microscopy, and neutron reflectivity experiments. Electrochemical impedance spectroscopy measurements upon cycling indicate that carbon coating also reduces the overall resistance as benchmarked against bare Si anodes. Galvanostatic cycling in half-cell studies revealed higher initial Coulombic efficiency and specific capacities with increasing carbon coating time. However, solid electrolyte interphase (SEI) investigations using XPS showed that the coated and uncoated samples have very similar characteristics, suggesting that the SEI may only play a minor role in enhancing the performance of Si@C. Full-cell evaluation of the Si electrodes was consistent with half-cell results relating to performance and SEI properties, further supporting the conclusion that electronic and ionic percolation, enabled by effective electrode manufacturing, are the dominant factors contributing to the favorable performance of Si@C.

25 ENERGY STORAGE↗

The Illumination of Thunderclouds by Lightning: 3. Retrieving Optical Source Altitude

Abstract Optical space‐based lightning sensors such as the Geostationary Lightning Mapper (GLM) detect and geolocate lightning by recording rapid changes in cloud top illumination. While lightning locations can be determined to within a pixel on the GLM imaging array, these instruments are not individually able to natively report lightning altitude. It has previously been shown that thunderclouds are illuminated differently based on the altitude of the optical source. In this study, we examine how altitude information can be extracted from the spatial distributions of GLM energy recorded from each optical pulse. We match GLM “groups” with Lightning Mapping Array (LMA) source data that accurately report the 3‐D positions of coincident Radio‐Frequency (RF) emitters. We then use machine learning methods to predict the mean LMA source altitudes matched to GLM groups using metrics from the optical data that describe the amplitude, breadth, and texture of the group spatial energy distribution. The resulting model can predict the LMA mean source altitude from GLM group data with a median absolute error of <1.5 km, which is sufficient to determine the location of the charge layer where the optical energy originated. This model is able to capture changes to the source altitude distribution in response to convective processes in the thunderstorm, and the GLM predictions can reveal the vertical structure of individual flashes ‐ enabling 3‐D flash geolocation with GLM for the first time. Future work will account for differences in thunderstorm charge/precipitation structures and viewing angle across the GLM Field of View.

54 ENVIRONMENTAL SCIENCES↗

One-Parameter Analytical Derivation in Modified Budyko Framework for Unsteady-State Streamflow Elasticity in Humid Catchments

The changes in climate and catchment properties have altered the hydrological processes significantly at different spatiotemporal scales around the world. In particular for finer time scales, changes in water storage, which has been commonly neglected for long-term temporal scales, may play an important role on hydrological processes. Nevertheless, few studies addressed this question in modifying the Budyko framework, with most of them remaining a steady-state assumption on catchment characteristics. Here, in this work, we derive new analytical formulas of unsteady-state streamflow elasticity in a modified Budyko framework, incorporating both storage change and one specific parameter for catchment properties. We study 78 humid catchments in the USA with simulation data of daily hydrological processes from a probability-distribution-based hydrological model (the modified HyMOD). As indicated by results, the annual storage change ratio is linearly correlated with the annual aridity index, and this relationship can be used to estimate elasticity coefficients with our formulas. The estimated elasticity coefficients perform well in simulating the annual streamflow with the power-law model. For different catchments, variability of the unsteady-state elasticity is higher than that of the steady-state elasticity. Unsteady-state streamflow coefficients show significant linear correlation with catchment properties, such as the average slope, average elevation, and catchment area. This study provides a new analytical approach to investigate the interannual stability of catchments with varying climate and catchment properties.

54 ENVIRONMENTAL SCIENCES↗