發(fā)布時間：2017-12-05 19:03

  本文關(guān)鍵詞：求解波動方程的高精度緊致顯式差分格式

  更多相關(guān)文章： 波動方程 Pade逼近 緊致格式 顯式差分 穩(wěn)定性

【摘要】：本文主要針對波動方程初邊值問題的有限差分法進行研究。首先,針對一維波動方程,利用泰勒級數(shù)展開公式與原方程推導(dǎo)建立了第一時間層上未知函數(shù)值的差分格式；之后,利用Pade逼近離散二階空間導(dǎo)數(shù)項在內(nèi)部節(jié)點上的導(dǎo)數(shù)值,時間方向采用中心差分公式進行離散,得到了一種時間二階、空間四階精度的緊致顯式差分格式,其截斷誤差為O(τ2+廳4);由于上述格式的時間精度與空間精度不相匹配,本文利用截斷誤差余項修正的方法,對時間離散進行了改進,改進后格式的截斷誤差為O(τ4+τ2h2+h4),即格式具有整體四階精度;然后,通過Fourier分析法分析了兩種格式的穩(wěn)定性,前者的穩(wěn)定性條件為|α|λ≤(2/3),后者的穩(wěn)定性條件為|α|λ∈[0,1]∪[2,3]。由于本文格式屬于顯式差分格式,只需進行一次追趕法求解和一次顯式遞推計算；最后通過數(shù)值實驗并與文獻中的格式進行對比,驗證了本文格式的精確性和穩(wěn)定性。接下來,將上述一維波動方程的兩種格式進行推廣,得到了兩種求解二維波動方程的高精度緊致顯式差分格式,一種格式的截斷誤差仍為O(τ2+h4),另一種格式的截斷誤差為O(τ4+τ2h2+h4)。利用Fourier分析法分析了兩種格式的穩(wěn)定性,前者的穩(wěn)定性條件為|α|λ≤1/3,后者的穩(wěn)定性條件為|α|λ∈[0,√2/2](?)[1,√6/2]。由于本文格式屬于顯式差分格式,只需進行兩次追趕法求解和一次顯式遞推計算,無需迭代。最后,通過數(shù)值實驗,將本文格式的計算結(jié)果與文獻中格式的計算結(jié)果進行了對比,驗證了本文格式的精確性和穩(wěn)定性。最后,將二維波動方程中具有四階精度的格式推廣到三維,推廣后的格式仍具有整體四階精度。通過Fourier分析法得出了該格式的穩(wěn)定性條件為|α|λ∈[0,√3/3](?)[√6/3,1]。由于本文格式屬于顯式差分格式,只需進行三次追趕法求解和一次顯式遞推計算,無需迭代。數(shù)值實驗結(jié)果驗證了本文格式的精確性和穩(wěn)定性。
【學(xué)位授予單位】：寧夏大學(xué)
【學(xué)位級別】：碩士
【學(xué)位授予年份】：2016
【分類號】：O241.82

【參考文獻】

中國期刊全文數(shù)據(jù)庫 前1條

1 田敏;羊丹平;;波動方程的重疊型區(qū)域分解并行有限差分算法[J];山東大學(xué)學(xué)報(理學(xué)版);2007年02期

，

本文編號：1255835