周?chē)?guó)霞，李 敏，陳豫眉

(西華師范大學(xué) 數(shù)學(xué)與信息學(xué)院，四川 南充 637009)

非定常Oseen方程的有限差分流線(xiàn)擴(kuò)散法

周?chē)?guó)霞，李 敏，陳豫眉

(西華師范大學(xué) 數(shù)學(xué)與信息學(xué)院，四川 南充 637009)

本文對(duì)非定常Oseen方程采用有限差分流線(xiàn)擴(kuò)散法，證明了此方程在FDSD格式下的穩(wěn)定性以及速度和壓力離散解的誤差估計(jì)。

非定常Oseen方程;有限差分流線(xiàn)擴(kuò)散法

粘性流體力學(xué)是流體力學(xué)中的一個(gè)重要分支。粘性流動(dòng)與生活密切相關(guān)，例如天氣的數(shù)值預(yù)報(bào)，空氣繞過(guò)汽車(chē)、地面建筑物和飛機(jī)的流動(dòng)，人造腎、人造心臟等與血液循環(huán)相關(guān)的問(wèn)題。不可壓縮粘性流動(dòng)的基本方程組為Navier-Stokes方程，它是不可壓縮流體的控制方程，對(duì)它的研究一直備受關(guān)注。但Navier-Stokes方程是偏微分方程組且對(duì)流項(xiàng)呈非線(xiàn)性，對(duì)其求解困難，且計(jì)算量大。而Oseen方程是Navier-Stokes方程的線(xiàn)性形式，它描述了一個(gè)在小雷諾數(shù)下的粘性不可壓縮流體的流動(dòng)，是研究Navier-Stokes方程的一個(gè)重要基礎(chǔ)。

在上世紀(jì)八十年代，Brooks和Hughes提出了流線(xiàn)擴(kuò)散法(SD)和流線(xiàn)迎風(fēng)Petrov-Galerkin法(SUPG)，此方法中增加的穩(wěn)定項(xiàng)只在流線(xiàn)方向起作用[1]。Johnson和Navert等人運(yùn)用流線(xiàn)擴(kuò)散法研究了含時(shí)間的線(xiàn)性問(wèn)題，獲得的解具有相應(yīng)的穩(wěn)定性和較高的精度[2-3]。而之前的有限元方法不具有這兩種特性，如Galerkin方法雖有高階精度但并不穩(wěn)定，甚至精確解的不光滑也會(huì)導(dǎo)致偽震蕩;又如人工粘性法(或迎風(fēng)法)包含了人工擴(kuò)散系數(shù)，雖然能得到穩(wěn)定性的解，但其精度不高。二十世紀(jì)九十年代，Saranen等人提出了在高雷諾數(shù)下非定常Navier-Stokes方程的流線(xiàn)擴(kuò)散有限元方法[4-5]。

1998年，孫澈等人提出了有限差分流線(xiàn)擴(kuò)散法，此方法既保留了傳統(tǒng)SD方法的本質(zhì)特征又簡(jiǎn)化了算法結(jié)構(gòu)，克服了SD方法計(jì)算量較大的不足[6-7]。對(duì)空間區(qū)域Ω∈Rd以及時(shí)間層t=tn到tn+1層，SD方法需要解決空間-時(shí)間區(qū)域Ω×[tn,tn+1]上的d+1維離散問(wèn)題，而FDSD方法只需解決Ω區(qū)域上的d維離散問(wèn)題。FDSD方法的計(jì)算規(guī)模和算法復(fù)雜度相當(dāng)于全離散的Galerkin有限元方法。

2001年和2002年，孫同軍和羊丹平先后運(yùn)用SD方法和FDSD方法研究了對(duì)流占優(yōu)的Sobolev方程[8-9]。2008年，駱艷和馮民富將壓力梯度投影和間斷有限元方法相結(jié)合，對(duì)可壓縮線(xiàn)性化Navier-Stokes方程提出了一種穩(wěn)定化間斷有限元格式[10]。2011年，劉程熙和朱登標(biāo)研究了瞬態(tài)Oseen方程的全離散局部投影穩(wěn)定化方法，給出了穩(wěn)定性證明，并使用一個(gè)具有正交性的插值得出誤差估計(jì)[11]。

本文安排如下: 第1部分給出了非定常Oseen方程的FDSD格式。第2部分給出此格式的穩(wěn)定性。第3部分給出速度和壓力相應(yīng)離散解的誤差估計(jì)[11]。

1 Oseen方程和FDSD格式

考慮非定常Oseen方程

(1)

其中Ω?R2為具有多邊形邊界?Ω的有界區(qū)域,u=u(x,t)∈R2為速度向量,p(x,t)∈R為壓力,f(x,t)為體力,u0(x)為初值向量函數(shù)。表示梯度,△表示Laplacian算子,ν>0為黏性系數(shù),α>0為標(biāo)量,流域a滿(mǎn)足·a=0。

令時(shí)間步長(zhǎng)Δt>0,N=[T/Δt],tn=nΔt(n=0,1,…,N),Θ={K}為區(qū)域Ω上的正則三角剖分.令0

定義

Qh={q∈H1(Ω)|q|K∈(P1(K)),?K∈Θ},

其中P1(K)表示K上的一次多項(xiàng)式。

1.1 FDSD格式

(2)

其中?n(x)?(x,tn),Rn=-()n為截?cái)嗾`差.

對(duì)方程(1)采用FDSD法:求(Un,Pn)∈Vh×Qh(n=1,2,…,N)滿(mǎn)足

=(fn,ν+δ1(a·ν+q)),?(ν,q)∈Vh×Qh,

(3)

U0=u0,

Un=0,在?Ω上,

其中δ1=ci(Δt)2,c1是獨(dú)立于h,ν,α的正常數(shù)。

2 FDSD格式的穩(wěn)定性

令v=Un,q=Pn則(3)為

=(fn,Un+δ1(a·Un+Pn))

(4)

分別估計(jì)上式中每一項(xiàng),結(jié)合δ1=c1(Δt)2,得

(a·Un+Pn,δ1(a·Un+Pn))+ν(Un,Un)

=δ1‖a·Un+Pn‖+ν‖Un‖2,

(αUn,δ1(a·Un+Pn))≤‖a·Un+Pn‖2+Pn‖2+α2δ1‖Un‖2

綜上所述,當(dāng)Δt適當(dāng)小且c1滿(mǎn)足0

(5)

記δ1=c1(Δt)2≤C,得

(fn,Un+δ1(a·Un+Pn))≤C‖fn‖2+c1‖Un‖2+‖a·Un+Pn‖2

(6)

由(5)和(6)得

(7)

將上式兩端同時(shí)乘以ζn-1,由0<ζn<ζn-1<1得

將上式進(jìn)行累加,結(jié)合ζn<ζi(i=1,2,…,n-1)得

當(dāng)Δt適當(dāng)小且0<4(c1Δt+α2δ1Δt)≤1/2時(shí),得

≤(1+16(c1Δt+α2δ1Δt))N

=(1+16(c1Δt+α2δ1Δt))T/Δt≤e16(c1+α2δ1)T,

由Gronwall不等式得

由上式可得下列結(jié)論。

定理1 當(dāng)Δt適當(dāng)小且c1滿(mǎn)足0

(8)

3 FDSD方法的誤差估計(jì)

+α‖un-Un‖2}≤C(h3+(Δt)2)。

(9)

令雙線(xiàn)性形式

則方程(2)和(3)等價(jià)于

(10)

和

(11)

由(10)和(11)得

(12)

當(dāng)Δt適當(dāng)小且0

(13)

估計(jì)(12)右端每一項(xiàng),得

|(Rn,δ1(a·ξn+γn))|≤‖a·ξn+γn‖2+4δ1‖Rn‖2,

對(duì)上式得

(diνηn,γn)=-(ηn,γn).

則

+α(ηn+ξn)+ν(ηn+ξn)+δ1(a·ηn+βn,a·ξn+γn)

+16δ1‖a·ηn+βn‖2+(16+16α2δ1)‖ηn‖2

綜上所述,得

(14)

由(13)和(14)得

+16δ1‖a·ηn+βn‖2+‖ηn‖2+Cδ1‖ηn‖2+A1‖ξn‖2+‖βn‖2

(15)

‖η‖W1,4(Ω)≤Ch‖u‖W2,4(Ω),‖η‖L∞(Ω)+h‖η‖L∞(Ω)≤Ch‖u‖W1,∞(Ω),

‖η‖+‖β‖+h(‖η‖+‖β‖)

≤Ch2(‖u‖2+‖p‖2),‖Rn‖2

‖a·ηn+βn‖2≤2‖a·ηn‖2+2‖βn‖2≤Ch2.

綜上所述,結(jié)合(15)得

由定理1,當(dāng)Δt適當(dāng)小且0

由上式,利用三角不等式可得

[1] Brooks A N,Hughes T J R.Streamline upwith/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations[J].Comput Meth Appl Mech Engrg,1982,32(1-3):199-259.

[2] Johnson C,Navert U.An analysis of some finite element methods for advection-diffusion problem[J].North-Holland Mathematics Studies, 1981,47:99-116.

[3] Johnson C,Navert U,PITKARANTA J.Finite element methods for linear hyperbolic problems[J].Comput Meth Appl Mech Engng,1985,45:285-312.

[4] Johnson C,Saranen J.Streamline diffusion problems for the incompressible Euler and Navier-Stokes equations[J].Math Comp,1986,47(175):1-18.

[5] Hansbo P,Szepessy A.A velocity-pressure streamline diffusion finite element method for the incompressible Navier-Stokes equations[J].Comput Meth Appl Mech Engrg,1990,84(2):175-192.

[6] Sun C,Shen H.The finite difference streamline diffusion method for time-dependent convection-diffusion equations[J].Numer Math J(English Series),1998,7(1):72-85.

[7] 張強(qiáng).不可壓縮N-S方程的差分流線(xiàn)擴(kuò)散法[J].計(jì)算數(shù)學(xué),2003,25(3):311-320.

[8] Sun T J.Finite difference streamline diffusion method for quasi-linear Sobolev equation with convection-dominated term[J].Numer Math J(Chinese Series),2001,23(4):340-356.

[9] Sun T J,Yang D P.The finite element difference streamline diffusion methods for Sobolev equation with convection-dominated term[J].Appl Math Comput,2002,125(2):325-345.

[10] 駱艷,馮民富.可壓縮Navier-Stokes方程的壓力梯度局部投影間斷有限元法[J].應(yīng)用數(shù)學(xué)和力學(xué)，2008,29(2):157-168.

[11] 劉程熙,朱登標(biāo).瞬態(tài)Oseen方程的全離散局部投影穩(wěn)定化方法[J].內(nèi)江師范學(xué)院學(xué)報(bào), 2011,4(26):1-6.

(責(zé)任編輯：張凱兵)

O241. 82

A

2095-4824(2013)06-0090-05

2013-09-22

國(guó)家自然科學(xué)基金項(xiàng)目(11271390)，西華師范大學(xué)基金項(xiàng)目(10B013)

周?chē)?guó)霞(1986- ),女,重慶人,西華師范大學(xué)數(shù)學(xué)與信息學(xué)院碩士。

李 敏(1986- ),女,山西朔州人,西華師范大學(xué)數(shù)學(xué)與信息學(xué)院碩士。

陳豫眉(1972- ),女,四川眉山人,西華師范大學(xué)數(shù)學(xué)與信息學(xué)院教授,博士。