求解非線性拋物組的新算法
——雙插值有限元

陳傳淼，胡宏伶

(湖南師范大學數學與計算機學院，中國 長沙 410081)

?

求解非線性拋物組的新算法
——雙插值有限元

陳傳淼*，胡宏伶

(湖南師范大學數學與計算機學院，中國 長沙 410081)

本文提出了雙插值有限元法求解一類非線性拋物組，它對未知函數和系數都采用了插值，于是某些常數矩陣可一次性計算好，每時間層組裝剛度矩陣很簡單.它是一種經濟格式.

非線性拋物組；有限元；雙插值；經濟格式

隨著科學技術的發展，非線性拋物問題出現在許多實驗領域，如高溫傳輸、核聚變、半導體、超導、石油開發、金融和圖像識別等.本文特別關注核聚變中的三溫計算問題.

在經典有限元法或有限體積法中，雖然有多種格式離散非線性拋物問題.但數值求解的主要困難是計算工作量極其巨大.

1)每個時間層的離散工作量巨大.

2)牛頓法求解非線性方程組需3～4次線性化(計算切矩陣).

3)必須計算數千數萬時間層.

三種困難交織在一起成為大規模求解的主要困難之一，因此發展高效算法有重要意義.

對一類具有散度型結構的非線性方程

ut-(a(u)u)+b(u)=f,u=0,x∈Γ,u(x,0)=ψ,x∈Ω.

(1)

直接插值u,w,b(u),稱為插值系統有限元(ICFEM)，它在每個時間層上的離散方程組可以一次性計算好，與tn無關，因此可節省離散時間，是一種經濟格式.Zlamal[1]在假設|uh|≤C時證明了最佳階O(h2);不要此假設，陳傳淼-Larsson-Zhang[2](1989)在分塊均勻網格三角形線元用超收斂技巧，首次證明幾乎最佳收斂性O(h2ln h).以后最佳收斂性及二次元的超收斂性也被研究.但是Kirchhoff變換只適合單個方程，對拋物組無效.

為求解非線性拋物組，我們提出雙插值有限元，即同時插值函數u及系數a(u),b(u),可一次性計算某些常數矩陣，每時間層組裝剛度矩陣很簡單.它是一種經濟格式.

1 雙插值有限元法(BIFE)

插值系數法的基本思想是：插值的函數與函數的插值在最佳階誤差意義不是等效的.

在d維單元σ上m個點xp,適當光滑函數u(x)及復合函數a(u)的l次插值分別為

容易證明最佳階誤差估計

‖Ina(u)-a(Inu)‖k,p,r≤C(u)hl+1-k,k=0,1.

(3)

因為‖u-Inu‖及‖a(u)-Ina(u)‖有這種估計.

我們討論擬線性拋物問題(也可以是方程組)

c(u)ut-Dj(a(u)Dju)+b(u)=f(t,x),w(0,x)=w0(x),w=0,x∈Γ.

(4)

其中系數a(u),c(u)為正定,b(u)為u光滑函數.

我們的新思想在于：取u為未知函數，同時引進a(u),b(u),c(u)作為相關的變量(雖然它們不是獨立變量).在每個單元σ直接作分片線性插值ui,aI(u),bI(u),cI(u).

為簡單起見，討論后向Euler-雙插值有限元Un∈Sh滿足

(5)

這里，v=Li,i=1,2,…,N,

是a,u,v的三線性型.

可以計算單元積分

Jσ(Li)=∫σ{aI(U)ULi+bI(U)Li}dx=

(6)

這里單元剛度矩陣系數和質量矩陣系數

關于i,j是對稱的，可以一次性計算好.于是單元測度矩陣可用簡單的乘法得到

類似地有單元質量矩陣Mσ(c),Mσ(b).對已給U，顯然Kσ(a),Mσ(c),Mσ(b)是對稱的.

組集所有單元積分并處理本質邊界條件，得到總體非線性方程組

Q(Un+1)≡M(cn)(Un-Un-1)+kK(an)Un+kM(bn)-kMFn=0.

(7)

這里M(cn),K(an)為對稱正定陣.由此看到，雙插值有限元法有以下3個優點:

1)像ICFEM一樣，矩陣Kijp,Mij可一次性算好，組裝Kij(U),Mij(U)很簡單；

2)對多未知變量方程組適用，只要系數a(u),b(u),c(u)與t,x無關即可.特別地，核聚變中的三溫方程正好可以變換為這種散度形式.

3)如何高效求解非線性方程組(7),是另一個重要問題[5]，將在其他論文討論.

[1] ZLAMAL M.A finite element method of the nonlinear heat equation[J].RAIRO Model Anal Numer,1980,(14):203-216.

[2] CHEN C M,LARSSON S,ZHANG N Y.Error estimates of optimal order for finite element methods with interpolated coefficients for the nonlinear heat equation[J].IMA J Numer Anal,1989,(9):507-524.

[3] 陳傳淼，黃云清.有限元高精度理論[M].長沙：湖南科技出版社,1995.

[4] 陳傳淼.有限元超收斂構造理論[M].長沙：湖南科技出版社,2001.

[5] HU H L,CHEN C M,PAN K J.Time extrapolation algorithm for parabolic problems[J].J Comput Math,2014,32(2):183-194.

(編輯 HWJ)

The Bi-Interpolation Finite Element Method —A New Algorithm for Solving Nonlinear Parabolic Systems

CHENChuan-miao*,HUHong-ling

(College of Mathematics and Computer Science,Hunan Normal University,Changsha 410081,China)

The bi-interpolation finite element method for solving nonlinear parabolic systems is proposed,in which both unknowns and their coefficients are interpolated,so some constant matrixes can be computed in one time,whereas at each time level the assembly stiffness matrixes are very simple.This is an economic scheme.

nonlinear parabolic systems; finte element method; bi-interpolation; economic scheme

10.7612/j.issn.1000-2537.2016.05.014

2016-06-30

國家自然科學基金資助項目(11301176,10771063)；湖南省自然科學基金資助項目(14JJ3070)

*通訊作者,E-mail：cmchen@hunnu.edu.cn

O241.82

A

1000-2537(2016)05-0081-02