skip to main content


Title: An Efficient Network Solver for Electromagnetic Transient Simulation of Power Systems Based on Hierarchical Inverse Computation and Modification
In both power system transient stability and electromagnetic transient (EMT) simulations, up to 90% of the computational time is devoted to solve the network equations, i.e., a set of linear equations. Traditional approaches are based on sparse LU factorization, which is inherently sequential. In this paper, EMT simulation is considered and an inverse-based network solution is proposed by a hierarchical method for computing and store the approximate inverse of the conductance matrix. The proposed method can also efficiently update the inverse by modifying only local sub-matrices to reflect changes in the network, e.g., loss of a line. Experiments on a series of simplified 179-bus Western Interconnection demonstrate the advantages of the proposed methods.  more » « less
Award ID(s):
1636772
NSF-PAR ID:
10381017
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
2021 North American Power Symposium (NAPS)
Page Range / eLocation ID:
01 to 06
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Electromigration (EM) becomes a major concern for VLSI circuits as the technology advances in the nanometer regime. With Korhonen equations, EM assessment for VLSI circuits remains challenged due to the increasing integrated density. VLSI multisegment interconnect trees can be naturally viewed as graphs. Based on this observation, we propose a new graph convolution network (GCN) model, which is called {\it EMGraph} considering both node and edge embedding features, to estimate the transient EM stress of interconnect trees. Compared with recently proposed generative adversarial network (GAN) based stress image-generation method, EMGraph model can learn more transferable knowledge to predict stress distributions on new graphs without retraining via inductive learning. Trained on the large dataset, the model shows less than 1.5% averaged error compared to the ground truth results and is orders of magnitude faster than both COMSOL and state-of-the-art method. It also achieves smaller model size, 4X accuracy and 14X speedup over the GAN-based method. 
    more » « less
  2. null (Ed.)
    The interconnection of distributed energy resources (DERs) in microgrids (MGs) operating in both islanded and grid-connected modes require coordinated control strategies. DERs are interfaced with voltage source inverters (VSIs) enabling interconnection. This paper proposes a load demand sharing scheme for the parallel operation of VSIs in an islanded voltage source inverter-based microgrid (VSI-MG). The ride-through capability of a heavily loaded VSI-MG, where some of the VSIs are fully loaded due to the occurrence of an event is investigated. In developing analytical equations to model the VSI, the concept of virtual synchronous machines (VSM) is applied to enable the VSI mimic the inertia effect of synchronous machines. A power frame transformation (PFT) that takes the line ratios of the MG network into account is also incorporated to yield satisfactory transient responses of both network frequency and bus voltages in the MG network. A Jacobian-based method is then developed to take into account the operational capacity of each VSI in the VSI-MG. The resulting amendable droop control constrains the VSIs within their power capabilities when an event occurs. Simulation results presented within demonstrate the effectiveness of the proposed procedure which has great potential to facilitate efforts in maintaining system reliability and resiliency. 
    more » « less
  3. Consider the scattering of a time-harmonic acoustic plane wave by a bounded elastic obstacle which is immersed in a homogeneous acoustic medium. This paper is concerned with an inverse acoustic-elastic interaction problem, which is to determine the location and shape of the elastic obstacle by using either the phased or phaseless far-field data. By introducing the Helmholtz decomposition, the model problem is reduced to a coupled boundary value problem of the Helmholtz equations. The jump relations are studied for the second derivatives of the single-layer potential in order to deduce the corresponding boundary integral equations. The well-posedness is discussed for the solution of the coupled boundary integral equations. An efficient and high order Nyström-type discretization method is proposed for the integral system. A numerical method of nonlinear integral equations is developed for the inverse problem. For the case of phaseless data, we show that the modulus of the far-field pattern is invariant under a translation of the obstacle. To break the translation invariance, an elastic reference ball technique is introduced. We prove that the inverse problem with phaseless far-field pattern has a unique solution under certain conditions. In addition, a numerical method of the reference ball technique based nonlinear integral equations is proposed for the phaseless inverse problem. Numerical experiments are presented to demonstrate the effectiveness and robustness of the proposed methods. 
    more » « less
  4. The inverse problem of recovery of a potential on a quantum tree graph from the Weyl matrix given at a number of points is considered. A method for its numerical solution is proposed. The overall approach is based on the leaf peeling method combined with Neumann series of Bessel functions (NSBF) representations for solutions of Sturm–Liouville equations. In each step, the solution of the arising inverse problems reduces to dealing with the NSBF coefficients. The leaf peeling method allows one to localize the general inverse problem to local problems on sheaves, while the approach based on the NSBF representations leads to splitting the local problems into two‐spectrum inverse problems on separate edges and reduces them to systems of linear algebraic equations for the NSBF coefficients. Moreover, the potential on each edge is recovered from the very first NSBF coefficient. The proposed method leads to an efficient numerical algorithm that is illustrated by numerical tests.

     
    more » « less
  5. The increase in penetration levels of inverter-based resources (IBRs) is changing the dynamic performance of power grids of different parts of the world. IBRs are now being more and more integrated into the grid at a single connection point as an IBR plant. Due to the complex nature and dynamicity of each inverter model, it is not realistic to build and analyze full complex models of each inverter in the IBR plant. Moreover, simulating a large plant including detailed models of all the IBRs would require high computing resources as well as a long simulation time. This has been the main issue addressed in the new IEEE Std 2800-2022. This paper proposes a novel approach to model an IBR plant, which can capture the transient nature at the plant level, detailed IBR control at the inverter level, interactions of multiple IBR groups in a plant structure, and a collector system connecting the IBRs to the grid. The IBRs in the plant use a voltage source inverter topology combined with a grid-connected filter. The control structure of the IBR includes a cascaded loop control where an inner current control and outer power control are designed in the dq-reference frame, and a closed-loop phase-locked loop is used for the grid synchronization. The mathematical study is conducted first to develop aggregated plant models considering different operating scenarios of active IBRs in an IBR plant. Then, an electromagnetic transient simulation (EMT) model of the plant is developed to investigate the plant’s dynamic performance under different operating scenarios. The performance of the aggregated plant model is compared with that of a detailed plant model to prove the effectiveness of the proposed strategy. The results show that the aggregated EMT simulation model provides almost the same result as the detailed model from the plant perspective while the running time/computation burden is much lower.

     
    more » « less