Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed analytics”

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 37 records · Page 2

Cutting Corners: Curvilinear-Surface-Based Gravity Models for Asteroids and Comets

We report that asteroids and comets are often highly irregular in shape and can contain large density inhomogenities. For such bodies, the simple point-mass gravity model does not capture the dynamics of the environment and this has driven a need for higher fidelity alternatives. A number of approaches have been developed with one of the most popular being the analytic polyhedral model. The analytic polyhedral model approximates the irregular shaped body as a polyhedron, typically consisting of triangular facets, with some assumed internal density distribution. Several analytic formulations have been derived for the gravitational fields of a constant-density polyhedron. A thorough history of the model is provided in Ref. [6] as well as a comparison of the implementations of Refs. [2–4]. The three versions provide near identical accuracy with the implementation of Werner requiring the fewest transcendental function evaluations suggesting its superior efficiency. Implementations with linear density contrasts and arbitrary polynomial density contrasts have also been developed.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear sparse Bayesian learning for physics-based models

This paper addresses the issue of overfitting while calibrating unknown parameters of over-parameterized physics-based models with noisy and incomplete observations. Here, a semi-analytical Bayesian framework of nonlinear sparse Bayesian learning (NSBL) is proposed to identify sparsity among model parameters during Bayesian inversion. NSBL offers significant advantages over machine learning algorithm of sparse Bayesian learning (SBL) for physics-based models, such as 1) the likelihood function or the posterior parameter distribution is not required to be Gaussian, and 2) prior parameter knowledge is incorporated into sparse learning (i.e. not all parameters are treated as questionable). NSBL employs the concept of automatic relevance determination (ARD) to facilitate sparsity among questionable parameters through parameterized prior distributions. The analytical tractability of NSBL is enabled by employing Gaussian ARD priors and by building a Gaussian mixture-model approximation of the posterior parameter distribution that excludes the contribution of ARD priors. Subsequently, type-II maximum likelihood is executed using Newton's method whereby the evidence and its gradient and Hessian information are computed in a semi-analytical fashion. We show numerically and analytically that SBL is a special case of NSBL for linear regression models. Subsequently, a linear regression example involving multimodality in both parameter posterior pdf and model evidence is considered to demonstrate the performance of NSBL in cases where SBL is inapplicable. Next, NSBL is applied to identify sparsity among the damping coefficients of a mass-spring-damper model of a shear building frame. These numerical studies demonstrate the robustness and efficiency of NSBL in alleviating overfitting during Bayesian inversion of nonlinear physics-based models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Performance Evaluation of Python Based Data Analytics Frameworks in Summit: Early Experiences

The explosion in the volumes of data generated from ever-larger simulation campaigns and experiments or observations necessitates competent tools for data wrangling and analysis). While the Oak Ridge Leadership Computing Facility (OLCF) provides a variety of tools to perform data wrangling and data analysis tasks, Python based tools often lack scalability, or the ability to fully exploit the computational capability of OLCF’s Summit supercomputer. NVIDIA RAPIDS and Dask offer a promising solution to accelerate and distribute data analytics workloads from personal computers to heterogeneous supercomputing systems. We discuss early performance evaluation results of RAPIDS and Dask on Summit to understand their capabilities, scalability, and limitations. Our evaluation includes a subset of RAPIDS libraries, i.e., cuDF, cuML, and cuGraph, and Chainer’s CuPy, and their multi-GPU variants when available.We also draw on the observed trends from the performance evaluation results to discuss best practices for maximizing performance.

Hernandez Arreguin, Benjamin↗

Pd/C promotes C–H bond activation and oxidation of p -hydroxybenzoate during hydrogenolysis of poplar

Hydrogenolysis of lignin generates a portfolio of products, the yields of which are generally calculated using a subset of phenolic monomers that are dependent on the lignin composition, product distribution, and analytical technique. Some lignins are naturally γ-acylated; poplar lignins, for example, have p-hydroxybenzoate groups on 1–15% of their syringyl subunits. Upon hydrogenolysis, it is generally assumed that the p-hydroxybenzoate is cleaved before the deacylated lignin is depolymerized. Hydrogenolysis of model γ-p-hydroxybenzoylated β-aryl ethers do not, however, produce the deacylated β-aryl ether intermediates, as was previously conjectured; products instead derive from palladium-assisted reactions on the cinnamyl p-hydroxybenzoates resulting in initial β-ether cleavage. The p-hydroxybenzoate moiety itself also undergoes carboxylate-assisted palladium-catalyzed C–H bond activation to form the 2,4-dihydroxybenzoate, that subsequently converts to the 2,4-dihydroxycyclohex-1-enoate. Furthermore, these details underscore previously unrecognized pathways and products that are key to understanding the different hydrogenolysis product distributions from naturally acylated lignins that are prevalent biomass-conversion feedstocks.

09 BIOMASS FUELS↗

Condensation and evaporation of boson stars

Abstract Axion-like particles, including the QCD axion, are well-motivated dark matter candidates. Numerical simulations have revealed coherent soliton configurations, also known as boson stars, in the centers of axion halos. We study evolution of axion solitons immersed into a gas of axion waves with Maxwellian velocity distribution. Combining analytical approach with controlled numerical simulations we find that heavy solitons grow by condensation of axions from the gas, while light solitons evaporate. We deduce the parametric dependence of the soliton growth/evaporation rate and show that it is proportional to the rate of the kinetic relaxation in the gas. The proportionality coefficient is controlled by the product of the soliton radius and the typical gas momentum or, equivalently, the ratio of the gas and soliton virial temperatures. We discuss the asymptotics of the rate when this parameter is large or small.

Physics↗

Dihadron correlations in small- x DIS at NLO: transverse momentum dependent fragmentation

We compute the inclusive dihadron cross-section in Deep Inelastic Scattering at next-to-leading order (NLO) and small x in the Color Glass Condensate. We focus on the kinematic limit where the hadrons are produced at forward rapidities (in the direction of the virtual photon) and back-to-back in the transverse plane. Our calculation demonstrates that the coefficient of the Sudakov double logarithm for this process is –$\frac{αs}{2π}$ [C F + $\frac{Nc}{2}$] instead of –$\frac{αsNc}{4π}$ when back-to-back jets are measured in the final state. To preserve the universality of the Sudakov soft factor associated with the Weizsäcker-Williams transverse momentum dependent (TMD) gluon distribution, we promote the collinear fragmentation functions into TMD fragmentation functions. We then perform the resummation of the Sudakov logarithms through Collins-Soper-Sterman evolution of the TMD fragmentation functions and the Weizsäcker-Williams TMD gluon distribution. Finally, analytic expressions are obtained for the NLO coefficient functions in the $\overline{MS}$-scheme. These results pave the way towards numerically calculating dihadron correlations at small x at the future Electron-Ion Collider with full NLO accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High-resolution particle-in-cell simulations of two-dimensional Bernstein–Greene–Kruskal modes

We present two-dimensional (2D) particle-in-cell (PIC) simulations of 2D Bernstein–Greene–Kruskal modes, which are exact nonlinear steady-state solutions of the Vlasov–Poisson equations, on a 2D plane perpendicular to a background magnetic field, with a cylindrically symmetric electric potential localized on the plane. PIC simulations are initialized using analytic electron distributions and electric potentials from the theory. We confirm the validity of such solutions using high-resolutions up to a 20482 grid. We show that the solutions are dynamically stable for a stronger background magnetic field, while keeping other parameters of the model fixed, but become unstable when the field strength is weaker than a certain value. When a mode becomes unstable, we observe that the instability begins with the excitation of azimuthal electrostatic waves that ends with a spiral pattern.

Physics↗

The importance of accounting for the Tolman correction to surface tension for nucleation and growth modeling of Fe clusters

Gibbs free energies of clusters are required for predictive modeling of cluster growth during condensation of a cooling vapor. Here, we present a straightforward method of calculating free energies of cluster formation using the data from molecular dynamics (MD) simulations. We apply this method to iron clusters having from 2 to 100 atoms. The energies obtained are verified by comparing to an MD-simulated equilibrium cluster size distribution in a sub-saturated vapor. We show that these free energies differ significantly from those obtained with a commonly used spherical cluster approximation, which relies on a surface tension coefficient of a flat surface, as it is used in the classical nucleation theory (CNT). We show that the spherical cluster approximation in CNT can be improved by using a cluster-size-dependent Tolman correction for the surface tension. The Tolman length and effective surface tension values were derived for iron clusters, and they significantly differ from the commonly used experimentally measured values. This improved approximation does not account for geometric magic number effects responsible for spikes and troughs in densities of neighbor cluster sizes. Nonetheless, it allows to more accurately model cluster formation from a cooling vapor. It better reproduces the condensation timeline, overall shape of the cluster size distribution, average cluster size, and the distribution width. In contrast, using a constant surface tension coefficient (as done in CNT) resulted in incorrect condensation dynamics and cluster size distributions. The analytical expression for cluster nucleation rate from CNT was updated to account for the size-dependence of cluster surface tension.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Flow annealed importance sampling bootstrap meets differentiable particle physics

High-energy physics requires the generation of large numbers of simulated data samples from complex but analytically tractable distributions called matrix elements. Surrogate models, such as normalizing flows, are gaining popularity for this task due to their computational efficiency. We adopt an approach based on flow annealed importance sampling bootstrap (FAB) that evaluates the differentiable target density during training and helps avoid the costly generation of training data in advance. We show that FAB reaches higher sampling efficiency with fewer target evaluations in high dimensions in comparison to other methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Machine Learning-Based Vulnerability Analysis for Cascading Failures of Integrated Power-Gas Systems

This article proposes a cascading failure simulation (CFS) method and a hybrid machine learning method for vulnerability analysis of integrated power-gas systems (IPGSs). The CFS method is designed to study the propagating process of cascading failures between the two systems, generating data for machine learning with initial states randomly sampled. The proposed method considers generator and gas well ramping, transmission line and gas pipeline tripping, island issue handling and load shedding strategies. Then, a hybrid machine learning model with a combined random forest (RF) classification and regression algorithms is proposed to investigate the impact of random initial states on the vulnerability metrics of IPGSs. Extensive case studies are carried out on three test IPGSs to verify the proposed models and algorithms. Simulation results show that the proposed models and algorithms can achieve high accuracy for the vulnerability analysis of IPGSs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Flattening the Radial Temperature Profile across the Transformational Challenge Reactor Core

The Transformational Challenge Reactor (TCR) program is demonstrating an agile development approach to advanced nuclear reactor design, which has traditionally utilized a linear design process. In leveraging artificial intelligence, additive manufacturing, advanced materials, and cutting-edge modeling and simulation, the TCR program aims to minimize the high cost and lengthy deployment timelines now standard in the nuclear industry. Within a relatively short period of time, a robust and mature advanced gas-cooled reactor was iteratively designed under the TCR program, using these cutting-edge technologies. The TCR is a 3 MWt gas-cooled microreactor fueled with uranium nitride (UN) tristructural isotropic (TRISO) fuel particles. Though manufactured via traditional means, these UN TRISO particles are loaded into additively manufactured silicon carbide (SiC) cans [4]. Once loaded with TRISO particles, the SiC cans are densified using a chemical vapor infiltration process. The additively manufactured SiC enables significantly more freedom in the design of the fuel form than could ever be achieved using traditionally manufactured SiC. The helium coolant, pressurized to 5 MPa, enters the core at 300°C and nominally exits it at 500°C. Typically, the most thermally limiting components in any reactor design are the fuel assemblies in the core center. To provide a wide thermal margin in these central fuel assemblies, the flow may be biased toward the center of the core to more effectively cool these fuel assemblies with more power deposition and flatten the core’s radial temperature distribution. An analytical fluid model of the TCR core was developed to explore methods for biasing the flow away from the cooler outer fuel assemblies and towards the hotter inner ones. Higher-fidelity models developed in STAR-CCM+ 2020.3.1, a computational fluid dynamics code, were then utilized to verify the analytical model’s findings.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Resilience Analytics for Interdependent Power and Water Distribution Systems

The occurrence of large-scale outages in power distribution systems (PDSs) caused, for instance, by extreme weather events has raised concerns on the cascading impact of the power outages on the resilience of other critical infrastructures, including water distribution systems (WDSs), whose healthy operation highly relies on the availability of power from PDSs. Here, this paper proposes an analytics model for quantifying the interdependence between the resilience of power and water distribution systems. The proposed model first performs spatio-temporal outage analysis on the PDS and WDS to determine the PDS power serving capability and the resulting cascading impact on the WDS water serving capability during the occurrence of disruptions. The resilience interdependence is then determined by a proposed set of six metrics calculated using the power and water serving capability profiles, which trace and quantify the cascading impacts of power outages in the operation of WDSs during the degradation and recovery stages. The proposed analytics model is implemented on the IEEE 33-bus power distribution system supplying power for pumping stations of a 15-node test WDS, considering multiple cases of water storage capacity, single and multiple outage scenarios, and availability of distributed generation. The numerical results demonstrate that the proposed metrics track the temporal and spatial interdependence between the resilience of power and water distribution systems, providing a mechanism to identify weak spots, and assess options to enhance the resilience of the critical water infrastructure.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Extremely Scalable Distributed Computation of Contour Trees via Pre-Simplification

Contour trees offer an abstract representation of the level set topology in scalar fields and are widely used in topological data analysis and visualization. However, applying contour trees to large-scale scientific datasets remains challenging due to scalability limitations. Recent developments in distributed hierarchical contour trees have addressed these challenges by enabling scalable computation across distributed systems. Building on these structures, advanced analytical tasks—such as volumetric branch decomposition and contour extraction—have been introduced to facilitate large-scale scientific analysis. Despite these advancements, such analytical tasks substantially increase memory usage, which hampers scalability. In this paper, we propose a pre-simplification strategy to significantly reduce the memory overhead associated with analytical tasks on distributed hierarchical contour trees. We demonstrate enhanced scalability through strong scaling experiments, constructing the largest known contour tree—comprising over half a trillion nodes with complex topology—in under 15 minutes on a dataset containing 550 billion elements.

Li, Mingzhe [University of Utah]↗

Metrics and Analytical Frameworks for Valuing Energy Efficiency and Distributed Energy Resources in the Built Environment: Preprint

This paper summarizes efforts to develop new—and enhance existing—analytical frameworks and metrics to quantify the value that grid-interactive efficient homes with solar (GEB-solar homes) can provide. Industry is working to characterize and understand these capabilities and benefits, but existing analytical frameworks for evaluating energy efficiency (EE) are often siloed from those that evaluate distributed energy resources (DERs). Five metrics were adapted from an extensive literature review and applied to case studies of a modeled home in Riverside, California: ramp up/down, cover factor demand/supply, and curtailable load. Eight different technology scenarios were analyzed using a more tightly connected suite of building-to-grid models (BEopt, REopt, ReEDS, and PLEXOS). Additionally, an initial version of the Cambium tool was developed, characterizing the marginal prices and emissions from NREL’s 2018 Standard Scenarios. These grid costs were extended to a time-varying proxy retail rate and applied as part of a new grid alignment metric. In the results, a more integrated combination of GEB-solar technologies led to a higher cover factor demand—the percentage of gross home load covered by on-site solar—however, a benchmark was required to determine what range of cover factor was “best” for given grid conditions. To that end, a grid alignment cost metric was applied to the case study scenarios. The average cost to serve the net load of the home decreased from a median of ~$0.24/kWh to ~$0.10/kWh when the most integrated technology scenario was optimized towards the grid pricing proxy versus the time-of-use (TOU) rate.

41 EE - Solar Energy Technologies Office (EE-4S)↗