Alternating direction implicit method

In numerical analysis, the Alternating Direction Implicit (ADI) method is a finite difference method for solving parabolic, hyperbolic and elliptic partial differential equations.^[1] It is most notably used to solve the problem of heat conduction or solving the diffusion equation in two or more dimensions. It is an example of an operator splitting method.^[2]

The traditional method for solving the heat conduction equation numerically is the Crank–Nicolson method. This method results in a very complicated set of equations in multiple dimensions, which are costly to solve. The advantage of the ADI method is that the equations that have to be solved in each step have a simpler structure and can be solved efficiently with the tridiagonal matrix algorithm.

The method

Stencil figure for the alternating direction implicit method in finite difference equations

Consider the linear diffusion equation in two dimensions,

${\partial u\over \partial t} = \left({\partial^2 u\over \partial x^2 } + {\partial^2 u\over \partial y^2 } \right) = ( u_{xx} + u_{yy} ) = \Delta u$

The implicit Crank–Nicolson method produces the following finite difference equation:

${u_{ij}^{n+1}-u_{ij}^n\over \Delta t} = {1 \over 2}\left(\delta_x^2+\delta_y^2\right) \left(u_{ij}^{n+1}+u_{ij}^n\right)$

where $\delta_p$ is the central difference operator for the p-coordinate. After performing a stability analysis, it can be shown that this method will be stable for any $\Delta t$.

A disadvantage of the Crank–Nicolson method is that the matrix in the above equation is banded with a band width that is generally quite large. This makes direct solution of the system of linear equations quite costly (although efficient approximate solutions exist, for example use of the conjugate gradient method preconditioned with incomplete Cholesky factorization).

The idea behind the ADI method is to split the finite difference equations into two, one with the x-derivative taken implicitly and the next with the y-derivative taken implicitly,

${u_{ij}^{n+1/2}-u_{ij}^n\over \Delta t/2} = \left(\delta_x^2 u_{ij}^{n+1/2}+\delta_y^2 u_{ij}^{n}\right)$

${u_{ij}^{n+1}-u_{ij}^{n+1/2}\over \Delta t/2} = \left(\delta_x^2 u_{ij}^{n+1/2}+\delta_y^2 u_{ij}^{n+1}\right).$

The system of equations involved is symmetric and tridiagonal (banded with bandwidth 3), and is typically solved using tridiagonal matrix algorithm.

It can be shown that this method is unconditionally stable and second order in time and space.^[3] There are more refined ADI methods such as the methods of Douglas,^[4] or the f-factor method^[5] which can be used for three or more dimensions.

References

1. Peaceman, D. W.; Rachford Jr., H. H. (1955), "The numerical solution of parabolic and elliptic differential equations", Journal of the Society for Industrial and Applied Mathematics 3 (1): 28–41, doi:10.1137/0103003, MR 0071874.
2. • Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). "Section 20.3.3. Operator Splitting Methods Generally". Numerical Recipes: The Art of Scientific Computing (3rd ed.). New York: Cambridge University Press. ISBN 978-0-521-88068-8.
3. Douglas, Jr., J. (1955), "On the numerical integration of u_xx+ u_yy= u_t by implicit methods", Journal of the Society of Industrial and Applied Mathematics 3: 42–65, MR 0071875.
4. Douglas Jr., Jim (1962), "Alternating direction methods for three space variables", Numerische Mathematik 4 (1): 41–63, doi:10.1007/BF01386295, ISSN 0029-599X.
5. Chang, M. J.; Chow, L. C.; Chang, W. S. (1991), "Improved alternating-direction implicit method for solving transient three-dimensional heat diffusion problems", Numerical Heat Transfer, Part B: Fundamentals 19 (1): 69–84, doi:10.1080/10407799108944957, ISSN 1040-7790.