A-stable and l-stable high-order integration methods for solving stiff differential equations
This paper describes a new A- and L-stable integration method for simulating the time-domain transient response of nonlinear circuits. The proposed method, which is based on the Obreshkov formula, can be made of arbitrary high order while maintaining the A-stability property. The new method allows for the adoption of higher order integration methods for the transient analysis of electronic circuits while enabling them to take larger step sizes without violating the stability, leading to faster simulations. The method can be run in an L-stable mode to handle circuits with extremely stiff equations. Necessary theoretical foundations, implementation details, error-control mechanisms, and computational results are presented.
|Keywords||A-stability, Circuit simulation, High-order integration methods, L-stability, Multiderivative methods, Numerical solution of differential equations (DEs), Stiff circuits|
|Journal||IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems|
Gad, E. (Emad), Nakhla, M.S, Achar, R, & Zhou, Y. (Yinghong). (2009). A-stable and l-stable high-order integration methods for solving stiff differential equations. IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, 28(9), 1359–1372. doi:10.1109/TCAD.2009.2024712