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

期刊菜单

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

文章导航

  • ●Abstract
  • ●Full-Text PDF
  • ●Full-Text HTML
  • ●Full-Text ePUB
  • ●Linked References
  • ●How to Cite this Article
International Journal of Fluid Dynamics 流体动力学, 2013, 1, 26-33
http://dx.doi.org/10.12677/ijfd.2013.12005 Published Online June 2013 (http://www.hanspub.org/journal/ijfd.html)
A Stabilized Crank-Nicolson Finite Volume Element
Formulation for Non-Stationary Stokes Equation*
Hong Li1, Zhihui Zhao1, Zhendong Luo2
1School of Mathematical Sciences, Inner Mongolia University, Huhhot
2School of Mathematics and Physics, North China Electric Power University, Beijing
Email: malhong@imu.edu.cn, zhaozhihui2328@sina.com, zhdluo@ncepu.edu.cn
Received: Mar. 18th, 2013; revised: Apr. 2nd, 2013; accepted: Apr. 18th, 2013
Copyright © 2013 Hong Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unre-
stricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract: A stabilized Crank-Nicolson (CN) finite volume element formulation with time second-order accuracy is
established for two-dimensional non-stationary Stokes equation. The error estimates of its numerical solutions are pro-
vided. Some numerical experiments are presented illustrating that the stabilized CN finite volume element formulation
with time second-order accuracy is far more advantageous than that with time first-order accuracy, thus validating that
the stabilized CN finite volume element formulation is feasible and efficient for finding the numerical solutions for two-
dimensional non-stationary Stokes equation.
Keywords: Stokes Equation; Stabilized Crank-Nicolson Finite Volume Element Formulation; Error Estimate
非定常 Stokes 方程的稳定化 CN 有限体积元格式*
李 宏1, 赵智慧 1, 罗振东 2
1内蒙古大学数学科学学院,呼和浩特
2华北电力大学数理学院,北京
Email: malhong@imu.edu.cn, zhaozhihui2328@sina.com, zhdluo@ncepu.edu.cn
收稿日期:2013 年3月18 日;修回日期:2013年4月2日; 录用日期:2013 年4月18 日
摘 要:建立二维非定常 Stokes 方程的时间二阶精度的稳定化 Crank-Nicolson (CN)有限体积元格式, 并给出其
稳定化 CN 有限体积元解的误差估计。数值实验说明时间二阶精度的稳定化CN 有限体积元格式比时间一阶精
度格式更优越, 从而表明稳定化 CN 有限体积元格式对于求解非定常Stokes 方程的数值解是有效可行的。
关键词:Stokes 方程;稳定化Crank-Nicolson 有限体积元格式;误差估计
1. 引言
设 是有界的连通凸多边形区域。考虑流体
动力学中的二维不可压的非定常Stokes方程
2
R
tp

uu f


,在 中 (1.1)

0,T
div 0u,在 中 (1.2)

0,T


,,
x
ytu

,在 上 (1.3)

0,T





,,0 ,
x
yxu

y

,在 中 (1.4)
其中是未知的流体速度向量, 是压力,
T是总体时间,

12
,T
uuup
1Re


,是 Reynolds 数, 是
已知的源函数,
Re f


,,
x
yt

和

,

x
y

分别是已知的边值
函数和初值函数。为了便于理论分析而且不失一般
性,不妨在下面的理论分析中假定

,,
x
yt

和


,
x
y

*资助信息:国家自然科学基金(批准号:11271127 和11061021)和
内蒙古自然科学基金(批准号:2012MS0106)。
Copyright © 2013 Hanspub
26
非定常 Stokes 方程的稳定化 CN 有限体积元格式
均为零向量。
非定常 Stokes 方程是流体力学的重要方程,已广
泛和成功地用于许多实际工程领域[1-5]。然而当实际问
题的计算域不规则或问题本身比较复杂时,要求出其
解析解是很困难的,有效的方法是求其数值解。由于
有限体积元法[6-8]能保持局部的质量或能量守恒,其精
度比有限差分法高而且能适应边界复杂的计算域,并
与有限元方法有同样精度但要比有限元方法便于计
算,因此有限体积元法被公认是最有效的数值计算方
法之一。有限体积元法也称为盒子法[8]或广义差分方
法[10],已经广泛应用于求解各种类型的偏微分方程(例
如,二阶椭圆方程、抛物型方程、Stokes 方程等)的数
值解[6-13]。但是求解非定常Stokes 方程的有限体积元
法的流速和压力的试验空间要求满足离散的 Babuska-
Brezzi (BB)条件[12],这对理论分析和实际应用都会产
生很多困难。为了避开离散 BB条件对 Stokes方程有
限体积元法的限制,一些基于两局部高斯积分或自由
参数的时间一阶精度的有限体积元格式已被提出
[13-15]。然而,据我们所知,到目前为止还没有非定常
Stokes 方程基于两局部高斯积分的时间二阶精度全离
散稳定化 Crank-Nicolson (CN)有限体积元格式及其误
差分析报道。因此,本文给出非定常 Stokes 方程的时
间二阶精度的全离散稳定化CN 有限体积元格式及误
差分析,并用数值例子说明时间二阶精度的稳定化
CN 有限体积元格式比时间一阶精度格式更优越,从
而表明用稳定化 CN 有限体积元格式去求解非定常
Stokes 方程的数值解是更有效可行。
本文的安排如下,第 2节建立非定常的Stokes方
程关于时间的 CN半离散格式。第3节直接从时间的
CN 半离散格式出发, 建立非定常的 Stokes 方程的基
于两局部高斯积分的全离散稳定化CN 有限体积元格
式,并做误差分析,这样可以避开关于空间变量的半
离散化 CN 有限体积元格式的讨论,使得论证过程更
简便。第 4节给出数值例子说明时间二阶精度的稳定
化CN 有限体积元格式比时间一阶精度格式更优越,
从而表明用稳定化 CN有限体积元格式去求解非定常
Stokes方程的数值解是更有效可行。第 5节是主要的
结论和讨论。
2. 关于时间半离散化 CN 格式
本文用到的 Sobolev 空间[16]都是标准的。设
 


2
12
0,;dUHM qLqxy
d0

 
。则 非定 常
Stokes 方程的变分形式为:
问题 1. 求


,pUM

u使得对于任意


0,tT
满足





 
,,,,,
tabp U

uvuvvfv v (2.1)


,0,bqq Mu (2.2)


,,0 ,xy

0u在 中 (2.3) 
其中 ,
 
,,,,,diabq

uvu vvv

v q


,

为


2
L

或

2
2
L

或 中的内积。

22
2
L

注意到双线性型


,bqv满足 BB条件[1-5]:


0
0
,
sup ,
U
bq qqM




v
v
v (2.4)
其中

是正常数。双线性型和 满足有界
性

,auv

,bqv



00
,,aU

,

 uvuv uv (2.5)


00
,,bqCqUqMvvv,

(2.6)
其中 C是正常数。双线性型 还满足正定性

,auv

2
0
,,a

U

vvv v (2.7)
归咎于(2.4)~(2.7),当 时,问题1的
解是存在唯一的[1-5]。

2
1
H
f
设 为正整数,时间步长为NkTN,n
tnk

,
为 在
n
uu


0,1, ,
n
tn N点关于时间的半离散化逼
近。则问题1关于时间 t的半离散化 CN格式为
问题 2. 求


,
nn
pUM

u

1, 2,,nN满足





12 11
2,, 2,
2,2, ,,
nn n
nnn
kakb p
kka


U

 
uv uvv
fv uvuvv
(2.8)


,0,
n
bqq Mu (2.9)
0

0u,在 中 (2.10)

其中


12
12
nn
t


ff。
利用[2,15]中类似的方法可以证明问题 2存在唯
一的解




,1,2
nn
pUMnN u,,满足

2
22
11
00
0
11
nn
niin
ii
kpCk



 



uuu f
2
0
i
(2.11)
Copyright © 2013 Hanspub 27
非定常 Stokes 方程的稳定化 CN 有限体积元格式
而且当问题1的解
 
2
12
,0,;0,;pHTUL TMu

时,有下面的误差估计




0
12 2
0
0
n
n
nn
nn
t
kt ptp


 


uu
uu Ck
(2.12)
其中是与k无关的常数。 C
3. 稳定化的全离散 CN 有限体积元格式及
误差分析
3.1. 有限体积元原理
为了获得问题2的有限体积元解,需要对计算域
做下面的三角形剖分和对偶剖分。 
设

h
K
 是具有最大直径


max h
K
K
hh

的
拟一致三角形剖分,其中
K
h是三角形 的直径
[4,5]。为了刻画有限体积元方法,我们还要引入基于
h
K
h

的对偶剖分 ,它的单元称为控制元,控制元的构造
与[10]中的构造相同。设
h


K
z是单元 的重心,将
h
K
K
z与三角形 K的各边中点连结,把K剖分成三个小
四边形
z
K
(如图 1的左图)。记

h
Z
K为K的顶点,则

h
hh
K
Z
ZK

为 的顶点集合。对于每个
h
h
zZ

,
由所有享有顶点 z的小四边形构成控制元
z
V(如图1
的右图)。所有的控制元覆盖

,构成 的重心型对
偶剖分 。用
h

h

h
Z
表示剖分 的顶点集合
h
h
Z
的内部顶
点集合。
重心型对偶剖分 称为是拟一致的,如果存在两
个不依赖于h的正数 和使得
h


1
C2
C
Z
k
K
Z
K
V
Z
Z
Z
Figure 1. Left chart is a triangle K partitioned into three
sub-regions Kz. Right chart is a sample region with dotted lines
indicating the corresponding control volume Vz
图1. 左图是将三角形 K剖分成三个子四边形 Kz。右图是用虚线连
成顶点为 Z的控制元 Vz的样本


22
12
,
z
zh
Chm VChV


其中


z
mV 表示
z
V的测度。若三角剖分 是拟一致
剖分,则相应的对偶剖分 也是拟一致的[10]。
h

h


定义 U和M的试验空间分别为






22
1
1
;,
;,
hh hh
K
hh hh
K
UCUPKK
MqMq PKK

 
 
vv
其中


1
PK是h
K

上的线性多项式构成的空间(可
由K的顶点唯一确定)。
设h

为U到上的插值算子。则由 Sobolev空
间的插值理论[4,5]知,当 时,有
h
U

2
2
Hu
2
2
π,0,
m
hmCh m

 uuu 1 (3.1)
其中 C及后面出现 C均表示与空间网格尺寸 h和时间
步长 k无关的正数,不同的出现可能不等。
速度的检验函数空间 取为
h
U

 





2
2
0
;,
0,
z
z
hhhz zhh
VV
zzh
UL PVV
VV



 



vvv
事实上, 可由下面的基函数
h
U

 

1, ,,
,0, ,,
z
z
h
z
xy V
x
yz
xy V


Z






张成。从而对于每个向量 有
hh
U
v
 
,
h
hhz
zZ
zxy





vv
对于 U

w,设是在上的插值,即 πh
wwh
U

 
π,
h
hz
zZ
zxy




ww
则由插值理论[4,5]有
1
0
πhCh

ww w
(3.2)
3.2. 稳定化全离散 CN 有限体积元格式
由于 hh
UUU


,所以需要双线性型


,auv 和


,bqv做修改。利用对偶元
z
V上的分部积分有
Copyright © 2013 Hanspub
28
非定常 Stokes 方程的稳定化 CN 有限体积元格式
Copyright © 2013 Hanspub 29
d

2
00 ,
hpppL


  (3.10)

dd dd
dd
zz
zh
zz
zh zh
VV
V
VV
VV
xy xy
s
xy




 
 





vu vu
vunvu
dd
h
(3.3)

1
10 ,
h
ppChppH


 (3.11)
 

,,
,
hhhhhh
hhhhhh
Dpq ppq
ppqq



 (3.12)
dd dd
ddiv
zz
zh
zz
zh zh
VV
V
VV
VV
pxy pxy
psp xy




 
 





vv
vnv (3.4)
则问题 2的稳定化全离散 CN有限体积元格式为:
问题 3. 求



,1,
nn
hhh h
pUMn N u2,,满足
当但 时,定义
hh
Uuh
U
v



,d
z
zh
hhhh h
V
V
as




uvv un (3.5)



 

121 1
2,π,,
2, ,2,π,
nnn
hhhhh hh
nnn
hhhhh
hh
kakb p
kka
U


 


uv uvv

h

f
vuvuv
v
(3.13)

,
z
zh
hhh hh
V
V
bpp s



vdvn
U
d
(3.6)




,,,
n
hhhhh h
bqDp qqMu (3.14)
利用 Green 格式可以证明下面两式成立[13-15]
0,
h

0u在 中 (3.15) 


,π,,,
hh hhhhhhh
aa
u vuvuv (3.7)
3.3. 全离散 CN 有限体积元解的误差分析


,π,, ,
hhhhhhhh hh
bpbpUp M
 vvv (3.8)
为了证明问题 3的解的存在唯一性及做误差分
析,需要引进下面两引理。
为了定义稳定化格式,引进基于两局部高斯积分
和自由参数的双线性型 引理 1 [13-15]. 有下面的对称性成立

,2 ,1
,
dd dd
h
hh
hh hh
KK
K
Dpq
pq xypq xy







 
,π,π,
,
hhh hhh
hh h
U



uv vu
uv (3.16)
对于任意 都有
 
20, 1
m
Hm u
其中表示在 K上精度为 i的
,
,,d
hh hhh
Ki
pq Mpqxy




,,π,
0, 1
mn
hhh h
mn
Ch
n



uvu vuv (3.17)
高斯积分,

,hh
g
xy pq是次数不超过 次的
多项式。于是,对于所有的检验函数 ,当

1, 2ii
hh
qM1i

时,试验函数必须是分片常数。再定义 投
影算子 如下
hh
pM

2
hh
LW
2
L
:

定义

12
0,π
hhhh

uuu,则 0
与0
等价,即
存在两个与时间步长 k及h无关的正常数 和 使
得
1
C2
C
 

2
,,
,
hhh
hh
qq qq
qLq W


 
,
(3.9)
1
00
0,
hh h
hh
CC
U


uu u
u
2
(3.18)
其中 表示相应于

2
h
WL h

的分片常数空间,即
在每个 是常数的空间。那么,投影算子
h
K h

有
下面的性质[4,5] 引理 2. 有下面上下确界不等式成立[13-15]







00
,
000
sup2 ,π,2,2, 2,
,,
hhh h
hhhhhh hhhhhhh
qUM
hh hhhh
kabqbp Dpqq
kppUM




 

 
vuv uvuvv
uu
非定常 Stokes 方程的稳定化 CN 有限体积元格式
归咎于引理 2,由 Lax-Milgram 定理[4,5]可得问题
3存在唯一的解 ,并
有


,1,
nn
hhh h
pUMn Nu2,,

2
21
00
0
1
2
12
0
1
n
nii
hhh
i
ni
i
kp
Ck





 





uuu
f
i
h
,2,,
(3.19)
利用 Stokes 方程的经典有限元方法[2-5]容易得到
下一引理。
引理 3. 设 是一个
Stokes 投影对,即对于

存在唯一的序
列对 满足

,:
hhhh
SQU MUM
,
nn
up UM


,1
nnhh
pUMn N
hh
SQu








,,,2,
,,, ,,
nn n
hh hhhh
nn hhhhhh
A
SQpqkDQpq
A
pq qU


uv
uv vM
h
(3.20)
其中







1
1
,,,
,2,2,
2,
nn hh
nn n
hh
nn
h
Apq
kabq bp



 


uv
uuvu v
uuv
,
并满足


00 0
nn n
hh
SQpC p uu
0
n
(3.21)
而且当 有

 
2
21
,
nn
up HH

2
22 4
00 0
1
n
nnnnii
hh h
i
SkpQpk S



uu uuCh
(3.22)
分析问题 3解的误差,还需要用到下面的离散
Gronwall 引理。
引理 4 (离散 Gronwall 引理[1,5,15]). 如果 和
及 是三个正的数列,且 是单调的,并满足
n
an
b
n
cn
c
1
00 0
0
;,
n
nn ni
i
abcabcan



 
1,0
那么有

exp, 0
nn n
abcnn

 
对于问题3的解,有下面的收敛性结果。
定理 5. 如果问题 I的

,
 
2
21
,pH Hu
而且


kOh充分小和 或,则
00
0
h
pp00
hh
pQp




1,2,N,
nn
hh h
pU ,
h
Mnu有下面的误差估计
 



12
00
0
1
22
n
nn
nhnhnh
i
kptpktt
Ck h



uuuu
n
(3.23)
证明. 在问题 2中取



,,
hh
qqvv

并与问题3相
减可得到下面的误差方程









11
1212
,,,2,
2π,π2π,π
2π,π,
,
nnnn n
hhhh hh
nn nn
h hhhhhhhhhhh
nn
hhhh
hhhh
AppqkDpq
k
qUM
 

 

 

uuv
uuvvuuvv
ffvv
v
(3.24)
令


1
,2,enn
h
Qp p
nnnnnn
hh
S
 Euuuuu n
h
。则
由(3.20)和误差方程(3.24)有






22
00
1
11
11
11212
1
2
,,
1
,,,
2
,,
,π,π
π,π,
1,π,π
2
nn
nnn nnn
hh
nnnnnnn n
hh
nnnnnn
hh
nnnnnnn
hhhh h
nnnnnnn
hhh hh
nnnnnn
hh
k
SkaS
ka
SkbQpp
kbp p
k







 

 

 


EE
uuEuuE
uuEuuE EE
uuEE
EuuEE
uuEEuuE
EEffEE




12121 1
11
1
π,π,
2
π,π,
nnnnnn
hh
nnnnnnnn
hh hhhh
k
kb e

 

 
ffEEEE
uuuuEEE
(3.25)
由引理 1、Hölder 不等式和 Cauchy 不等式,当


kOh时,有


1212
22
212412
10 1
π,π
4
nnnn
hh
nn n
k
k
Ckh Ckh





ffEE
fE fE
0
n

(3.26)
由逆估计定理[4,5]和Taylor 公式,有
Copyright © 2013 Hanspub
30
非定常 Stokes 方程的稳定化 CN 有限体积元格式




1
11
21
10
2
2
33
00
22
31
00
22
332
0
22
31
00
π,π
4
4
nnnnnn
hh hhhh
nn n
hh
nn
hh
nn
ntLH
nn
Ch
ChChSS
k
Ch
ChCh k
k
Ch



 




 

 
 

 
uuuuEE
uuE
Eu
EE
Eu
EE
1n
u
(3.27)
由(3.9)和(3.14) 及Hölder 不等式和Cauchy不等式
有


11
22
11
00
,,
2
,
2
44
nnn nn n
hh
nnn
h
nn nn
hh
be eeee
eeee
ee ee
n
h







 
 
 
E
(3.28)
结合(3.25)~(3.28)得

1
22 2
11
00 0
22
412 32
1
222
11
000
44
nnnnn
h
ntLH
nnn
kk
keee
CkhCh k
Ch Ch
2
0
n
h
e







 


EE
fu
EEE

(3.29)
对上式从 1到n求和,当 h充分小使得 12Ch

和
或 时,有
00
0
h
pp0
hh
pQp0
22
00
1
12
4
0
0
n
nin
h
i
ni
i
kke
Ch Ch








EE
E
2
0
n
e
(3.30)
对上式用离散Gronwall 引理 4可得到

22
00
1
44
exp
n
nin
h
i
kke
ChChn Ch





EE 2
0
n
e
(3.31)
注意到 2
2
11
nn
ii
ii
aan



 及abab 由
上式


12
00
2
0
1
nnn n
hh h
nii
hh
i
Ske
kS Ch


 
  

uu
uu
若0
nnn
hh
eQpp

,则 ,则由(3.22)
有
n
h
Qp pn
h
12 122
0
nnnn
hh
kppkpQp Ch
0
(3.33)
若0
nnn
hh
eQpp

,则 22
00
0
nn
h
ee

。因
此存在常数 01


使得
00
n
h
ee

n
。于是,有


00 0
1nn nnn
h
ee eee


0
h


(3.34)
取

2
1


 ,由(3.32)和(3.34)及引理 3得
12 2
001
1
n
nnnnii
hhh
i
kpp kCh



uuuu (3.35)
把(3.33)和(3.35)与结合(2.12)得(3.23)。定理 5证
毕。
4. 数值实验
下面给出数值例子,验证非定常 Stokes 的时间二
阶精度的稳定化 CN 有限体积元格式比时间一阶精度
格式更优越,从而表明用稳定化 CN有限体积元格式
去求解非定常Stokes 方程的数值解是更有效可行。
计算域取为











1, 20,10,11, 0 
(如图 2所示),源函数

0f,Re = 1000,2T

,初
始函数和边值函数取为:左入口

,,
x
yt

=

(x, y) =




0, 1yy,右出口为




0,2,, ,t xy

xy ,其
他的边界的速度均 。 0
将计算域











1, 20,10,11, 0 剖分
为边长 0.01xy

 小正方形,然后在同一方向连结
其对角线将每个小正方形剖分成为两个小三角形构
成20.01h的三角形剖分 。对偶剖分取为重心
h

e
(3.32) Figure 2. Computational field
图2. 计算域
Copyright © 2013 Hanspub 31
非定常 Stokes 方程的稳定化 CN 有限体积元格式
型对偶剖分 。时间步长
h

0.01k

。
用[15]中的时间一阶精度的有限体积元格式算出
t = 2 时的速度和压力分别解画在图3和图 5中。
利用时间二阶精度的稳定化CN 有限体积元格式
问题 3算出 t = 2时的速度和压力分别解画在图 4和图
6。
Figure 3. Time first-order accuracy finite volume element solution
of velocity at t = 2
图3. 在t = 2处的时间一阶精度的速度有限体积元解
Figure 4. Time second-order accuracy CN finite volume element
solution of velocity at t = 2
图4. 在t = 2处时间二阶精度的速度 CN有限体积元解
Figure 5. Time first-order accuracy finite volume element solution
of pressure at t = 2
图5. 在t = 2处的时间一阶精度的压力有限体积元解
Figure 6. Time second-order accuracy CN finite volume element
solution of pressure at t = 2
图6. 在t = 2处时间二阶精度的压力 CN有限体积元解
比较图 3和图 4以及图5和图6可看出,时间二
阶精度的稳定化 CN有限体积元格式的速度(见图4)
和压力(见图6)的数值结果明显比时间一阶精度的有
限体积元格式的速度(见图 3)和压力(见图 5)的数值结
果要好。这也说明了用稳定化 CN有限体积元格式问
题3求解非定常 Stokes 方程的数值解更有效可行。
在用时间一阶精度的有限体积元格式[15]求非定
常Stokes 方程的数值解时,为了得到 的精度,
时间步长取为

4
10O

4
10

。这样,当用时间一阶精度的有限
体积元格式[15]计算 t = 2 的解时,要计算 步。而
用稳定化的 CN有限体积元降格式问题 3只要计算
200步,计算步数减少 100 倍,从而能减少计算过程
中截断误差的积累。因此利用稳定化 CN 有限体积元
格式问题 3求解非定常 Stokes 方程的数值解可以极大
地减少计算量,从而能减少计算过程中截断误差的积
累。所以非定常Stokes 方程的稳定化 CN 有限体积元
格式问题 3比时间一阶精度的有限体积元格式[15]更优
越。
4
210
5. 结论和讨论
本文建立了二维非定常 Stokes 方程的时间二阶精
度的稳定化 CN 有限体积元格式、分析了稳定化 CN
有限体积元解的误差。并用数值实验了二维非定常
Stokes 方程的时间二阶精度的稳定化 CN 有限体积元
格式在求解二维非定常Stokes 方程的数值解时,比时
间一阶精度的有限体积元格式 更有效。
]15[
特别是,利用时间二阶精度的稳定化 CN 有限体
积元格式问题3求解非定常 Stokes 方程的数值解时,
至少要比时间一阶精度的有限体积元降格式节省100
Copyright © 2013 Hanspub
32
非定常 Stokes 方程的稳定化 CN 有限体积元格式
Copyright © 2013 Hanspub 33
[7] E. Suli. Convergence of finite volume schemes for Poisson’s
equation on no-nuniform meshes. SIAM Journal on Numerical
Analysis, 1991, 28 (5): 1419-1430.
倍的计算量(当 时),从而能减少计算过程中截
断误差的积累,提高数值解的计算精度。因此,非定
常Stokes 方程的时间二阶精度的稳定化CN 有限体积
元格式问题 3比时间一阶精度的有限体积元格式[15]
更优越。
0.01k
[8] W. P. Jones, K. R. Menziest. Analysis of the cell-centered finite
volume method for the diffusion equation. Journal of Computa-
tional Physics, 2000, 165: 45-68.
[9] R. E. Bank, D. J. Rose. Some error estimates for the box methods.
SIAM Journal on Numerical Analysis, 1987, 24(4): 777-787.
[10] R. H. Li, Z. Y. Chen and W. Wu. Generalized difference methods
for differential equations-numerical analysis of finite volume
methods. New York: Marcel Dekker Inc., 2000.
[11] P. Chatzipantelidis, R. D. Lazarrov and V. Thomee. Error esti-
mates for a finite volume element method for parabolic equa-
tions in convex in polygonal domains. Numerical Methods for
Partial Differential Equations, 2004, 20: 650-674.
参考文献 (References)
[1] V. Girault, P. A. Raviart. Finite element methods for Navier-
Stokes equations: Theory and algorithms. Berlin Heidelberg:
Springer-Verlag, 1986. [12] X. Ye. On the relation between finite volume and finite element
methods applied to the Stokes equations. Numerical Methods for
Partial Differential Equations, 2001, 17: 440-453.
[2] J. G. Heywood, R. Rannacher. Finite element approximation of
the non-stationary Navier-Stokes problem part IV: Error analysis
for second-order time discretization. SIAM Journal on Numeri-
cal Analysis, 1990, 27(2): 353-384.
[13] M. Yang, H. L. Song. A post processing finite volume element
method for time-dependent Stokes equations. Applied Numerical
Mathematics, 2009, 59(8): 1922-1932.
[3] F. Brezzi, J. Douglas, Jr. Stabilized mixed method for the Stokes
problem. Numerische Mathematik, 1988, 53(1-2): 225-235. [14] J. Li, Z. X. Chen. A new stabilized finite volume method for the
stationary Stokes equations. Advances in Computational Mathe-
matics, 2009, 30(2): 141-152.
[4] F. Brezzi, M. Fortin. Mixed and hybrid finite element methods.
New York: Springer-Verlag, 1991. [15] 安静, 孙萍, 罗振东, 黄晓鸣. 非定常 Stokes方程的稳定化全
离散有限体积元格式[J]. 计算数学, 2011, 33(2): 213-224.
[5] 罗振东. 混合有限元法基础及其应用[M]. 北京: 科学出版社,
2006. [16] R. A. Adams. Sobolev Space [M]. New York: Academic Press,
1975.
[6] Z. Cai, S. McCormick. On the accuracy of the finite volume ele-
ment method for diffusion equations on composite grid. SIAM
Journal on Numerical Analysis, 1990, 27(3): 636-655.

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