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.

SourceResearch connectionFrom yearThrough year

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

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Auto-adaptive Harris corner detection algorithm based on entropy-improved block processing

Extracting well distributed control points (CPs) is a very challenging task for remote sensing image registration, particularly for large high-resolution images over heterogeneous landscape. Based on image analysis such as edge detection, corner detection, and information theory, a new CP detection approach is proposed to select high- quality, evenly distributed CPs. The Entropy-Block-Based variant of the Harris Corner Detector (EBB-HCD) is achieved by dividing the image into blocks and by allocating the number of CP's based upon the entropy of each block. While the block-based strategy improves the CP balance problem, a factor calculated from entropy avoids overdetection. We conducted a comparison study utilizing the well-known Harris Corner Detector (HCD) and an implementation of the Block-Based Harris Corner Detector (BB-HCD). Experimental results indicate that using EBB-HCD to find the CPs improves the overall alignment accuracy during registration compared with HCD or BB-HCD.

Sun, Yihang↗

DOE OSTI2017

ch-placement

libch-placement is a modular consistent hashing library. It can be used by a distributed storage system to access multiple hashing algorithms, distance metrics, and virtual node settings using a consistent API.

Carns, Phil↗

DOE OSTI2021

Elga v.1

SAND2021-4635 O Elga is a software system that provides distributed and elastic computation for temporal graph algorithms. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Gabert, Kasimir↗

DOE OSTI2026

Dynamic Transmission Line Switching Amid Wildfire-Prone Weather Under Decision-Dependent Uncertainty

During dry and windy seasons, environmental conditions significantly increase the risk of wildfires, exposing power grids to disruptions caused by transmission line failures. Wildfire propagation exacerbates grid vulnerability, potentially leading to prolonged power outages. To address this challenge, we propose a multistage optimization model that dynamically adjusts transmission grid topology in response to wildfire propagation, aiming to develop an optimal response policy. By accounting for decision-dependent uncertainty, where line survival probabilities depend on usage, we employ distributionally robust optimization to model uncertainty in line survival distributions. We adapt the stochastic nested decomposition algorithm and derive a deterministic upper bound for its finite convergence. To enhance computational efficiency, we exploit the Lagrangian dual problem structure for a faster generation of Lagrangian cuts. Using realistic data from the California transmission grid, we demonstrate the superior performance of dynamic response policies against two-stage alternatives through a comprehensive case study. In addition, after solving the multistage formulation, we construct easy-to-implement policies that significantly reduce computational burden while maintaining good performance in real-time deployment. History: Accepted by Russell Bent, Area Editor for Network Optimization: Algorithms and Applications. Funding: This work was supported by the U.S. Department of Energy, Office of Electricity [Grant DE-AC02-05CH11231]. The work of R. Jiang was supported in part by the U.S. National Science Foundation, Division of Electrical, Communications and Cyber Systems [Grant ECCS-1845980] and the U.S. Air Force Office of Scientific Research [Grant FA9550-23-1-0323]. Supplemental Material: The software that supports the findings of this study is available within the paper and its Supplemental Information ( https://pubsonline.informs.org/doi/suppl/10.1287/ijoc.2025.1210 ) as well as from the IJOC GitHub software repository ( https://github.com/INFORMSJoC/2025.1210 ). The complete IJOC Software and Data Repository is available at https://informsjoc.github.io/ .

Estrada-Garcia, Juan-Alberto↗

DOE OSTI2014

On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability

The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the electromagnetic PIC algorithm to analyze the origin of these instabilities. We rigorously derive the faithful 3D numerical dispersion of the PIC algorithm, and then specialize to the Yee FDTD scheme. In particular, we account for the manner in which the PIC algorithm updates and samples the fields and distribution function. Temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme are also explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical 1D modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Advanced Offshore Wind Turbine/Foundation Concept for the Great Lakes

This project investigated a conceptual 2-bladed rotor wind turbine design and assessed its feasibility for installation in the Great Lakes. The levelized cost of energy was used for this purpose. A location in Lake Erie near the coast of Cleveland, Ohio was selected as the application site. The loading environment was defined using wind and wave data collected at a weather station in Lake Erie near Cleveland. In addition, the probability distributions of the annual significant wave height and wind speed were determined. A model of the dependence of the above two quantities was also developed and used in the study of wind turbine system loads. Loads from ice floes and ridges were also included.The NREL 5 MW 3-bladed rotor wind turbine concept was used as the baseline design. The proposed turbine design employs variable pitch blade control with tip-brakes and a teeter mechanism. The rotor diameter, rated power and the tower dimensions were selected to closely match those of the NREL 5 MW wind turbine.A semi-floating gravity base foundation was designed for this project primarily to adapt to regional logistical constraints to transport and install the gravity base foundation. This foundation consists of, from bottom to top, a base plate, a buoyancy chamber, a taper zone, a column (with ice cone), and a service platform. A compound upward-downward ice cone was selected to secure the foundation from moving because of ice impact.The turbine loads analysis was based on International ElectroTechnical Committee (IEC) Standard 61400-1, Class III winds. The NREL software FAST was the primary computational tool used in this study to determine all design load cases. An initial set of studies of the dynamics of wind turbines using Automatic Dynamic Analysis of Mechanical Systems (ADAMS) demonstrated that FAST and ADAMS load predictions were comparable. Because of its relative simplicity and short run times, FAST was selected for this study. For ice load calculations, a method was developed and implemented in FAST to extend its capability for ice load modeling.Both upwind and downwind 2-bladed rotor wind turbine designs were developed and studied. The new rotor blade uses a new twist angle distribution design and a new pitch control algorithm compared with the baseline model. The coning and tilt angles were selected for both the upwind and downwind configurations to maximize the annual energy production. The risk of blade-tower impact is greater for the downwind design, particularly under a power grid fault; however, this risk was effectively reduced by adjusting the tilt angle for the downwind configuration.

17 WIND ENERGY↗

DOE OSTI2017

Partnership for Edge Physics (EPSI), University of Texas Final Report

Simulations of tokamak plasmas require a number of inputs whose values are uncertain. The effects of these input uncertainties on the reliability of model predictions is of great importance when validating predictions by comparison to experimental observations, and when using the predictions for design and operation of devices. However, high fidelity simulation of tokamak plasmas, particular those aimed at characterization of the edge plasma physics, are computationally expensive, so lower cost surrogates are required to enable practical uncertainty estimates. Two surrogate modeling techniques have been explored in the context of tokamak plasma simulations using the XGC family of plasma simulation codes. The first is a response surface surrogate, and the second is an augmented surrogate relying on scenario extrapolation. In addition, to reduce the costs of the XGC simulations, a particle resampling algorithm was developed, which allows marker particle distributions to be adjusted to maintain optimal importance sampling. This means that the total number of particles in and therefore the cost of a simulation can be reduced while maintaining the same accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Fracture Diagnostics Using Low Frequency Electromagnetic Induction and Electrically Conductive Proppants

The monitoring and diagnostics of induced fractures are important for the real-time performance evaluation of hydraulic fracturing operations. Previous electromagnetic (EM) based studies show that single backbone tri-axial induction logging tools are promising candidates for real-time monitoring and diagnosis of fractures in open-hole wells. To support the development of field deployable tools, the concept must be tested in experiments, in a controllable environment, before it is tested under field-like conditions. To this end, we have developed numerical tools which can simulate any wellbore environment while logging hydraulic fractures with the induction tool. We have designed and built a prototype induction tool and performed two sets of tests to compare with numerical simulation results. The computational and experimental setup consists of tri-axial transmitter- and receiver-coils in co-axial, co-planar and cross-polarized configurations. Both lab and shallow earth measurements are shown to be in good agreement with simulations for all examined cases. The average relative and maximum discrepancies of the measured signals from the simulated ones were lower than 3% and 10%, respectively. With the prototype tool, strong signals sensitive to the fracture’s surface area and dip-angle were measured in co-axial coil configuration, while weaker signals sensitive to the fracture’s aspect ratio were observed in the co-planar configuration. Cross-polarized signals are also shown to be strong and sensitive to the fracture’s dip. Lastly, we resolved the detectable components of the signal tensor to obtain parameters of simplified fracture geometries. The inversion algorithm, a derivative free directional search model, uses an objective function defined as a combination of co-axial and cross-polarized signals from different tool spacings to provide a well behaved global minimum. The robustness of inversion algorithm is tested on synthetic data for single cluster fractures in homogeneous and heterogeneous background electrical conductivity. All the effective model parameters for different cases, fracture conductivity (width), propped area (length) and dip-angle, are shown to be recovered with sufficient accuracy. We also evaluated the effect of multiple neighboring fractures and suggested a multi-cluster inversion algorithm that can accurately recover the proppant distribution in each fracture. Based on these results we provide a downhole tool design, with detailed specifications, that can effectively recover far-field proppant distribution in hydraulic fractures.

02 PETROLEUM↗

DOE OSTI2023

Integrated Design of Ultradurable, Low CO 2 Alternative Binder Systems via Machine Learning

This ARPA-E project developed a machine learning tool to use in formulation design of cementitious binders for concrete having 50% less embodied CO 2 and possessing twice the durability compared to concrete based on ordinary portland cement (OPC) binders. The technical focus was on limestone/calcined clay cement (LC3), the leading replacement for OPC. Here, hierarchical machine learning (HML) was used to model the flowability, set time, strength, and durability of LC3 concrete. This methodology identifies latent variables derived from domain knowledge and empirical models that develop an accurate model for a response surface from small datasets. For the flowability metric, particle packing was a dominant factor, while strength and durability were both strongly determined by the fraction of metakaolin and the water:solids ratio. Under constraints of water:binder ratio, material performance metrics, embodied CO 2 , and cost per tonne of OPC, multi-objective optimization was used to design binders parameterized by the mineral composition replacing OPC, particle size distributions, and water:solids ratio. The trained algorithm was able to predict multiple mixes met these performance criteria, and experimental testing validated the predictions. The model demonstrated here is relevant for North America, where pure kaolin deposits are found broadly. The approach is being taken forward into commercial application by Ansatz AI, a materials informatics company founded by PI Washburn and co-PI Poczos. Through collaborations with the cement and concrete industry, and funding from SBIR programs, a commercial software will be developed in future research.

36 MATERIALS SCIENCE↗

DOE OSTI2017

A framework to combine low- and high-resolution spectroscopy for the atmospheres of transiting exoplanets

Current observations of the atmospheres of close-in exoplanets are predominantly obtained with two techniques: low-resolution spectroscopy with space telescopes and high-resolution spectroscopy from the ground. Although the observables delivered by the two methods are in principle highly complementary, no attempt has ever been made to combine them, perhaps due to the different modeling approaches that are typically used in their interpretation. Here, we present the first combined analysis of previously published dayside spectra of the exoplanet HD 209458 b obtained at low resolution with HST/Wide Field Camera 3 (WFC3) and Spitzer/IRAC and at high resolution with VLT/CRIRES. By utilizing a novel retrieval algorithm capable of computing the joint probability distribution of low- and high-resolution spectra, we obtain tight constraints on the chemical composition of the planet’s atmosphere. In contrast to the WFC3 data, we do not confidently detect H{sub 2}O at high spectral resolution. The retrieved water abundance from the combined analysis deviates by 1.9σ from the expectations for a solar-composition atmosphere in chemical equilibrium. Measured relative molecular abundances of CO and H{sub 2}O strongly favor an oxygen-rich atmosphere (C/O < 1 at 3.5σ) for the planet when compared to equilibrium calculations including O rainout. From the abundances of the seven molecular species included in this study we constrain the planet metallicity to 0.1–1.0× the stellar value (1σ). This study opens the way to coordinated exoplanet surveys between the flagship ground- and space-based facilities, which ultimately will be crucial for characterizing potentially habitable planets.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Characterization of Fresh Nuclear Fuel using Time-Of-Flight Neutrons

Neutrons offer excellent characterization capabilities for nuclear fuels and can provide non-destructive insight on micro structure, chemistry and crystallographic phase. Measurements are possible at reactors or pulsed sources, but the inherent energy resolution implicit at pulsed sources offers unique opportunities. Recent advances in 2D neutron detection enabled new capabilities of material characterization using neutron imaging and tomography for nuclear materials. Moreover, time-of-flight (TOF) neutron measurements can enable isotope specific tomography by deriving isotope densities from transmission spectra exhibiting neutron resonances, which are recorded for each pixel. Tomographic algorithms allow 3D reconstructions of the spatial distribution of isotopes, density and features in typical fuel pellet samples. The ability to select radiographs for different energies allows characterizing samples which are opaque for thermal neutrons, such as enriched fuels or fuels containing plutonium, by selecting epi-thermal neutron energies and by knowledge of the total neutron cross sections of the constituents, isotopic enrichment levels can be determined. Nuclear cross sections are typically not calculated directly and databases such as the Evaluated Nuclear Data File (ENDF) database provide this information with cross sections determined by standardized experimental measurements, such as neutron transmission measurements for individual isotopes. Utilizing the information in the ENDF database, the 'inverse' problem can be solved by assuming cross section parameters to be known parameters and areal densities can be determined. This work presents recent advances at the Los Alamos Neutron Science Center (LANSCE) facility at Los Alamos National Laboratory (LANL), which provides the infrastructure to handle radioactive samples, in the development of fuel characterization with TOF neutrons. Examples include recent work on prototype accident tolerant nuclear fuels (ATF). The potential application of these techniques for nondestructive examination of irradiated fuel is under consideration. Non-destructive access to materials features for specimens embedded in containers, particularly for hazardous samples, such as a radioactive nuclear fuel, is indispensable to characterize and ultimately license new fuel forms. Energy-resolve neutron imaging is a new tool for this research field. The results reported here present a new imaging technique, setting a milestone in the capabilities of measuring non-destructively absolute quantities of physical properties. The enabling physics behind this new technique are nuclear resonances. We have demonstrated how absolute areal densities can be extracted from radiographs utilizing pulsed neutron sources. This technique provides information that is otherwise difficult and costly to obtain. In the case of nuclear fuels, the results from energy-dispersive neutron imaging may guide destructive examination of irradiated fuel. This summery does not describe recent analysis, which is ongoing. Particularly, the implementation of calibration standards has improved the determination of enrichment levels without additional calibration. The analysis of the samples was performed with full tomographic reconstructions for isotope specific data and presents a first approach in absolute isotope specific 3D reconstruction. This work was predominantly funded by and received substantial support from the Fuel Cycle R and D program at Los Alamos National Laboratory. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2015

Statistical behaviour of adaptive multilevel splitting algorithms in simple models

Adaptive multilevel splitting algorithms have been introduced rather recently for estimating tail distributions in a fast and efficient way. In particular, they can be used for computing the so-called reactive trajectories corresponding to direct transitions from one metastable state to another. The algorithm is based on successive selection–mutation steps performed on the system in a controlled way. It has two intrinsic parameters, the number of particles/trajectories and the reaction coordinate used for discriminating good or bad trajectories. We investigate first the convergence in law of the algorithm as a function of the timestep for several simple stochastic models. Second, we consider the average duration of reactive trajectories for which no theoretical predictions exist. The most important aspect of this work concerns some systems with two degrees of freedom. They are studied in detail as a function of the reaction coordinate in the asymptotic regime where the number of trajectories goes to infinity. We show that during phase transitions, the statistics of the algorithm deviate significatively from known theoretical results when using non-optimal reaction coordinates. In this case, the variance of the algorithm is peaking at the transition and the convergence of the algorithm can be much slower than the usual expected central limit behaviour. The duration of trajectories is affected as well. Moreover, reactive trajectories do not correspond to the most probable ones. Such behaviour disappears when using the optimal reaction coordinate called committor as predicted by the theory. We finally investigate a three-state Markov chain which reproduces this phenomenon and show logarithmic convergence of the trajectory durations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Which Combinations of Environmental Conditions and Microphysical Parameter Values Produce a Given Orographic Precipitation Distribution?

This study applies an idealized modeling framework, alongside a Bayesian Markov chain Monte Carlo (MCMC) algorithm, to explore which combinations of upstream environmental conditions and cloud microphysical parameter values can produce a particular precipitation distribution over an idealized two-dimensional, bell-shaped mountain. Simulations focus on orographic precipitation produced when an atmospheric river interacts with topography. MCMC-based analysis reveals that different combinations of parameter values produce a similar precipitation distribution, with the most influential parameters being relative humidity (RH), horizontal wind speed (U), surface potential temperature (θ sfc ), and the snow fall speed coefficient (A s ). RH, U, and A s exhibit interdependence: changes in one or more of these factors can be mitigated by compensating changes in the other(s) to produce similar orographic precipitation rates. The results also indicate that the parameter sensitivities and relationships can vary for spatial subregions and given different environmental conditions. In particular, high θ sfc values are more likely to produce the target precipitation rate and spatial distribution, and thus the ensemble of simulations shows a preference for liquid precipitation at the surface. The results presented here highlight the complexity of orographic precipitation controls, and have implications for flood and water management, observational efforts, and climate change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Detection of Periodic Forced Oscillations in Power Systems

An algorithm for the detection and frequency estimation of forced oscillations in power systems is proposed. The method operates by comparing the periodogram of synchrophasor measurements to a detection threshold. This threshold is established by deriving a general expression for the distribution of the periodogram and is related to the algorithm's probabilities of false alarm and detection. Unlike classic detection algorithms designed for use with white Gaussian noise, the proposed algorithm uses a detection threshold that varies with frequency to account for the colored nature of synchrophasor measurements. Further, a detection scheme based on multiple segments of data is proposed to improve the algorithm's performance as a monitoring tool in the online environment. A design approach that helps to ensure that the best-available probability of detection from any one detection segment is constantly increasing with the duration of the forced oscillation is also developed. Results from application of the detection algorithm to simulated and measured power system data suggest that the algorithm provides the expected detection performance and can be used to detect forced oscillations in practical monitoring of power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2015

Simultaneous beam sampling and aperture shape optimization for SPORT

Purpose: Station parameter optimized radiation therapy (SPORT) was recently proposed to fully utilize the technical capability of emerging digital linear accelerators, in which the station parameters of a delivery system, such as aperture shape and weight, couch position/angle, gantry/collimator angle, can be optimized simultaneously. SPORT promises to deliver remarkable radiation dose distributions in an efficient manner, yet there exists no optimization algorithm for its implementation. The purpose of this work is to develop an algorithm to simultaneously optimize the beam sampling and aperture shapes. Methods: The authors build a mathematical model with the fundamental station point parameters as the decision variables. To solve the resulting large-scale optimization problem, the authors devise an effective algorithm by integrating three advanced optimization techniques: column generation, subgradient method, and pattern search. Column generation adds the most beneficial stations sequentially until the plan quality improvement saturates and provides a good starting point for the subsequent optimization. It also adds the new stations during the algorithm if beneficial. For each update resulted from column generation, the subgradient method improves the selected stations locally by reshaping the apertures and updating the beam angles toward a descent subgradient direction. The algorithm continues to improve the selected stations locally and globally by a pattern search algorithm to explore the part of search space not reachable by the subgradient method. By combining these three techniques together, all plausible combinations of station parameters are searched efficiently to yield the optimal solution. Results: A SPORT optimization framework with seamlessly integration of three complementary algorithms, column generation, subgradient method, and pattern search, was established. The proposed technique was applied to two previously treated clinical cases: a head and neck and a prostate case. It significantly improved the target conformality and at the same time critical structure sparing compared with conventional intensity modulated radiation therapy (IMRT). In the head and neck case, for example, the average PTV coverage D99% for two PTVs, cord and brainstem max doses, and right parotid gland mean dose were improved, respectively, by about 7%, 37%, 12%, and 16%. Conclusions: The proposed method automatically determines the number of the stations required to generate a satisfactory plan and optimizes simultaneously the involved station parameters, leading to improved quality of the resultant treatment plans as compared with the conventional IMRT plans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2017

Stochastic level-set method for shape optimisation

We present a new method for stochastic shape optimisation of engineering structures. The method generalises an existing deterministic scheme, in which the structure is represented and evolved by a level-set method coupled with mathematical programming. The stochastic element of the algorithm is built on the methods of statistical mechanics and is designed so that the system explores a Boltzmann–Gibbs distribution of structures. In non-convex optimisation problems, the deterministic algorithm can get trapped in local optima: the stochastic generalisation enables sampling of multiple local optima, which aids the search for the globally-optimal structure. The method is demonstrated for several simple geometrical problems, and a proof-of-principle calculation is shown for a simple engineering structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Exactness of OPF Relaxation on Three-phase Radial Networks with Delta Connections

Simulations have shown that while semi-definite relaxations of AC optimal power flow (AC-OPF) on three-phase radial networks with only wye connections tend to be exact, the presence of delta connections seem to render them inexact. This paper shows that such inexactness originates from the non-uniqueness of relaxation solutions and numerical errors amplified by the non-uniqueness. This finding motivates two algorithms to recover the exact solution of AC-OPF in unbalanced distribution networks featuring both wye and delta connections. In simulations using IEEE 13, 37 and 123-bus systems, the proposed algorithms provide exact optimal solutions up to numerical precision.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2019

Online Optimization as a Feedback Controller: Stability and Tracking

This paper develops and analyzes feedback-based online optimization methods to regulate the output of a linear time-invariant (LTI) dynamical system to the optimal solution of a time-varying convex optimization problem. The design of the algorithm is based on continuous-time primal-dual dynamics, properly modified to incorporate feedback from the LTI dynamical system, applied to a proximal augmented Lagrangian function. The resultant closed-loop algorithm tracks the solution of the time-varying optimization problem without requiring knowledge of (time-varying) disturbances in the dynamical system. The analysis leverages integral quadratic constraints to provide linear matrix inequality (LMI) conditions that guarantee global exponential stability and bounded tracking error. Analytical results show that, under a sufficient time-scale separation between the dynamics of the LTI dynamical system and the algorithm, the LMI conditions can be always satisfied. The paper further proposes a modified algorithm that can track an approximate solution trajectory of the constrained optimization problem under less restrictive assumptions. As an illustrative example, the proposed algorithms are showcased for power transmission systems, to compress the time scales between secondary and tertiary control, and allow to simultaneously power re-balancing and tracking of DC optimal power flow points.

24 POWER TRANSMISSION AND DISTRIBUTION↗