Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed algorithms”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Towards FAIR Workflows for Federated Experimental Sciences

A de-centralized, peer-to-peer AI metadata framework is demonstrated which can enable end-to-end metadata & lineage tracking for distributed Machine Learning pipelines spanning edge, High Performance Computing, and cloud environments. With a specific example of end-to-end microscopy algorithm and datasets, the proposed method shows how to enable reproducibility, audit trail, provenance of metadata artifacts. The emerging needs of automation in experimental sciences, ML-centric workflows, and FAIR metadata management across federated compute environments is addressed.

machine learning↗
DOE OSTI2017

Distributed Estimation of All the Eigenvalues and Eigenvectors of Matrices Associated With Strongly Connected Digraphs

This letter considers the problem of estimating all the eigenvalues and eigenvectors of an irreducible matrix, corresponding to a strongly connected digraph, in the absence of knowledge on the global network topology. To this end, in this work we propose a unified distributed strategy performed by each node in the network and relies only on the local information. The key idea is to transform the nonlinear problem of computing both the eigenvalues and eigenvectors of an irreducible matrix into a linear one. Specifically, we first transform distributively the irreducible matrix into a nonsingular irreducible matrix. Each node in the network then estimates in a distributed fashion the inverse of the nonsingular matrix by solving a set of linear equations based on a consensus-type algorithm. The eigenvalues and the corresponding eigenvectors are finally computed by exploiting the relations between the eigenvalues and eigenvectors of both the inverse and the original irreducible matrices. A numerical example is provided to demonstrate the effectiveness of the proposed distributed strategy.

14 SOLAR ENERGY↗
DOE OSTI2018

Dynamic Power Distribution System Management With a Locally Connected Communication Network

Coordinated optimization and control of distribution-level assets can enable a reliable and optimal integration of massive amount of distributed energy resources (DERs) and facilitate distribution system management (DSM). Accordingly, the objective is to coordinate the power injection at the DERs to maintain certain quantities across the network, e.g., voltage magnitude, line flows, or line losses, to be close to a desired profile. By and large, the performance of the DSM algorithms has been challenged by two factors: i) the possibly non-strongly connected communication network over DERs that hinders the coordination; ii) the dynamics of the real system caused by the DERs with heterogeneous capabilities, time-varying operating conditions, and real-time measurement mismatches. In this paper, we investigate the modeling and algorithm design and analysis with the consideration of these two factors. In particular, a game theoretic characterization is first proposed to account for a locally connected communication network over DERs, along with the analysis of the existence and uniqueness of the Nash equilibrium (NE) therein. To achieve the equilibrium in a distributed fashion, a projected-gradient-based asynchronous DSM algorithm is then advocated. The algorithm performance, including the convergence speed and the tracking error, is analytically guaranteed under the dynamic setting. In conclusion, extensive numerical tests on both synthetic and realistic cases corroborate the analytical results derived.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

On Distributed PV Hosting Capacity Estimation, Sensitivity Study, and Improvement

This paper first studies the estimated distributed PV hosting capacities of seventeen utility distribution feeders using the Monte Carlo simulation based stochastic analysis, and then analyzes the sensitivity of PV hosting capacity to both feeder and photovoltaic system characteristics. Furthermore, an active distribution network management approach is proposed to maximize PV hosting capacity by optimally switching capacitors, adjusting voltage regulator taps, managing controllable branch switches and controlling smart PV inverters. The approach is formulated as a mixed-integer nonlinear optimization problem and a genetic algorithm is developed to obtain the solution. Multiple simulation cases are studied and the effectiveness of the proposed approach on increasing PV hosting capacity is demonstrated.

14 SOLAR ENERGY↗
DOE OSTI2011

Characterization and monitoring of subsurface processes using parallel computing and electrical resistivity imaging

This newsletter discusses recent advancement in subsurface resistivity characterization and monitoring capabilities. The BC Cribs field desiccation treatability test resistivity monitoring data is use an example to demonstrate near-real time 3D subsurface imaging capabilities. Electrical resistivity tomography (ERT) is a method of imaging the electrical resistivity distribution of the subsurface. An ERT data collection system consists of an array of electrodes, deployed on the ground surface or within boreholes, that are connected to a control unit which can access each electrode independently (Figure 1). A single measurement is collected by injecting current across a pair of current injection electrodes (source and sink), and measuring the resulting potential generated across a pair of potential measurement electrodes (positive and negative). An ERT data set is generated by collecting many such measurements using strategically selected current and potential electrode pairs. This data set is then processed using an inversion algorithm, which reconstructs an estimate (or image) of the electrical conductivity (i.e. the inverse of resistivity) distribution that gave rise to the measured data.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2012

SLUG-STOCHASTICALLY LIGHTING UP GALAXIES. I. METHODS AND VALIDATING TESTS

The effects of stochasticity on the luminosities of stellar populations are an often neglected but crucial element for understanding populations in the low-mass or the low star formation rate regime. To address this issue, we present SLUG, a new code to 'Stochastically Light Up Galaxies'. SLUG synthesizes stellar populations using a Monte Carlo technique that properly treats stochastic sampling including the effects of clustering, the stellar initial mass function, star formation history, stellar evolution, and cluster disruption. This code produces many useful outputs, including (1) catalogs of star clusters and their properties such as their stellar initial mass distributions and their photometric properties in a variety of filters, (2) two dimensional histograms of color-magnitude diagrams of every star in the simulation, and (3) the photometric properties of field stars and the integrated photometry of the entire simulated galaxy. After presenting the SLUG algorithm in detail, we validate the code through comparisons with STARBURST99 in the well-sampled regime, and with observed photometry of Milky Way clusters. Finally, we demonstrate SLUG's capabilities by presenting outputs in the stochastic regime. SLUG is publicly distributed through the Web site http://sites.google.com/site/runslug/.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Adaptive PID Gain Scheduling Control for Hydropower Turbine Using Neural CDE and Stochastic Distribution Shaping

This paper introduces a gain-scheduling PID controller design strategy for hydroturbine frequency control mode. This scheme first uses real data to learn the nonlinear dynamics of the hydroturbine using neural controlled differential equations and then perturbs the obtained nonlinear system at different equilibrium points, based on which a static output feedback adaptive dynamic programming algorithm is then used to optimize the PID gains for each equilibrium point. Moreover, a continuous-time version of stochastic distribution control is proposed to further fine-tune the optimized PID gains. Finally, the controller is obtained by implementing linear interpolation between the optimized PID control gains. The simulation results show that the proposed gain-scheduling PID controller can control a larger range of operation points compared with the given fixed PID controller and the baseline method. Compared with the given fixed PID controller, the proposed gain-scheduling PID controller can regulate hydroturbine frequency against disturbances induced by power-load variation with over 50% less overshoot for some operation points.

13 HYDRO ENERGY↗
DOE OSTI2016

Applying the new gamma ray imager diagnostic to measurements of runaway electron Bremsstrahlung radiation in the DIII-D Tokamak (invited)

A new gamma ray imager (GRI) is developed to probe the electron distribution function with 2D spatial resolution during runaway electron (RE) experiments at the DIII-D tokamak. The diagnostic is sensitive to 0.5–100 MeV gamma rays, allowing characterization of the RE distribution function evolution during RE growth and dissipation. The GRI consists of a lead “pinhole camera” mounted on the DIII-D midplane with 123 honeycombed tangential chords 20 cm wide that span the vessel interior. Up to 30 bismuth germanate (BGO) scintillation detectors capture RE bremsstrahlung radiation for Pulse Height Analysis (PHA) capable of discriminating up to 20 000 pulses per second. Digital signal processing routines combining shaping filters are performed during PHA to reject noise and record gamma ray energy. The GRI setup and PHA algorithms will be described and initial data from experiments will be presented. A synthetic diagnostic is developed to generate the gamma ray spectrum of a GRI channel given the plasma information and a prescribed distribution function. Magnetic reconstructions of the plasma are used to calculate the angle between every GRI sightline and orient and discriminate gamma rays emitted by a field-aligned RE distribution function.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-E-T-478: Sliding Window Multi-Criteria IMRT Optimization

Purpose: To demonstrate a method for what-you-see-is-what-you-get multi-criteria Pareto surface navigation for step and shoot IMRT treatment planning. Methods: We show mathematically how multiple sliding window treatment plans can be averaged to yield a single plan whose dose distribution is the dosimetric average of the averaged plans. This is incorporated into the Pareto surface navigation based approach to treatment planning in such a way that as the user navigates the surface, the plans he/she is viewing are ready to be delivered (i.e. there is no extra ‘segment the plans’ step that often leads to unacceptable plan degradation in step and shoot Pareto surface navigation). We also describe how the technique can be applied to VMAT. Briefly, sliding window VMAT plans are created such that MLC leaves paint out fluence maps every 15 degrees or so. These fluence map leaf trajectories are averaged in the same way the static beam IMRT ones are. Results: We show mathematically that fluence maps are exactly averaged using our leaf sweep averaging algorithm. Leaf transmission and output factor corrections effects, which are ignored in this work, can lead to small errors in terms of the dose distributions not being exactly averaged even though the fluence maps are. However, our demonstrations show that the dose distributions are almost exactly averaged as well. We demonstrate the technique both for IMRT and VMAT. Conclusions: By turning to sliding window delivery, we show that the problem of losing plan fidelity during the conversion of an idealized fluence map plan into a deliverable plan is remedied. This will allow for multicriteria optimization that avoids the pitfall that the planning has to be redone after the conversion into MLC segments due to plan quality decline. David Craft partially funded by RaySearch Laboratories.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Assessment of uncertainty in the numerical simulation of solar irradiance over inclined PV panels: New algorithms using measurements and modeling tools

Development of accurate transposition models to simulate plane-of-array (POA) irradiance from horizontal measurements or simulations is a complex process mainly because of the anisotropic distribution of diffuse solar radiation in the atmosphere. The limited availability of reliable POA measurements at large temporal and spatial scales leads to difficulties in the comprehensive evaluation of transposition models. This paper proposes new algorithms to assess the uncertainty of transposition models using both surface-based observations and modeling tools. We reviewed the analytical derivation of POA irradiance and the approximation of isotropic diffuse radiation that simplifies the computation. Two transposition models are evaluated against the computation by the rigorous analytical solution. We proposed a new algorithm to evaluate transposition models using the clear-sky measurements at the National Renewable Energy Laboratory's (NREL's) Solar Radiation Research Laboratory (SRRL) and a radiative transfer model that integrates diffuse radiances of various sky-viewing angles. We found that the radiative transfer model and a transposition model based on empirical regressions are superior to the isotropic models when compared to measurements. We further compared the radiative transfer model to the transposition models under an extensive range of idealized conditions. Our results suggest that the empirical transposition model has slightly higher cloudy-sky POA irradiance than the radiative transfer model, but performs better than the isotropic models under clear-sky conditions. Significantly smaller POA irradiances computed by the transposition models are observed when the photovoltaics (PV) panel deviates from the azimuthal direction of the sun. The new algorithms developed in the current study have opened the door to a more comprehensive evaluation of transposition models for various atmospheric conditions and solar and PV orientations.

14 SOLAR ENERGY↗
DOE OSTI2021

Communication-Avoiding and Memory-Constrained Sparse Matrix-Matrix Multiplication at Extreme Scale

Sparse matrix-matrix multiplication (SpGEMM) is a widely used kernel in various graph, scientific computing and machine learning algorithms. In this paper, we consider SpGEMMs performed on hundreds of thousands of processors generating trillions of nonzeros in the output matrix. Distributed SpGEMM at this extreme scale faces two key challenges: (1) high communication cost and (2) inadequate memory to generate the output. Furthermore, we address these challenges with an integrated communication-avoiding and memory-constrained SpGEMM algorithm that scales to 262,144 cores (more than 1 million hardware threads) and can multiply sparse matrices of any size as long as inputs and a fraction of output fit in the aggregated memory. As we go from 16,384 cores to 262,144 cores on a Cray XC40 supercomputer, the new SpGEMM algorithm runs 10x faster when multiplying large-scale protein-similarity matrices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

A Survey of Recent Advances in Transmission Network Switching

Optimal Transmission Switching for electric power networks is a research topic of great interest for both researchers and electric power engineers. Transmission switching is commonly considered a corrective mechanism tool for overloaded networks; however in recent years it has been sought as a means to reduce power loss in the network. Ultimately the objective of this analysis is to provide smart tools for power and economic dispatch planning, allowing utilities and engineers to utilize existing infrastructure as a means of automatic control. The optimization problem is a complex multi-objective search including criteria such as system power loss, equality/non-equality constraints, islanding management, and radial topology configuration in the case of distribution networks. Many optimization methods have been adopted for the transmission switching problem in past years using various heuristic search algorithms. The scope of methodology extends from tried and true methods such as simulated annealing to more modern methods like neural network optimization. Solution benchmarks under consideration commonly include solution state accuracy, computational burden and solution time. The objective of this presentation is to provide a comprehensive account of development in transmission switching methodologies. Benchmarking and testing strategies will be discussed as well as application of the Transmission Expansion Planning Tool in finding transmission switching solutions. Limited Discrepancy Search (LOS) and Probe Discrepancy Search (PDS) Algorithms are compared with more classic optimization strategies, as seen in defining literature and publications in this research field.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Load-Flow in Multiphase Distribution Networks: Existence, Uniqueness, Non-Singularity, and Linear Models

This paper considers unbalanced multiphase distribution systems with generic topology and different load models, and extends the Z-bus iterative load-flow algorithm based on a fixed-point interpretation of the AC load-flow equations. Explicit conditions for existence and uniqueness of load-flow solutions are presented. These conditions also guarantee convergence of the load-flow algorithm to the unique solution. The proposed methodology is applicable to generic systems featuring (i) wye connections; (ii) ungrounded delta connections; (iii) a combination of wye-connected and delta-connected sources/loads; and, (iv) a combination of line-to-line and line-to-grounded-neutral devices at the secondary of distribution transformers. Further, a sufficient condition for the non-singularity of the load-flow Jacobian is proposed. Lastly, linear load-flow models are derived, and their approximation accuracy is analyzed. Theoretical results are corroborated through experiments on IEEE test feeders.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Graphene K‐Tape Meshes for Densely Distributed Human Motion Monitoring

Abstract Wearable sensors that measure parameters associated with physical activity and bodily motions have been regarded as an indispensable tool for assessing personal wellness. Recent advances in nanocomposite strain sensors have been successfully used for monitoring skin strains and other strain‐derived physiological parameters. This study complements the existing body of work and presents a flexible, self‐adhering, fabric‐based wearable sensor for measuring skin strains and human motions. Graphene nanosheet thin films are directly spray‐coated onto kinesiology tape (K‐Tape) to obtain a self‐adhering strain sensor of high sensitivity and linearity. Their stable sensing performance and high repeatability are verified, while human subject tests confirm that they adequately capture muscle engagement during functional movements. In addition, densely distributed strain monitoring can be achieved using an electrical impedance tomography measurement approach and algorithm. Spatial strain sensing is successfully demonstrated using a single strip of piezoresistive K‐Tape, as well as when a K‐Tape network or mesh pattern is formed. The results show promise for using graphene K‐Tape meshes to measure how specific major muscle groups engage during different physical activities.

Lin, Yun‐An↗
DOE OSTI2024

Data-driven nonlocal model for fragmentation in the crushing of solids

A technique is proposed for reproducing particle size distributions in three-dimensional simulations of the crushing and comminution of solid materials. The method is designed to produce realistic distributions over a wide range of loading conditions, especially for small fragments. In contrast to most existing methods, the new model does not explicitly treat the small-scale process of fracture. Instead, it uses measured fragment distributions from laboratory tests as the basic material property that is incorporated into the algorithm, providing a data-driven approach. The algorithm is implemented within a nonlocal peridynamic solver, which simulates the underlying continuum mechanics and contact interactions between fragments after they are formed. Finally, the technique is illustrated in reproducing fragmentation data from drop weight testing on sandstone samples.

58 GEOSCIENCES↗
DOE OSTI2019

Turbulent flow and vortex characteristics in a blocked subchannel of a helically wrapped rod bundle

The study of flow and heat transfer through a bundle of rods is important, as this type of configuration is widely used in many engineering applications, for example, in heat exchangers, clusters of support structures, steam generators, and nuclear power reactors. This work experimentally investigated the flow-field characteristics of rod bundles with helically wrapped wires (wire-wrapped rod assemblies) proposed for advanced burner reactors. Subchannel blockage is one of the most interesting accident scenarios for liquid metal fast reactors as different mechanisms can lead to a partial or total subchannel blockage. It is of paramount important to understand the thermal-hydraulic behavior of the flow in the vicinity of such blockages. The flow mixing characteristics in two interior subchannels located near the core of the bundle, one without and one with the presence of a blockage, were experimentally investigated. The velocity fields within the spatial gaps of two subchannels were obtained by applying the matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) technique for Reynolds numbers of 4000 and 17,000. The first- and second-order statistics, including the mean velocity, root-mean-square fluctuating velocity, and Reynolds stress, were computed from the obtained TR-PIV velocity vector fields. In the blocked sub-channel, the TR-PIV results revealed a large recirculation flow region in the vicinity of the blockage that was created by the interaction between inflows from neighboring subchannels. Spectral analysis was performed and the vortex shedding frequency indicated by the Strouhal number was found to be St = 0.16 for the two Reynolds numbers investigated. The flow characteristics of the interior subchannels were investigated via the two-point cross-correlation of fluctuating velocities. In addition, proper orthogonal decomposition (POD) analysis was applied to the instantaneous vorticity fields to extract the statistically dominant flow structures. Lastly, the statistical characteristics of the vortex, such as the vortex populations, spatial distributions, and vortex strengths, were acquired by combining POD analysis and the vortex identification algorithm.

42 ENGINEERING↗
DOE OSTI2023

Long–short-term memory encoder–decoder with regularized hidden dynamics for fault detection in industrial processes

The ability of recurrent neural networks (RNN) to model nonlinear dynamics of high dimensional process data has enabled data-driven RNN-based fault detection algorithms. Previous studies have focused on detecting faults by identifying the discrepancies in data distribution between the faulty and normal data, as reflected in prediction errors generated by RNN models. However, in industrial processes, variations in data distribution can also result from changes in normal control setpoints and compensatory control adjustments in response to disturbances, making it hard to differentiate between normal and faulty conditions. This paper proposes a fault detection method utilizing a long short-term memory (LSTM) encoder–decoder structure with regularized hidden dynamics and reversible instance normalization (RevIN) to compactly represent high-dimensional measurements for effective monitoring. During training, the hidden states of the model are regularized to form a low-dimensional latent space representation of the original multivariate time series data. As a result, the prediction errors of the latent states can be used to monitor the abnormal dynamic variations, while the reconstruction errors of the measured variables are used to monitor the abnormal static variations. Furthermore, the proposed indices can reflect operating conditions, even when the distribution of test data changes, which helps distinguish faults from normal adjustments and disturbances that controllers can settle. Here, data from numerical simulation and the Tennessee Eastman process are used to illustrate the effectiveness of the proposed fault detection method.

42 ENGINEERING↗
DOE OSTI2025

Medium-induced modification of groomed and ungroomed jet mass and angularities in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The ALICE Collaboration presents a new suite of jet substructure measurements in Pb–Pb and pp collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV. These measurements provide access to the internal structure of jets via the momentum and angle of their constituents, probing how the quark–gluon plasma modifies jets, an effect known as jet quenching. Jet grooming additionally removes soft wide-angle radiation to enhance perturbative accuracy and reduce experimental uncertainties. We report the groomed and ungroomed jet mass m jet and jet angularities λ$^{κ}_{α}$ using κ = 1 and α > 0. Charged-particle jets are reconstructed at midrapidity using the anti-k T algorithm with resolution parameter R = 0.2. A narrowing of the jet mass and angularity distributions in Pb–Pb collisions with respect to pp is observed and is enhanced for groomed results, confirming modification of the jet core. By using consistent jet definitions and kinematic cuts between the mass and angularities for the first time, previous inconsistencies in the interpretation of quenching measurements are resolved, rectifying a hurdle for understanding how jet quenching arises from first principles and highlighting the importance of a well-controlled baseline. These results are compared with a variety of theoretical models of jet quenching, providing constraints on jet energy-loss mechanisms in the quark–gluon plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗