A high-order ADI method for parabolic problems with variable coefficients

A high-order compact alternating direction implicit (ADI) method is proposed for solving two-dimensional (2D) parabolic problems with variable coefficients. The computational problem is reduced to sequence one-dimensional problems which makes the computation cost-effective. The method is easily extendable to multi-dimensional problems. Various numerical tests are performed to test its high-order accuracy and efficiency, and to compare it with the standard second-order Peaceman–Rachford ADI method. The method has been applied to obtain the numerical solutions of the lid-driven cavity flow problem governed by the 2D incompressible Navier–Stokes equations using the stream function-vorticity formulation. The solutions obtained agree well with other results in the literature.

[1]  Graham F. Carey,et al.  Extension of high‐order compact schemes to time‐dependent problems , 2001 .

[2]  Murli M. Gupta High accuracy solutions of incompressible Navier-Stokes equations , 1991 .

[3]  W. Spotz High-Order Compact Finite Difference Schemes for Computational Mechanics , 1995 .

[4]  Alexander Ostermann,et al.  Dimension splitting for evolution equations , 2008, Numerische Mathematik.

[5]  Irfan Altas,et al.  A high accuracy defect-correction multigrid method for the steady incompressible Navier-Stokes equations , 1994 .

[6]  Alain Rigal High order difference schemes for unsteady one-dimensional diffusion-convection problems , 1994 .

[7]  Murli M. Gupta,et al.  A single cell high order scheme for the convection‐diffusion equation with variable coefficients , 1984 .

[8]  Jun Zhangy Multigrid Method and Fourth Order Compact Diierence Scheme for 2d Poisson Equation with Unequal Meshsize Discretization , 2001 .

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

[10]  B. J. Noye,et al.  Finite difference methods for solving the two‐dimensional advection–diffusion equation , 1989 .

[11]  G. Carey,et al.  High‐order compact scheme for the steady stream‐function vorticity equations , 1995 .

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

[13]  Jun Zhang Multigrid Method and Fourth-Order Compact Scheme for 2D Poisson Equation with Unequal Mesh-Size Discretization , 2002 .

[14]  J. Geiser,et al.  Higher-Order Difference and Higher-Order Splitting Methods for 2D Parabolic Problems with Mixed Derivatives , 2007 .

[15]  Hai-wei Sun,et al.  High order compact scheme with multigrid local mesh refinement procedure for convection diffusion problems , 2002 .

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

[17]  R. K. Mohanty,et al.  A fourth order difference method for the one‐dimensional general quasilinear parabolic partial differential equation , 1990 .

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

[20]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

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

[22]  Schémas compacts d'ordre élevé: application aux problèmes bidimensionnels de diffusion-convection instationnaire I , 1999 .

[23]  S. Dennis,et al.  Compact h4 finite-difference approximations to operators of Navier-Stokes type , 1989 .

[24]  Samir Karaa An accurate LOD scheme for two-dimensional parabolic problems , 2005, Appl. Math. Comput..

[25]  R. K. Mohanty,et al.  High order difference schemes for the system of two space second order nonlinear hyperbolic equations with variable coefficients , 1996 .

[26]  Robert J. MacKinnon,et al.  Differential‐equation‐based representation of truncation errors for accurate numerical simulation , 1991 .

[27]  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 .

[28]  B. J. Noye,et al.  A third-order semi-implicit finite difference method for solving the one-dimensional convection-diffusion equation , 1988 .

[29]  W. Spotz Accuracy and performance of numerical wall boundary conditions for steady, 2D, incompressible streamfunction vorticity , 1998 .

[30]  E. Erturk,et al.  Fourth‐order compact formulation of Navier–Stokes equations and driven cavity flow at high Reynolds numbers , 2004, ArXiv.