應(yīng)用于偏微分方程的增量未知元和類小波增量未知元方法.pdf

1、提出不同類型的增量未知元用于構(gòu)造有限差分數(shù)值格式。本文主要考慮增量未知元以下方面。 首先，通過增量未知元方法建立適合三維偏微分方程的清晰矩陣框架，在三維情形提出增量未知元的多層格式，通過數(shù)值試驗，證實包括求解泊松方程，增量未知元的分層預(yù)處理形式可以用于更一般的方程。 其次，將增量未知元方法和一些現(xiàn)代迭代方法相結(jié)合，如MR，GCR，Orthomin(k)，Bi-CGSTAB，HSS，BTSS等，并求解由三維對流擴散方程多層

2、離散生成的非對稱正定線性系統(tǒng)。通過理論分析，估計出與靜態(tài)對流擴散方程相關(guān)的增量未知元矩陣條件數(shù)，以及MR，GCR，Orthomin(k)迭代法求解數(shù)值解所需的迭代步數(shù)。給出數(shù)值試驗，通過兩層逼近證明增量未知元是一個好的預(yù)處理子，尤其在結(jié)合一些迭代方法時。 第三，在二維情形，使用兩層增量未知元的θ格式求解依賴時間的對流擴散方程。證實關(guān)于二階增量未知元的慣性流形多重網(wǎng)格算法是收斂的。增量未知元θ格式在0≤θ<1/2時是條件穩(wěn)定的，在

3、1/2≤θ≤1時是無條件穩(wěn)定的。在時間層每一步迭代使用GMRES方法求解依賴時間的對流擴散方程，數(shù)值試驗說明增量未知元方法能控制計算中帶來的擾動。 第四，提出一個修正的，基于一邊差分的Crank-Nicolson格式，用于求解二維依賴時間的對流占優(yōu)擴散方程。二階一邊差分逼近用于對流項的離散，二階中心差分逼近用于擴散項的離散。該數(shù)值格式是相容的，無條件穩(wěn)定的。對全離散格式，推導(dǎo)了先驗誤差估計。使用增量未知元預(yù)處理子之后，數(shù)值結(jié)果證

4、實了修正的Crank-Nicolson格式的穩(wěn)定性和有效性。 接著，本文提出一個隱式差分格式求解上述類似問題。對常系數(shù)的對流占優(yōu)擴散方程而言，基于一邊有限差分和中心差分的隱格式是相容的并且無條件穩(wěn)定，對應(yīng)L2誤差估計被推導(dǎo)。而且Burgers類型的方程使用類似半隱格式，其中只是對非線性項線性化，其時間上采用顯式，這種半隱格式是條件穩(wěn)定的。通過增量未知元方法的預(yù)處理技術(shù)，這種隱格式或半隱格式是一種有效格式，它能避免數(shù)值擾動，在穩(wěn)定

5、性和精度上表現(xiàn)和理論結(jié)果一致。 同時，提出一種新的類小波增量未知元，它具有一些好的性質(zhì)。一方面，近似解空間能被類小波增量未知元分解成兩個L2正交子空間，這能使耦合系統(tǒng)自動消掉一些項，簡化計算。另一方面，從類小波增量未知元到節(jié)點未知元的變化矩陣是正交矩陣，更方便計算。當(dāng)應(yīng)用于慣性流形多網(wǎng)格算法和非線性伽略金方法時，可證明類小波增量未知元方法的收斂性。 最后，提出二維和三維類小波增量未知元方法，用于求解二維和三維具有多項式非