Magnus series expansion method for solving nonhomogeneous stiff systems of ordinary differential equations
Keywords:
Linear differential equations, Magnus Series Expansion Method, Stiff Systems, Geometric Integration, Lie group method.Abstract
In this paper, Magnus Series Expansion Method which is based on Lie Groups and Lie algebras is proposed with different orders to solve nonhomogeneous stiff systems of ordinary differential equations. Using multivariate Gaussian quadrature, fourth (MG4) and sixth (MG6) order method are presented. Then, it is applied to nonhomogeneous stiff systems using different step sizes and stiffness ratios. In addition, approximate and exact solutions are demonstrated with figures in detail. Moreover, absolute errors are illustrated with detailed tables.
References
Aminikhah, H. & Hemmatnezhad, M. 2011. An effective modification of the homotopy perturbation method for stiff systems of ordinary differential equations. Applied Mathematics Letters, 24:1502-1508.
Blanes, S. & Ponsoda, E. 2012. Time-averaging and exponential integrators for non-homogeneous linear IVPs and BVPs. Applied Numerical Mathematics, 62: 875-894.
Brockett, R. W. 1976. Volterra series and geometric control theory. Automatica 12:167-176.
Bui, T.D. & Bui, T.R. 1979. Numerical methods for extremely stiff systems of ordinary differential equations. Appl. Math. Modell. 3:355–358.
Butcher, J.C. 2003. Numerical Methods for Ordinary Differential Equations. Wiley, New York.
Chen, K.T. 1957. Integration of paths, geometric invariants, and a generalized Baker-Hausdorff formula. Annals Math. 67:164-178.
Curtiss, C. F. & Hirschfelder, J. O. 1952. Integration of Stiff Equations. Proc. Nat. Acad. Sci. 38:235-243.
Dahlquist, G. 1963. A Special Stability Problem for Linear Multi-Step. BIT Numer. Math. 3:27-43.
Fatunla, S.V. 1978. An implicit two-point numerical integration formula for linear and nonlinear stiff systems of ordinary differential equations. Mathematics of Computation, 32:1-11.
Flaherty, J.E. & OMalley, R.E. 1997. The numerical solution of boundary value problems for stiff differential equations. Math. Comput. 31:66–93.
Gear, C.W. 1971. Automatic Integration of Ordinary Differential Equations. Communications of the ACM 14:176-179.
Hojjati, G., Rahimi Ardabili, M.Y. & Hosseini, S.M. 2004. A-EBDF: an adaptive method for numerical solution of stiff systems of ODEs. Mathematics and Computers in Simulation 66:33–41.
Iserles, A., Marthinsen, A. & Norsett, S. P. 1999. On the implementation of the method of Magnus series for linear differential equations. BIT, 39:281-304.
Iserles, A., Munthe-Kaas, H. Z., Norsett, S. P. & Zanna, A. 2000. Lie-group methods. Acta Numer, 9: 215-365.
Iserles, A. & Norsett, S.P. 1997. On the solution of linear differential equations in Lie groups. Technical Report 1997/NA3. Department of Applied Mathematics and Theoretical Physics, University of Cambridge, England.
Iserles, A. & Norsett, S. P. 1999. On the solution of linear differential equations in Lie groups. Phil. Trans. R. Soc. A 357: 983-1019.
Iserles, A., Norsett, S.P. & Rasmussen, A.F. 1998. Time-symmetry and high-order Magnus methods. Tech. Rep. 1998/NA06, DAMTP, Cambridge, England.
Ixaru, L.Gr., Vanden Berghe, G. & De Meyer, H. 2000. Frequency evaluation in exponential fitting multistep algorithms for ODEs. J. Comput. Appl. Math. 140:423–434.
Lee, H.C., Chen, C.K. & Hung, C.I. 2002. A modified group-preserving scheme for solving the initial value problems of stiff ordinary differential equations. Applied Mathematics and Computation 133:445-459.
Magnus, W. 1954. On the exponential solution of differential equations for a linear operator. Comm. Pure and Appl. Math. 7:639-673.
Wilcox, R. 1972. Bounds for approximate solutions to the operator differential equation Y ̇(t)=M(t)Y(t); applications to Magnus expansions an to u ̈+[1+f(t) ]u=0. J. Math. Analyt. Applic. 39:92-111.
