二維熱傳導方程的交替方向隱格式

2022-06-03 23:27:04 字數 586 閱讀 9849

考慮二維熱傳導方程的初邊值問題

(1)(2)(3)

其中,且當時有

取空間步長,時間步長,其中,為整數。記,,,,,,,,。此外記

記為上的網格函式,引人如下記號:

,,,易知

記則有即是等價的。

若令則為上的網格函式。

差分格式的建立

定義上的網格函式

其中在點處考慮微分方程(1),有

4)由泰勒展開式,有

將以上三式代入(4),得

上式又可進一步寫為

5)其中

由初邊值條件(1)和(2)有

6)7)

在(5)~(7)略去小量項,並用代替,可得如下差分格式8)9)

10)它是乙個2層隱式差分格式。差分格式的截斷誤差為。

差分格式的求解

差分格式(8)可以寫成

(11)

其中令則(11)等價於

12)當第層上的值已知時,由(11)求出過渡層變數的值:對任意固定的,取邊值條件

13)求解(12)得到

當已求出時,由求出第層上的值

:對任意固定的,取邊界條件

14)求解得到

熱傳導方程的差分格式

熱傳導方程的差分格式上機實習報告 二零一四年五月 分別用向前差分格式 向後差分格式 六點對稱格式,求解下列問題 在時刻的數值解,並與解析解。1差分格式形式 設空間步長,時間步長,網比.1 向前差分格式 向前差分格式,即 1 其中,以表示網比。1 式可改寫成如下 此格式為顯格式。其矩陣表示式如下 2 ...

ug轉二維cad的方法

方法1 關於ug二維圖與autocad圖紙的轉換最佳方法以下給出乙個較佳方案供大家參考!我平時都是這樣做的 有必要做這種轉換時 發現效果很好!希望能起到拋磚引玉的作用在用ugtodxf 轉換複雜的裝配圖形時,為避免因有關引數設定不當而引起多餘線條和檢視丟失等問題,在轉換過程中,採用cgm資料格式過度...

關於二維陣列和指向指標的指標

以前一直有種誤解 二維陣列的是陣列的陣列,所以陣列的首位址是指向第乙個元素指標,而這個元素又是乙個陣列,所以把陣列首位址理解為指向指標的指標。如int a 3 2 以前一直認為a是乙個指向int指標的指標,即是乙個int 最近發現這是錯的。如果int p a 編譯就會報錯。如果強制轉換int p i...