Search NASA⌕ Search

SEARCH · Search NASA

Results for “edge computing”

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

Using Machine Learning for Quantum Annealing Accuracy Prediction

Quantum annealers, such as the device built by D-Wave Systems, Inc., offer a way to compute solutions of NP-hard problems that can be expressed in Ising or quadratic unconstrained binary optimization (QUBO) form. Although such solutions are typically of very high quality, problem instances are usually not solved to optimality due to imperfections of the current generations quantum annealers. In this contribution, we aim to understand some of the factors contributing to the hardness of a problem instance, and to use machine learning models to predict the accuracy of the D-Wave 2000Q annealer for solving specific problems. We focus on the maximum clique problem, a classic NP-hard problem with important applications in network analysis, bioinformatics, and computational chemistry. By training a machine learning classification model on basic problem characteristics such as the number of edges in the graph, or annealing parameters, such as the D-Wave’s chain strength, we are able to rank certain features in the order of their contribution to the solution hardness, and present a simple decision tree which allows to predict whether a problem will be solvable to optimality with the D-Wave 2000Q. We extend these results by training a machine learning regression model that predicts the clique size found by D-Wave.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Fractional chiral hinge insulator

We propose and study a wave function describing an interacting three-dimensional fractional chiral hinge insulator (FCHI) constructed by Gutzwiller projection of two noninteracting second-order topological insulators with chiral hinge modes at half filling. We use large-scale variational Monte Carlo computations to characterize the model states via the entanglement entropy and charge-spin fluctuations. We show that the FCHI possesses fractional chiral hinge modes characterized by a central charge c = 1 and Luttinger parameter K = 1/2, like the edge modes of a Laughlin 1/2 state. The bulk and surface topology is characterized by the topological entanglement entropy (TEE) correction to the area law. While our computations indicate a vanishing bulk TEE, we show that the gapped surfaces host an unconventional two-dimensional topological phase. In a clear departure from the physics of a Laughlin 1/2 state, we find a TEE per surface compatible with (ln √2 )/2, half that of a Laughlin 1/2 state. This value cannot be obtained from topological quantum field theory for purely two-dimensional systems. For the sake of completeness, we also investigate the topological degeneracy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Multiobjective Dynamic Optimization for Optimal Load-Following of Natural Gas Combined Cycle Power Plants under Stress Constraints

With the increasing penetration of intermittent renewable energy sources into the electric grid, there is an associated need for conventional thermal power plants that are designed to operate at base-load conditions to cycle their load rapidly and frequently and operate under low-load conditions. Fast load-following operations lead to plant efficiency loss and the reduction of equipment life. A high-fidelity dynamic model of a natural gas combined cycle power plant, with rigorous equipment level submodels, is developed to capture the plant transient performance and the thermomechanical stress evolution at the high-pressure drum is computed to assess the drum life consumption. Stress evolution is modeled for the location of the edge at the drum/downcomer junction that experiences high circumferential stress amplitude during the fast load-following operation, leading to higher fatigue damage than locations considered under existing design standards. A dynamic optimization problem is solved for optimal load-following operation. Depending on the value of the stress constraint and the desired ramp rate, satisfying the desired stress constraint may be infeasible without relaxing the ramp rate. A multiobjective dynamic optimization problem is solved using a lexicographic approach that minimizes ramp rate relaxation and maximizes efficiency while satisfying stress constraints, avoiding spraying to saturation, and maintaining main steam and reheat steam temperature within bounds. It was observed that the optimal ramp rate can be highly nonlinear as opposed to the industry-standard linear ramp rate. Nonlinear optimal ramp-rate profiles not only help to avoid stress constraints that may be unavoidable by using the linear profile but also result in higher efficiency than the linear profile. Furthermore, the study shows that there is a strong tradeoff between the relaxation in the plant ramp rate and the time-average thermal efficiency.

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

Effect of aspect ratio on plasma response to resonant magnetic perturbations in tokamak devices

A systematic numerical study is carried out, computing and comparing the plasma response to the resonant magnetic perturbation (RMP) field, applied for controlling edge localized modes (ELMs), in a series of tokamak plasmas with varying aspect ratio and utilizing the MARS-F code. The aspect ratio is scanned either by varying the plasma major radius at a fixed minor radius or by varying the latter while fixing the former. Both approaches yield similar results when compared in terms of quantities with proper normalizations. In general, a non-monotonic dependence of the resonant response field (normalized by the vacuum counterpart) near the plasma edge is found with varying aspect ratio, indicating that a given ELM control coil current configuration strongly favors plasmas with a certain aspect ratio. This optimal aspect ratio, on the other hand, depends on the toroidal as well as poloidal (i.e., coil phasing) spectra of the applied RMP field. The equilibrium (edge) safety factor, the plasma shape, and the plasma toroidal flow are all fixed to ensure that the effects identified here are predominantly due to the plasma aspect ratio.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Topological data analysis of task-based fMRI data from experiments on schizophrenia

We use methods from computational algebraic topology to study functional brain networks, in which nodes represent brain regions and weighted edges represent similarity of fMRI time series from each region. With these tools, which allow one to characterize topological invariants such as loops in high-dimensional data, we are able to gain understanding into low-dimensional structures in networks in a way that complements traditional approaches based on pairwise interactions. In the present paper, we analyze networks constructed from task-based fMRI data from schizophrenia patients, healthy controls, and healthy siblings of schizophrenia patients using persistent homology, which allows us to explore the persistence of topological structures such as loops at different scales in the networks. We use persistence landscapes, persistence images, and Betti curves to create output summaries from our persistent-homology calculations, and we study the persistence landscapes and images using k-means clustering and community detection. Based on our analysis of persistence landscapes, we find that the members of the sibling cohort have topological features (specifically, their 1-dimensional loops) that are distinct from the other two cohorts. From the persistence images, we are able to distinguish all three subject groups and to determine the brain regions in the loops (with four or more edges) that allow us to make these distinctions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

A Scalable, Distribution Network-Aware, Customer Privacy-Preserving Framework for Operation of Virtual Power Plants

This poster presents a hierarchical control framework for a virtual power plant that leverages behind-the-meter resources for grid services while maintaining customer privacy during setpoint disaggregation. Unlike many existing approaches, the virtual power plant model uses a hierarchical control strategy and an iterative approach to determine the optimal set point dis-aggregation without direct load control while maintaining system-level power flow and voltage constraints. The proposed approach is numerically validated on a synthetic distribution feeder in San Francisco, demonstrating the ability of the framework to provide privacy-preserving virtual power plant services.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

ELM control optimization for various ITER scenarios based on linear and quasi-linear figures of merit

For the purpose of better understanding type-I edge localized mode (ELM) control in ITER with resonant magnetic perturbation (RMP) fields, the plasma response to RMP is computed by a resistive full magneto-hydrodynamic model in toroidal geometry. Five scenarios designed for ITER are considered, ranging from the pre-nuclear to nuclear phases. The plasma response to RMP is quantified by the plasma surface displacement near the X-point of the divertor plasma and at the outboard mid-plane. The optimal coil configurations between two high- Q deuterium-tritium (DT) scenarios (at the same plasma current of 15 MA and the same magnetic field of 5.3 T but different fusion gains, Q = 5 and 10) are predicted to be similar. For the other ITER scenarios with similar edge safety factor q 95 ~ 3 to that of the baseline scenario, the optimal coil phasing is also similar. The optimization results are different for a half-current full-field (7.5 MA/5.3 T) scenario, largely due to the difference in q 95 . The RMP coil currents are also optimized to tailor the core vs edge toroidal torques exerted by the 3D RMP fields on the plasma column. Torque optimization, with various objective functions proposed in the study, is useful for minimizing the side effects of RMP on the plasma core flow in ITER, while still maintaining the ELM control capability. Full utilization of three rows of ELM control coils in ITER is found to be essential to ensure both flexibility and robustness of ELM control, in terms of both linear and quasilinear plasma responses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Keeping LAMMPS cutting edge

Since its inception 30 years ago, LAMMPS has grown to be a world-class molecular dynamics code and a cornerstone of computational materials science research. This project aimed to keep LAMMPS at the forefront of molecular dynamics simulations by adapting LAMMPS to the latest developments in machine learning technology and hardware. Initially, the project set out to provide a unified implementation of active learning for efficient training data generation in LAMMPS, but the research trajectory pivoted to address more immediate and impactful opportunities. On the hardware side, recent record-breaking molecular dynamics simulations were developed on the Cerebras wafer-scale AI chip, and this project has developed an interface between LAMMPS and the hardware-specific molecular dynamics code to accelerate and simplify development and user adoption. On the software side, PyTorch’s Ahead-of-Time (AOT) compilation features promised increased performance for state-of-the-art equivariant neural network potentials, and this project laid the groundwork for their adoption in LAMMPS, resulting in a nearly 20x acceleration in extreme cases. Combined with a comprehensive benchmark study of LAMMPS across all current exascale systems, this project has reinforced LAMMPS’s role as a versatile, high-performance tool for current and future materials science applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Energy-Efficient and Resilient Infrastructure: Simulation, Validation, and Installation

Advanced, high-performance computing at the National Renewable Energy Laboratory (NREL) has enabled access to vast data resources with cutting-edge software techniques to understand, design, plan for, and maintain energy-efficient and resilient infrastructure. We have focused on cities and airports, but the technology we have developed will easily translate to seaports, inland ports, military installations, or other complex and large-scale energy-intensive systems. We can digitally simulate and explore current and future scenarios to make datadriven decisions for optimizing advanced energy systems, transportation and building operations, infrastructure planning and expansion, and battery storage to guide short- and long-term investments, electrification strategies, and integration of new technologies.

Athena↗
DOE OSTI2018

An approach to structural analysis using Moore-Shannon network reliability

We develop a formalism for approaching many questions dealing with network structure. The formalism is based on Shannon-Moore network reliability and on the Birnbaum definition of importance that uses it. The methods we suggest based on this formalism apply to a wide variety of dynamics on very general inter- action networks. They are designed to answer specific questions about how network structure affects particular aspects of the dynamics. We first introduce the concepts that make up the formalism, explain basic methods of using these concepts, and demonstrate applications to problems such as edge ranking.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

An Upper Bound on the Number of Edges of a Graph Whose kth Power Has a Connected Complement

We say that a graph is k-wide if for any partition of its vertex set into two subsets, one can choose vertices at distance at least k in these subsets (i.e., the complement of the kth power of this graph is connected). We say that a graph is k-mono-wide if for any partition of its vertex set into two subsets, one can choose vertices at distance exactly k in these subsets.We prove that the complement of a 3-wide graph on n vertices has at least 3n − 7 edges, and the complement of a 3-mono-wide graph on n vertices has at least 3n − 8 edges. We construct infinite series of graphs for which these bounds are attained.We also prove an asymptotically tight bound for the case k ≥ 4: the complement of a k-wide graph contains at least (n − 2k)(2k − 4[log{sub 2}k] − 1) edges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Community Project for Accelerator Science and Simulation-3

The focus of the Rensselaer efforts for COMPASS-3 was the development of parallel mesh adaptation procedures for high-order unstructured meshes where elements must be curved to properly approximate the domain geometry. The procedures developed were integrated into the ACE3P accelerator modeling code. Areas of development included: High order curved mesh representations that can provide the orders of geometric approximation needed. A curved mesh adaptation procedure that can match an adaptively defined anisotropic mesh metric field. In-memory integration of these procedures with the Omega3P code to provide an effective parallel adaptive simulation work flow. Improved dynamic load balancing tools that account for the degrees of freedom associated with the various order mesh entities (vertices, edges, faces and regions).

42 ENGINEERING↗
DOE OSTI2017

Model Order Reduction Algorithm for Estimating the Absorption Spectrum

The ab initio description of the spectral interior of the absorption spectrum poses both a theoretical and computational challenge for modern electronic structure theory. Due to the often spectrally dense character of this domain in the quantum propagator’s eigenspectrum for medium-to-large sized systems, traditional approaches based on the partial diagonalization of the propagator often encounter oscillatory and stagnating convergence. Electronic structure methods which solve the molecular response problem through the solution of spectrally shifted linear systems, such as the complex polarization propagator, offer an alternative approach which is agnostic to the underlying spectral density or domain location. This generality comes at a seemingly high computational cost associated with solving a large linear system for each spectral shift in some discretization of the spectral domain of interest. In this work, we present a novel, adaptive solution to this high computational overhead based on model order reduction techniques via interpolation. Model order reduction reduces the computational complexity of mathematical models and is ubiquitous in the simulation of dynamical systems and control theory. The efficiency and effectiveness of the proposed algorithm in the ab initio prediction of X-ray absorption spectra is demonstrated using a test set of challenging water clusters which are spectrally dense in the neighborhood of the oxygen K-edge. On the basis of a single, user defined tolerance we automatically determine the order of the reduced models and approximate the absorption spectrum up to the given tolerance. We also illustrate that, for the systems studied, the automatically determined model order increases logarithmically with the problem dimension, compared to a linear increase of the number of eigenvalues within the energy window. Furthermore, we observed that the computational cost of the proposed algorithm only scales quadratically with respect to the problem dimension.

IDREAM, K-Edge, theoretical, computational, electr↗
DOE OSTI2020

Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination

We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.

Lopez, Graham↗
DOE OSTI2017

An improved hyperbolic embedding algorithm

Because hyperbolic space has properties that make it amenable to graph representations, there is significant interest in scalable hyperbolic-space embedding methods. These embeddings enable constant-time approximation of shortest-path distances, and so are significantly more efficient than full shortest-path computations. In this article, we improve on existing landmark-based hyperbolic embedding algorithms for large-scale graphs. Whereas previous methods compute the embedding by using the derivative-free Nelder–Mead simplex optimization method, our approach uses the limited-memory BFGS (LBFGS) method, which is quasi-Newton optimization, with analytic gradients. Our method is not only significantly faster but also produces higher-quality embeddings. Moreover, we are able to include the hyperbolic curvature as a variable in the optimization. We compare our hyperbolic embedding method implementation in Python (called Hypy) against the best publicly available software, Rigel. Our method is an order of magnitude faster and shows significant improvements in the accuracy of the shortest-path distance calculations. Furthermore, tests are performed on a variety of real-world networks, and we show the scalability of our method by embedding a graph with 1.8 billion edges and 65 million nodes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Modal Substructuring of Geometrically Nonlinear Finite Element Models with Interface Reduction

Substructuring methods have been widely used in structural dynamics to divide large, complicated finite element models into smaller substructures. For linear systems, many methods have been developed to reduce the subcomponents down to a low order set of equations using a special set of component modes, and these are then assembled to approximate the dynamics of a large scale model. In this paper, a substructuring approach is developed for coupling geometrically nonlinear structures, where each subcomponent is drastically reduced to a low order set of nonlinear equations using a truncated set of fixedinterface and characteristic constraint modes. The method used to extract the coefficients of the nonlinear reduced order model (NLROM) is non-intrusive in that it does not require any modification to the commercial FEA code, but computes the NLROM from the results of several nonlinear static analyses. The NLROMs are then assembled to approximate the nonlinear differential equations of the global assembly. The method is demonstrated on the coupling of two geometrically nonlinear plates with simple supports at all edges. The plates are joined at a continuous interface through the rotational degrees-of-freedom (DOF), and the nonlinear normal modes (NNMs) of the assembled equations are computed to validate the models. The proposed substructuring approach reduces a 12,861 DOF nonlinear finite element model down to only 23 DOF, while still accurately reproducing the first three NNMs of the full order model.

42 ENGINEERING↗
DOE OSTI2023

Organo‐Functionalized Lacunary Double Cubane‐Type Oxometallates: Synthesis, Structure, and Properties of [(M II Cl) 2 (V IV O) 2 {((HOCH 2 CH 2 )(H)N(CH 2 CH 2 O))(HN(CH 2 CH 2 O) 2 )} 2 ] (M=Co, Zn)

Abstract Organofunctionalized tetranuclear clusters [(M II Cl) 2 (V IV O) 2 {((HOCH 2 CH 2 )(H)N(CH 2 CH 2 O))(HN(CH 2 CH 2 O) 2 )} 2 ] (1, M=Co,2: M=Zn) containing an unprecedented oxometallacyclic {M 2 V 2 Cl 2 N 4 O 8 } (M=Co, Zn) framework have been prepared by solvothermal reactions. The new oxo‐alkoxide compounds were fully characterized by spectroscopic methods, magnetic susceptibility measurement, DFT and ab initio computational methods, and complete single‐crystal X‐ray diffraction structure analysis. The isostructural clusters are formed of edge‐sharing octahedral {VO 5 N} and trigonal bipyramidal {MO 3 NCl} units. Diethanolamine ligates the bimetallic lacunary double cubane core of1and2in an unusual two‐mode fashion, unobserved previously. In the crystalline state, the clusters of1and2are joined by hydrogen bonds to form a three‐dimensional network structure. Magnetic susceptibility data indicate weakly antiferromagnetic interactions between the vanadium centers [J iso (V IV −V IV )=−5.4(1); −3.9(2) cm −1 ], and inequivalent antiferromagnetic interactions between the cobalt and vanadium centers [J iso (V IV −Co II )=−12.6 and −7.5 cm −1 ] contained in1.

Chemistry↗