王潔

（臺(tái)州學(xué)院 數(shù)學(xué)與信息工程學(xué)院，浙江 臨海 317000）

二維熱傳導(dǎo)方程的緊ADI法

王潔

（臺(tái)州學(xué)院 數(shù)學(xué)與信息工程學(xué)院，浙江 臨海 317000）

對(duì)二維熱傳導(dǎo)方程進(jìn)行緊交替方向有限差分，該方法在空間方向上具有四階精度，在時(shí)間方向上具有二階精度.證明了當(dāng)rx，ry≥1/6時(shí)該有限差分解收斂于連續(xù)解.數(shù)值例子驗(yàn)證了該有限差分法具有高階精度.

二維熱傳導(dǎo)方程；緊ADI法；高階精度；收斂性

1 引言

熱傳導(dǎo)方程在物理學(xué)、生物學(xué)中有很重要的作用,已受到廣泛關(guān)注．工程技術(shù)中許多重要問題如微波的熱處理、自燃、地下水的傳輸、擴(kuò)散物質(zhì)濃度、電纜的傳輸?shù)葐栴}都可用熱傳導(dǎo)方程來描述.

通常用有限差分法求解依賴于時(shí)間的熱傳導(dǎo)方程，如顯格式、隱格式及C-R格式.對(duì)于高維熱傳導(dǎo)方程可采用加權(quán)差分格式、交替方向顯格式、交替方向隱格式（ADI）及局部一維（LOD）法.對(duì)于二維熱傳導(dǎo)方程數(shù)值解法出現(xiàn)了許多研究成果,但它們主要進(jìn)行的是解的存在性與穩(wěn)定性分析［1］-［4］，很少有收斂性方面的工作.

廖文遠(yuǎn)等人［5］提出了一種緊交替方向隱式有限差分方法，即在時(shí)間方向上用Crank-Nicolson方法近似，在空間方向上用Padé方法逼近．這種方法具有緊的性質(zhì)，并且在時(shí)間方向上具有二階精度，在空間方向上具有四階精度．優(yōu)點(diǎn)是將一個(gè)二維問題降為兩個(gè)一維問題，這將大大降低求解難度．本文利用緊ADI方法對(duì)二維熱傳導(dǎo)方程進(jìn)行有限差分，并把它寫成矩陣的形式,證明了當(dāng)rx，ry≥1/6時(shí)有限差分解收斂于連續(xù)解，并給出數(shù)值例子驗(yàn)證。

2 緊ADI有限差分格式

考慮如下二維熱傳導(dǎo)方程：

其中Ω=（0，1）×（0，1），u（x，y，t）表示溫度，f（x，y，t）為源項(xiàng)，φ和ψ是已知函數(shù)，且具有充分的光滑性，a為導(dǎo)溫系數(shù)或熱擴(kuò)散系數(shù).

將Ω分別在x和y方向上進(jìn)行分割，步長分別為hx和hy．將（0，T]進(jìn)行分割，步長為τ．記My=1/hx，My=

1/hy，Mt=T/τ．離散點(diǎn)記作．為方便起見，我們用數(shù)對(duì)（i，j）表示離散點(diǎn)（xi，yj），記：

定義標(biāo)準(zhǔn)的中心差商算子：

在時(shí)間方向上用Crank-Nicolson格式，在空間方向上用Padé逼近，我們有如下有限差分逼近［5］：

此處，ui，j，n表示在點(diǎn)（xi，yj，tn）處的u逼近．因?yàn)椋?）左邊部分分別只含有三點(diǎn)中心差分算子，所以它將兩維問題（4）降低為兩個(gè)一維問題，故得到一個(gè)ADI算法．即第一步是沿x方向利用Thomas算法逐行求解三對(duì)角方程，得到vi，j，n-1，然后沿y方向逐列求解三對(duì)角方程得到解ui，j，n，由于每步只需求解一組三對(duì)角方程，且兩個(gè)方向交替變換，故通稱這類格式為交替方向隱格式.

定義如下的列向量：

我們也定義如下（Mx-1）階的對(duì)稱三對(duì)角矩陣：

和（My-1）階對(duì)稱三對(duì)角矩陣：

則（5）式可以寫成如下矩陣形式：

這里對(duì)每個(gè)n和j，U0=UM，n=0，Gj，n和Ri，n是與邊界函數(shù)相關(guān)的兩個(gè)向量。

為了將（10）寫成更緊湊的形式，定義如下向量：

向量Un和有相同的分量，向量Un是先排x方向，后排y方向，而向量是先排y方向，后排x方向,它們的唯一的區(qū)別是分量的次序.

記M=（Mx-1）（My-1），引入以下M階塊矩陣：

于是，（10）可寫成

3 收斂性

定理2.1.設(shè)Un為方程（1）的精確解，Wn為方程（13）的數(shù)值解.當(dāng)rx，ry≥1/6時(shí)，Wn收斂于Un，n=0，1，…，Mt．

這就證明了定理結(jié)論．

4 數(shù)值例子

下面通過具體數(shù)值例子來分析前面格式的收斂性.在0≤x，y≤1內(nèi)考察精確解為u（x，y，t）=e-2π2tsin （πx）sin（πy）的初邊值問題，初始條件由精確解給出.表1給出了當(dāng)r=1/2時(shí)，在（1/2，1/2）點(diǎn)處精確解和緊ADI方法數(shù)值解的比較，由此看出，緊ADI方法得到的數(shù)值解誤差很小，數(shù)值解收斂于精確解，收斂階數(shù)是4階.比［4］中方法得到的誤差要小.

表1 r=1/2時(shí)的緊ADI法最大誤差和收斂階Tab.1 Maximum error and convergence rate for compact ADI method when r=1/2

［1］葉其孝，李正元.反應(yīng)擴(kuò)散引論［M］.北京：科學(xué)出版社，1994.

［2］戴嘉尊，邱建賢.微分方程數(shù)值解法［M］.南京：東南大學(xué)出版社，2002.

［3］張文生.科學(xué)計(jì)算中的偏微分方程［M］.北京：高等教育出版社，2006.

［4］葛永斌.高維熱傳導(dǎo)方程的高精度交替方向隱格式［J］.上海理工大學(xué)學(xué)報(bào),2007,29（1）：56-58.

［5］Wenyuan Liao,Jianping Zhu,and Abdul Q.M.Khaliq.An Efficient High-Order Algorithm for Solving System of Reaction-Diffusion Equations［J］.Numerc Methods Partial Differential Eq,2002,18：340-354.

［6］Yuan-Ming Wang，Jie Wang．A higher-order compact ADI method with monotone iterative procedure for systems of reaction-diffusion equations［J］．Computers and Mathematics with Applications，2011，62：2434-2451．

［7］A.Berman，R.Plemmons，Nonnegative Matrix in the Mathematical Science，Academic Press，New York，1979．

Compact ADI Method for Solving Heat Equation with Two Dimensions

WANG Jie

（School of Mathematics and Information Engineering,Taizhou University,Linhai 317000,China）

A high order algorithm is discussed in the paper. This algorithm has the accuracy of four-order in space and two-order in time.The finite difference scheme can be solved by ADI method.The convergence of finite difference solution to the continuous solution is proved when it is between 1/6 and 5/6.A numerical example is shown the high accuracy of the finite difference in the last section.

heat equation with two dimensions;compact ADI method;high order accuracy;convergence

10.13853/j.cnki.issn.1672-3708.2014.03.003

2014-04-17

臺(tái)州學(xué)院青年項(xiàng)目（2013QN09）

王 潔（1982- ），女，浙江仙居人，講師，碩士，主要從事計(jì)算數(shù)學(xué)研究。