Study of a Numerical Integral Interpolation Method for Electromagnetic Transient Simulations

In the fixed time-step electromagnetic transient (EMT)-type program, an interpolation process is applied to deal with switching events. The interpolation method frequently reduces the algorithm's accuracy when dealing with power electronics. In this study, we use the Butcher tableau to analyze the defects of linear interpolation. Then, based on the theories of Runge-Kutta integration, we propose two three-stage diagonally implicit Runge-Kutta (3S-DIRK) algorithms combined with the trapezoidal rule (TR) and backward Euler (BE), respectively, with TR-3S-DIRK and BE2-3S-DIRK for the interpolation and synchronization processes. The proposed numerical integral interpolation scheme has second-order accuracy and does not produce spurious oscillations due to the size change in the time step. The proposed method is compared with the critical damping adjustment method (CDA) and the trapezoidal method, showing that it does not produce spurious numerical oscillations or first-order errors.