Next Article in Journal
Strong Convergence of a System of Generalized Mixed Equilibrium Problem, Split Variational Inclusion Problem and Fixed Point Problem in Banach Spaces
Previous Article in Journal
An Integrated Model for Demand Forecasting and Train Stop Planning for High-Speed Rail
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fast Calculation Method for Transient Response of Transmission Line Based on Chebyshev Pseudospectral–Two-Step Three-Order Boundary Value Coupled Method

1
China Three Gorges University College of Electrical Engineering & New Energy, Yichang 443002, China
2
Hubei Provincial Collaborative Innovation Center for New Energy Microgrid, China Three Gorges University, Yichang 443002, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(5), 721; https://doi.org/10.3390/sym11050721
Submission received: 25 April 2019 / Revised: 19 May 2019 / Accepted: 23 May 2019 / Published: 27 May 2019

Abstract

:
A Chebyshev pseudospectral–two-step three-order boundary value coupled method is proposed and presented for handling the issue associated with complicated calculation, low precision, and poor stability in the process of transient response of transmission line. The first order differential equation in time domain is obtained via dispersing the telegraph equation in space domain by utilizing the pseudospectral method (PSM) based on Chebyshev polynomial. Then the two-step three-order boundary value method (BVM3) is presented and employed to resolve the obtained differential equation, so the numerical solution of the space discrete points can be obtained. Furthermore, the Chebyshev pseudospectral–two-step three-order boundary value coupled method (PSM-BVM3) is presented and compared with the Chebyshev pseudospectral–two-step two order boundary value coupled method (PSM-BVM2), the pseudospectral–differential quadrature method (PSM-DQM), and the pseudospectral method–trapezoid rule (PSM-TR) to validate the feasibility of the new proposed method. Theoretical analysis and numerical simulation reveal that the proposed Chebyshev PSM-BVM3 has a higher performance than the conventional method. For the proposed Chebyshev PSM-BVM3, the higher precision, efficiency, and numerical stability can be obtained and achieved only with fewer discrete points in the space domain, which is suitable for solving the transient response of transmission line. The proposed PSM-BVM3 can improve the drawback of numerical instability of the PSM and can also improve the disadvantage of the BVM as it is not easy to change the latter’s timestep size.

1. Introduction

As the main carrier of information data of transmission line response, the research of transmission line response on excitation signal is significantly crucial in the electromagnetic transient analysis of transmission line. It is investigated that the transmission frequency of signal in high-speed and large-scale lumped circuit will gradually extend to high frequency and even ultra-high frequency [1,2]. Therefore, the issues of delay, distortion, and crosstalk caused by excitation signal passing through transmission line is becoming increasingly serious, which will directly influence the accuracy of transmission line response. As a result, in order to fully and accurately reflect the transient electromagnetic effect of the excitation signal passing through the transmission line, it is of great necessity to analyze the transient response process of transmission line.
Scholars have proposed many numerical methods to resolve the issues in transient response of transmission line, which mainly include the time domain analysis method and frequency domain analysis method. These methods mainly focused on solving the measurement-efficiency and precision problems of non-uniform multi-conductor transmission lines. Among them, Kokkinos et al. have proposed a finite difference time domain method (FDTD) that is simple and easy to implement and employ in practical applications but the calculation timestep of this method is constrained by the stability condition, which resulted in low computation efficiency [3,4]. Jiang et al. have proposed a precise integration method (PIM) that employed in the numerical calculation of telegraph equations in transient response process of transmission line [5,6]. PIM has the performance of good numerical stability and high precision as the calculation of numerical matrix exponent increases with the increase of system dimension. Feng et al. have proposed an improved FDTD (IFDTD) method that is has the higher calculation accuracy than the previous FDTD [7]. By using a time domain method that is similar to Maxwell’s equations, the modeling problem of transmission lines is resolved, and the accuracy of time domain solution is improved remarkably. However, as the calculation time step is quite small and the overall computation is accordingly increasing, it is hard to employ in practical applications. Xu et al. have applied the differential quadrature method (DQM) to transient response of transmission lines and the accuracy can be guaranteed due to its global approximation but with a complex calculation [8]. Hagci et al. have proposed a fast Fourier transform (FFT) method for calculating transient response of transmission line [9]. By using a high-order time-based function to approximate the transmission line parameters, it can preferably overcome electromagnetic interference and have good compatibility while bein susceptible to numerical dispersion and performing with deficiency. Thielen et al. have employed the convolution technique (CT) which shows that the instantaneous response of transmission line is calculated quickly, does not need to solve a large number of matrix equations, and has certain advantages in computation efficiency [10]. However, the stability is limited and numerical oscillation occurs easily. In general, the existing methods cannot take into full account stability, computation accuracy, and efficiency. Therefore, it is necessary to research the coupled algorithm with good comprehensive performance to analyze the transient response of transmission line.
Javidi et al. have proposed the pseudospectral method (PSM), which is also called spectral collocation method [11,12,13]. It is a numerical method widely used to resolve partial differential equation. The principle is to define a set of orthogonal functions as the basis function in the computation region, and then the basis function is used to approximate the resolved variable as its spectral approximation. Because PSM requires fewer discrete nodes to obtain higher approximation precision, Arribas et al. have successfully applied it in many fields. However, there are still some shortcomings for the PSM as the singularity of conventional PSM at the boundary will lead to numerical instability during the practical calculation [14,15,16]. Besides, most of the problems in engineering applications are non-periodic and the solution area is irregular, which will influence the calculation accuracy.
Truhlar et al. have proposed the boundary value method (BVM), Dahlquist et al. have studied the two-step boundary method, and Noor et al. have investigated the convergence and stability of the three and four step boundary value method. The results have shown BVM can control global error and has a high accuracy and great advantages in stability [17,18,19,20,21,22]. However, BVM still has some disadvantages such as its heavy and deficient calculation and it cannot change the step size. A synthesis of the transient response calculation approach proposed in the literature, depending on the type of analysis, is presented in Table 1.
By coupling the Chebyshev pseudospectral method and the boundary value method, a Chebyshev pseudospectral–two-step three-order boundary value coupled method is proposed and presented to resolve the transient response process of transmission line. Firstly, the first order differential equation in time domain is obtained by using the Chebyshev pseudospectral method to discrete telegraph equation in space domain, and then the equation is resolved by using the two-step three-order boundary value method, so the time domain solution of each discrete node can be obtained. Compared with the same level of pseudospectral–differential quadrature method and the Chebyshev pseudospectral–two-step two order boundary value coupled method, the simulation result shows that the coupled method proposed in this paper can simulate the transient response process of transmission line more accurately, and has high precision and efficiency. Moreover, it is unconditionally stable in the time domain, is no longer constrained by Courant conditions, and the spectral precision is convergent in space domain.

2. Chebyshev Pseudospectral Method

The p-order Chebyshev differential matrix D ˜ N ( p ) is defined as:
G = D ˜ N ( p ) U
where: G = [ g ( x ˜ 0 ) , , g ( x ˜ N ) ] T , D ˜ N ( p ) = ( d ˜ i j ( p ) ) ( N + 1 ) × ( N + 1 ) , U = [ u ( x ˜ 0 ) , , u ( x ˜ N ) ] T , x ˜ i = cos ( i π / N ) , i , j ( 0 , N ) .
Babolian et al. have defined the formula of any element of the one-order Chebyshev differential matrix when p = 1 , as shown in Formula (2) [23,24,25,26,27]. And the two-order Chebyshev differential matrix can be directly obtained as shown in Formula (3).
{ d ˜ 00 ( 1 ) = d ˜ N N ( 1 ) = ( 2 N 2 + 1 ) / 6 d ˜ i j ( 1 ) = γ i 2 γ j ( 1 ) i + j sin ( ( i + j ) π 2 N ) sin ( ( i j ) π 2 N ) , i j d ˜ i j ( 1 ) = 1 2 cos ( i π N ) ( 1 + cot 2 ( i π N ) ) , i = j , i 0 , N γ k = { 2 , k = 0 , N 1 , k ( 1 , N 1 )
D ˜ N ( 2 ) = D ˜ N ( 1 ) D ˜ N ( 1 )
As can be known that the Chebyshev point x i in the general interval [ a , b ] can be obtained by simple linear transformation as shown in Formula (4). Where, i ( 0 , N ) .
x i = a + b a 2 ( x ˜ i + 1 )
Accordingly, the Chebyshev differential matrix D N ( p ) corresponding to Chebyshev point x i in any interval [ a , b ] can be given as:
D N ( p ) = ( 2 b a ) p D ˜ N ( p )
For simplicity, the above method is called Chebyshev pseudospectral method in the following text.

3. Two-Step Three-Order Boundary Value Method

Consider the initial value problem of the following first-order ordinary differential equation:
{ y ˙ ( t ) = f ( t , y ( t ) ) , t ( t 0 , t f ) y ( t 0 ) = y 0 ,
The two-step boundary value method is used to solve the first-order ordinary differential equation as shown in Formula (6). In addition, the calculation form of the two-step two-order boundary value method is shown in Formula (7).
y m + 1 ( 1 + β ) y m + β y m 1 = h 2 [ ( 3 β ) f m ( 1 + β ) f m 1 ] , m ( 1 , M 1 )
In the same way, Dahlquist et al. [19] gave the calculation equation of the two-step three-order boundary value method and is shown in Formula (8).
y m + 1 ( 1 + θ ) y m + θ y m 1 = h 12 [ ( 5 + θ ) f m + 1 + 8 ( 1 θ ) f m ( 1 + 5 θ ) f m 1 ] , m ( 1 , M 1 )
In Formulas (7) and (8), h is the length of the time-integral step, and M is the number of time intervals. Where, h t m + 1 t m = ( t f t 0 ) / M , f m + i = f ( t m + i , y m + i ) , i = 1 , 0 , 1 , y m + i = y ( t m + i ) , t m + i = t m + i × h .
In Formulas (7) and (8), if we attempt to resolve the unknown points y 1 , y 2 , , y M 1 , the two boundary values y 0 , y M must be obtained. For the initial value problem (6), only y 0 is known. Therefore, the additional equation should be introduced to solve the unknown point y M .
y M = y M 1 + h 2 ( f M + f M 1 )
Combining Formulas (7) or (8) and (9) to resolve the initial value question (6), we can get the value of each time discrete point y k , k ( 1 , M ) .

4. Calculation Method for Transient Response of Transmission Line Based on Chebyshev Pseudospectral–Two-Step Three-Order Boundary Value Coupled Method

Topology Model of a single uniform transmission line is shown in Figure 1. In Figure 1, R 0 , L 0 , G 0 , and C 0 are the resistance, inductance, conductance, and capacitance per unit length of the single uniform transmission line, respectively.
The famous telegraph equations known to describe voltage and current variations on a uniform transmission line can be written as:
2 ω ( z , t ) t 2 + k 0 ω ( z , t ) t = a 0 2 2 ω ( z , t ) z 2 + b 0 ω ( z , t )
where, a 0 = 1 L 0 C 0 , b 0 = R 0 G 0 L 0 C 0 , k 0 = G 0 L 1 + C 0 R 1 C 0 L 1 .
The new variable ω ( z , t ) represents the voltage or current of the transmission line. The definition domain is ( z , t ) [ c , d ] × [ 0 , T ] , and the initial conditions are:
ω ( z , 0 ) = υ 1 ( z ) , z [ c , d ]
ω t ( z , 0 ) = υ 2 ( z ) , z [ c , d ]
The Dirichlet boundary conditions are:
ω ( c , t ) = h 1 ( t ) , t [ 0 , T ]
ω ( d , t ) = h 2 ( t ) , t [ 0 , T ]
By multiplying e k 0 t / 2 in both sides of Equation (10), the equation can be simplified as:
2 ϖ ( z , t ) t 2 = a 0 2 2 ϖ ( z , t ) z 2 + ( b 0 + 1 4 k 0 2 ) ϖ ( z , t )
In Formula (15), ϖ ( z , t ) = e k 0 t / 2 ω ( z , t ) , the definition domain of the new variable is not changing, and the initial conditions are:
ϖ ( z , 0 ) = υ 1 ( z ) , z [ c , d ]
ϖ t ( z , 0 ) = υ 2 ( z ) + 1 2 k 0 υ 1 ( z ) , z [ c , d ]
The Dirichlet boundary condition is changed to:
ϖ ( c , t ) = e k 0 t / 2 h 1 ( t ) , t [ 0 , T ]
ϖ ( d , t ) = e k 0 t / 2 h 2 ( t ) , t [ 0 , T ]
By using Chebyshev pseudospectral method, the ϖ z z in Formula (15) can be dispersed as formula as:
ϖ z z ( z i , t ) = j = 0 N d i j ( 2 ) ϖ ( z j , t ) , i ( 0 , N )
where, d i j ( 2 ) = ( 2 d c ) 2 d ˜ i j ( 2 ) , z i = c + d c 2 ( x ˜ i + 1 ) . Substituting Formula (20) into Formula (15), we can get:
d 2 ϖ d t 2 = [ a 0 2 D N ( 2 ) + ( b 0 + 1 4 k 0 2 ) I N + 1 ] ϖ
where, ϖ = [ ϖ ( z i , t ) ] T , i ( 0 , N ) . I N + 1 is a N + 1 dimensional identity matrix. Define the following vector ζ , and let ϑ = d ϖ d t :
ζ = [ ϑ ϖ ]
Then the following initial value problem can be obtained:
{ ζ ˙ = Δ ζ , t ( 0 , T ] ζ ( t = 0 ) = ζ 0
in Formula (23),
Δ = ( 0 P Γ 0 )
in Formula (24), P = a 0 2 D N ( 2 ) + ( b 0 + 1 4 k 0 2 ) I N + 1 , and Γ is the N + 1 dimensional identity matrix.
The first-order ordinary differential initial value model for transient response calculation of transmission line is obtained. By using the two-step three-order boundary value method as the basic method, and using the implicit trapezoid formula as the last point method, the global discrete solution of Formula (23) in the time domain can be carried out, and the following equations can be obtained:
J 0 ζ = Ζ
in Equation (25), J 0 is a constant coefficient matrix, q = 2 ( N + 1 ) , i ( 1 , M ) . The form of each variable is:
ζ [ ξ 1 T , ξ 2 T , , ξ M T ] T , ξ i R q × 1
Ζ [ z 1 T , z 2 T , , z M T ] T , z i R q × 1
J 0 = [ J 1 J 5 0 J 3 J 1 J 5 J 3 J 1 J 5 0 J 4 J 2 ]
in Formula (28):
J 1 = [ ( 1 + θ ) I q + 2 ( 1 θ ) 3 h Δ ]
J 2 = ( I q 1 2 h Δ )
J 3 = [ θ I q + ( 1 + 5 θ ) 12 h Δ ]
J 4 = ( I q + 1 2 h Δ )
J 5 = [ I q ( 5 + θ ) 12 h Δ ]
where, I q is a q-order identity matrix.
As can be seen from the above derivation, if Gauss elimination is directly used to solve Equation (25), the trig decomposition of a q × M dimensional matrix will be involved. In this way, the calculated amount is about 2 q 3 M 3 / 3 times multiplication and division without considering the sparsity of J 0 . However, because J 0 is a large dimensional block tridiagonal matrix, the catch-up method for a block tridiagonal matrix equation is applied to solve this equation [28]. In this case, the calculated amount can be reduced to 2 M q 3 / 3 + 2 M q 2 . As a contrast, the s-step s-order differential quadrature method based on network uniformity and fast numerical algorithm based on V-conversion are used to solve the initial value problem (23). Danesh et al. have presented the calculation process and calculation amount of this method in detail, which is not repeated in this paper [29]. Through comparison, the advantages of the proposed algorithm in calculation precision and efficiency can be further illustrated. The final point method is the same, and only J 0 should be adjusted accordingly.

5. Simulation Test

For different values of parameters in Equation (10), two examples are given to validate the improved performance of the proposed algorithm. The Chebyshev pseudospectral–two-step two-order boundary value coupled method (PSM-BVM2), the Chebyshev pseudospectral–two-step three-order boundary value coupled method (PSM-BVM3), and the Chebyshev pseudospectral–two-step differential quadrature method (PSM-DQM2) are adopted to solve the initial value problem (23), to validate the performance of the proposed method. All three methods are approximated dispersing by the Chebyshev pseudospectral method in space domain and adopting the same step size h in time domain. The simulation software is Matlab7.14, and the hardware platform is CPU Intel core i7-8700k 3.70 GHz.

5.1. Test for Example 1

In Equation (10), a 0 = 2 , b 0 = 3 π 2 , k 0 = 8 π , c = 0 , d = 1 . The initial value and Dirichlet boundary condition are determined by analytic solution of ω ( z , t ) = e π t sin ( π z ) . Here the PSM-BVM2 ( β = 6 ) , PSM-BVM3 ( θ = 3 ) , and PSM-DQM2 are used to resolve the initial value problem (23). The spatial discrete points of the Chebyshev pseudospectral method are N . The absolute error (AE) and relative error (RE) were compared.
In order to test the computational accuracy of the proposed algorithm, different simulation times are selected using the analytic solution as the benchmark. The numerical errors of the above algorithms are compared and analyzed. As shown in Table 2 and Table 3, the results of the PSM-BVM2, PSM-BVM3, and PSM-DQM2 are given.
From Table 2 and Table 3, it can be seen that the measurement error of the PSM-BVM3 and the PSM-BVM2 are less than that of the PSM-DQM2, which can prove that the proposed PSM-BVM algorithm exhibits improved performance. Moreover, the measurement error of the PSM-BVM3 is three times less than that of the PSM-DQM2 and the PSM-BVM2. In conclusion, according to the result, the PSM-BVM3 has higher accuracy than the other two algorithms of the PSM-BVM2 and the PSM-DQM2.
To validate the computation efficiency of the proposed algorithm, the CPU consumption time of the above three methods was tested under different spatial discrete points, as shown in Table 4. In Table 4, T = 2   s , h = 0.1   ms .
From Table 4, it can be seen that the average CPU consumption time is 0.9745 s for the PSM-BVM2, 1.560 s for the PSM-BVM3, and 5.564 s for the PSM-DQM2, which indicated that the PSM-BVM3 and PSM-BVM2 have higher computation efficiency. Moreover, the runtime of the PSM-BVM3 is larger than the PSM-BVM2. But considering the measurement accuracy and the runtime cost, the PSM-BVM3 exhibits the most improved performance than the PSM-DQM2 and the PSM-BVM2, which can prove that the higher runtime required by PSM-BVM3 (see Table 4) is justified by an improvement in the quality of the solutions.

5.2. Test for Example 2

In Equation (10), a 0 = 1 / π , b 0 = 4 , k 0 = 4 , c = 0 , d = 1 . In the same way, the initial value and Dirichlet boundary condition are determined by analytic solution of ω ( z , t ) = e t cosh ( π z ) . The number of spatial discrete points is N = 14 . h = 0.001   s , T = 15   s . Here the PSM-BVM2 ( β = 6 ) and PSM-BVM3 ( θ = 3 ) are used to solve, by using the analytic solution of (10) as the benchmark. The absolute error of the two methods e r r ( z , t ) = | ω ˜ ( z , t ) ω ( z , t ) | are tracked, where ω ˜ is the numerical solution. The error curves are shown in Figure 2 and Figure 3. The numerical results of PSM-BVM3 ω ¯ ( z , t ) and analytical solutions of example 2 are shown in Figure 4 and Figure 5. Only part of sample points are intercepted in Figure 2, Figure 3, Figure 4 and Figure 5.
From Figure 2 and Figure 3 we can see that both the PSM-BVM2 and PSM-BVM3 can well track the transient response process of analog transmission lines, and the calculation errors are generally decreasing with the change of time, so the algorithm has good convergence. In the time domain, the PSM-BVM3 has better calculation accuracy than the PSM-BVM2.
From Figure 4 and Figure 5, the numerical results of the PSM-BVM3 are almost the same as the real analytical solutions, which fully shows that the algorithm is unconditionally stable in the time domain and the spectral precision is convergent in the space domain.
In conclusion, the example 1 shows that the PSM-BVM2 and PSM-BVM3 are more accurate and efficient than the PSM-DQM2, and example 2 shows that the numerical stability of the proposed algorithm can be guaranteed after long time integration process. Therefore, compared with the PSM-DQM2, PSM-BVM2, and PSM-BVM3 have higher precision and efficiency, and good numerical stability, so they can simulate the transient response process of transmission line for a long time.

5.3. Calculation of Transient Response of Nondestructive Transmission Line

In order to validate the feasibility and effectiveness of the proposed coupled method in practical cases, the example of the transient response of transmission line is carried out, in which the comparison between the PSM-BV3, PSM-BVM2 and the pseudospectral method–trapezoid rule (PSM-TR) is presented.
The simulation circuit diagram of nondestructive transmission line is shown in Figure 6. And the parameters of the transmission are set as:
The excitation source is a step voltage source e i = 100 ε ( t )   kV where ε ( t ) is a unit step function. The internal resistance of source is R S = 100   Ω . The transmission line length is l = 400   km , and the distribution parameters of transmission line are L 0 = 2.5 × 10 7 H / m , C 0 = 1.0 × 10 10 F / m . The maximum transmission speed in the transmission line is v = 2 × 10 8 m / s . The load parameters are: resistance R L = 200   Ω , inductance L L = 300   mH , capacitance C L = 3000   pF . The number of spatial discrete points of the PSM is N = 7 . The time integral step length of the PSM-BVM3, PSM-BVM2 and PSM-TR are h 1 = 8.0 × 10 8 s , h 2 = 8.0 × 10 8 s and h 3 = 1.0 × 10 8 s , respectively. The calculation results of end terminal load voltage of nondestructive transmission line are shown in Figure 7, in which (a) is the result of the PSM-BVM3, (b) is the result of the PSM-BVM2 and (c) is the result of the PSM-TR.
It can be obtained from Figure 7 that the calculation results of steady state voltage of the three numerical methods are consistent, and the value is about 130 kV. The transient load voltage calculated by the PSM-BVM3 did not generate numerical oscillation due to the abrupt change of excitation signal, while for the PSM-BVM2 a little numerical oscillation appears in the calculation result and for the PSM-TR a large number of high frequency components are included in the calculation result. The reason is the time integral step of the PSM-BVM3 is large, and error accumulation can be well avoided, so numerical oscillation can be effectively avoided. Moreover, the time step of the PSM-BVM3 is larger than that of the PSM-TR, but the PSM-BVM3 is better than the PSM-TR in transient response of transmission line, which proves the improved effectiveness of the proposed coupled method.
The simulation results show that the proposed new method exhibits improved performance of faster decay rate to reduce the truncation error, which can effectively suppress the numerical oscillation phenomenon and reduce the transient calculation error of the transmission line caused by the oscillation. Besides, the proposed coupled method can compensate the drawback of the pseudospectral method which results in numerical instability caused by the singularity at the boundary and can also improve the calculation accuracy. Moreover, it can make improvements regarding the disadvantage of the boundary value method in which it is not easy to change the step size and the computational efficiency of which is not high. In conclusion, the method has the advantages of high computational precision, high computational efficiency, and numerical stability and is easy to expand in the practical application of transmission line problems in complex systems.

6. Conclusions

A Chebyshev pseudospectral–two-step three-order boundary value coupled method to resolve the transient response of transmission line is proposed in this study. In the given Chebyshev grid points, the telegraph equation is dispersed in space domain by the Chebyshev pseudospectral method, and the first-order linear differential equation is obtained. Moreover, such equation can be obtained by using the two-step three-order boundary value method to discrete the differential equation in the time domain. To avoid dimension disaster, the catch-up method of block tridiagonal is used to solve this equation. The numerical examples and simulation example prove that:
(1) The numerical results show that the coupled method in this paper is more accurate, efficient, and stable than the conventional PSM-DQM and PSM-BVM2 in the time domain. Moreover, the proposed coupled method is unconditionally stable in the time domain, and the spectral precision is convergent in the space domain. Furthermore, through the simulation of the lossless transmission line, it is proved that the proposed coupled method can better simulate the transient response of transmission line and can be applied in the practical application.
(2) The proposed coupled method exhibits improved performance of faster decay rate to reduce the truncation error, which can effectively suppress the numerical oscillation phenomenon and reduce the transient calculation error of the transmission line caused by the oscillation.
(3) The proposed coupled method can compensate the drawback of the pseudospectral method which results in numerical instability caused by the singularity at the boundary and can also improve the calculation accuracy. Moreover, it can improve on the disadvantage of the boundary value method, in which it is not easy to change the step size and the computational efficiency of which is not high.
Theoretical analysis and numerical simulation reveal that the proposed Chebyshev PSM-BVM3 has higher computation accuracy, efficiency, and better numerical stability. The coupled method is suitable for solving the transient response of the transmission line in the power system. Besides, applying the couple method, the issue of transmission oscillation and inaccurate signal transmission caused by delay, distortion, and crosstalk on the transmission line can be resolved. Furthermore, by employing the proposed method, the excitation signal in the power system can be accurately transmitted to the load. Therefore, the calculation benefit problem caused by the measurement inaccuracy can be improved, and the economic loss in the power system can be reduced.
Transient electromagnetic response on the transmission line can be coupled to the secondary and protection device, which will cause the secondary and protection device work abnormaly. Therefore, in order to accurately predict the electromagnetic transient process of transmission line, it is significant to analyze the transient response process of the transmission line. Thus, the proposed Chebyshev PSM-BVM3 method for transmission line transient response is significant in electromagnetic compatibility research.

Author Contributions

Z.L. provided the main idea for this paper, designed the research of the proposed method, and wrote the paper. T.H., Y.T., T.Z., and Z.L. conducted the test data collection and designed the experiment.

Funding

This work was supported by the Key Project of Science and Technology Research Plan of Education Department of Hubei (Grant Number D20181204).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Li, Z.; Tao, Y.; Abu-Siada, A.; Masoum, M.A.; Li, Z.; Xu, Y.; Zhao, X. A New Vibration Testing Platform for Electronic Current Transformers. IEEE Trans. Instrum. Meas. 2019, 68, 704–712. [Google Scholar] [CrossRef]
  2. Lee, S. Dynamics of Trapped Solitary Waves for the Forced KdV Equation. Symmetry 2018, 10, 129. [Google Scholar] [CrossRef]
  3. Kokkinos, T.; Sarris, C.D.; Eleftheriades, G.V. Periodic finite-difference time-domain analysis of loaded transmission-linenegative-refractive-index metamaterials. IEEE Trans. Microw. Theory Tech. 2005, 53, 1488–1495. [Google Scholar] [CrossRef]
  4. Zhang, Y. Finite-difference time-domain analysis of integrated ceramic ball grid array package antenna for highly integrated wireless transceivers. IEEE Trans. Antennas Propagat. 2004, 52, 435–442. [Google Scholar] [CrossRef]
  5. Jiang, L.; Chen, Z.; Mao, J. On the numerical stability of the Precise Integration Time-Domain (PITD) method. IEEE Microw. Wirel. Compon. Lett. 2007, 17, 471–473. [Google Scholar] [CrossRef]
  6. Tang, M.; Mao, J.; Gao, C.; Chen, B.; Chen, H. A precise Time-Step Integration method for transient analysis of lossy nonuniform transmission lines. IEEE Trans. Electromagn. Compat. 2008, 50, 166–174. [Google Scholar] [CrossRef]
  7. Feng, X.; Chen, L.; Wei, H.; Wu, K. Fast and accurate transient analysis of buried wires and its applications. IEEE Trans. Electromagn. Compat. 2014, 56, 188–199. [Google Scholar] [CrossRef]
  8. Xu, Q.; Li, Z.; Wang, J. Modeling of transmission lines by the differential quadrature method. IEEE Microw. Guided Wave Lett. 1999, 9, 145–147. [Google Scholar] [CrossRef]
  9. Bagci, H.; Yilmaz, A.E.; Michielssen, E. An FFT-accelerated time domain multiconductor transmission line simulator. IEEE Trans. Electromagn. Compat. 2010, 52, 199–214. [Google Scholar] [CrossRef]
  10. Thielen, B.L.A.; Vandenbosch, G.A.E. Fast transmission line coupling calculation using a convolution technique. IEEE Trans. Electromagn. Compat. 2001, 43, 11–17. [Google Scholar] [CrossRef]
  11. Wang, T.; Guo, B. Composite generalized Laguerre–Legendre pseudospectral method for Fokker-Planck equation in an infinite channel. Appl. Numer. Math. 2008, 58, 1448–1466. [Google Scholar] [CrossRef]
  12. Ghoshi, S. Relative Effects of Asymmetry and Wall Slip on the Stability of Plane Channel Flow. Fluids 2017, 2, 66. [Google Scholar] [CrossRef]
  13. Izadkhah, M.M.; Saberi, J.; Toutounian, F. An extension of the Gegenbauer pseudospectral method for the time fractional Fokker-Planck equation. Math. Method. Appl. Sci. 2018, 41, 1301–1315. [Google Scholar] [CrossRef]
  14. Arribas, D.G.; Rivo, M.S.; Arnedo, M.S. Optimization of path-constrained systems using pseudospectral methods applied to aircraft trajectory planning. IFIC Pap. 2015, 48, 192–197. [Google Scholar] [CrossRef]
  15. Argyros, I.K.; Shakhno, S.; Yarmola, H. Two-Step Solver for Nonlinear Equations. Symmetry 2019, 11, 128. [Google Scholar] [CrossRef]
  16. Zhang, Y.; Tan, Z. On the existence of solutions to the Navier-Stokes-Poisson equations of a two-dimensional compressible flow. Math. Method Appl. Sci. 2011, 30, 305–329. [Google Scholar] [CrossRef]
  17. Truhlar, D.G. Finite difference boundary value method for solving one-dimensional eigenvalue equations. J. Comput. Phys. 2015, 10, 123–132. [Google Scholar] [CrossRef]
  18. Zhu, B.; Zeng, Q.; Chen, Y.; Zhao, Y.; Liu, S. A dual-input high step-up DC/DC converter with ZVT auxiliary circuit. IEEE Trans. Energy Convers. 2019, 34, 161–169. [Google Scholar] [CrossRef]
  19. Chen, Z.; Liu, C.; Simos, T.E. New three–stages symmetric two step method with improved properties for second order initial/boundary value problems. J. Math. Chem. 2018, 56, 2591–2616. [Google Scholar] [CrossRef]
  20. Jagannadha Rao, G.V.V.; Padhan, S.K.; Postolache, M. Application of Fixed Point Results on Rational F-Contraction Mappings to Solve Boundary Value Problems. Symmetry 2019, 11, 70. [Google Scholar] [CrossRef]
  21. Wattanasakulpong, N.; Pornpeerakeat, S.; Chaikittiratana, A. Chebyshev Collocation Solutions for Vibration Analysis of Circular Cylindrical Shells with Arbitrary Boundary Conditions. Int. J. Struct. Stab. Dyn. 2017, 17, 1750020. [Google Scholar] [CrossRef]
  22. Noor, M.A.; Al-Said, E.A. Finite-Difference Method for a System of Third-Order Boundary-Value Problems. J. Optim. Theory Appl. 2002, 112, 627–637. [Google Scholar] [CrossRef]
  23. Babolian, E.; Fattahzadeh, F. Numerical solution of differential equations by using Chebyshev wavelet operational matrix of integration. Appl. Math. Comput. 2007, 188, 417–426. [Google Scholar] [CrossRef]
  24. Qiu, L.; Yu, Y.; Xiong, Q.; Deng, C.; Cao, Q.; Han, X.; Li, L. Analysis of electromagnetic force and deformation behavior in electromagnetic tube expansion with concave coil based on finite element method. IEEE Trans. Appl. Supercond. 2018, 28. [Google Scholar] [CrossRef]
  25. Doha, E.H.; Abd-Elhameed, W.M.; Youssri, Y.H. Second kind Chebyshev operational matrix algorithm for solving differential equations of Lane–Emden type. New Astron. 2013, 24, 113–117. [Google Scholar] [CrossRef]
  26. Xiong, Q.; Tang, H.; Wang, M.; Huang, H.; Qiu, L.; Yu, K.; Chen, Q. Design and implementation of tube bulging by an attractive electromagnetic force. J. Mater. Process. Tech. 2019. [Google Scholar] [CrossRef]
  27. Koksa, M.E. An operator-difference method for telegraph equations arising in transmission lines. Discret. Dyn. Nat. Soc. 2011, 2011, 561015. [Google Scholar] [CrossRef]
  28. Yang, N.; Huang, Y.; Hou, D.; Liu, S.; Ye, D.; Dong, B.; Fan, Y. Adaptive Nonparametric Kernel Density Estimation Approach for Joint Probability Density Function Modeling of Multiple Wind Farms. Energies 2019, 12, 1356. [Google Scholar] [CrossRef]
  29. Danesh, M.; Farajpour, A.; Mohammadi, M. Axial vibration analysis of a tapered nanorod based on nonlocal elasticity theory and differential quadrature method. Mech. Res. Commun. 2012, 39, 23–27. [Google Scholar] [CrossRef]
Figure 1. Topological model of a single uniform transmission line.
Figure 1. Topological model of a single uniform transmission line.
Symmetry 11 00721 g001
Figure 2. Three-dimensional curve of absolute error of the PSM-BVM3.
Figure 2. Three-dimensional curve of absolute error of the PSM-BVM3.
Symmetry 11 00721 g002
Figure 3. Three-dimensional curve of absolute error of the PSM-BVM2.
Figure 3. Three-dimensional curve of absolute error of the PSM-BVM2.
Symmetry 11 00721 g003
Figure 4. Numerical results for the PSM-BVM3.
Figure 4. Numerical results for the PSM-BVM3.
Symmetry 11 00721 g004
Figure 5. Analytic solutions of example 2.
Figure 5. Analytic solutions of example 2.
Symmetry 11 00721 g005
Figure 6. Nondestructive transmission line circuit.
Figure 6. Nondestructive transmission line circuit.
Symmetry 11 00721 g006
Figure 7. The calculation results of end terminal load voltage of nondestructive transmission line. (a) PSM-BVM3; (b) PSM-BVM2; (c) PSM-TR.
Figure 7. The calculation results of end terminal load voltage of nondestructive transmission line. (a) PSM-BVM3; (b) PSM-BVM2; (c) PSM-TR.
Symmetry 11 00721 g007
Table 1. Synthesis of transient response calculation approach. Finite difference time domain method (FDTD); precise integration method (PIM); improved FDTD (IFDTD); differential quadrature method (DQM); fast Fourier transform (FFT); convolution technique (CT); pseudo spectral method (PSM); boundary value method (BVM).
Table 1. Synthesis of transient response calculation approach. Finite difference time domain method (FDTD); precise integration method (PIM); improved FDTD (IFDTD); differential quadrature method (DQM); fast Fourier transform (FFT); convolution technique (CT); pseudo spectral method (PSM); boundary value method (BVM).
Transient Response Calculation ApproachAdvantagesDrawbacks
FDTDSimple and easy to implement and employ in practical applicationCalculation timestep limitation, low computation efficiency
PIMGood numerical stability, high precisionCalculation of numerical matrix exponent increase due to increase of system dimension
IFDTDHigher calculation accuracy than conventional FDTDCalculation burden due to its small calculation timestep
DQMAccuracy can be guaranteed due to its global approximationComplex calculation
FFTCan preferably reduce electromagnetic interference and has good compatibilitySusceptible to numerical dispersion and performed with deficiency
CTGood computation efficiencyStability limitation, numerical oscillation
PSMRequires fewer discrete nodes to obtain higher precisionSingularity at the boundary, non-periodic
BVMHigh accuracy and great stabilityHeavy and deficient calculation
Table 2. The calculation error of the Chebyshev pseudospectral–two-step two-order boundary value coupled method (PSM-BVM2, the Chebyshev pseudospectral–two-step three-order boundary value coupled method (PSM-BVM3), and the Chebyshev pseudospectral–two-step differential quadrature method (PSM-DQM2). ( N = 8 , h = 0.1   ms ) .
Table 2. The calculation error of the Chebyshev pseudospectral–two-step two-order boundary value coupled method (PSM-BVM2, the Chebyshev pseudospectral–two-step three-order boundary value coupled method (PSM-BVM3), and the Chebyshev pseudospectral–two-step differential quadrature method (PSM-DQM2). ( N = 8 , h = 0.1   ms ) .
T (s)PSM-BVM2 ( × 10 8 ) PSM-BVM3 ( × 10 11 ) PSM-DQM2 ( × 10 8 )
AEREAEREAERE
0.54.36708.43413.13426.23795.29169.0227
11.51478.79791.15786.45691.93409.6047
20.12388.85080.10586.57360.16059.6879
Table 3. The calculation error of PSM-BVM3, PSM-BVM2 and PSM-DQM2. ( N = 12 , h = 0.1   ms ) .
Table 3. The calculation error of PSM-BVM3, PSM-BVM2 and PSM-DQM2. ( N = 12 , h = 0.1   ms ) .
T (s)PSM-BVM2 ( × 10 8 ) PSM-BVM3 ( × 10 11 ) PSM-DQM2 ( × 10 8 )
AEREAEREAERE
0.52.81294.13651.25152.93447.41306.2612
11.91385.36620.98213.59073.08207.6619
20.10175.79010.10483.68080.26647.8543
Table 4. Comparison of CPU consumption time between three algorithms.
Table 4. Comparison of CPU consumption time between three algorithms.
N PSM-BVM2PSM-BVM3PSM-DQM2
80.936 s1.172 s5.386 s
121.013 s1.948 s5.742 s

Share and Cite

MDPI and ACS Style

Li, Z.; Hu, T.; Tao, Y.; Zhang, T.; Li, Z. Fast Calculation Method for Transient Response of Transmission Line Based on Chebyshev Pseudospectral–Two-Step Three-Order Boundary Value Coupled Method. Symmetry 2019, 11, 721. https://doi.org/10.3390/sym11050721

AMA Style

Li Z, Hu T, Tao Y, Zhang T, Li Z. Fast Calculation Method for Transient Response of Transmission Line Based on Chebyshev Pseudospectral–Two-Step Three-Order Boundary Value Coupled Method. Symmetry. 2019; 11(5):721. https://doi.org/10.3390/sym11050721

Chicago/Turabian Style

Li, Zhenhua, Tinghe Hu, Yuan Tao, Tao Zhang, and Zhenxing Li. 2019. "Fast Calculation Method for Transient Response of Transmission Line Based on Chebyshev Pseudospectral–Two-Step Three-Order Boundary Value Coupled Method" Symmetry 11, no. 5: 721. https://doi.org/10.3390/sym11050721

APA Style

Li, Z., Hu, T., Tao, Y., Zhang, T., & Li, Z. (2019). Fast Calculation Method for Transient Response of Transmission Line Based on Chebyshev Pseudospectral–Two-Step Three-Order Boundary Value Coupled Method. Symmetry, 11(5), 721. https://doi.org/10.3390/sym11050721

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