Interpolation for power electronic circuit simulation revisited with matrix exponential and dense outputs

Abstract With a high penetration of power electronic equipment, transient simulation for power electronic circuit has been a main challenge for performance improvement of the electromagnetic transient simulation tools. In this paper, two new solvers for the matrix exponential-based simulation method are proposed based on the dense output formula and Pade approximation of different orders. The proposed solvers are implemented with the optimal combination of numerical integration method and interpolation method. Both solvers are L-stable. One has 3rd-order accuracy which is more accurate than existing simulation tools in power electronic simulation. The other solver achieves 1st-order accuracy with a lower precision numerical integration method than the same-order algorithms and is appealing from computation speed perspective. Numerical studies including TCR circuit and two types of VSC-HVDC systems are conducted to demonstrate the effectiveness of the proposed solvers.

[1]  C. Brezinski,et al.  A Taste of Padé Approximation , 1995, Acta Numerica.

[2]  Jean Mahseredjian,et al.  A Comparison of Numerical Integration Methods and Discontinuity Treatment for EMT Simulations , 2018, 2018 Power Systems Computation Conference (PSCC).

[3]  Hermann W. Dommel,et al.  Digital Computer Solution of Electromagnetic Transients in Single-and Multiphase Networks , 1969 .

[4]  Akihiro Ametani Numerical analysis of power system transients and dynamics , 2015 .

[5]  Awad H. Al-Mohy,et al.  Computing the Action of the Matrix Exponential, with an Application to Exponential Integrators , 2011, SIAM J. Sci. Comput..

[6]  Johan Driesen,et al.  On the Numerical Accuracy of Electromagnetic Transient Simulation With Power Electronics , 2018, IEEE Transactions on Power Delivery.

[7]  Mohamed Abdelrahman,et al.  Practical frequency-domain characterization of switching converters under PWM via Floquet theory , 2019, Electric Power Systems Research.

[8]  Jean Mahseredjian,et al.  Safe Operation of DFIG-Based Wind Parks in Series-Compensated Systems , 2018, IEEE Transactions on Power Delivery.

[9]  N. Higham The Scaling and Squaring Method for the Matrix Exponential Revisited , 2005, SIAM J. Matrix Anal. Appl..

[10]  G. Irwin,et al.  The implementation and effectiveness of linear interpolation within digital simulation , 1997 .

[11]  Jean Mahseredjian,et al.  On a new approach for the simulation of transients in power systems , 2007 .

[12]  Jianzhong Wu,et al.  Accurate dense output formula for exponential integrators using the scaling and squaring method , 2015, Appl. Math. Lett..

[13]  Wilsun Xu,et al.  Algorithms for the accounting of multiple switching events in digital simulation of power-electronic systems , 2005, IEEE Transactions on Power Delivery.

[14]  Taku Noda,et al.  Numerical integration by the 2-stage diagonally implicit Runge-Kutta method for electromagnetic transient simulations , 2009, 2009 IEEE Power & Energy Society General Meeting.

[15]  Peng Li,et al.  Exponential integration algorithm for large-scale wind farm simulation with Krylov subspace acceleration , 2019, Applied Energy.

[16]  C. Dufour,et al.  Two Variable Time-Step Algorithms for Simulation of Transients , 2019, 2019 IEEE Milan PowerTech.

[17]  Jianzhong Wu,et al.  Multiscale Simulation of Power System Transients Based on the Matrix Exponential Function , 2017, IEEE Transactions on Power Systems.

[18]  Johan Driesen,et al.  Accurate second-order interpolation for power electronic circuit simulation , 2017, 2017 IEEE 18th Workshop on Control and Modeling for Power Electronics (COMPEL).

[19]  Jos Arrillaga,et al.  Power Systems Electromagnetic Transients Simulation , 2002 .

[20]  Jean Mahseredjian,et al.  A combined state-space nodal method for the simulation of power system transients , 2011, 2011 IEEE Power and Energy Society General Meeting.