復(fù)習(xí)有限差分法基礎(chǔ)市公開課特等獎市賽課微課一等獎?wù)n件

鑄造過程計算機模擬講義

復(fù)習(xí)-有限差分法

冒國兵

1第1頁主要內(nèi)容1、差分原理及迫近誤差2、差分方程，截斷誤差和相容性3、收斂性與穩(wěn)定性4、Lax等價定理2第2頁第一節(jié)差分原理及迫近誤差/差分原理向前差分

（1-2）向后差分（1-3）中心差分（1-4）1．差分原理3第3頁第一節(jié)差分原理及迫近誤差/差分原理對一階差分再作一階差分，所得到稱為二階差分，記為。以向前差分為例，有

（1-5）4第4頁第一節(jié)差分原理及迫近誤差/差分原理

依這類推，任何階差分都可由其低一階差分再作一階差分得到。比如n階前差分為

（1-6）5第5頁第一節(jié)差分原理及迫近誤差/差分原理函數(shù)差分與自變量差分之比，即為差商。

一階向前差商為

一階向后差商為

（1-7）（1-8）6第6頁第一節(jié)差分原理及迫近誤差/差分原理一階中心差商為或（1-9）（1-10）7第7頁第一節(jié)差分原理及迫近誤差/差分原理二階差商多取中心式，即

當(dāng)然，在一些情況下也可取向前或向后二階差商。

（1-11）8第8頁第一節(jié)差分原理及迫近誤差/差分原理多元函數(shù)f(x,y,…)差分與差商也能夠類推。如一階向前差商為（1-12）（1-13）9第9頁第一節(jié)差分原理及迫近誤差/迫近誤差差商與導(dǎo)數(shù)之間誤差表明差商迫近導(dǎo)數(shù)程度，稱為迫近誤差。由函數(shù)Taylor展開，能夠得到迫近誤差相對于自變量差分（增量）量級，稱為用差商代替導(dǎo)數(shù)精度，簡稱為差商精度。（1-15）2．迫近誤差10第10頁第一節(jié)差分原理及迫近誤差/迫近誤差一階向后差商也含有一階精度。（1-16）同理：11第11頁第一節(jié)差分原理及迫近誤差/迫近誤差將與Taylor展開式相減可得

可見一階中心差商含有二階精度。（1-17）12第12頁第一節(jié)差分原理及迫近誤差/迫近誤差這說明二階中心差商精度也為二階。

（1-18）將與Taylor展開式相加可得13第13頁第二節(jié)差分方程、截斷誤差和相容性/差分方程假如將微分方程中導(dǎo)數(shù)用對應(yīng)差商近似代替，就可得到有限形式差分方程?，F(xiàn)以對流方程為例，列出對應(yīng)差分方程。（2-1）14第14頁圖2-1差分網(wǎng)格第二節(jié)差分方程、截斷誤差和相容性/差分方程15第15頁若時間導(dǎo)數(shù)用一階向前差商近似代替，即

空間導(dǎo)數(shù)用一階中心差商近似代替，即則在點對流方程就可近似地寫作（2-2）（2-3）（2-4）第二節(jié)差分方程、截斷誤差和相容性/差分方程16第16頁第二節(jié)差分方程、截斷誤差和相容性/截斷誤差按照前面關(guān)于迫近誤差分析知道，用時間向前差商代替時間導(dǎo)數(shù)時誤差為,用空間中心差商代替空間導(dǎo)數(shù)時誤差為因而對流方程與對應(yīng)差分方程間也存在一個誤差：（2-5）17第17頁第二節(jié)差分方程、截斷誤差和相容性/截斷誤差一個與時間相關(guān)物理問題，還須給定初始條件，從而形成完整初值問題。對流方程初值問題為：對應(yīng)差分格式：（2-7）（2-8）18第18頁第二節(jié)差分方程、截斷誤差和相容性/截斷誤差FTCS格式（2-9）FTFS格式（2-10）19第19頁第二節(jié)差分方程、截斷誤差和相容性/截斷誤差（2-11）FTBS格式20第20頁第二節(jié)差分方程、截斷誤差和相容性/截斷誤差(a)FTCS(b)FTFS(c)FTBS圖2-2差分格式21第21頁第二節(jié)差分方程、截斷誤差和相容性/截斷誤差FTCS格式截斷誤差為

FTFS和FTBS格式截斷誤差為（2-12）（2-13）3種格式對都有一階精度。22第22頁第二節(jié)差分方程、截斷誤差和相容性/相容性假如當(dāng)時，差分方程截斷誤差某種范數(shù)也趨近于零，即則表明從截斷誤差角度來看，此差分方程是能用來迫近微分方程，通常稱這么差分方程和對應(yīng)微分方程相容。（2-17）23第23頁

第二節(jié)差分方程、截斷誤差和相容性/相容性只有方程相容，定解條件也相容，即和整個問題才相容。

（2-21）24第24頁第三節(jié)收斂性與穩(wěn)定性/收斂性當(dāng)自變量步長趨于零時，要求差分格式（準(zhǔn)確）解趨于微分方程定解問題（準(zhǔn)確）解。我們稱這種是否趨于微分方程定解問題解情況為差分格式收斂性。，也是微分問題定解區(qū)域上一固定點，設(shè)差分格式在此點解為,對應(yīng)微分問題解為，二者之差為（3-1）更明確地說，對差分網(wǎng)格上任意結(jié)點稱為離散化誤差。25第25頁第三節(jié)收斂性與穩(wěn)定性/收斂性假如當(dāng)、時，離散化誤差某種范數(shù)趨近于零，即則說明此差分格式是收斂，即此差分格式解收斂于對應(yīng)微分問題解。（3-2）26第26頁第三節(jié)收斂性與穩(wěn)定性/穩(wěn)定性首先介紹一下差分格式依賴區(qū)間、決定區(qū)域和影響區(qū)域。還是以初值問題

為例。先看FTCS格式，如圖3-1（a）,p點解依賴于初值線AB段上全部結(jié)點初值，故稱AB段上全部結(jié)點為p點依賴區(qū)間。

（3-17）(a)FTCS(b)FTFS(c)FTBS

圖3-1差分格式依賴區(qū)間27第27頁第三節(jié)收斂性與穩(wěn)定性/穩(wěn)定性(a)FTCS(b)FTFS(c)FTBS

圖3-1差分格式依賴區(qū)間又，三角形pAB區(qū)域內(nèi)任一結(jié)點依賴區(qū)間都包含在AB之內(nèi)，即該區(qū)域內(nèi)任一結(jié)點上解都由AB段上一些結(jié)點初值所決定，而與AB以外結(jié)點初值無關(guān)

，故稱此三角形區(qū)域為AB區(qū)間所決定區(qū)域。28第28頁第三節(jié)收斂性與穩(wěn)定性/穩(wěn)定性全部受p點函數(shù)值影響結(jié)點總和為p點影響區(qū)域，如圖3-2中陰影所表示區(qū)域。FTCS格式

(b)FTFS格式

(c)FTBS格式圖3-2 差分格式影響區(qū)域29第29頁第三節(jié)收斂性與穩(wěn)定性/穩(wěn)定性比如微分問題

其解為零，即若用FTBS格式計算，且計算中不產(chǎn)生任何誤差，則結(jié)果也是零，即（3-18）（3-19）30第30頁第三節(jié)收斂性與穩(wěn)定性/穩(wěn)定性假設(shè)在第k層上第j點，因為計算誤差得到不妨設(shè)k=0,j=0,

即相當(dāng)于FTBS格式寫成（3-20）31第31頁第三節(jié)收斂性與穩(wěn)定性/穩(wěn)定性現(xiàn)分別取和2，列表計算以下：（1）400001161438141163000018383818020000141214001000012120000000010000i-4-3-2-101234n32第32頁第三節(jié)收斂性與穩(wěn)定性/穩(wěn)定性（2）n40000000013000000010200000010010000010000000010000i-4-3-2-10123433第33頁第三節(jié)收斂性與穩(wěn)定性/穩(wěn)定性（3）n400001-824-321630000-16-1280200001-440010000-120000000010000I-4-3-2-10123434第34頁第三節(jié)收斂性與穩(wěn)定性/穩(wěn)定性這個例子首先顯示了該格式影響區(qū)域，另首先還顯示了當(dāng)值不一樣時，計算誤差所產(chǎn)生影響在數(shù)值上有很大不一樣。當(dāng)