一類四階方程基于降階格式的譜 Galerkin 逼近及誤差估計

打開文本圖片集

摘要：本文針對一類四階方程提出了一種基于降階格式的有效譜 Galerkin 逼近.首先，引入一個輔助函數(shù)，將四階方程化為兩個耦合的二階方程，并推導(dǎo)了它們的弱形式及其離散格式.其次，利用 Lax-Milgram 引理和非一致帶權(quán) Sobolev 空間中正交投影算子的逼近性質(zhì)，嚴(yán)格地證明了弱解和逼近解的存在唯一性及它們之間的誤差估計.最后，通過一些數(shù)值算例，數(shù)值結(jié)果表明該算法是收斂和高精度的.

關(guān)鍵詞：四階方程；降階格式；譜 Galerkin 逼近；誤差估計

中圖分類號： O241   文獻(xiàn)標(biāo)識碼：A   文章編號：1009-3583（2024）-0081-04

Spectral Galerkin Approximation and Error Estimates Based on Reduced Order Scheme for A Class of Fourth Order Equations

WANG Yuan-lu ， JIANG Jian-tao

（School of Mathematical Science， Guizhou Normal University， Guiyang 550025， China）

Abstract： In this paper， we propose a spectral Galerkin approximation and error estimates based on reduced order scheme for a class of fourth order equations. Firstly， by introducing a auxiliary function， we transform the original problems to two coupled second order equa- tions， and their weak form and corresponding discrete format are also derived. Secondly， by using Lax-Milgram lemma and the approxi- mation properties of orthogonal projection operators in non-uniform weighted Sobolev spaces， we strictly prove the existence and uni- queness of weak solution and approximate solution and as well the error estimate. At the end， we conduct some numerical experiments， which show that the algorithm is convergent and high accurate.

keywords： fourth order equation; reduced order scheme; spectral Galerkin approximation; error estimation

四階方程在流體力學(xué)、物理學(xué)和生物學(xué)等多個領(lǐng)域有著廣泛應(yīng)用，許多復(fù)雜的非線性問題[1]的計算最終都?xì)w結(jié)于反復(fù)求解一個四階方程問題，如Cahn- Hiliard 方程[1-4]、擴(kuò)展的 Fisher-Kolmogorov 方程[5-8]等.因此，提出一種有效求解四階方程的高精度數(shù)值方法是非常有意義的.

到目前為止，已有很多數(shù)值方法求解不同邊界條件下四階方程問題，主要包括有限元法[9-12]、譜方法及其他高階數(shù)值方法[13-19].對于有限元方法，區(qū)域剖分和對連續(xù)可微的有限元空間的要求將產(chǎn)生大量的自由度，要獲得高精度的數(shù)值解將需要很多計算時間和內(nèi)存容量.

眾所周知，譜方法是一種具有譜精度的高階數(shù)值方法[19-20]，是求解微分方程的一種重要的數(shù)值方法.據(jù)文獻(xiàn)所知，很少提到關(guān)于四階方程基于降階格式的譜方法.因此，本文針對一類四階方程提出了一種基于降階格式的有效譜 Galerkin 逼近.該方法首先是通過引入一個輔助函數(shù)，將四階方程化為兩個耦合的二階方程，并推導(dǎo)它們的弱形式及其離散格式.其次，利用 Lax-Milgram 引理和正交投影算子的逼近性質(zhì)，嚴(yán)格地證明弱解和逼近解的存在唯一性及它們之間的誤差估計.最后，給出了一些數(shù)值算例，通過數(shù)值結(jié)果表明算法是收斂和高精度的.

1降階格式及其變分形式

作為一個模型，本文考慮下面的四階問題

令 .則方程（1）可化為下面兩個耦合的二階問題

引入通常的Sobolev 空間： ； ， 相應(yīng)的內(nèi)積和范數(shù)分別為，其中I=（-1，1），則（2）、（3）的弱形式為：找， 使得

2 誤差估計

3算法的有效實(shí)現(xiàn)

4數(shù)值實(shí)驗

為了測試該算法的有效性，本文將給出2個數(shù)值算例，并在MATLAB 2019b 平臺上進(jìn)行編程計算.

例1：令 ，顯然滿足邊界條件. 將代入（3）可得，再將代入（2）可得f.對不同的 N，精確解，和逼近解 N，N 在三種模下的誤差分別在表1和表2中列出，方程（2）和（3）的精確解和逼近解的擬合圖及誤差圖分別在圖1和圖2中給出.

從表1和表2可知，當(dāng)N≥25時，逼近解 N，N 達(dá)到大約10-14的精度.另外， 從圖1和圖2中可以看出該算法是收斂的和高精度的.

例2：令 ，顯然滿足邊界條件。（剩余3842字）