GLOBAL ASYMPTOTIC BEHAVIOR OF ITERATIVE IMPLICIT SCHEMES

The global asymptotic nonlinear behavior of some standard iterative procedures in solving nonlinear systems of algebraic equations arising from four implicit linear multistep methods (LMMs) in discretizing three models of 2×2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed using the theory of dynamical systems. The iterative procedures include simple iteration and full and modified Newton iterations. The results are compared with standard Runge-Kutta explicit methods, a noniterative implicit procedure, and the Newton method of solving the steady part of the ODEs. Studies showed that aside from exhibiting spurious asymptotes, all of the four implicit LMMs can change the type and stability of the steady states of the differential equations (DEs). They also exhibit a drastic distortion but less shrinkage of the basin of attraction of the true solution than standard nonLMM explicit methods. The simple iteration procedure exhibits behavior which is similar to standard nonLMM explicit methods except that spurious steady-state numerical solutions cannot occur. The numerical basins of attraction of the noniterative implicit procedure mimic more closely the basins of attraction of the DEs and are more efficient than the three iterative implicit procedures for the four implicit LMMs. Contrary to popular belief, the initial data using the Newton method of solving the steady part of the DEs may not have to be close to the exact steady state for convergence. These results can be used as an explanation for possible causes and cures of slow convergence and nonconvergence of steady-state numerical solutions when using an implicit LMM time-dependent approach in computational fluid dynamics.

[1]  Arieh Iserles,et al.  A unified approach to spurious solutions introduced by time discretization. Part I: basic theory , 1991 .

[2]  H. Flashner,et al.  A computational approach for studying domains of attraction for non-linear systems , 1988 .

[3]  H. C. Yee,et al.  DYNAMICAL APPROACH STUDY OF SPURIOUS STEADY-STATE NUMERICAL SOLUTIONS OF NONLINEAR DIFFERENTIAL EQUATIONS II. GLOBAL ASYMPTOTIC BEHAVIOR OF TIME DISCRETIZATIONS ∗ , 1995 .

[4]  Luca Dieci,et al.  Some Stability Aspects of Schemes for the Adaptive Integration of Siff Initial Value Problems , 1991, SIAM J. Sci. Comput..

[5]  P. Holmes,et al.  Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields , 1983, Applied Mathematical Sciences.

[6]  A. Iserles.  Stability and Dynamics of Numerical Methods for Nonlinear Ordinary Differential Equations , 1990 .

[7]  Christian P. Ullrich,et al.  Computer Arithmetic and Self-Validating Numerical Methods , 1990, Notes and reports in mathematics in science and engineering.

[8]  H. C. Yee,et al.  On the Dynamics of Some Grid Adaption Schemes , 1994 .

[9]  H. C. Yee,et al.  On spurious asymptotic numerical solutions of explicit Runge-Kutta methods , 1992 .

[10]  J. Hale,et al.  Dynamics and Bifurcations , 1991 .

[11]  E. Adams.  Periodic Solutions: Enclosure, Verification, and Applications , 1990, Computer Arithmetic and Self-Validating Numerical Methods.

[12]  R. S. Guttalu,et al.  A Method of Analyzing Generalized Cell Mappings , 1982 .

[13]  Andrew M. Stuart.  Linear Instability Implies Spurious Periodic Solutions , 1989 .

[14]  H. C. Yee,et al.  Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. I. The dynamics of time discretization and its implications for algorithm development in computational fluid dynamics☆ , 1991 .