許 超,周家全,唐啟立

(1.洛陽理工學(xué)院 數(shù)理部,河南 洛陽 471023;2.湘潭大學(xué) 數(shù)學(xué)與計(jì)算科學(xué)統(tǒng)計(jì)學(xué)院,湖南 湘潭 411105)

血吸蟲廣泛分布于亞洲、非洲及拉丁美洲的70多個(gè)國家和地區(qū),血吸蟲病是人或哺乳動(dòng)物感染血吸蟲所引起的一種疾病,其傳播環(huán)節(jié)多、流行因素復(fù)雜,不僅嚴(yán)重危害人體健康,同時(shí)對家畜也會(huì)造成極大的危害而影響農(nóng)業(yè)和畜牧業(yè)的發(fā)展,從而嚴(yán)重影響疫區(qū)經(jīng)濟(jì)發(fā)展.因此,血吸蟲病的防治是世界上重要的公共衛(wèi)生問題之一.

1 有限元構(gòu)造

有限元空間Vh定義如下

??K,?K∈Th},

有限元空間Vh上插值算子Ih定義如下:Ιh:H2(Ω)→Vh,Ιh|Kv=IKv,?v∈H2(Ω),滿足

此外,論文中Wm,p(Ω)為標(biāo)準(zhǔn)的Sobolev空間,且Wm,2(Ω)=Hm(Ω),Hm(Ω)上的范數(shù)和半范數(shù)分別為||·||m和|·|m.當(dāng)m=0時(shí),記H0(Ω)=L2(Ω),L2(Ω)上的范數(shù)為||·||0.

2 血吸蟲病數(shù)學(xué)模型及其非協(xié)調(diào)有限元格式構(gòu)造

論文考慮如下人、牛共患的血吸蟲病數(shù)學(xué)模型[4]

(1)

由文獻(xiàn)[8]知,F(u,v)對u,v的偏導(dǎo)數(shù)直到二階有界,且0≤F(u,v)<1.在易受傳染人口總數(shù)不變的情況下,F(u,v)關(guān)于u,v滿足Lipschitz條件,即存在常數(shù)L1,L2,使得

|F(u,v)|≤L1|u1-u2|+L2|v1-v2|.

(2)

問題(2)的非協(xié)調(diào)有限元半離散逼近格式為:求(uh,vh)∈Vh×Vh,使得對任意φh∈Vh,有

(3)

為了進(jìn)行誤差分析,給出以下重要引理.

((u-Ihu),φh)=0,

(4)

(5)

(6)

其中:n為單元邊界?K的單位外法向量，這里及下文中出現(xiàn)的C均為與h無關(guān)的正常數(shù),且不同地方可取不同值.

3 收斂性分析

首先,給出血吸蟲病數(shù)學(xué)模型的L2模的最優(yōu)誤差估計(jì).

||u-uh||0+||v-vh||0≤Ch2(|u|2+|v|2+m0),

(7)

證明令

u-uh=(u-Ιhu)+(Ιhu-uh)?η+ξ,v-vh=(v-Ιhv)+(Ιhv-vh)?ρ+θ.

由類似文獻(xiàn)[12]的方法可得

||η||0+||ηt||0+h||η||h≤Ch2(|u|2+|ut|2),

(8)

||ρ||0+||ρt||0+h||ρ||h≤Ch2(|v|2+|vt|2).

(9)

由(2),(3)式和引理1的(4)式,可得如下誤差方程

(ξt,φh)+d1(ξ,φh)=k1(F(u,v)-F(uh,vh),φh)-

(10)

(θt,φh)+d2(θ,φh)=k2(F(u,v)-F(uh,vh),φh)-

(11)

在(10)式中取φh=ξ,可得

(ξt,ξ)+d1(ξ,ξ)

由Cauchy-Schwarz不等式、引理1的(6)式和Young不等式,可得

(||η||0+||ξ||0)||ξ||0+||ηt||0||ξ||0+h2|u|3||ξ||h),

又由于F(u,v)關(guān)于u,v滿足Lipschitz條件,有

||F(u,v)-F(uh,vh)||0≤L1||u-uh||0+L2||v-vh||0≤

L1(||η||0+||ξ||0)+L2(||ρ||0+||θ||0).

進(jìn)而可得

(12)

類似上述過程,在(11)式中取φh=θ,有

(13)

將(12),(13)式相加,再由(8)，(9)式可得

注意到ξ(0)=θ(0)=0,上式兩端從0到t取積分,再利用Gronwall引理,有

(14)

由(8),(9)和(14)式及三角不等式,即可得(7)式.

注1如果在上述過程中對邊界的估計(jì)利用引理1的(5)式,相比于(14)式只能得到如下結(jié)果

(15)

而無法得到L2模的最優(yōu)誤差估計(jì)(7).

下面給出H1模的最優(yōu)誤差估計(jì).

||u-uh||h+||v-vh||h≤Ch(|u|2+|v|2+m1)，

(16)

證明在(10)式取φh=ξt,由Cauchy-Schwarz不等式、Young不等式和引理1的(5)式,有

即

(17)

同理,在(11)式取φh=θt,有

(18)

將(17),(18)式相加,并利用(8)，(9)式，有

注意到ξ(0)=θ(0)=0,上式兩端從0到t取積分,再利用引理1的(5)式、Cauchy-Schwarz不等式、Young不等式、Gronwall引理以及(15)式,有

(19)

由(8),(9)和(19)式及三角不等式,即可得到(16)式.

完全類似于定理2的證明,利用引理1的(6)，(14)式,通過精細(xì)估計(jì)可以得到下面的超逼近結(jié)果.

||Ιhu-uh||h+||Ιhv-vh||h≤Ch(|u|3+|v|3+m0+m2),

(20)

進(jìn)一步地,類似于文獻(xiàn)[10],構(gòu)造滿足如下性質(zhì)的插值后處理算子I2h

I2hIhu=I2hu, ||I2hu-u||h≤Ch2|u|3, ?v∈H3(Ω)，

||I2hv||h≤C||v||h,?v∈Vh，

(21)

并由類似于文獻(xiàn)[10]的方法可得如下超收斂結(jié)果.

定理4在定理3假設(shè)下,有

||u-I2huh||h+||v-I2hvh||h≤Ch2(|u|3+|v|3+m0+m2).

(22)

注2引理1中的(4)，(6) 式是論文定理1，3和4成立的關(guān)鍵所在.可以證明正方形網(wǎng)格剖分下的旋轉(zhuǎn)Q1元[18]滿足上述引理,因而論文的結(jié)論對此非協(xié)調(diào)元也成立.