The rate of convergence to steady state of parabolic ADI solvers is analyzed in terms of the L2-norms of the residuals. The analysis allows one to predict the number of iterations necessary for convergence as function of the Courant number, λ. A simple modification of existing ADI codes is devised. It improves the convergence rate substantially and is insensitive to the Courant number in a large range of λ.
Bibliographical noteFunding Information:
* Research was supported in part by the U.S. Air Force Office of Scientific Research under Grant AFOSR-80-0249 and in part by the NASA Cooperative Agreement NCCl-45. + Research was supported in part by the U.S. Army Research and Standardization Group (Europe) under Contract DAJA 38-80-C-0032 and in part by NASA Contract NAS l-15810 while the author was in residence at ICASE, NASA Langley Research Center, Hampton, Virginia 23665.