设为首页 加入收藏 期刊导航 网站地图
  • 首页
  • 期刊
    • 数学与物理
    • 地球与环境
    • 信息通讯
    • 经济与管理
    • 生命科学
    • 工程技术
    • 医药卫生
    • 人文社科
    • 化学与材料
  • 会议
  • 合作
  • 新闻
  • 我们
  • 招聘
  • 千人智库
  • 我要投搞
  • 办刊

期刊菜单

  • ●领域
  • ●编委
  • ●投稿须知
  • ●最新文章
  • ●检索
  • ●投稿

文章导航

  • ●Abstract
  • ●Full-Text PDF
  • ●Full-Text HTML
  • ●Full-Text ePUB
  • ●Linked References
  • ●How to Cite this Article
Pure Mathematics 理论数学, 2011, 1, 73-79
http://dx.doi.org/10.12677/pm.2011.12016 Published Online July 2011 (http://www.hanspub.org/journal/pm/)
Copyright © 2011 Hanspub PM
H1-Galerkin Mixed Element Method for the Coupling
Nonlinear Parabolic Partial Equations*
Jinfeng Wang1, Yang Liu2, Hong Li2, Xiaoyu Li3
1School of Statistics and Mat h e matics, Inner Mongolia Finance and E conomics College, Hohhot
2School of Mathematical Sciences, Inner Mongolia University, H oh ho t
3College of Sciences, Inner Mongolia University of Technology, Hohh ot
Email: w45j85f@163.com; mathliuyang@yahoo.cn; smslh@imu.edu.cn
Received: Apr. 1st, 2011; revised: Apr. 13th, 2011; accepted: Apr. 15th, 2011.
Abstract: An H1-Galerkin mixed finite element method is discussed for the coupling nonlinear parabolic par-
tial equations. Semidiscrete and fully discrete schemes and optimal error estimates of the scalar unk nown and
its gradient are derived for problems in one space dimension, and it does not require the LBB consistency
condition. Finally, a numerical example is presented to illustrate the effectiveness of the proposed method.
Keywords: Coup ling Nonlinear Parabolic Partial Equations; H1-Galerkin Mixed Element Method; Backward
Euler’s Method; Optimal Error Estimates
耦合非线性抛物方程组的 H1-Galerkin 混合元方法*
王金凤 1,刘 洋2,李 宏2,李晓瑜 3
1内蒙古财经学院统计与数学学院,呼和浩特
2内蒙古大学数学科学学院,呼和浩特
3内蒙古工业大学理学院,呼和浩特
Email: w45j85f@1 63 .com; mathliuyang@yahoo.cn; smslh@imu.edu.cn
收稿日期:2011 年4月1日;修回日期:2011 年4月13 日;录用日期:2011 年4月15 日
摘 要:利用 H1-Galerkin 混合有限元方法讨论耦合非线性抛物方程组,得到一维情形的半离散和全离
散格式和未知存量函数和它的梯度的最优收敛阶误差估计,而且不用验证 LBB 相容性条件。最后,通
过数值例子验证了该算法的可行性。
关键词:耦合非线性抛物方程组;H1-Galerkin 混合元方法;向后欧拉方法;最优阶误差估计
1. 引言
考虑耦合非线性抛物方程组的初边值问题


 
 
112
223
00
,,,
,,,
,0,,0,,
,0, ,0,
t
t
uBukukuvfxtxt J
vBvkuvkvgxtxt J
u xtv xtxtJ
uxux vxvx x
 
 





(1)
这里常数表示变量 之间的相互作用,
123
,,kkk

,,
,uv
 
00
,, ,
耦合非线性抛物型方程组(1)描述了两类物质的
反应扩散过程及相互作用。在环境污染源扩散,液体
渗透,生物群体繁殖扩散等领域有广泛的应用。因此
对此类方程的数值研究是非常有价值的。对于发展方
程的数值方法的研究主要有差分方法、有限体积方法、
连续有限元方法、间断有限元方法、传统的混合元方
法等。其中传统混合元方法一直以来备受关注,得到
了很好的发展,如[1,2]中分别利用了混合有限元方法
研究了正则长波方程和双曲方程等。但传统的混合元
方法需要满足 LBB 相容性条件,限制了有限元逼近空
间的选取。为了克服这些不利因素,Pani[3]于1998 年
提出了 H1-Galerkin 混合有限元方法,该方法不同于传
统的混合元方法,该方法可以使逼近有限元空间 和
具有不同次数的多项式,而且H1-Galerkin 混合有
h
V
h
W

f
xtgx vx
1,2, 3n

xt u

是已知函数,
是具有 Lipschitz 连续边界  的凸有界区
域,
,
n
R

0,
J
T为时间区域,BB 表示热传导系数
(或反应扩散系数)。
12
,0
*基金项目:国家自然科学基金(11061021);内蒙古自治区高校科学
研究基金(NJ10006);内蒙古大学青年科学基金(ND0702)资助。
王金凤 等耦合非线性抛物方程组的1混合元方法
74 | H-Galerkin
限元方法不需验证LBB相容性条件。随后Pani 等人
将此方法进一步推广到了一些发展型积分微分方程
[4-6]、双曲型和伪双曲型方程[7,8]和Sobolev 方程[9]等。
但是对于耦合非线性抛物型方程组的研究还没有见
到,本文给出该非线性问题的 H1-Galerkin 混合有限元
方法的半离散和全离散格式的最优收敛阶误差估计,
并通过数值算例验证了算法的可行性。注意在本文中
作估计时所有的C都是与空间网格参数h和时间步长
无关的正常数。 t
2. H1-Galerkin 混合有限元格式
本文考虑问题(1)的一维情形

 

  
 
112
223
00
,,,
,,,
0,1,0,1,0,
,0,( ,0),
txx
txx
uBu kukuvfxtxtJ
vBvkuvkvgxtxtJ
ututvtvttJ
uxux vxvxx
 
 





(2)
这里
 

0,1 ,0,
J
T。
首先给出定义记号及引理, 内积用

2
L


,

表
示,




11
0|01 0HvHv v

,传统的Sobolev 空
间记作 ,1
,mp
Wp

,相应的范数为 ,
|| ||mp

,特
别当 2p

时,把 记作
,2m
W2
H
,其范数简记为 || ||m

。
引理[4] 若

于




0, ,0,tt T定义的可积函数,则
有积分不等式
 
2
00 0
dd d
tt
2
s
sC s



 s
(3)
下面对问题(2)应用 H1-Galerkin 混合有限元方法。
令,
x
x
qu v


可将原问题(2)化为一阶系统:


112
223
,
,
,
xx
tx
tx
uqv
uBqkukuvfxt
vB kuvkvgxt







(4)
为了对系统(4)应用 H1-Galerkin混合有限元方法,考
虑如下弱形式:求


,;,uvq



:0,T11
0
H
H
使得:





 
 


 



1
0
1
0
1
112
1
22 3
,,,,
,,,,
,,, ,,,,
,, ,,,,
xx x
xx x
txx x
txx x
uq Ha
vHb
qBq kqkqvufHc
Bkqvukg H


 
 






,d
(5)
对于(5c)和(5d)是通过分部积分并由 Dirichlet边界条 件




0,1, 0
tt
utut
和




0,1,0
tt
vtvt而得。
设 和分别是
h
Vh
W1
0
H
和的有限维子空间,对
1
H
1p

及正整数 有文献[7 ]中的逼近性质。 ,kr
则半离散 H1-Galerkin 混合有限元格式为:求


,;,
hhh h
uvq

:[0, ]h
TVW
h
,
h
d
使得:


 

 

 

112
22 3
,,,,
,,,,
,,, ,,,,
,, ,,,,
hhhh h
xx xh
hhhhh
xx xh
hhhhhhhh hhhhh
txx xh
hhhhhh hhhhhhh
txx x
uq Va
vVb
qBqkqkqvufWc
Bkqvukg W


 










(6)
为了得到误差估计首先定义 的投影,求
满足
,uv
,uv
hh h
V


,0, ,0,,
hhhhhh
x
xxxxxh
uu vvV




(7)
定义 ,q

的椭圆投影,求 满足
,
hh h
qW







,0,, 0,,
hhhhhh h
A
qq AW


 (8)
其中






,, ,
xx
A
ww
 
w,

是保证 A的1
H

正定,即


2
01
,||||,
设,, ,
hh h
uu vvqqh



 
 
0, 1j

,参见文
献[10]有:对于

1
11
kj
tt
jkjk
Chu u





 

,


0,
h
uC

u
(9)
1
11
kj
tt
jk
jk
Chv v





 

,


0,
h
vC

v
(10)
11
11
,
rj rj
t
jr t
j
r
Ch qCh q

 


 (11)
1
A
www w

H
。且易知


,A

是
有界的,且 00

为常数。 11
11
,
rj rj
tt
jr
j
r
Ch Ch

 



 (12)
Copyright © 2011 Hanspub PM
王金凤 等耦合非线性抛物方程组的1混合元方法 75
| H-Galerkin
3. 半离散格式误差估计
令hhhh
uuuu u u



 ;
hhhh
vvvvv v



 ;
hhhh
qqqqqq



 
 ;
hhhh



  

由式(5)~(8)可以得到误差方程为:

 








 


112
1
22 3
2
,,,,,
,,,, ,
,, ,,
,,,,
,, ,
,,,,
hhhh
xxx xh
hhhh
xxx xh
hhhhhhh
txx
hhh
th
hh hhhhh
txx
hhh
th
Va
Vb
Bk kqvqvuu
BWc
Bkqvqvuuk
BWd
 

 


,
h
h

 







 


 




(13)
定理 1若,, ,
 
00
hh
uu
 
00
hh
vv(0) (0)
hh
qq




0
hh

0,则

22 2
22
2min1, 1
2
11
10
1d
t
kr
hh t
kr
k
uu huuChuuqs





 






22 2
22
2min1, 1
2
11
10
1d
t
kr
hh t
kr
k
vv hvvChvvs






 






22 22
2min1, 1
2
11
0
1d
t
rk
hh
rr
qq hqqChqs





 






22 22
2min1, 1
2
11
0
1d
t
rk
hh
rr
hChq
 

 s


 








222 22min1,1
2
0 0
11
dd
t t
rk
hhh h
qqh qqsChs
 


 



其中
22 2222 2
111 1
111
tt t
kkr r
kk r
uu vv qq

 

 2
1
t
r

证明:因为,,


和

已知,故只需估计 ,,


和

。在(13a)和(13b)中分别取 h


和h


,利用
Cauchy-Schwarz 不等式,我们有
 
,
xx
CC
又因为 1
0
,
H

,所以由Poincare 不等式,得


,CC




  (15)
在(13c)中令h



,并使用(15),有



  (14)



2
22
10,
,, h
tt
BACC qvv





0,
0, 0,
hhhh h
Cvq qCu uCu

 


 


  
222
0, 0, 0,
hh
t
CCq CvCu

 

 
0, (C)



 (16)
在(13d)中令 h


,并使用(15),同理有




222
20,
,,
tt
BA CCq



 







0,
0, 0,
hh
CvCu C




 (17)
将(16)和(17)两式相加,得
 




22
2222 2222
12
d,,
dtt
BA BA CC
t



  


(18)
对上式两端关于时间从0到t积分,即得
Copyright © 2011 Hanspub PM
王金凤 等耦合非线性抛物方程组的1混合元方法
76 | H-Galerkin
 



22
2222 2222
12
00 0
,,d d
tt t
tt
BABAsCs Cs
  
 
 

d (19)
对上式应用 Gronwall 引理及积分不等式(3),可得



2222 222222
11
00
dd
tt
tt
s
Cs
 
 
 (20)
在(13c, d )中分别令,
hh
tt


,并使用(15),类似于上面的推导过程,有
 




22 22
22 2222
012
d,,
d
tt tt
BA BACC
t

 

 


 (21)
对上式两端关于时间从0到t积分,利用Gronwall 引理,可得
 


22 22222
012
0 0
d, ,
t t
tt tt
2
d
s
BA BACs
 






(22)
联合(14),(15),(20)和(22),有


2222 22222
10d
t
tt
CC
 
 
s (23)
和


2222 22222
10d
t
tt
CC
 
 
s (24)
联合(20),(22)~(24),(9)~(12)和三角不等式 ,即得
定理结论。
4. 全离散格式误差估计
对时间区间


0,T进行剖分,为了使用向后欧拉公
式,记 ,
n
tnttTM

 ,其中
M
为任意的正整数。
进一步令


nn
t

和

1nn

n
tt



 ,其中


10,C

T。则得到系统(5)的全离散向后欧拉格式
为:存在一系列

使得
0
,,,
nn n
n
UVQZ 
M
n




 

 

 

112
22 3
,,,,
,,,,
,,, ,,,,
,, ,,,,
nhnh h
xx xh
nhnh h
xx xh
nhnhnhnn nnhnhh
txx xh
nhnhnn nnhnhnhh
txx xh
UQ Va
VZ Vb
QBQkQkQVUZ fWc
,
Z
BZkQVUZkZgWd


 
 




 

 


(25)
为了得到误差估计,令
 



nhhnn
nnnn
ut Uututut Un


 
 ;










nhhnn
nnnn
vtVvtvtv tVn



  
 ;
 



nhhnn
nnnn
qtQqtq tq tQn

 
 ;












nhhn
nnnn
tZtttZ nn

 


 

在 处使用(5),(7)~(8),(25)可以得到误差方程为:
n
tt
 

 

 


 

112
1
22 3
,,,,
,,,,
,, ,,
,,,,
,, ,
,
nhn nhh
xxxh
nhn nhh
xxx h
nhnhn nhnnnn nnnnh
txx
nnhnhh
th
n hn hnnnnnnnn hnn h
txx
nnh
t
Va
Vb
Bk kqvQVuUZ
BWc
BkqvQVuUZk
 

 

,

 

 
 

 
 
 


2,, ,
nhhh
BWd
 













(26)
这里
 
,
nn
tn tntn tn
qt qttt
 
  。
Copyright © 2011 Hanspub PM
王金凤 等耦合非线性抛物方程组的1混合元方法 77
| H-Galerkin
0

定理 2 设 ,则对于1
 
00
00,0
hh
QqZ



J
M

,有
11
JJJJ JJ JJJJJJ
qQZuUvVhuUhvV



  
11111 121 2
min1, 1|| ||||||||||
kkrkkrr
rk tt t
LH LH LHLH LHLH LHLH
Chu uv vqq

 

1r
t



 


 
22 22
tt t
LL LL
tq


 


证明:在(26a,b)中分别取 hn


和hn



,
易得




,
nnnnn
xx
CC
n

 
  (27)
再由 1
0
,
nn
H

,所以由Poincare 不等式,得




,
nnnnn
CC
在(31c,d)中分别令 ,
hnhn

 

,并将两式相
加,同时使用
 
22
11
,,,
22
nnn nnn
tttt




和Young-不等式,可得
n

 
  (28)

2222
12
11
222222222
11
22
nnnn
tt
nnnnnnnnnn
tt
BB
C


 
  2
(29)
这里
11
22
22
1
1d, d
nn
nn
tt
nn
tt tt
s
r
tt
qs Ctq
t



 

,
11
22
22
1
1d, d
nn
nn
tt
nn
tt
r
tt
tt
s
Ct s
t
 


 

上式两端关于 n = 1到J,

1
J
M求和,可得





2222 22222
00
111
11
12
JJ
JJ nnnnnn
nn
CtBtC t


  

2








12
222 2
2
21
00
0
dd
JJ
J
tt
rnn
t ttt ttn
hq stqst






 2

0
(30)
选择 ,当 时1,使用离散
Gronwall 引理,并注意到
0
t0
0tt Ct
00
0, 0


,应用三角不
等式得定理结论。
5. 数值算例
为了验证所提出算法的可行性,时空区域都取为


0, 1
,0ux
,对(2)中参数取值为 1,初值
,,此时
123 12
,,,,kkkBB
 
,0 sinvxx
 
sin x




,
f
xt

和

22ex


sin
tt
exsin





sin
t
,
g
xt ex


,容易验证精确解


sin
t
e


2
 x


,uxt

。采用分片线性基函数,时间离散
采用欧拉向后差分。表 1中分别给出了变量

vxt

,sin
t
e



x
h
uu

和
的空间 时间最大模的误差估计结果和 收敛
阶。表 2分别给出了变量
h
vv2
L
h
qq

和h


的空间 时
间最大模的误差估计结果和收敛阶。从数据上可以看
出,该数值方法对非线性耦合方程组问题是有效的。
图1~4中分别给出了在时间和空间步长分别
2
L
t


和
0.0125 0.025h

时的四个数值变量 和,,
hhh
uvq h

的表面图像,并且在图 5~8中分别给出了在同样剖分
下的真解 ,,,uvq

和数值解 uv,,,
hhh
qh

的对比图像。
从图像上直观的看到了该数值算法的可行性。
Table 1.


2
LL
-errors and order of convergence
表1.


2
LL
-误差和收敛阶
2ht



2
h
L
L
uu
Order

2
h
L
L

vv Order
1
40 1.03422E-03- 9.12765E-04 -
1
80 4.48740E-041.20459 5.48770E-04 0.73404
1
160 2.07280E-041.11430 2.97770E-04 0.88200
1
320 9.93750E-051.06063 1.54770E-04 0.94407
1
640 4.86210E-051.03130 7.88600E-05 0.97276
Copyright © 2011 Hanspub PM
王金凤 等耦合非线性抛物方程组的1混合元方法
78 | H-Galerkin
Table 2.


2
LL
-errors and order of convergence
表2.


2
LL
-误差和收敛阶
2ht

2
h
L
L
qq 
 Order

2
h
L
L


 Order
1
40 2.37388E-03 - 3.87404E-03 -
1
80 1.19331E-03 0.99228 1.97945E-03 0.96874
1
160 5.98310E-04 0.99600 1.00032E-03 0.98464
1
320 2.99580E-04 0.99795 5.02810E-04 0.99238
1
640 1.49900E-04 0.99894 2.52070E-04 0.99619
Figure 1. The surface sho ws the value of uh
图1. 数值解 uh的表面图
Figure 2. The surface sho ws the value of vh Figure 3. The surface shows the value of qh
图2. 数值解 vh的表面图 图3. 数值解 qh的表面图
Figure 4. The surface sho ws the value of

h Figure 5. Comparison between uh and u
图4. 数值解 h

的表面图 图5. 数值解 uh和真解 u的对比单线图
Copyright © 2011 Hanspub PM
王金凤 等耦合非线性抛物方程组的1混合元方法 79
| H-Galerkin
Figure 6. Comparison between vh and v Figure 7. Comparison between qh and q
图6. 数值解 vh和真解 v的对比单线图 图7. 数值解 qh和真解 q的对比单线图
Figure 8. Comparison between h

and

图8. 数值解 h

和真解

的对比单线图
参考文献 (References)
[1] Z. D. Luo, R. X. Liu. Mixed finite element analysis and numeri-
cal solitary solution for the RLW equation. SIAM Journal on
Numerical Analysis, 1998, 36(1): 89-104.
[2] Y. P. Chen, Y. Q. Huang. The superconvergence of mixed finite
element methods for nonlinear hyperbolic equations. Communi-
cations in Nonlinear Science and Numerical Simulation, 1998,
3(3): 155-158.
[3] A. K. Pani. An H1-Galerkin mixed finite element method for
parabolic partial differential equations. SIAM Journal on Nu-
merical Analysis, 1998, 35(2) : 712-727.
[4] A. K. Pani, G. Fairweather. H1-Galerkin mixed finite element
methods for parabolic partial integro-differential equations. IMA
Journal of Numerical Analysis, 2002, 22(2): 231-252.
[5] D. Y. Shi, H. H. Wang. An H1-Galerkin nonconforming mixed
finite element method for integro-differential equation of para-
bolic type. Journal of Mathematical Research and Exposition,
2009, 29(5): 871-881.
[6] 王瑞文. 双曲型积分微分方程的 H1-Galerkin 混合有限元方法
误差估计[J]. 计算数学, 2006, 28(1): 19-30.
[7] A. K. Pani, R. K. Sinha, and A. K. Otta. An H1-Galerkin mixed
method for second order hyperbolic equations. International
Journal of Numerical Analysis and Modeling, 2004, 1(2):
111-129.
[8] Y. Liu, H. Li. H1-Galerkin mixed finite element methods for
pseudo-hyperbolic equations. Applied Mathematics and Compu-
tation, 2009, 212(2): 446-457.
[9] 郭玲, 陈焕贞. Sobolev方程的 H1-Galerkin混合有限元方法[J].
系统科学与数学, 2006, 26(3): 301-314.
[10] M. F. Wheeler. A priori L2-error estimates for Galerkin ap-
proximations to parabolic differential equation. SIAM Journal on
Numerical Analysis, 1973, 10(4): 723-749.
Copyright © 2011 Hanspub PM

版权所有:汉斯出版社 (Hans Publishers) Copyright © 2012 Hans Publishers Inc. All rights reserved.