A high-order Padé ADI method for unsteady convection-diffusion equations

A high-order alternating direction implicit (ADI) method for computations of unsteady convection-diffusion equations is proposed. By using fourth-order Pade schemes for spatial derivatives, the present scheme is fourth-order accurate in space and second-order accurate in time. The solution procedure consists of a number of tridiagonal matrix operations which make the computation cost effective. The method is unconditionally stable, and shows higher accuracy and better phase and amplitude error characteristics than the standard second-order ADI method [D.W. Peaceman, H.H. Rachford Jr., The numerical solution of parabolic and elliptic differential equations, Journal of the Society of Industrial and Applied Mathematics 3 (1959) 28-41] and the fourth-order ADI scheme of Karaa and Zhang [High order ADI method for solving unsteady convection-diffusion problem, Journal of Computational Physics 198 (2004) 1-9].

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

[2]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .

[3]  P. Moin.  Fundamentals of Engineering Numerical Analysis , 2001 .

[4]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[5]  Willem Hundsdorfer,et al.  Stability and convergence of the Peaceman-Rachford ADI method for initial-boundary value problems , 1987 .

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

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

[8]  J. Nordström,et al.  Summation by Parts Operators for Finite Difference Approximations of Second-Derivatives with Variable Coefficients , 2004, Journal of Scientific Computing.

[9]  Haecheon Choi,et al.  Turbulent Drag Reduction: Studies of Feedback Control and Flow Over Riblets , 1993 .

[10]  A fourth order ADI method for semidiscrete parabolic equations , 1983 .

[11]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

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

[13]  Meng Wang,et al.  A Computational Methodology for Large-Eddy Simulation of Tip-Clearance Flows , 2003 .

[14]  R. J. MacKinnon,et al.  Analysis of material interface discontinuities and superconvergent fluxes in finite difference theory , 1988 .

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

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

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

[18]  Miguel R. Visbal,et al.  High-Order-Accurate Methods for Complex Unsteady Subsonic Flows , 1999 .

[19]  D. Gottlieb,et al.  A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy , 1999 .

[20]  A comparative study of ADI splitting methods for parabolic equations in two space dimensions , 1984 .

[21]  P. Moin,et al.  Suitability of upwind-biased finite difference schemes for large-eddy simulation of turbulent flows , 1997 .

[22]  N. S. Barnett,et al.  Private communication , 1969 .