交替方向隱式法

數值分析中,交替方向隱式法(Alternating direction implicit method)是有限差分法的一種,用於求解拋物線型偏微分方程橢圓型偏微分方程[1]。特別適用於求解二維及更高維度的熱傳導方程擴散方程

求解熱傳導方程在傳統上使用Crank-Nicolson方法,該方法較為耗時。ADI的優點在於,每一迭代步中,所求解的方程具有更為簡單的結構,因此更易於求解。

方法

考慮二維擴散方程,

 

隱式Crank-Nicolson方法將給出以下有限差分方程:

 

其中, 是關於坐標方向p上的中心差分算符。通過穩定性分析可以證明該方法對於任意 都表現穩定。

但是,Crank-Nicolson方法的缺點在於,上述方程中的帶狀矩陣分佈過寬,這使得求解方程相當耗時。

ADI方法的思想在於將一個有限差分方程分割為兩個,一個在x方向上隱式求導,另一個在y方向上隱式求導。

 
 

這樣,該方程系統涉及一個對稱陣和一個三角矩陣,可以用三對角陣的求解算法進行計算。

可以證明,二維條件下該方法無條件穩定[2]

在此基礎上擴展有更多的ADI方法,如Douglas[3],f-factor方法[4],可用於求解三維及更高維的問題。

參考文獻

  1. ^ Peaceman, D. W.; Rachford Jr., H. H., The numerical solution of parabolic and elliptic differential equations, Journal of the Society for Industrial and Applied Mathematics, 1955, 3: 28–41, MR0071874 .
  2. ^ Douglas, Jr., J., On the numerical integration of uxx+ uyy= utt by implicit methods, Journal of the Society of Industrial and Applied Mathematics, 1955, 3: 42–65, MR0071875 .
  3. ^ Douglas Jr., Jim, Alternating direction methods for three space variables, Numerische Mathematik, 1962, 4: 41–63, ISSN 0029-599X, doi:10.1007/BF01386295 .
  4. ^ Chang, M. J.; Chow, L. C.; Chang, W. S., Improved alternating-direction implicit method for solving transient three-dimensional heat diffusion problems, Numerical Heat Transfer, Part B: Fundamentals, 1991, 19 (1): 69–84, ISSN 1040-7790, doi:10.1080/10407799108944957 .