Maximum norm error bounds of ADI and compact ADI methods for solving parabolic equations

Alternating direction implicit (ADI) schemes are computationally efficient and widely utilized for numerical approximation of the multidimensional parabolic equations. By using the discrete energy method, it is shown that the ADI solution is unconditionally convergent with the convergence order of two in the maximum norm. Considering an asymptotic expansion of the difference solution, we obtain a fourth-order, in both time and space, approximation by one Richardson extrapolation. Extension of our technique to the higher-order compact ADI schemes also yields the maximum norm error estimate of the discrete solution. And by one extrapolation, we obtain a sixth order accurate approximation when the time step is proportional to the squares of the spatial size. An numerical example is presented to support our theoretical results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010

[1]  Jianping Zhu,et al.  An efficient high‐order algorithm for solving systems of reaction‐diffusion equations , 2002 .

[2]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[3]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

[4]  Samir Karaa.  A high-order compact ADI method for solving three-dimensional unsteady convection-diffusion problems , 2006 .

[5]  Graeme Fairweather,et al.  Improved forms of the alternating direction methods of Douglas, Peaceman, and Rachford for solving parabolic and elliptic equations , 1964 .

[6]  Noel J. Walkington,et al.  An ADI Method for Hysteretic Reaction-Diffusion Systems , 1997 .

[7]  J. J. Douglas.  Alternating direction methods for three space variables , 1962 .

[8]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947 .

[9]  Jiten C. Kalita,et al.  A class of higher order compact schemes for the unsteady two‐dimensional convection–diffusion equation with variable convection coefficients , 2002 .

[10]  N. N. Yanenko,et al.  The Method of Fractional Steps , 1971 .

[11]  Zhi-zhong Sun,et al.  An unconditionally stable and O(τ2 + h4) order L∞ convergent difference scheme for linear parabolic equations with variable coefficients , 2001 .

[12]  Weizhong Dai,et al.  A second‐order ADI scheme for three‐dimensional parabolic differential equations , 1998 .

[13]  Jun Zhang,et al.  High order ADI method for solving unsteady convection-diffusion problems , 2004 .

[14]  J. Douglas,et al.  A general formulation of alternating direction methods , 1964 .

[15]  B. Fornberg,et al.  A compact fourth‐order finite difference scheme for the steady incompressible Navier‐Stokes equations , 1995 .

[16]  Weizhong Dai,et al.  Compact ADI method for solving parabolic differential equations , 2002 .

[17]  Graeme Fairweather,et al.  A New Computational Procedure for A.D.I. Methods , 1967 .