若干發(fā)展方程的譜方法和譜元法.pdf

1、譜方法是求解偏微分方程的重要數值方法。它的主要優(yōu)點是高精度，這使得該方法能夠與有限差分、有限元一起而成為偏微分方程的三大數值方法之一。它的缺點是不能靈活地適應復雜的計算區(qū)域，從而阻礙了它的廣泛應用和發(fā)展。解決的辦法之一就是將問題所在的區(qū)域分解成若干子區(qū)域，在每個子區(qū)域上使用譜方法。這種技巧通常叫做譜元法。 當針對具體的問題建立了譜(元)方法的數值格式以后，格式的誤差分析就非常重要。關于譜方法和譜元法的誤差分析已經有大量的工作，但

2、豐滿的最優(yōu)的收斂階估計并不很多，尤其對于非線性問題，好的結果更少。本文就致力于若干個發(fā)展方程的譜(元)方法，討論其誤差分析，即格式的數值穩(wěn)定性和收斂性(收斂階)。特別關注的是非線性問題以及收斂階的最優(yōu)估計。 本文的主要工作為： 首先對五階KdV方程建立了Legendre-Petrov-Galerkin譜方法，該方法是三階KdV方程建立的Legendre-Petrov．Galerkin譜方法[64]的自然推廣。證明了該格式

3、的數值穩(wěn)定性以及收斂性。結論表明，收斂階是最優(yōu)的。 其次，討論了一維情形的對流-耗散方程的譜元法。證明了該方法的收斂性與穩(wěn)定性，并給出了收斂階的很好的估計。另外，由于譜元法中區(qū)域的分解，算法的并行化顯得更為重要。根據算法的特點，描述了數值實驗中并行化的過程，數值例子也驗證了方法的有效性。 再次，對于二維的線性Schrdinger方程，空間上采用譜元法，時間上采用Crank-Nicolson離散得到的全離散格式，證明了格式

4、在L<'2>以及H<'1>意義下的穩(wěn)定性和收斂性，得到了最優(yōu)的收斂階估計。為了計算上的簡化，利用算子分裂，將交替方向方法應用于格式設計，并證明了交替方向的譜元方法在L<'2>以及H<'1>意義下的穩(wěn)定性和收斂性。同樣得到的收斂階也是最優(yōu)的。 接著，對一類非經典非線性拋物型方程建立了Legendre譜方法的Galerkin離散格式，推導了格式的穩(wěn)定性，得到了收斂階估計．同時也給出了Chebyshev擬譜格式和Chebyshev-L

5、egendre擬譜格式，進行了數值實驗．在這類方程數值求解的研究中，嘗試應用譜方法。理論結果和數值實驗表明，這種應用是很成功的。 最后，對于無界區(qū)域上的Burgers方程，給出了一種穩(wěn)定化的Hermite譜方法。該方法首次直接使用Hermite多項式作為基函數來逼近問題的解，解決了格式的數值穩(wěn)定性問題，并得到了格式最優(yōu)的收斂階估計，數值結果驗證了理論的正確性。 接著討論了擬譜方法以及混合的Hermite譜方法。在整個誤差