Two-Level Decentralized-Centralized Control of Distributed Energy Resources in Grid-Interactive Efficient Buildings
Explore the source record for details and available documents.
Huo, Xiang↗
SEARCH · Engineering Papers
Results for “Grid 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.
RLC4CLR (Reinforcement Learning Controller for Critical Load Restoration Problems)
RLC4CLR demonstrates using a reinforcement learning controller (RLC) to solve a critical load restoration (CLR) problem, which improves the grid resilience after a substation outage event. RLC4CLR consists of two parts. (1) RL environment: This environment encapsulates the CLR problem to be solved and provides interfacing functions to follow the standard OpenAI Gym format. A power system simulator, i.e., OpenDSS, is included to provide the power flow solution. Controller inputs and outputs (RL state and action) as well as the reward are defined in this environment as well. In summary, the RL environment is the problem formulation from which the RL agent can learn. (2) RL training script: The training script enables the RL agent to learn its control policy by interacting with the RL environment. For RL training, an open-sourced RL library, i.e., RLlib, is leveraged which is based on a distributed computing framework (Ray). The training script is designed to be able to be run on both local machine or the NREL HPC system. Other components of RLC4CLR include input data, e.g., grid model (standard IEEE test feeders), and other files used for results analysis.
Zhang, Xiangyu↗
Identification of topology changes in power grids using phasor measurements
Phasor measurement units (PMUs) are increasingly important for monitoring the state of an electrical power grid and quickly detecting topology changes caused by events such as lines going down or large loads being dropped. Phasors are complex‐valued measurements of voltage and current at various points of generation and consumption. If a line goes down or a load is removed, power flows change throughout the grid according to known physical laws, and the probability distribution of phasor measurements changes accordingly. This paper develops a method to estimate the current topology of a power grid from phasor measurements and considers the design goal of placing PMUs at strategic points in a distribution system to achieve good sensitivity to single‐line outages. From a vector of phasor measurements, probabilities are computed corresponding to the scenario that all power lines are operational and to alternate scenarios in which each line goes down individually. These probabilities are functions of the joint distributions of phasor measurements under each possible scenario, obtained through Monte Carlo simulations with random load profiles. We use log‐spline densities to estimate marginal distributions of phasor measurements and fold these into a multivariate Gaussian copula to capture important correlations. Sensitivity to outages varies according to which line goes down and where PMUs are placed on the grid. A greedy search algorithm is demonstrated for placing PMUs at locations that provide good sensitivity to single‐line outages. Published 2014. This article is a U.S. Government work and is in the public domain in the USA.
Wiel, Scott Vander↗
SU-E-I-52: Effect of Various X-Ray Beam Qualities On the Exposure Index
Purpose: The exposure index (EI) proposed by the International Electrotechnical Commission (IEC) 62494-1 is expected to be utilized as a standard dose index by every manufacturer. The IEC recommended the usage of RQA5 for the EI. However, X-ray beam qualities, particularly in clinical practices, vary depending on the examination objects and exposure conditions, including usage of anti-scatter grids. We investigated the effects of the X-ray beam qualities other than RQA5 on the EI. Methods: The Xray beam qualities of RQA3, 5, 7, and 9 in IEC 61267 Ed. 1.0 were adopted in a computed radiography system. A uniform exposure without objects was performed to measure the exposure indicators (S values) and air kerma (K). The relational equations between the S values and K were derived for the determination of the EI values. The EI values for RQA3, 7, and 9 were compared to those for RQA5 at the fixed S values of 100, 200, 400, and 600. Finally, the half-value layers (HVLs) using four grids (ratio 6:1, 8:1, 10:1, and 12:1) for the RQA5 X-ray were compared to those with RQA3–9. Results: The EI values for RQA3, 7, and 9 were up to 35.3%, 11.8%, and 38.7% higher, respectively, than that for RQA5 at the S value of 600. The HVLs without grids and with various grids for RQA5 were 6.85 mm Al. and in the range of 6.94–7.29 mm Al. (ΔHVL: up to 0.44 mm Al.), respectively. This variation in the HVLs with grids was smaller than that observed for RQA3–9 (ΔHVL: 2.0–7.5 mm Al.). Conclusion: Although the usage of grids may not greatly affect the EI, the X-ray beam quality for the determination of the EI cannot be ignored in the clinical evaluation of the dose index.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Regional Earthquake Ground Motion Simulations for Southern California With EQSIM: Insights From the 2008 Chino Hills, 2024 Highland Park, and 2021 Carson Earthquakes
This study presents physics-based, 3D simulations using the EQSIM framework for several earthquakes in the Los Angeles region. The primary objective was to assess the ability of deterministic physics-based ground motion simulations to reproduce the observed motions from historical events. The selected events included the mathematical equation M w 5.4 2008 Chino Hills, the mathematical equation M w 4.4 2024 Highland Park, and the mathematical equation M w 4.3 2021 Carson events. The simulated motions were evaluated by comparing the recorded and simulated seismograms, as well as the Fourier amplitude spectra, across multiple seismic stations. The SCEC 3D velocity model, CVM-S4.26.M01, was used to represent the regional geology, and ground motion simulations were carried out with a resolution of up to 5 Hz. The results indicate that the simulated motions captured the recorded motions up to approximately 4 Hz. While careful iterations regarding source parameters and corner frequencies were required, and, for the case of the Highland Park event, some of the near-source stations had relatively low accuracy, the present study established a positive step toward the utilization of physics-based simulations in practical applications. The computational efficiencies exhibited by EQSIM, especially on GPU clusters, further supported this assertion, as wall-clock times of simulations involving more than 10 billion grid points were as low as mathematical equation minutes. This permits ensemble simulations for a considered scenario event so that modeling uncertainties (e.g., source and geology) can be bracketed.
EQSIM↗
Towards Adaptive Mesh Refinement for the Spectral Element Solver Nek5000
When performing computational fluid dynamics (CFD) simulations of complex flows, the a priori knowledge of the flow physics and the location of the dominant flow features are usually unknown. For this reason, the development of adaptive remeshing techniques is crucial for large-scale computational problems. Some work has been made recently to provide Nek5000 with adaptive mesh refinement (AMR) capabilities in order to facilitate the generation of the grid and push forward the limit in terms of problem size and complexity [10]. Nek5000 is an open-source, highly scalable and portable code based on the spectral element method (SEM) [4], which offers minimal dissipation and dispersion, high accuracy and exponential convergence. It is aimed at solving direct and large-eddy simulations of turbulent incompressible or low Mach-number flows with heat transfer and species transport. The approach chosen for adapting the mesh is the h-refinement method, where elements are split locally, which requires the relaxation of the conforming grid constraint currently imposed by Nek5000. Other challenges include the implementation of an efficient management of the grid as refinement is applied, the development of tools to localize the critical flow regions via error estimators and the extension of the current preconditioning strategy to non-conforming grids. In this paper, we present a new procedure to setup an algebraic multigrid solver used as part of the preconditioner for the pressure equation.
Offermans, N.↗
Enhanced Oblique Decision Tree Enabled Policy Extraction for Deep Reinforcement Learning in Power System Emergency Control
Deep reinforcement learning (DRL) algorithms have successfully solved many challenging problems in various power system control scenarios. However, their decision-making process is usually regarded as black-boxes. Furthermore, how DRL models interact with human intelligence remains an open problem. Thus, this paper proposes a policy extraction framework to extract a complex DRL model into an explainable policy. This framework includes three parts: 1) DRL training and data generation. We train an agent for a specific control task and generate data, which contains the control policy of the agent. 2) Policy extraction. We propose an information gain rate based weighted oblique decision tree (IGR-WODT) for DRL policy extraction. 3) Policy evaluation. We define three metrics to evaluate the performance of the proposed approach. A case study for the under-voltage load shedding problem shows that the IGR-WODT presents a performance enhancement compared with DRL, weighted oblique decision tree, and univariate decision tree. The proposed policy extraction method could provide an intuitive explanation of the neural network decision-making process to the dispatchers when making final decisions on power grid operation. Also, the resulted rule-based controller could replace the deep neural network-based controller in many field edge devices with limited computing resources, providing comparable performance.
deep reinforcement learning↗
Numerical stability of relativistic beam multidimensional PIC simulations employing the Esirkepov algorithm
Rapidly growing numerical instabilities routinely occur in multidimensional particle-in-cell computer simulations of plasma-based particle accelerators, astrophysical phenomena, and relativistic charged particle beams. Reducing instability growth to acceptable levels has necessitated higher resolution grids, high-order field solvers, current filtering, etc. except for certain ratios of the time step to the axial cell size, for which numerical growth rates and saturation levels are reduced substantially. This paper derives and solves the cold beam dispersion relation for numerical instabilities in multidimensional, relativistic, electromagnetic particle-in-cell programs employing either the standard or the Cole–Karkkainnen finite difference field solver on a staggered mesh and the common Esirkepov current-gathering algorithm. Good overall agreement is achieved with previously reported results of the WARP code. In particular, the existence of select time steps for which instabilities are minimized is explained. Additionally, an alternative field interpolation algorithm is proposed for which instabilities are almost completely eliminated for a particular time step in ultra-relativistic simulations.
43 PARTICLE ACCELERATORS↗
Amplitude flux, probability flux, and gauge invariance in the finite volume scheme for the Schrödinger equation
The time-dependent Schrödinger equation can be put in a probability conserving, gauge invariant form, on arbitrary structured grids via finite volume discretization. The gauge terms in the discrete system cancel with a portion of the amplitude flux to produce abbreviated flux functions. The resulting time translation operator is strictly unitary, and is compatible with an efficient operator splitting scheme that allows for multi-dimensional simulation with complex grid geometries. Moreover, the abbreviated amplitude flux is necessary to the construction of a conservative probability current. This construction turns out to be important when computing Bohmian trajectories in multi-dimensions. Bohmian trajectories are useful in the interpretation of quantum mechanical phenomena such as tunneling ionization, and provide a bridge between quantum and classical regimes.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A WENO-solver combined with adaptive momentum discretization for the Wigner transport equation and its application to resonant tunneling diodes
We present a novel numerical scheme for the deterministic solution of the Wigner transport equation, especially suited to deal with situations in which strong quantum effects are present. The unique feature of the algorithm is the expansion of the Wigner function in local basis functions, similar to finite element or finite volume methods. This procedure yields a discretization of the pseudo-differential operator that conserves the particle density on arbitrarily chosen grids. The high flexibility in refining the grid spacing together with the weighted essentially non-oscillatory (WENO) scheme for the advection term allows for an accurate and well-resolved simulation of the phase space dynamics. A resonant tunneling diode is considered as test case and a detailed convergence study is given by comparing the results to a non-equilibrium Green's functions calculation. The impact of the considered domain size and of the grid spacing is analyzed. The obtained convergence of the results towards a quasi-exact agreement of the steady state Wigner and Green's functions computations demonstrates the accuracy of the scheme, as well as the high flexibility to adjust to different physical situations.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Shock interaction with deformable particles using a constrained interface reinitialization scheme
In this paper, we present axisymmetric numerical simulations of shock propagation in nitromethane over an aluminum particle for post-shock pressures up to 10 GPa. We use the Mie-Gruneisen equation of state to describe both the medium and the particle. The numerical method is a finite-volume based solver on a Cartesian grid, that allows for multi-material interfaces and shocks, and uses a novel constrained reinitialization scheme to precisely preserve particle mass and volume. We compute the unsteady inviscid drag coefficient as a function of time, and show that when normalized by post-shock conditions, the maximum drag coefficient decreases with increasing post-shock pressure. We also compute the mass-averaged particle pressure and show that the observed oscillations inside the particle are on the particle-acoustic time scale. Finally, we present simplified point-particle models that can be used for macroscale simulations. In the Appendix, we extend the isothermal or isentropic assumption concerning the point-force models to non-ideal equations of state, thus justifying their use for the current problem.
Sridharan, P. (ORCID:0000000287902423)↗
Modeling and Simulation of Power-Electronic Inverters in Analog Electronic Circuit Simulators
This paper demonstrates how equivalent-circuit representations of grid-following power-electronic inverters can be realized within a SPICE-based development environment using common circuit components and VerilogA code. This facilitates computationally lean simulations of inverter networks leveraging the strengths of SPICE in large-scale simulations. We validate the approach with time-domain simulations for a modified version of the IEEE standard 118-bus system modeled in Virtuoso (a SPICE-based solver). Simulation results are compared-focusing on accuracy and computation speed-with results from a commonly used power-electronics simulation package. We note a significant decrease in simulation time with comparable signal resolution when simulating the network in SPICE.
Billmeyer, Ryan↗
Model-Free Control of Grid-Interactive Efficient Buildings Under Communication Time Delays
Grid-interactive efficient buildings (GEBs) have recently been used to enhance the reliability and stability of the electric grid through demand response (DR) programs. However, most existing DR control strategies require accurate modeling of the various building thermostatically controlled loads (TCLs) and are computationally expensive. To address these challenges, a model-free control (MFC)-based strategy has recently been introduced for coordinating and controlling GEBs. MFC is a data-enabled control strategy that is computationally efficient and does not require the analytical models of the various building equipment. In this paper, we numerically investigate the impact of communication time delays on the performance of MFC in maintaining the TCLs' temperatures within the desired comfort levels while meeting the assigned power allocation constraint.
Telsang, Bhagyashri [University of Tennessee, Knox↗
SU-E-T-423: Fast Photon Convolution Calculation with a 3D-Ideal Kernel On the GPU
Purpose: The calculation time is a trade-off for improving the accuracy of convolution dose calculation with fine calculation spacing of the KERMA kernel. We investigated to accelerate the convolution calculation using an ideal kernel on the Graphic Processing Units (GPU). Methods: The calculation was performed on the AMD graphics hardware of Dual FirePro D700 and our algorithm was implemented using the Aparapi that convert Java bytecode to OpenCL. The process of dose calculation was separated with the TERMA and KERMA steps. The dose deposited at the coordinate (x, y, z) was determined in the process. In the dose calculation running on the central processing unit (CPU) of Intel Xeon E5, the calculation loops were performed for all calculation points. On the GPU computation, all of the calculation processes for the points were sent to the GPU and the multi-thread computation was done. In this study, the dose calculation was performed in a water equivalent homogeneous phantom with 150{sup 3} voxels (2 mm calculation grid) and the calculation speed on the GPU to that on the CPU and the accuracy of PDD were compared. Results: The calculation time for the GPU and the CPU were 3.3 sec and 4.4 hour, respectively. The calculation speed for the GPU was 4800 times faster than that for the CPU. The PDD curve for the GPU was perfectly matched to that for the CPU. Conclusion: The convolution calculation with the ideal kernel on the GPU was clinically acceptable for time and may be more accurate in an inhomogeneous region. Intensity modulated arc therapy needs dose calculations for different gantry angles at many control points. Thus, it would be more practical that the kernel uses a coarse spacing technique if the calculation is faster while keeping the similar accuracy to a current treatment planning system.
60 APPLIED LIFE SCIENCES↗
A review of privacy in energy applications
As the distribution system continues to experience an increase in distributed energy resource (DER) and electric vehicle (EV) penetration, so does the need for new solutions that can help grid operators manage and leverage their capabilities. This will undoubtedly lead to new operational schemes and business opportunities that will transform the traditional consumer into a prosumer who will be more actively engaged in grid operations. Although the field is still under active development, many of the potential use cases presented in literature or industry are built upon edge computing, two-way communications, and other innovative computational constructs to attain their goals. However, at their core, many use cases assume a great level of data access to aid with the decision-making process, an assumption that may need to be revised to ensure fair and equitable operational processes are maintained. This may be particularly true as edge resources are predicted to participate in retail-side, many-to-many, or peer-to-peer markets and thus may lead to financial impacts if data access considerations are ignored. The need to revise data access mechanisms can be further justified by the introduction of new participants into the operational process, who do not have the same level of trust, nor the incentives to focus on energy delivery as their primary objective. At the same time, more consumers are becoming aware of their own data, and the potential impacts of its abuse. To help solution developers better understand these risks, this report has been developed to offer an initial introduction to the topic of privacy. This is achieved by 1) Highlighting the need for privacy-aware solutions; 2) Encouraging system designers to be inquisitive about the status quo; 3) Documenting the existing threat space; 4) Presenting and evaluating tools that may be helpful towards enabling better privacy postures; and 5) Making recommendations to encourage the adoption of better practices. From a technical perspective, the report focuses on evaluating two potential techniques by applying them to the Transactive Energy Space. Based on the obtained results, it can be established that differential privacy (DP) methods may have limited applicability when highly correlated, time-series data records need to be protected. However, DP may be a powerful tool when it is used to aggregate and analyze mid-size and large-size data sets in a more traditional statistical environment. The second tool under evaluation is threshold cryptography, which can guarantee complete secrecy (and thus privacy) but requires the establishment of key management procedures and dedicated communication channels for key coordination. Therefore, due to its increased computational overhead, the use of threshold cryptography must be weighted using a cost/benefit analysis on a per-application basis.
24 POWER TRANSMISSION AND DISTRIBUTION↗
SymPix: A spherical grid for efficient sampling of rotationally invariant operators
We present SymPix, a special-purpose spherical grid optimized for efficiently sampling rotationally invariant linear operators. This grid is conceptually similar to the Gauss–Legendre (GL) grid, aligning sample points with iso-latitude rings located on Legendre polynomial zeros. Unlike the GL grid, however, the number of grid points per ring varies as a function of latitude, avoiding expensive oversampling near the poles and ensuring nearly equal sky area per grid point. The ratio between the number of grid points in two neighboring rings is required to be a low-order rational number (3, 2, 1, 4/3, 5/4, or 6/5) to maintain a high degree of symmetries. Our main motivation for this grid is to solve linear systems using multi-grid methods, and to construct efficient preconditioners through pixel-space sampling of the linear operator in question. As a benchmark and representative example, we compute a preconditioner for a linear system that involves the operator D-hat + B-hat {sup T}N{sup −1} B-hat , where B-hat and D-hat may be described as both local and rotationally invariant operators, and N is diagonal in the pixel domain. For a bandwidth limit of ℓ{sub max} = 3000, we find that our new SymPix implementation yields average speed-ups of 360 and 23 for B-hat {sup T}N{sup −1} B-hat and D-hat , respectively, compared with the previous state-of-the-art implementation.
79 ASTRONOMY AND ASTROPHYSICS↗
Enhancing the Representation of Subgrid Land Surface Characteristics in Land Surface Models
Land surface heterogeneity has long been recognized and increasingly incorporated in the land surface modelling. In most existing land surface models, the spatial variability of surface cover is represented as subgrid composition of multiple surface cover types. In this study, we developed a new subgrid classification method (SGC) that accounts for the topographic variability of the vegetation cover. Each model grid cell was represented with a number of elevation classes and each elevation class was further described by a number of vegetation types. The numbers of elevation classes and vegetation types were variable and optimized for each model grid so that the spatial variability of both elevation and vegetation can be reasonably explained given a pre-determined total number of classes. The subgrid structure of the Community Land Model (CLM) was used as an example to illustrate the newly developed method in this study. With similar computational burden as the current subgrid vegetation representation in CLM, the new method is able to explain at least 80% of the total subgrid PFTs and greatly reduced the variations of elevation within each subgrid class compared to the baseline method where a single elevation class is assigned to each subgrid PFT. The new method was also evaluated against two other subgrid methods (SGC1 and SGC2) that assigned fixed numbers of elevation and vegetation classes for each model grid with different perspectives of surface cover classification. Implemented at five model resolutions (0.1°, 0.25°, 0.5°, 1.0° and 2.0°) with three maximum-allowed total number of classes N_class of 24, 18 and 12 representing different computational burdens over the North America (NA) continent, the new method showed variable performances compared to the SGC1 and SGC2 methods. However, the advantage of the SGC method over the other two methods clearly emerged at coarser model resolutions and with moderate computational intensity (N_class = 18) as it explained the most PFTs and elevation variability among the three subgrid methods. Spatially, the SGC method explained more elevation variability in topography-complex areas and more vegetation variability in flat areas. Furthermore, the variability of both elevation and vegetation explained by the new method was more spatially homogeneous regardless of the model resolutions and computational burdens. The SGC method will be implemented in CLM over the NA continent to assess its impacts on simulating land surface processes.
characteristics↗
Analysis of Electric Power Board of Chattanooga Smart Grid Investment
The deployment of a Smart Grid within the EPB electrical system has brought significant value to EPB customers and to EPB operations. The Smart Grid implementation has permitted the capability to automate many of the existing systems, increasing the overall efficiency and decreasing round trip time for many required operation functions. This has both reduced the cost of operations and increased the reliability of customer service. This boost in performance can be primarily linked to two main components: the communication network and the availability of data. The existence of a fiber communication backbone has opened numerous opportunities for EPB. Large data sets from remote devices can be continuously polled and used to compute numerous metrics. The speed of communication has provided a pathway to sophisticated distribution automation that would not otherwise be possible. The data have been shown to provide a much broader view of the actual electrical system that was not previously available, opening many new opportunities in data analytics. This report documents many of the existing systems and operational features of the EPB system. Many of these systems have seen significant improvement with the deployment of Smart Grid technologies. The research behind this report found that in some cases, the value was directly computable or estimable for these enhancements. Currently, several ongoing research topics are under investigation within EPB. Through a grant from TVA, EPB is currently installing a 1.3MW solar photovoltaic system, which will be the first-generation asset owned by EPB.
29 ENERGY PLANNING, POLICY, AND ECONOMY↗