Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed method”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

2-D seismic wave propagation using the distributional finite-difference method: further developments and potential for global seismology

SUMMARY We present a time-domain distributional finite-difference scheme based on the Lebedev staggered grid for the numerical simulation of wave propagation in acoustic and elastic media. The central aspect of the proposed method is the representation of the stresses and displacements with different sets of B-splines functions organized according to the staggered grid. The distributional finite-difference approach allows domain-decomposition, heterogeneity of the medium, curvilinear mesh, anisotropy, non-conformal interfaces, discontinuous grid and fluid–solid interfaces. Numerical examples show that the proposed scheme is suitable to model wave propagation through the Earth, where sharp interfaces separate large, relatively homogeneous layers. A few domains or elements are sufficient to represent the Earth’s internal structure without relying on advanced meshing techniques. We compare seismograms obtained with the proposed scheme and the spectral element method, and we show that our approach offers superior accuracy, reduced memory usage, and comparable efficiency.

Geochemistry & Geophysics↗
DOE OSTI2022

Progress in calculation of the fourth Mellin moment of the pion light-cone distribution amplitude using the HOPE method

The pion light-cone distribution amplitude (LCDA) is a central non-perturbative object of interest for the calculation of high-energy exclusive processes in quantum chromodynamics. This article describes the progress in the lattice QCD calculation of the fourth Mellin moment of the pion LCDA using a heavy-quark operator product expansion (HOPE).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Blockchain-Enabled Cyber-Secure Microgrid Control Using Consensus Algorithm

This paper proposes a distributed control method based on a consensus algorithm for distributed energy resources (DERs) using blockchain as a secure communication medium for cyber resilience. Each DER communicates with a local blockchain server that is maintained by the distributed ledger technology to allow for securely sharing local measurements among neighbouring assets to achieve the global control objectives, i.e., voltage and frequency regulation as well as accurate power sharing among the DERs, including collective grid-forming capability. To prove that the distributed control can retain system stability under a blockchain-induced variable communication delay, Lyapunov function-based stability analysis is carried out. This paper demonstrates the concept on an 11-bus test case developed in MATLAB Simulink, which has been modified from the IEEE 9-bus test case, to study dynamic operations of the five inverter-based DERs working with a blockchain-induced variable delay. The results validate the superior performance of the proposed control method, compared to heavily compromised operations of the other test cases without the security measure, affected by heavy communication delays and communication interruptions.

blockchains↗
DOE OSTI2025

End-to-end microgrid protection using distributed data-driven methods

This paper introduces an end-to-end microgrid protection framework that offers real-time system monitoring, fault-related decision making, and circuit breaker control. This is achieved through the design of distributed data-driven techniques based on the support vector machine method, where each relay is responsible for distributed data collection, fault detection, fault localization, and fault isolation. Local communication is established among neighboring relays, fostering cooperative fault localization and isolation. This decentralized design not only reduces the computational and communication requirements but also enables the adaptability of each relay under varying operational dynamics. The proposed end-to-end protection framework was validated using MATLAB/Simulink simulations on a 100% renewable microgrid, achieving an accuracy of 93.1% with response time of 0.0523 s, in protecting against a range of fault scenarios that are characterized by various types, locations, impedances, load conditions, photovoltaic power levels, and microgrid operating modes.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Advancing electrochemical impedance analysis through innovations in the distribution of relaxation times method

Electrochemical impedance spectroscopy (EIS) is a key tool across various scientific disciplines, including energy sciences, chemistry, and biology, enabling the analysis of electrochemical systems. However, conventional methods for interpreting EIS data are often complex and model dependent. The distribution of relaxation times (DRT) offers a non-parametric approach that simplifies the interpretation process by providing a timescale interpretation of EIS data. This article provides a comprehensive review of current methods for DRT inversion. Additionally, a survey of practitioners highlights key challenges in the field. Here, the findings underscore the need for standardized DRT analysis and benchmarks, as well as the development of automated analysis tools. These advancements would improve the usability and interpretability of EIS data. Ultimately, implementing these improvements could not only propel the field forward but also expand the application of DRT in scientific research by making it accessible to a broader range of researchers, including those without specialized expertise in programming or statistics.

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

Using distributed power electronics-based devices to improve the voltage and frequency stability of distribution systems

Methods include, in response to a line frequency variation of a power grid, adjusting a voltage setpoint of a voltage regulator coupled to the power grid at a grid edge to maintain a voltage at the grid edge, wherein the adjusting the regulated voltage setpoint is configured to reduce the line frequency variation to stabilize the line frequency of the power grid. Apparatus include a voltage regulator configured to couple to a power grid at a grid edge and to maintain a voltage at the grid edge, wherein the voltage regulator is further configured to adjust a voltage setpoint of a voltage regulator in response to a line frequency variation of the power grid to reduce the line frequency variation and stabilize the line frequency of the power grid.

Du, Wei↗
DOE OSTI2020

Sequential ensemble transform for Bayesian inverse problems

In this work, we present the Sequential Ensemble Transform (SET) method, an approach for generating approximate samples from a Bayesian posterior distribution. The method explores the posterior distribution by solving a sequence of discrete optimal transport problems to produce a series of transport plans which map prior samples to posterior samples. We prove that the sequence of Dirac mixture distributions produced by the SET method converges weakly to the true posterior as the sample size approaches infinity. Furthermore, our numerical results indicate that, when compared to standard Sequential Monte Carlo (SMC) methods, the SET approach is more robust to the choice of Markov mutation kernels and requires less computational efforts to reach a similar accuracy when used to explore complex posterior distributions. Finally, we describe adaptive schemes that allow to completely automate the use of the SET method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

DLMP of Competitive Markets in Active Distribution Networks: Models, Solutions, Applications, and Visions

Traditionally, the electric distribution system operates with uniform energy prices across all system nodes. However, as the adoption of distributed energy resources (DERs) propels a shift from passive to active distribution network (ADN) operation, a distribution-level electricity market has been proposed to manage new complexities efficiently. In addition, distribution locational marginal price (DLMP) has been established in the literature as the primary pricing mechanism. The DLMP inherits the LMP concept in the transmission-level wholesale market but incorporates characteristics of the distribution system, such as high $R/X$ ratios and power losses, system imbalance, and voltage regulation needs. The DLMP provides a solution that can be essential for competitive market operation in future distribution systems. This article first provides an overview of the current distribution-level market architectures and their early implementations. Next, the general clearing model, model relaxations, and DLMP formulation are comprehensively reviewed. The state-of-the-art solution methods for distribution market clearing are summarized and categorized into centralized, distributed, and decentralized methods. Then, DLMP applications for the operation and planning of DERs and distribution system operators (DSOs) are discussed in detail. Finally, visions of future research directions and possible barriers and challenges are presented.

42 ENGINEERING↗
DOE OSTI2025

Analysis of Niobium Electropolishing Using a Generalized Distribution of Relaxation Times Method

Using electrochemical impedance spectroscopy, we have devised a method of sensing the microscopic surface conditions on the surface of niobium as it is undergoing an electrochemical polishing (EP) treatment. The method uses electrochemical impedance spectroscopy (EIS) to gather information on the surface state of the electrode without disrupting the polishing reaction. The EIS data is analyzed using a so-called distribution of relaxation times (DRT) method. Using DRT, the EIS data can be deconvolved into discrete relaxation time peaks without any a priori knowledge of the electrode dynamics. By analyzing the relaxation time peaks, we are able to distinguish two distinct modes of the EP reaction. As the polishing voltage is increased, the electrode transitions from the low voltage EP mode, characterized by a single relaxation time peaks, to the high voltage EP mode, characterized by two relaxation time peaks. We theorize that this second peak is caused by the formation of an oxide layer on the electrode. We also find that this oxide induced peak transitions from to a negative relaxation time, which is indicative of a blocking electrode process. By analyzing EPed samples, we show that samples polished in the low voltage mode have significantly higher surface roughness due to grain etching and faceting. We find that the surface roughness of the samples only improves when the oxide film peak is present and in the negative relaxation time region. This shows that EIS combined with DRT analysis can be used to predict etching on EPed Nb. This method can also be performed before or during the EP, which could allow for adjustment of polishing parameters to guarantee a smooth cavity surface finish.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Method to determine the electron–ion temperature relaxation rate from test particle distributions

Here a method to calculate the electron–ion energy exchange rate from the forces on and velocities of an ensemble of test particles is described. The essence of the method is that the energy exchange rate is related to the covariance between the distribution of velocities of test particles and the forces acting on them. The energy exchange rate is related to the electron–ion temperature relaxation rate in the limit of a low test particle speed. A proof of concept is conducted with first-principles molecular dynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Fast and Scalable Genetic Algorithm-Based Approach for Planning of Microgrids in Distribution Networks

As a result of climate change, extreme weather events are occurring more frequently and with increasing impact. This trend poses a significant challenge for distribution utilities and system operators to ensure that there is uninterrupted power supply to critical loads in their networks; thus, the level of proactive preparation of the distribution system to be able to handle severe impacts of extreme weather events represents the system's resilience. One method that distribution system planners can use to prepare for future extreme events is to plan multiple microgrids which can use local generation as much as possible to supply critical loads. But partitioning an existing distribution system such that multiple feasible islands are planned and which are capable of supporting critical loads is still challenging for distribution systems - first, because of the size of the network graph partitioning problem and, second, because of the difficulty in properly formulating the desired attributes of such islands or microgrids. Therefore, this paper presents a genetic algorithm based approach that facilitates incorporating multiple objectives for grid partitioning by formulating two types of problems - node allocation and edge elimination - and it considers multiple topological and resilience-enhancing objectives. The performance of the proposed genetic algorithm-based approach is numerically evaluated on multiple test systems as well as on a real distribution feeder in Colorado, United States.

genetic algorithm↗
DOE OSTI2021

AC Power Flow Based DLMP Calculation and Decomposition Method to Smooth Power Fluctuation of Distributed Renewable Energy Sources

As the penetration of renewable energy sources increases, the growing renewable power variability brings ramping issues to power systems. Meanwhile, the development of distributed energy resources (DERs) makes the distribution systems to provide both energy and ancillary services. To incentivise individual resources and customers to alleviate ramping issues on the demand side, a two-stage distribution locational marginal price (DLMP) calculation and decomposition method is developed to formulate the marginal power ramping price for DERs. In the first stage of the proposed method, a distribution system operator market scheduling model based on AC optimal power flow is designed to estimate the optimal operating point of the distribution system. Subsequently, the voltage and power flow constraints are linearised in stage two to calculate DLMP. Finally, based on the Lagrange function and sensitivity factors, DLMP is decomposed to the marginal costs for active/reactive power, voltage management, power loss and power variability. Case studies demonstrate that the proposed model can effectively smooth the power fluctuation and reduce the ramping flexibility requirements of distribution systems.

AC optimal power flow↗
DOE OSTI2023

A Sparse Distributed Gigascale Resolution Material Point Method

In this paper, we present a four-layer distributed simulation system and its adaptation to the Material Point Method (MPM). The system is built upon a performance portable C++ programming model targeting major High-Performance-Computing (HPC) platforms. A key ingredient of our system is a hierarchical block-tile-cell sparse grid data structure that is distributable to an arbitrary number of Message Passing Interface (MPI) ranks. We additionally propose strategies for efficient dynamic load balance optimization to maximize the efficiency of MPI tasks. Our simulation pipeline can easily switch among backend programming models, including OpenMP and CUDA, and can be effortlessly dispatched onto supercomputers and the cloud. Finally, we construct benchmark experiments and ablation studies on supercomputers and consumer workstations in a local network to evaluate the scalability and load balancing criteria. We demonstrate massively parallel, highly scalable, and gigascale resolution MPM simulations of up to 1.01 billion particles for less than 323.25 seconds per frame with 8 OpenSSH-connected workstations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Embedded, Real-Time, and Distributed Traveling Wave Fault Location Method Using Graph Convolutional Neural Networks

This work proposes and develops an implementation of a fault location method to provide a fast and resilient protection scheme for power distribution systems. The method analyzes the transient dynamics of traveling waves (TWs) to generate features using the discrete wavelet transform (DWT), which are then used to train several graph convolutional network (GCN) models. Faults are simulated in the IEEE 34-node system, which is divided into three protection zones (PZs). The goal is to identify the PZ in which the fault occurs. The GCN models create a distributed protection scheme, as all nodes are able to retrieve a prediction. Given that message-passing between nodes occurs both during training and in the execution of the model, the resiliency of such schemes to communication losses was analyzed and demonstrated. One of the models, which only uses voltage measurements, was implemented on a Texas Instruments F28379D development board. The execution times were monitored to assess the speed of the protection scheme. It is shown that the proposed method can be executed in approximately a millisecond, which is comparable to existing TW protection in the transmission system. For experimental purposes, a DWT-based detection method is employed. A design of a setup to playback TWs using two development boards is also addressed.

Jiménez-Aparicio, Miguel (ORCID:000000016864461X)↗
DOE OSTI2021

Multi-Area Model-Free State Estimation via Distributed Tensor Decomposition

This paper proposes a model-free method for distribution system state estimation based on tensor completion using canonical polyadic decomposition. In particular, we consider a setting where the network is divided into multiple areas. The measured physical quantities at buses located in the same area are processed by an area controller. A third-order tensor is constructed to collect these measured quantities. The measurements are analyzed locally to recover the full state information of the network. A closed-form iterative algorithm based on the alternating direction method of multipliers is developed to obtain the low-rank factors of the whole network state tensor where information exchange happens only between neighboring areas. To demonstrate the efficacy of the developed algorithm, numerical simulations are carried out using an IEEE test system.

alternating direction method of multipliers↗
DOE OSTI2025

Local structure of zinc–indium–tin oxide films via grazing-incidence x-ray pair-distribution functions and theoretical methods

A detailed experimental and theoretical study on the local (r ≤ 4.5 Å) atomic structure of amorphous and crystalline zinc–indium–tin oxide (ZITO) thin films using grazing-incidence x-ray Pair-Distribution Functions (PDFs), ab initio Molecular Dynamics (MD), and Empirical Potential Structure Refinement (EPSR) Monte Carlo simulations is presented. High-energy synchrotron x rays, a two-dimensional detector, and different incident angles were used to probe the depth uniformity of five (ZnO) 0.15 (In 2 O 3 ) 0.70 (SnO 2 ) 0.15 films that were deposited via pulsed-laser deposition at growth temperatures (T G ) ranging from 25 to 300 °C. Films deposited at T G ≤ 150 °C were amorphous. The partially crystalline (T G = 200 °C) and fully crystalline (T G = 300 °C) films were highly textured. Both crystalline and amorphous structures were investigated using ab initio MD and EPSR Monte Carlo simulations. The density of the amorphous films determined from the experimental data agreed with MD calculations. Coordination numbers, bond lengths, and distortion for metal–oxygen and for both the edge- and corner-shared In–metal shells up to 4.5 Å obtained from PDF analysis closely agreed with MD and EPSR simulations. There is a pronounced decrease in the edge- and corner-shared In–Zn distances arising from the shorter Zn–O bond length, Zn–O tetrahedral coordination, and In–O–Zn angle in amorphous ZITO compared to its crystalline counterpart. A maximum in electrical mobility was observed for the amorphous film just before crystallization occurred. While the peak is broad, consistent with nearly unchanged overall cation–oxygen coordination in the amorphous films, ESPR results indicate that the tetrahedral coordination follows the conductivity trend.

Grazing Incidence X-ray↗
DOE OSTI2021

Review of Intrusion Detection Methods and Tools for Distributed Energy Resources

Recent trends in the growth of distributed energy resources (DER) in the electric grid and newfound malware frameworks that target internet of things (IoT) devices is driving an urgent need for more reliable and effective methods for intrusion detection and prevention. Cybersecurity intrusion detection systems (IDSs) are responsible for detecting threats by monitoring and analyzing network data, which can originate either from networking equipment or end-devices. Creating intrusion detection systems for PV/DER networks is a challenging undertaking because of the diversity of the attack types and intermittency and variability in the data. Distinguishing malicious events from other sources of anomalies or system faults is particularly difficult. New approaches are needed that not only sense anomalies in the power system but also determine causational factors for the detected events. In this report, a range of IDS approaches were summarized along with their pros and cons. Using the review of IDS approaches and subsequent gap analysis for application to DER systems, a preliminary hybrid IDS approach to protect PV/DER communications is formed in the conclusion of this report to inform ongoing and future research regarding the cybersecurity and resilience enhancement of DER systems.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Absolute mass calibration of fission product distributions measured with the E-$v$ method

The mass calibration of fission product distributions measured with the energy–velocity (E-v) method constitutes a technical challenge. The energy loss of the fission fragments in the various dead layers of the spectrometer and other sources of pulse-height defects of the energy detectors are a significant source of systematic uncertainty in the mass calculations used to determine the measured fission product yields. Here, in this study, the absolute calibration of the fission mass distributions is accomplished by measuring prompt γ-rays in coincidence with the fission fragments. This allows for the direct calibration of the mass spectra and eliminates any dependence on potentially complicated energy corrections. The first test measurement was performed at the Los Alamos Neutron Science Center employing the SPectrometer for Ion DEtermination in fission Research (SPIDER). SPIDER is a 2E-2v spectrometer designed for measuring independent fission product yields from neutron-induced fission. In this test, the single-arm SPIDER (E-v) system and an array of 252 Cf sources were used. The single-arm system consisted of two time pick-off detectors for measuring the time-of-flight of the fission fragments and a double-sided silicon strip detector (DSSD) for measuring the kinetic energy. Characteristic γ-rays from fission fragments were detected using three high-purity germanium (HPGe) detectors. For the mass calibration, γ–mass coincidence events from twelve product isotopes were used. The measured FPYs from 252 Cf spontaneous fission were found to be in excellent agreement with the evaluated data after applying the absolute mass calibration. From the γ–mass coincidence events, the mass resolution of the system was also extracted. An average mass resolution of ~1.4 AMU (FWHM) for the light fragments and ~2.4 AMU (FWHM) for the heavy fragments were found. This was the first in-situ calibration and characterization of the SPIDER spectrometer, which paves the way for high-quality FPY measurements with this instrument.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗