淺談求解四階拋物方程的幾種方法

蔣超　周甜　胥康

摘 ? 要：高階拋物型方程問題在工程應用和科學研究中占據著重要地位。本文主要針對四階拋物型方程混合問題，先后給出了Crank-Nicolson隱式差分格式、Saulyev非對稱差分隱格式和一種三層顯式差分格式。簡要闡述了這三種方法的構造過程，并通過給出其穩定性和截斷誤差分析，表明這三種方法在求解四階拋物型方程時是穩定可靠的。

關鍵詞：四階拋物方程 ?隱式差分格式 ?顯示差分格式

中圖分類號：O175.26 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻標識碼：A ? ? ? ? ? ? ? ? ? ? ? ?文章編號：1674-098X（2019）06（b）-0045-02

1 ?問題背景

在科學和工程應用中，拋物型方程[1-2]占據著重要地位，構造可靠的差分格式求解拋物型方程一直是研究重點。雖然目前針對二階拋物方程的研究頗多，然而對高階拋物問題的研究同樣意義重大。本文針對四階拋物型方程[3-5]，給出并分析了3種可靠的差分格式。

2 ?Crank-Nicolson隱格式求解四階拋物方程

考慮如下的拋物型方程

（1）

先取空間步長，時間步長為τ，對網格進行剖分，網格比。記號是函數u（x，t）在點處的數值逼近。邊界條件可處理為：

。用關于t的向前向后差商和關于x的二階中心差商近似代替和，得到兩種古典格式：

（2）

（3）

將上述（2）式和（3）式相加得到Crank-Nicolson隱格式：

由文獻可知其具局部截斷誤差為，且絕對穩定。

3 ?Saulyev非對稱差分隱格式求解四階拋物方程

針對問題（1），結合中值定理可以構造出如下的Saulyev非對稱差分隱格式：

（4）

（5）

（6）

（7）

通過Fouier方法可分析（4）、（5）在時穩定，（6）、（7）絕對穩定，且局部截斷誤差均為。

4 ?三層顯式差分格式[5，6]求解四階拋物方程

為了構造三層顯式格式，需要添加耗散項得到方程。

其中為待定參數。類似于Crank-Nicolson隱格式的構造過程，用關于t的向前向后差商及二階中心差商近似代替和，關于x的二階中心差商近似代替，得到如下格式：。

式中α為權數，和為中心差商。文獻證明了其是無條件穩定的，且當時，局部截斷誤差均為，當α=1/2時，局部截斷誤差均為。

5 ?結語

四階拋物型方程問題在工程應用和科學研究中占據著重要地位，本文針對四階混合問題，先后給出了Crank-Nicolson隱式差分格式、Saulyev非對稱差分隱格式和三層顯式差分格式。根據其構造過程，及穩定性與局部截斷誤差分析，可見這三種方法在求解四階拋物方程時是可靠的。

參考文獻

[1] N.Khan， M.S.Hashmi， S.Iqbal， T.Mahmood. Optimal homotopy asymptotic method for solving Volterra integral equation of first kind[J]. Alexandria Engineering Journal， 2014（53）： 751-755.

[2] CUI Xiao-peng，SHAN Shuang-rong.High-Accurate and Two-Layer Difference Schemes for Solving Four-Order Parabolic Equation[J].Journal of Huaqiao University（Natural Science），2011，32（6）：710.

[3] ZHANG Xing，SHAN Shuang-rong.Explicit Difference Scheme of High Accuracy for Solving Four-Order Parabolic Equation[J].Journal of Huaqiao University（Natural Science），2010，31（6）：703-705.

[4] 詹涌強，張傳林.解拋物型方程的一族高精度隱式差分格式[J].應用數學和力學， 2014， 35（7）： 790-797.

[5] Shan Shuangrong.Difference Schemes of High Accuracy for SolvingParabolic Equation of Four Order[J].Journal of Huaqiao University（Natural Science），2003，24（1）：11-15.

[6] Jin Xianghua，Zeng Wenping.Several New Difference Schemes for Solving Fourth Order Parabolic Equation[J].Journal of Huaqiao University（Natural Science），2006，27（1）：238.