Matrix exponential based electromagnetic transients simulation algorithm with Krylov subspace approximation and accurate dense output

The ever-increasing scale of system interconnection and diverse dynamic time constants of power equipment pose challenges to power system electromagnetic transients simulation algorithm. In this paper, an advanced numerical integration method based on matrix exponential operator is proposed. Krylov subspace approximation is utilized to effectively reduce the problem dimension, and an accurate dense output formula is established to enable the method adopting larger simulation step size and revealing higher frequency dynamic details simultaneously. Accuracy, numerical stability and computation efficiency is demonstrated through numerical case studies.