Next Article in Journal
Analysis of the Effect of Parameter Variation on a Dynamic Cost Function for Distributed Energy Resources: A DER-CAM Case Study
Previous Article in Journal
Numerical Simulation of the Percolation Threshold in Non-Overlapping Ellipsoid Composites: Toward Bottom-Up Approach for Carbon Based Electromagnetic Components Realization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extended Critical Damping Adjustment Method for Electromagnetic Transients Simulation in Power Systems

1
College of Electrical Engineering & New Energy, China Three Gorges University, Yichang 443002, China
2
China Shanghai Administration Bureau of State Grid Operation Company, Shanghai 200000, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2018, 8(6), 883; https://doi.org/10.3390/app8060883
Submission received: 13 April 2018 / Revised: 16 May 2018 / Accepted: 22 May 2018 / Published: 28 May 2018

Abstract

:
Numerical oscillation is a common phenomenon in the current electromagnetic transient digital simulation in power systems. From the origin of the numerical oscillation, we proposed an extended critical damping adjustment method (ECDA) to solve the numerical oscillation problem in this paper, after inspired by the critical damping adjustment method (CDA). ECDA was composed of a class of new low-order, L-stable linear multistep and implicit trapezoidal integration methods. ECDA used two-half-step, three-step or four-steps linear multistep method to perform numerical integrations when experienced with sudden changes, and the implicit trapezoidal integration method was still used under normal conditions. The numerical simulation results of the typical electromagnetic transient examples show that the ECDA could completely eliminate the numerical oscillation, and it had higher computational accuracy than the CDA.

Graphical Abstract

1. Introduction

The phenomenon of numerical oscillations is a common problem in the electromagnetic transient simulation in power systems [1], which is due to the reason that the numerical integration algorithm in power-system transient simulation calculation is mainly an implicit trapezoidal integration method [1,2]. However, the implicit trapezoidal integration method is not L-stable and does not have the ability to eliminate numerical oscillations [3]. Therefore, the traditional power-system electromagnetic transient simulation software, such as EMTP 3.0, use critical damping adjustment (CDA) to eliminate numerical oscillations [4,5].
The CDA method is essentially a combined integration algorithm. The core step of the CDA method is to use two-half-step backward Euler method for numerical integrations to eliminate numerical oscillations [6,7]. The backward Euler method is a one-step, low-order, L-stable numerical integration algorithm, which eliminates numerical oscillations from its source.
In addition, other schemes for solving the problem of numerical oscillations indirectly use different technical approaches to approximating non-state variables for obtaining their “real” values [1,5]. For example, PSCAD version 4.0 uses a half-step interpolation algorithm to eliminate numerical oscillations, but this method is only suitable for offline simulations [8].
Based on the state space method in electromagnetic transient simulations [9,10], we proposed an extended critical damping adjustment method (ECDA), in this paper, to solve the problem of numerical oscillations directly from the origin of the numerical oscillations. The core idea of the ECDA is similar to that of CDA. The difference between these two methods lies in the use of a new class of low-order linear multistep formulae (L2MF) to avoid the numerical oscillations. The absolute value of the main coefficient of the local truncation error of a new four-step L2MF using Chebyshev grid points was smaller than that of the backward Euler method. The pure L2MF method and the implicit trapezoidal integration method were used to form the ECDA to solve numerical examples. The simulation results show that the L2MF method can effectively avoid the numerical oscillations, and also indirectly show that the ECDA has a better simulation effect than the CDA in theory, and hence has the potential applications in EMTP-type programs.

2. A New Class of Linear Multistep Formulae

2.1. Construction of a New Class of Linear Multistep Formulae by Classic Differential Quadrature

The initial-value problem of an ordinary differential equation was shown in the following:
{ x ˙ ( t ) = f ( t , x ) x ( t = 0 ) = x 0
where t is time and x is a one-dimensional variable; f ( t , x ) is a function of the time t and the variable x ; and x 0 is the initial value of an unknown state variable.
One of main methods for solving the initial-value problem, as indicated in Equation (1), is the use of classical linear multistep formulae (LMF) [11], and one kind of the LMF is the following backward differentiation formula (BDF), which was described as following:
i = 0 s α i x n + i = h f n + s
where s is the number of steps; α i is integration coefficient; h denotes the integration step; x n + i x ( t n + i ) is the approximate value of a state variable at a time point t n + i ; and f n + s = f ( t n + s , x n + s ) .
Obviously, when s takes one, Equation (2) is the first-order backward Euler method that is A-stable and L-stable. Based on the traditional BDF, new linear multistep methods were constructed as follows:
i = 0 s α ˜ i x n + i = h β s f n + s
where α ˜ i and β s are integration coefficients of new LMF, respectively.
In order to obtain the unknown integration coefficients of new LMF, we similarly used s-stage s-order time-domain differential quadrature method (DQM) [12,13,14] to solve Equation (1). For any time-step ( h = t n + 1 t n ), we have the following equation:
G [ x ˜ n + 1 x ˜ n + s ] + G 0 x n = h [ f ( t ˜ 1 , x ˜ n + 1 ) f ( t ˜ s , x ˜ n + s ) ]
where c i = ( t ˜ i t n ) / h , i ( 0 , s ) is the non-equidistant grid point, which is distributed over the regularization interval [ 0 , 1 ] ; t ˜ 0 = t n is a time-starting point; t ˜ s = t n + 1 is a time-end point; t ˜ i = t n + c i h , i ( 1 , s ) are internal points; and x ˜ n + i x ( t ˜ i ) is the approximate value of a state variable at a time-internal point t ˜ i .
G = [ g 11 g 1 s g s 1 g s s ] ,   G 0 = [ g 10 g s 0 ]
where G is called the weighted coefficient matrix of DQM, and the calculation formulae of each element in G is given in Reference [15] G 0 = G e , where e is the column vector with each element equaling to one.
In Equation (4), the grid point distribution has to be a non-equidistant grid, including Chebyshev grid, Chebyshev-Gauss-Lobatto grid, and Legendre grid [16]. In this paper, we used the DQM with Chebyshev grid points [16] for constructing the new LMF. Suppose:
α ˜ i = g s i / g s s , i ( 0 , s ) ,   β s = i = 1 s i α ˜ i
where α ˜ s = 1 .
As an example, the new three-step LMF using Chebyshev grid points was shown below:
15 4 2 9 h f n + 3 = 1 9 x n + 4 2 4 9 x n + 1 4 2 + 4 9 x n + 2 + x n + 3
The new four-step LMF using Chebyshev grid points is as follows:
108 32 3 57 h f n + 4 = 1 19 x n + 16 3 32 57 x n + 1 + 4 57 x n + 2 16 3 + 32 57 x n + 3 + x n + 4
In addition, the backward Euler method is obtained when s takes one in Equation (3) regardless of the kind of the used grid. Similarly, when s takes two, we get the second-order Gear method.

2.2. The Compatibility and Convergence of the L2MF

The new linear multistep methods have practical applications only if they satisfy compatibility and convergence. Therefore, according to Equation (3), the first-characteristic polynomial and the second-characteristic polynomial of z were as follows [17]:
ρ ( z ) = i = 0 s α ˜ i z i ,   σ ( z ) = β s z s
where ρ ( 1 ) = 0 and ρ ( 1 ) = σ ( 1 ) . Therefore, the new linear multistep methods satisfy the compatibility. Suppose ρ ( z ) = 0 . A unique root modulus is | z 1 | = 1 , and the rest of the roots’ modulus satisfies | z i | < 1 , i ( 2 , s ) . Therefore, it also meets the root condition. Since the new linear multistep formulae satisfy both the compatibility and the root condition, the new linear multistep formulae are convergent according to the Dahlquist equivalence theorem [17,18].

2.3. Stability of the L2MF

For the stability analysis of the new linear multistep method, the following standard test equations can be considered:
{ d x d t = λ x x ( 0 ) = 1
where Re ( λ ) < 0 .
According to Equation (3), the corresponding characteristic equation was obtained:
ρ ( z ) q σ ( z ) = 0 ,   q = λ h
where the stability domains of different s values in Equation (3), supposing, were shown in Figure 1 and Figure 2. In Figure 1 and Figure 2, the outside regions of boundary curves were the stability domain of the new LMF. Obviously, the numerical stability regions of the new linear multistep methods contained the entire left-half complex plane, so the new LMF were A-stable.
According to the Dahlquist barriers theorem [19], new two-step LMF were second-order methods, and other LMF were first-order methods. The main coefficient of the local truncation error for the three-step LMF and the four-step LMFs are 0.6683 and 0.3009 , respectively.
In addition, Equation (9) was substituted into Equation (3); hence the following equation was obtained:
x n + s x n + s 1 = α ˜ s 1 1 β s λ h + i = 0 s 2 α ˜ i x n + i x n + s 1 1 β s λ h
where | x n + s / x n + s 1 | 0 if Re ( λ h ) . Therefore, the new linear multistep formulae were L-stable [19].
In summary, the new linear multistep methods are low-order methods. So it is named as low-order LMF (L2MF) method. In numerical stability, the L2MF method was both A-stable and L-stable. Therefore, compared with the A-stable second-order implicit trapezoidal integration method, the L2MF method has strong damping characteristics, which can effectively suppress the numerical oscillations, and hence is suitable for electromagnetic transient simulations in power systems.

3. Extended Critical Damping Adjustment Scheme

From the above analysis, it can be seen that the L2MF method belongs to the numerical methods with low-order and strong numerical stability. For this purpose, it is necessary to combine the implicit trapezoidal integration method, which is widely used in electromagnetic transient calculation, to calculate the implicit trapezoidal integration method under normal circumstances, and to switch to L2MF only when the system is mutated, such as switch operation and line fault. In practice, we can use the same method of auxiliary detection mutation as the CDA method and the algorithm switching strategy. Because the absolute value of the local truncation error of the four-step L2MF method is smaller than that of the backward Euler method (the absolute value of the main coefficient of local truncation error is 0.5), it was recommend that a four-step L2MF method and an implicit trapezoid integration method are used to form that combined integration algorithm for the electromagnetic transient simulation in power systems.
In the new combination integration scheme, the implicit trapezoidal integration method had the ability of self-starting computation, and the L2MF method did not have the self-starting calculation ability, which needed a startup algorithm for independent calculation. The starting numerical method used in this paper was explicit Euler method (forward Euler method). Figure 3 is a flow chart of electromagnetic transient simulation in power systems on a basis of extended CDA method (ECDA).
The combined integration algorithm proposed in this paper was a class of algorithms, and the combination algorithm was the CDA method when s takes one. Therefore, the new combination algorithm proposed in this paper was also named as the extended CDA method.

4. Simulation Examples

4.1. Simulation of Half-Wave Rectifier Circuit

The example is half-wave rectifier circuit. As shown in Figure 4, the excitation source was an ideal sinusoidal voltage source E ( t ) = 100 2 sin ( 100 π t ) . The value of the DC voltage source is 20 V . The first-order differential equation for the circuit is as follows:
d i d t = 1 L ( E ( t ) d ( t ) 20 R i )
where i is the current, R is the resistance, and L is the inductance.
A three-step L2MF method and an implicit trapezoidal integration method were used to solve Equation (13), and the integration step in this case was set to h = 0.1 ms. The simulation results are shown in Figure 5.
As shown in Figure 5, since the three-step L2MF method was L-stable and had a strong damping characteristic, which could suppress the numerical oscillations of non-state variables when the diode D turned on. In contrast, the implicit trapezoidal integration method used in the style of EMTP produces a serious numerical oscillation when the network changed suddenly. It is worth mentioning that the implicit trapezoidal integration method directly used for solving Equation (13) did not produce numerical oscillations.

4.2. Simulation of High-Voltage Transmission Line without Load Being Switched-In Suddenly

This example is an electromagnetic transient simulation of a single-phase middle-long transmission line with distributed parameters. The initial phase of the sinusoidal voltage source was φ 0 = 90 ° . The other parameters of the transmission line were given in Figure 6.
The mathematical model of the electromagnetic transient process for the transmission line was described by the telegraph equation as follows:
{ v ( y , t ) y = l 0 i ( y , t ) t + r 0 i ( y , t ) i ( y , t ) y = c 0 v ( y , t ) t + g 0 v ( y , t )
where y and t are the spatial variables of the length and the time variable, respectively; v ( y , t ) and i ( y , t ) denote the changes of the voltage and the current with respect to the space and the time, respectively; l 0 , c 0 , r 0 and g 0 are the unit parameters of inductance, capacitance, resistance, and susceptance, respectively.
In this paper, we adopt the time-domain method to settle the telegraph Equation (14). First, in order to obtain the voltage and current variables along the transmission line, we decompose the whole line into N sub-intervals uniformly, as is shown in Figure 7. Supposing the total length of the transmission line is l , we define: the length of each small interval is Δ y = l / N . Let v k denote the voltage at the point, where y = y k = k Δ y , k ( 0 , N ) ; i k , k ( 1 , N ) is the current from the point y k 1 to the point y k ; and i 0 and i N + 1 are the currents of two terminals in the line, respectively.
Furthermore, we apply the fourth-order accuracy interpolation formula to do spatial discretization of the telegraph equation [20]:
θ 1 λ j ( y ) y | j = k 1 + θ 2 λ j ( y ) y | j = k + θ 1 λ j ( y ) y | j = k + 1 = λ k λ k + 1 Δ y
where λ denotes ν or i ; θ 1 and θ 2 are the constant coefficients decided by the interpolation formula; and θ 1 = 1 / 24 , θ 2 = 1 2 θ 2 .
Substituting Equation (15) into the telegraph equation (Equation (14)), we obtained the spatial discretization form of the telegraph equation as follows:
{ v k v k + 1 Δ y = θ 1 ( l 0 i ˙ k 1 + r 0 i k 1 ) + θ 2 ( l 0 i ˙ k + r 0 i k ) + θ 1 ( l 0 i ˙ k + 1 + r 0 i k + 1 ) , k ( 1 , N 2 ) i k i k + 1 Δ y = θ 1 ( c 0 v ˙ k 1 + g 0 v k 1 ) + θ 2 ( c 0 v ˙ k + g 0 v k ) + θ 1 ( c 0 v ˙ k + 1 + g 0 v k + 1 ) , k ( 1 , N 1 )
For the head-end and the tail-end of the line, we used the second-order accuracy interpolation formula [20] to do spatial discretization as follows:
{ v 0 v 1 Δ y = θ 3 ( l 0 i ˙ 1 + r 0 i 1 ) + θ 1 ( l 0 i ˙ 2 + r 0 i 2 ) i 0 i 1 Δ y = θ 4 ( c 0 v ˙ 0 + g 0 v 0 ) + θ 1 ( c 0 v ˙ 1 + g 0 v 1 ) ,
and
{ v N 1 v N Δ y = θ 3 ( l 0 i ˙ N + r 0 i N ) + θ 1 ( l 0 i ˙ N 1 + r 0 i N 1 ) i N i N + 1 Δ y = θ 4 ( c 0 v ˙ N + g 0 v N ) + θ 1 ( c 0 v ˙ N 1 + g 0 v N 1 )
where θ 3 = 1 θ 1 and θ 4 = 1 / 2 θ 1 .
Considering the boundary conditions of both ends of the transmission line, we obtained the following equations by combining Equations (16)–(18):
K x ˙ + ρ x + ω ( t ) = 0
where K , ρ R ( 2 N + 1 ) × ( 2 N + 1 ) is the constant coefficient matrix; ω ( t ) is the input excitation source for the transmission line; and where:
x = [ ν 0 ν 1 ν N i 1 i N ] T
ω ( t ) = [ i 0 , 0 , , 0 , i N + 1 , 0 , , 0 ] T
K = [ a 4 a 1 0 a 1 a 2 a 2 a 1 a 1 a 4 b 3 b 1 b 1 b 2 b 1 0 b 1 b 3 ]
ρ = [ β 4 β 1 1 β 1 β 2 1 1 β 2 β 1 1 1 β 1 β 4 1 1 1 φ 3 φ 1 1 1 φ 1 φ 2 φ 1 1 1 φ 1 φ 3 ]
where
a j = c 0 Δ y θ j ,   b j = l 0 Δ y θ j ,   β j = g 0 Δ y θ j ,   φ j = r 0 Δ y θ j ,   ( j = 1 , 2 , 3 , 4 )
The four-step L2MF method and the implicit trapezoidal integration method were used to solve Equation (19), and the integration step of the simulation test was h = 1 μ s for both L2MF and the implicit trapezoidal integration method. Moreover, the transmission line was uniformly divided into N = 50 sections. The simulation results are shown in Figure 8, and Figure 9 is an enlargement of Figure 8.
As shown in Figure 8 and Figure 9, the calculation results of the four-step L2MF method were basically consistent with the second-order the implicit trapezoidal integration method. In the moment of the circuit closing, the implicit trapezoidal integration method had a continuous numerical oscillation due to the lack of L-stability. On the contrary, the four-step L2MF method was L-stable, and therefore had strong damping property, which effectively inhibited the numerical oscillations.

4.3. Very Fast Transient Overvoltage Calculation

In this case, we simulate a situation when the disconnector of 550 kV gas insulated switchgear (GIS) acts on the no-load bus. During this process, very fast transient overvoltage (VFTO) is caused by slow contact action and poor arc-extinguishing performance of switch [21,22]. Figure 10 is a simplified equivalent circuit diagram for this process. In Figure 10, e ( t ) is the AC power supply; L t and C t are used to simulate the lumped inductance and capacitance parameters of the transformer and other equipment. DS is the isolation switch, and its mathematical model can be established by a nonlinear time-varying resistance. C r is the lumped ground capacitance with initial voltage to simulate the residual charge; L1 and L2 are import and export short bus, respectively. All parameters of this example in Figure 10 were taken from Reference [23].
The mathematical model of the VFTO simulation was established by using the telegraph equation for the short bus L 1 and L 2 . The short bus, L1, was uniformly divided into N1 = 20 sections, and N2 = 7 of the short bus L2 was selected. The following ordinary differential equations were obtained by the interpolation formulae [20]:
x ˙ = A x + δ ( t )
As the simulation time-step of VFTO was nanosecond scaled, it could be considered that A was the constant matrix in a time-step and δ ( t ) was the incentive source of VFTO.
The four-step L2MF and the implicit trapezoidal integration method were used to solve Equation (16). The integration step of the implicit trapezoidal integration method was h = 1 ns and a smaller integration step of L2MF was h = 0.1 ns . The simulation results in this case are shown in Figure 11.
In this simulation test, the simulation time-step of VFTO was nanosecond scaled, which accurately simulated the change rule of the open-circuit voltage u ( t ) at the end of the line. Neither of these methods produced numerical oscillations, and the waveforms of VFTO for both the methods were almost the same.

4.4. RC Series Circuit Simulation Calculation

This example is the series of a RC circuit, as shown in Figure 12. Figure 12a is a basic linear RC series test circuit. Figure 12b is the excitation voltage source waveform. The exact solution of the capacitor voltage in the circuit was given as follows:
u C ( t ) = 1 Δ t { [ t τ ( 1 e t / τ ) ] [ t Δ t τ ( 1 e ( t Δ t ) / τ ) ] u ( t Δ t ) }
where Δ t = 0.1 ms; a time constant, τ = R C , τ is 0.1 ms; and u ( t Δ t ) is a delayed unit-step signal.
In this paper, the four-step L2MF and the backward Euler method were used to solve the voltages at both the ends of the capacitor. The calculated time-step of the two algorithms was both h = 0.1 ms. The absolute error curve of the calculation results solved with the use of the two algorithms is shown in Figure 13. As shown in Figure 13, the calculation error of the four-step L2MF was smaller than that of the backward Euler method. This is because, although both of them were first-order, the absolute value of the local phase error of the four-step L2MF method was smaller than that of the backward Euler method. This example also indirectly verifies that the calculation effect of the ECDA combined with the four-step L2MF and the implicit trapezoidal integration method was better than that of the CDA.

5. Conclusions

Based on the traditional LMF, a class of low-order LMF was derived by the classic differential quadrature method. The L2MF method had L-stability, and the absolute value of the basic coefficient of the truncation error of four-step L2MF was smaller than that of the backward Euler method. Inspired by the CDA method, we proposed ECDA for electromagnetic transient simulation, which used the implicit trapezoidal integration method to perform integration operation under normal circumstances and switched the integration algorithm to the four-step L2MF method only when the network changed. The simulation results of the simple L2MF method and the implicit trapezoidal integration method showed that the L2MF method could avoid the numerical oscillation completely, and the ECDA had a better simulation effect than the CDA. The next step will apply the ECDA in EMTP-type programs.

Author Contributions

Y.W. proposed the methods and derived the equations; Y.W. and Q.L. performed the numerical simulations; L.Z., J.Y., X.L. and H.R. analyzed the simulation results and given the conclusions; L.Z., J.Y. and B.C. discussed the numerical results and provided useful suggestions on the overall organization of the paper; Y.W. and Q.L. wrote the paper.

Acknowledgments

The authors are grateful to the reviewers and editors for their useful comments and suggestions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ji, Y.; Xu, J.; Xu, Y.; Zhao, C. Design of damping circuit for depressing numerical oscillation in interpolation-free simulation environment. Electr. Power Autom. Equip. 2016, 36, 143–149. [Google Scholar]
  2. Chi, Y.; Shu, D.; Zhang, C.; Wei, L.; Jiang, Q. A parameterized rational-fraction fitting algorithm for numerical oscillation suppression in weak AC grid. Power Syst. Technol. 2016, 40, 1792–1796. [Google Scholar]
  3. Yang, M.; Wang, F. 2-stage 3-order diagonally implicit Runge-Kutta method for electromagnetic transient calculation. Power Syst. Protect. Control 2017, 45, 68–73. [Google Scholar]
  4. Noda, T.; Takenaka, K.; Inoue, T. Numerical integration by the 2-stage diagonally implicit Runge-Kutta method for electromagnetic transient simulations. IEEE Trans. Power Deliv. 2009, 24, 390–399. [Google Scholar] [CrossRef]
  5. Shang, Y.; Yu, Y.M.; Zou, Z.Y.; Li, F.J.; Chen, B. An advanced method of non-linear circuit eliminating numerical oscillation in electromagnetic transient simulation. Power Syst. Protect. Control 2011, 39, 142–146. [Google Scholar]
  6. Lin, J.; Marti, J.R. Implementation of the CDA procedure in the EMTP. IEEE Trans. Power Syst. 1990, 5, 394–402. [Google Scholar]
  7. Marti, J.R.; Lin, J. Suppression of numerical oscillation in the EMTP. IEEE Trans. Power Syst. 1989, 4, 739–747. [Google Scholar] [CrossRef]
  8. Watson, N.; Arrillaga, J. Power Systems Electromagnetic Transients Simulation; IET Press: London, UK, 2003; pp. 221–230. [Google Scholar]
  9. Wang, C.; Li, P.; Wang, L. Progresses on algorithm of electromagnetic transient simulation for electric power system. Autom. Electr. Power Syst. 2009, 33, 97–103. [Google Scholar]
  10. Ji, F.; Wei, X.; Wu, X.; Liu, J.; Xu, W. State space method to analyze the electromagnetic transient of linear switching circuit. Proc. CSEE 2016, 36, 6028–6037. [Google Scholar]
  11. Gear, C.W. Numerical Initial Value Problems in Ordinary Differential Equations; Prentice Hall: Upper Saddle River, NJ, USA, 1971. [Google Scholar]
  12. Bellman, R.; Casti, J. Differential quadrature and long-term integration. J. Math. Anal. Appl. 1971, 34, 235–238. [Google Scholar] [CrossRef]
  13. Wang, F.; Liao, X.; Xie, X. Stability Analysis and Order Improvement for Time Domain Differential Quadrature Method. Adv. Appl. Math. Mech. 2016, 8, 128–144. [Google Scholar]
  14. Wang, F.; Wang, Y. A Coupled Pseudospectral-Differential Quadrature Method for a Class of Hyperbolic Telegraph Equations. Math. Probl. Eng. 2017, 2017, 9013826. [Google Scholar] [CrossRef]
  15. Quan, J.R.; Chang, C.T. New insights in solving distributed system equations by the quadrature methods, Part I: Analysis. Comput. Chem. Eng. 1989, 13, 779–788. [Google Scholar] [CrossRef]
  16. Wang, F.; Liao, X.; Xie, X. Characteristics of the Differential Quadrature Method and Its Improvement. Math. Probl. Eng. 2015, 2015, 657494. [Google Scholar] [CrossRef]
  17. Brugnano, L.; Trigiante, D. Boundary value methods: The third way between linear multistep and Runge-Kutta methods. Comput. Math. Appl. 1998, 36, 269–284. [Google Scholar] [CrossRef]
  18. Liu, D.; Yang, D. A linear 3-step formula of 3-order with a large absolutely stable region. J. Panzhihua Univ. 2008, 6, 66–69. [Google Scholar]
  19. Yuan, Z.; Fei, J.; Liu, D. Numerical Solution of the Initial Value Problem of Stiff Ordinary Differential Equations; Science Press: Beijing, China, 1987. [Google Scholar]
  20. Cangellaris, A.C.; Pasha, S.; Prince, J.L.; Celik, M. A new discrete transmission line model for passive model order reduction and macromodeling of high-speed interconnections. IEEE Trans. Adv. Packag. 1999, 22, 356–364. [Google Scholar] [CrossRef]
  21. Duan, S.F.; Zhao, L.; Li, Z.B.; Zhan, H.M.; Li, C.R. Experimental and simulation study on statistical characteristics of VFTO and VFTC in GIS. Power Syst. Technol. 2015, 39, 3634–3640. [Google Scholar]
  22. Du, W.; Li, Y.G.; Liu, J.M. A new VFTO whole process simulation method based on time domain finite element. Proc. CSEE 2014, 34, 3004–3012. [Google Scholar]
  23. Wang, F.; Yang, M. Fast electromagnetic transient simulation for over-voltages of transmission line by high order Radau method and V-transformation. IET Gener. Transm. Distrib. 2016, 10, 3639–3645. [Google Scholar] [CrossRef]
Figure 1. Stability region of new linear multistep formulae (LMF) for even numbers of s.
Figure 1. Stability region of new linear multistep formulae (LMF) for even numbers of s.
Applsci 08 00883 g001
Figure 2. Stability region of new LMF for odd numbers of s.
Figure 2. Stability region of new LMF for odd numbers of s.
Applsci 08 00883 g002
Figure 3. Flow chart of electromagnetic transient simulation in power systems on a basis of extended critical damping adjustment method (ECDA).
Figure 3. Flow chart of electromagnetic transient simulation in power systems on a basis of extended critical damping adjustment method (ECDA).
Applsci 08 00883 g003
Figure 4. Half-wave rectifier circuit. D represents the diode, which is an ideal switching device with voltage described as d(t). u R L ( t ) denotes a load voltage.
Figure 4. Half-wave rectifier circuit. D represents the diode, which is an ideal switching device with voltage described as d(t). u R L ( t ) denotes a load voltage.
Applsci 08 00883 g004
Figure 5. Load voltage calculated by the three-step low-order linear multistep formulae (L2MF) method and the implicit trapezoidal integration method(TR).
Figure 5. Load voltage calculated by the three-step low-order linear multistep formulae (L2MF) method and the implicit trapezoidal integration method(TR).
Applsci 08 00883 g005
Figure 6. Schematic diagram of a single-phase high-voltage overhead transmission line without load being switched-in suddenly.
Figure 6. Schematic diagram of a single-phase high-voltage overhead transmission line without load being switched-in suddenly.
Applsci 08 00883 g006
Figure 7. A piecewise model of a transmission line.
Figure 7. A piecewise model of a transmission line.
Applsci 08 00883 g007
Figure 8. Terminal voltage of the transmission line calculated by the four-step L2MF method and the implicit trapezoidal integration method.
Figure 8. Terminal voltage of the transmission line calculated by the four-step L2MF method and the implicit trapezoidal integration method.
Applsci 08 00883 g008
Figure 9. An enlargement of Figure 8.
Figure 9. An enlargement of Figure 8.
Applsci 08 00883 g009
Figure 10. A simplified GIS (gas insulated switchgear) for very fast transient overvoltage (VFTO) simulation.
Figure 10. A simplified GIS (gas insulated switchgear) for very fast transient overvoltage (VFTO) simulation.
Applsci 08 00883 g010
Figure 11. VFTO calculated by the four-step L2MF method and the TR method.
Figure 11. VFTO calculated by the four-step L2MF method and the TR method.
Applsci 08 00883 g011
Figure 12. RC (resistance and capacitor) series test circuit and its excitation source. (a) RC series test circuit; and (b) excitation source waveform in the RC test circuit.
Figure 12. RC (resistance and capacitor) series test circuit and its excitation source. (a) RC series test circuit; and (b) excitation source waveform in the RC test circuit.
Applsci 08 00883 g012
Figure 13. Computation errors of the L2MF and backward Euler algorithms as a function of time.
Figure 13. Computation errors of the L2MF and backward Euler algorithms as a function of time.
Applsci 08 00883 g013

Share and Cite

MDPI and ACS Style

Li, Q.; Wang, Y.; Rao, H.; Zhang, L.; Ye, J.; Lei, X.; Chen, B. Extended Critical Damping Adjustment Method for Electromagnetic Transients Simulation in Power Systems. Appl. Sci. 2018, 8, 883. https://doi.org/10.3390/app8060883

AMA Style

Li Q, Wang Y, Rao H, Zhang L, Ye J, Lei X, Chen B. Extended Critical Damping Adjustment Method for Electromagnetic Transients Simulation in Power Systems. Applied Sciences. 2018; 8(6):883. https://doi.org/10.3390/app8060883

Chicago/Turabian Style

Li, Qiang, Yong Wang, Huaxing Rao, Lei Zhang, Jing Ye, Xiaolin Lei, and Bingwen Chen. 2018. "Extended Critical Damping Adjustment Method for Electromagnetic Transients Simulation in Power Systems" Applied Sciences 8, no. 6: 883. https://doi.org/10.3390/app8060883

APA Style

Li, Q., Wang, Y., Rao, H., Zhang, L., Ye, J., Lei, X., & Chen, B. (2018). Extended Critical Damping Adjustment Method for Electromagnetic Transients Simulation in Power Systems. Applied Sciences, 8(6), 883. https://doi.org/10.3390/app8060883

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop