Electromagnetic Transients Research Articles

Multimessenger emission from tidal waves in neutron star oceans

ABSTRACT Neutron stars in astrophysical binary systems represent exciting sources for multimessenger astrophysics. A potential source of electromagnetic transients from compact binary systems is the neutron star ocean, the external fluid layer encasing a neutron star. We present a groundwork study into tidal waves in neutron star oceans and their consequences. Specifically, we investigate how oscillation modes in neutron star oceans can be tidally excited during compact binary inspirals and parabolic encounters. We find that neutron star oceans can sustain tidal waves with frequencies between 0.01 and 20 Hz. Our results suggest that tidally resonant neutron star ocean waves may serve as a never-before studied source of precursor electromagnetic emission prior to neutron star–black hole and binary neutron star mergers. If accompanied by electromagnetic flares, tidally resonant neutron star ocean waves, whose energy budget can reach 1046 erg, may serve as early warning signs (≳1 min before merger) for compact binary mergers. Similarly, excited ocean tidal waves will coincide with neutron star parabolic encounters. Depending on the neutron star ocean model and a flare emission scenario, tidally resonant ocean flares may be detectable by Fermi and Nuclear Spectroscopic Telescope Array (NuSTAR) out to ≳100 Mpc with detection rates as high as ∼7 yr−1 for binary neutron stars and ∼0.6 yr−1 for neutron star–black hole binaries. Observations of emission from neutron star ocean tidal waves along with gravitational waves will provide insight into the equation of state at the neutron star surface, the composition of neutron star oceans and crusts, and neutron star geophysics.

Analysis of Cross-Bonded Cables Using Accurate Model Parameters

This paper presents a thorough investigation of frequency and electromagnetic transient (EMT) responses of cross-bonded cables. The study adopts newly developed Line/ Cable Data (LCD) Module in an existing EMT-type simulation software. The input impedances evaluated by positive sequence and various mode energizations on a typical cross-bonded cable (CBC) in frequency domain are investigated.In addition to cross-bonded, the solid-bonded cable (SBC) and core-transposed cable (CTC) are analyzed for characteristics of input impedance. The mixed propagation modes of CBC and CTC are studied. The influence of bonding methods and cable models on transient voltages is investigated. The results simulated by the new LCD are compared with the Classical Parameter Calculation (CPC) method and validated by recently proposed extended transmission line (ETL) approach.

Full-Wave 3-D Transient Analysis With Method of Moments and Numerical Laplace Transform Including Resistive Non-Linear Elements

Methods based on both time ( <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">t</i> )-domain and frequency ( <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">f</i> )-domain have been used for electromagnetic transient (EMT) analysis using both circuit-theory and field-theory approaches. Methods based on circuit theory are fast and convenient but may become inapplicable to many engineering problems due to simplifications such as the transverse electromagnetic assumptions. Moreover, <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">t</i> -domain-based techniques rely on approximate methods to account for frequency dependence, which is an essential characteristic of electric conductors, equipment and soil environment. In this paper, <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">f</i> -domain solutions obtained from the method of moments (MoM) discretization of the electric field integral equation (EFIE) are converted to <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">t-</i> domain using the numerical Laplace transform (NLT) by means of post-processing. This circumvents the need to formulate the EFIE and MoM based on the complex frequency which would be required for their application to conventional NLT. Hence, existing MoM implementations can be used to perform a full-wave 3-D EMT analysis for a wide range of power system applications. Examples include energizations of 3-D power system networks with fast and non-vanishing excitations such as step functions, as well as modeling non-linear elements using the piecewise linear approximation. Results from experimental measurements, finite-difference <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">t</i> -domain method, and EMT-type software confirm the validity of the proposed method for power system transients in the range of microseconds down to nanoseconds.

Detailed EMT-Type Load Modeling for Power System Dynamic and Harmonic Studies

The aim of this paper is to develop a detailed model and investigate its importance for capturing load behavior in both dynamic and harmonic simulations. A detailed electromagnetic transient (EMT) type model of a mixed residential-commercial load is developed using a component-based technique, and its dynamic characteristic is validated against existing field measurements. The model is subsequently used to study the harmonic impedance of an existing underground distribution network. Errors caused by model inadequacies are also studied in both dynamic and harmonic case studies.

Fault current control of MMC in HVDC-connected offshore wind farm: A coordinated perspective with current differential protection

The modular multilevel converter (MMC) based high-voltage dc (HVDC)-connected offshore wind farm (OWF) is power electronic converter dominated power system, where the conventional current differential protection may not operate effectively. Based on the basic operation principle of the current differential relay, this paper points out that the phase difference (ϕdiff) of fault currents from MMC and OWF should be limited within a specific range, i.e. (-ϕdiffmax, ϕdiffmax), to guarantee the reliable tripping of the current differential relay. Yet, the highly control-dependent fault currents from the MMC and OWF could yield an arbitrary ϕdiff that might beyond (-ϕdiffmax, ϕdiffmax), which leads to the possible malfunction of the current differential relay. To tackle this challenge, the coordinated control method of MMC is proposed to align the phase angle of its fault current with that of the OWF, such that ϕdiff ∈ (-ϕdiffmax, ϕdiffmax) can be always guaranteed, so as the reliable operation of the current differential relay. Finally, the proposed coordinated control is verified by the electromagnetic transient (EMT) simulations in PSCAD and the real-time digital simulator (RTDS).

Improvement of distance protection algorithms based on phasor measurements

Distance protection of overhead lines is widely used and is continually being improved. However, the distance protection cannot work stably and is blocked in several electromechanical transients. Additionally, a non-linear arc can affect the stability of the distance protection in case of a short-circuit. The authors consider ways to improve distance protection during such transients. In particular, this paper proposes the use of synchronized phasor measurement technology to improve the efficiency of distance protection. This paper considers a new approach for the analysis of electromechanical and electromagnetic transients in a power system using process synchrophasors. As examples, cases of a short-circuit in a line with oneway and two-way powers are considered. The authors consider the influence of a non-linear arc at the fault point and consider the possibility of the coincidence of the electromechanical and electromagnetic transients. The result is a new algorithm for operating distance protection, an improvement in the parametric method for determining the location of line damage. The calculations performed and mathematical modeling confirm the effectiveness of the proposed algorithm.

Evaluation of Inherent Damping Introduced by Thyristor Controlled Series Compensators

Series compensation of transmission lines increase the risk of Sub-Synchronous Resonance (SSR) problems. Use of Thyristor Controlled Series Compensators (TCSC) in place of Fixed Series Capacitors (FSC) has benefited in terms of damping SSR. SSR Damping with TCSC can be achieved with supplementary controls or through its inherent damping capability. Use of discrete time models for small signal stability analysis of TCSC is well established. This paper uses a simplified Dynamic Phasor (DP) model of the TCSC for an in-depth evaluation of the inherent SSR damping capability through small signal stability analysis. Source of inherent damping in open loop operation and the representation of it in the DP model is elaborated and the small signal model is validated against Electro Magnetic Transient (EMT) simulations. Effect of inductive-resistive nature of the TCSC on Induction Generator Effect (IGE) and Torsional Interactions (TI) is shown through small signal stability assessment. It is shown that the network resonant frequency exists only in a narrow range of firing angles of the TCSC and thus IGE and TI are avoided in most part of its operating range. IEEE first benchmark system for SSR studies is used to demonstrate the inherent damping capability of a TCSC.

SSR Stable Wind Speed Range Quantification for DFIG-Based Wind Power Conversion System Considering Frequency Coupling

The subsynchronous resonance (SSR) characteristics of doubly-fed induction generation-based wind power conversion system (DFIG-WPCS) varies with the wind speed. Obtaining the small signal SSR stable wind speed range is urgently important for wind farm dispatch planning. However, there are few methods to quantify the small signal SSR stable wind speed range of DFIG-WPCS. Moreover, the inherent frequency coupling when DFIG-WPCS is integrated to modular multilevel converter-high voltage direct current (MMC-HVDC) may lead to inaccurate stability assessment results. In this article, the frequency coupling mechanism between WPCS and MMC is revealed by analyzing the law of harmonic interaction in the interconnected system. The wind speed-frequency two-variable equivalent admittance model of DFIG-WPCS is established by taking the frequency coupling account. Relying on the two-variable open-loop transfer function, the multiple phase margin contours plot approach is proposed to quantify the small signal SSR stable wind speed range and damping properties. The effectiveness of two-variable admittance and analysis method are validated against electromagnetic transient (EMT) testbed with MATLAB /SimPowerSystems. The analysis results can provide a reference for the early planning and stable operation of DFIG-WPCS.

Analysis of Power Quality for Distribution Networks Using Active Compensator

This paper concentrates on compensating the power quality issues which have been increased in day-to-day life due to the enormous usage of loads with power electronic control. One such solution is compensating devices like Pension Protection Fund (PPF), Active power filter (APF), hybrid power filter (HPF), etc., which are used to overcome Power Quality (PQ) issues. The proposed method used here is an active compensator called unified power quality conditioner (UPQC) which is a combination of shunt and series type active filter connected via a common DC link. The primary objective is to investigate the behavior of the compensators in the distribution networks. The performance of two configurations of UPQC, Right Shunt UPQC (RS-UPQC) and Left Shunt UPQC (LS-UPQC) are tested in the distribution networks under various load conditions by connecting them at the source side of harmonic generation using a specially constructed transformer called inductively filtered converter transformer which adopts special wiring scheme at the secondary side. PSCAD (Power Systems Computer Aided Design)/EMTDC (Electromagnetic Transients with DC Analysis) software is used to model the compensators connected to the nonlinear load. Both RS-UPQC and LS-UPQC are tested at the distribution side of the supply system with Hysteresis current controller for shunt and Sinusoidal pulse with modulation controller for series at various locations of power system network and their results are compared.

Faster-Than-Real-Time Hardware Emulation of Extensive Contingencies for Dynamic Security Analysis of Large-Scale Integrated AC/DC Grid

The rapid expansion of modern power systems has brought a tremendous computational challenge to dynamic security analysis (DSA) tools which consequently need to process extensive contingencies. In this work, hardware emulation is investigated to accelerate the DSA solution of a large-scale AC/DC system deployed on the field-programmable gate arrays (FPGAs) faster-than-real-time (FTRT) execution. Electromagnetic transient (EMT) modeling of the DC grid is conducted since the fast converter dynamics require a small time-step for accuracy; in contrast, the transient stability (TS) simulation is applicable to the AC grid which tolerates a much larger step size. To coordinate the 2 different types of simulation, an interface based on dynamic voltage injection is proposed to integrate the AC and DC grids, in addition to maintaining a low hardware latency. An emulation platform consisting of multiple FPGA boards is established so that with a proper allocation it has a sufficient capacity to accommodate the system under study which has 6 ACTIVSg 500-bus systems interconnected by a 6-terminal DC grid. The efficacy of the proposed FTRT hardware emulation platform is demonstrated by 2 case studies with more than 5500 contingencies analyzed in total, where an FTRT ratio of more than 208 is achieved for the hybrid AC/DC grid, while it is over 277 times for a single 500-bus system. Furthermore, the FTRT dynamic emulation results, including the security indices, are validated by the simulation tool DSATools/TSAT.

Faster-Than-Real-Time Hardware Emulation of Transients and Dynamics of a Grid of Microgrids

Enhanced environmental standards are leading to an increasing proportion of microgrids (MGs) being integrated with renewable energy resources in modern power systems, which brings new challenges to simulate such a complex system. In this work, comprehensive modeling of a grid of microgrids for faster-than-real-time (FTRT) emulation is proposed, which can be utilized in the energy control center for contingencies analysis and dynamic security assessment. Electromagnetic transient (EMT) modeling is applied to the microgrid in order to reflect the detailed device processes of the converter and renewable energy sources, while the AC grid utilizes the transient stability modeling to reduce the computational burden and obtain a high acceleration value over real-time execution. Consequently, a dynamic power injection interface is proposed for the coexistence of the two simulation types. The reconfigurability and parallelism of the field-programmable gate arrays (FPGAs) enable the whole system to be executed in FTRT mode with 51 times acceleration over real-time. Meanwhile, three case studies are emulated and the results are validated by the off-line simulation tool Matlab/Simulink®.

Data-Driven Interarea Oscillation Analysis for a 100% IBR-Penetrated Power Grid

In this research, a 100% inverter-based resource (IBR)-penetrated bulk power system (BPS) is examined for possible interarea oscillations. A testbed of a BPS constructed in PSCAD, with both grid-forming inverter (GFM)-based IBRs and grid-following inverter (GFL)-based IBRs, is utilized for simulation. For linear analysis, a data-driven method is used to obtain the <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$dq$ </tex-math></inline-formula> admittance models of the IBRs. This information, combined with our knowledge of the network admittance, leads to eigenvalue analysis of the power grid. Furthermore, for the dominant 3-Hz mode, extended frequency-domain mode shape analysis, along with subgroup effect analysis, and network decomposition analysis are conducted. The analysis results can precisely tell the influencing factors of the oscillation mode. This mode is identified as an interarea oscillation mode. It has a different characteristics compared to the traditional electromechanical interarea mode. Specifically, the network decomposition analysis provides new insights on interarea oscillation modes. The analysis results are all confirmed by the electromagnetic transient (EMT) simulation results.

Small-signal stability analysis of the asymmetrical DC operation in HVDC networks

• Stability of interconnection of monopolar and bipolar HVDC stations is analysed. • A small-signal stability model for asymmetric DC operation studies is developed. • The model is validated with EMT dynamic simulation. • The findings are useful for systematic design of the HVDC grounding systems. • A controller is designed to enhance stability during asymmetric DC operation. The increase in High Voltage Direct Current (HVDC) projects favors their future interconnection to form regional networks, but this interconnection can be hindered by the lack of standardization in station's configuration and converter's differing grounding schemes. Whilst other studies have studied the effect of Direct Current (DC) asymmetrical operation, they are limited in their scope. This paper proposes a more realistic perspective, considering stations with different topologies, permanent DC asymmetries and the influence of the grounding systems. This more complex point of departure requires the re-examination of the modeling methods, from time-domain to frequency-domain. This paper proposes a small-signal stability model of the system, which is validated with the electromagnetic transient (EMT) dynamic simulations using PowerFactory. This model provides guidelines for designing the grounding impedance and allows the proper design of a specific controller to dampen the asymmetrical DC operation mode.

A New Type of Weak Grid IBR Oscillations

This letter presents the discovery of a new type oscillations associated with inverter-based resources (IBRs) in weak grids. While typical weak grid oscillations become worse when the grid strength decreases and the power exporting level increases, this new type of oscillations appear when the grid becomes weaker and the power exporting level decreases. Electromagnetic transient (EMT) simulation results are presented, along with eigenvalue loci for a varying power. The notable signatures of such oscillations are also remarked. In further analysis, a feedback system consisting of two blocks is derived with one block associated with the phase-locked-loop (PLL) and the other associated with the rest of the system. This feedback system quantitatively demonstrates why insufficient damping in the PLL may be an issue for stability. Open-loop analysis shows that a larger phase lag in angle tracking in the low-frequency range due to high power exporting level at the weak grid condition is beneficial to the open-loop poles introduced by the PLL. And this is the reason why such oscillations may appear when power is ramping down.

A Generalized Fault Location Approach for Any Type of Power Network

This paper presents a generalized fault location algorithm that can locate faults for any type of power network configuration. It enhances previous work presented by the authors by reducing the estimation steps but keeping its robustness and redundancy against traveling wave attenuation by power transformers, bad data, voltage signal errors and the impact of Inverter Based Power Sources (IBPSs). This is accomplished by forming arrays containing the network parameters and time-of-arrivals (ToAs) for fault signals captured by Digital Fault Recorders (DFRs), which are used as inputs to the algorithm to compute the accurate fault localization. The IEEE 39-bus system is used to validate the algorithm under different fault conditions using Electromagnetic Transients Program and Matlab as simulation tools.

Time-Domain Protection Scheme for Microgrids With Aggregated Inverter-Based Distributed Energy Resources

This paper proposes a time-domain-based protection scheme for radial and loop microgrid systems with inverter-based resources (IBRs), such as solar photovoltaic (PV) and type-4 wind turbines. The protection scheme is designed to function during both grid-interconnected and grid-isolated modes. The proposed scheme provides an ultra-high-speed sub-cycle directional element aided with low bandwidth communication between relays. Like directional comparison schemes, relays identify whether faults are in-zone or out-zone. The directional element is based on time-domain superimposed quantities and Park’s transformation algorithms. More specifically, the element calculates the superimposed positive-sequence direct component of transient energy during faults. Superimposed voltage and current quantities are calculated using delta filters and decoupled double synchronous reference frame (DDSRF) filters. The proposed filtering method improves the reliability of the superimposed directional element when IBRs are present. The protection scheme is evaluated on a modified IEEE 34-bus distribution system simulated using an electromagnetic transients program.

МОДЕЛИРОВАНИЕ АСИНХРОННОГО ЭЛЕКТРОПРИВОДА ДЛЯ ИССЛЕДОВАНИЯ ПРОЦЕССА ПУСКА

A model of an asynchronous motor is proposed, which makes it possible to investigate the process of motor acceleration, taking into account changes in the resistance of the rotor circuit, consumed current and voltage in the power supply network. A method for calculating the parameters of the equivalent circuit of an asynchronous motor according to passport data is proposed, the calculation takes into account changes in the resistance of the rotor circuit during start-up. The electric motor model is based on a static mechanical characteristic with a linear dependence of the rotor resistance on slip. The model is designed to study the electromechanical parameters of the electric drive during the start-up process, taking into account the parameters of the power supply circuit. The analysis of existing functions for modeling the mechanical characteristics of an induction motor was carried out and a formula was chosen for determining the electromagnetic torque from the parameters of a simplified L-shaped equivalent circuit without taking into account electromagnetic transients. Therefore, the proposed model based on the static mechanical characteristic is applicable for mechanisms with a long start. The parameters of the equivalent circuit of an asynchronous motor are determined according to the passport data. The dependence of the rotor circuit resistance on slip is described by a linear function of the first order. The results of calculating the mechanical characteristic according to the compiled model do not differ from the parameters at four known characteristic points corresponding to the passport data of the induction motor. Because when calculating the no-load current and the active resistances of the stator and rotor windings, some assumptions were made, it should be noted that the proposed method is preferable for calculating high-power motors.

Experimental and Analytical Studies on Lightning Surge Response of the Quadruple-Circuit Transmission Line

The compact design of a 500 kV quadruple-circuit transmission line can effectively reduce the line corridor area, but the height of the tower also increases, increasing the probability of suffering a lightning strike. The 500 kV quadruple-circuit transmission lines carry more energy, and because of this, lightning strikes that cause power line trips are more likely to result in large-scale power outages. Therefore, it is necessary to make an accurate assessment of the lightning performance of 500 kV quadruple-circuit transmission lines. First, simulated lightning-striking experiments were carried out on a scaled 500 kV quadruple-circuit transmission line in the laboratory, where transient voltages and currents were measured. Second, a numerical model was established with the Finite-Difference Time Domain (FDTD) method, which was then verified with the experimental results. Third, lightning surge responses of a 500 kV quadruple-circuit transmission line under near-real facility conditions are estimated with the verified FDTD model. In the simulations, influencing factors, such as the rise time of injecting current, the velocity of return-stroke current and the terrains, were taken into consideration, but not in previous lightning surge analysis with the Electromagnetic Transients Program (EMTP). Results show that insulator voltages on the same tower cross-arm are nearly identical, although the length of the cross-arm is large enough. Furthermore, it is found that the rise time and the lightning current velocity have great effects on the lightning surge response, and the terrains are less impactful but not negligible. Therefore, these factors should be considered carefully where higher accuracy lightning protection design is necessary.

Modeling for large-scale offshore wind farm using multi-thread parallel computing

Large-scale Offshore Wind Farms (OWFs) face difficulty when carrying out microsecond-level Electro-Magnetic Transient (EMT) simulations due to the large number of Wind Turbines (WTs) and electrical nodes. This paper provides an OWF modeling method considering the electrical collection network based on the integrated equivalent modeling method of single WT, which can eliminate the internal nodes and reduce the order of admittance matrix. Meanwhile, the multi-thread parallel computing is integrated into the whole solution procedure, which achieves further simulation speedup. Moreover, simulations are performed in PSCAD/EMTDC to compare the model accuracy and time efficiency of the proposed method with the detailed model. The simulation results show that the proposed method can accurately simulate the dynamic characteristics of large-scale OWFs with large speedup factors.

An Improved High-Resistance Fault Detection Method in DC Microgrid Based on Orthogonal Wavelet Decomposition

High-resistance faults in direct current (DC) microgrids are small and thus difficult to detect. Such faults may be “invisible” in that grid operation continues for a considerable time, which damages the grid. It is essential to detect and remove high-resistance faults; we present a detection method herein. First, the transient DC current during the fault is subjected to hierarchical wavelet decomposition to identify high-resistance faults accurately and sensitively; the wavelet coefficients are detected using the singular value decomposition (SVD) method. The SVD valve can denoise the dc microgrid fault current, which eliminates the influence of converter switching frequency and background noise effectively. Power system computer-aided design (PSCAD)/electromagnetic transients including direct current (EMTDC)-based simulations showed that our method successfully identified high-resistance faults.