Improvement of numerical stability for the computation of transients in lines and cables

This paper discusses numerical stability problems of a frequency-dependent transmission-line and cable modeling approach used for electromagnetic transient analysis. Time-domain numerical errors due to the discrete computation of convolution integrals can be estimated in terms of transfer function parameters for a given line or cable model. Based on this estimation, a methodology for the improvement of numerical stability is presented. The numerical advantages of the new method are supported by demonstrations and comparisons with existing models. The method presented in this paper is applicable to power cables and transmission lines.