文章编号:
1674-7070(2010
)05-0477-04
基于有限体积
HWENO
格式的一维溃坝流模拟
摘要
基于浅水波方程组建立一维渍坝
JA
模型,并给出数值模拟结果.其中,空间
离散采用
HWENO
(Hermit Weighted
Es-
盼
ntially
Non-Oscillat01Y)
格式,时问离散
采用四步
TVD
(Total
飞
V
怆咄
a
缸
nat
m
喀
g
引)
Rung
悖
e
←
ι-
1<
占
K
豆
Qω~l
凶川
u
川
1
川
tta
方法,模拟堤坝渍决时
洪
7
水扎
7
演突进过程.模拟结果表明:较采用
WENO
格式所得数值解更精确;同时,相
比
WENO
格式的相应算法,该算法解决
一维渍坝流问题能更有效地减弱振荡,
对间断具有更高的分辨率.
关键词
;责坝流;有限体积法;
HWENO
格式
中图分类号。
175.2
文献标志码
A
收稿日期
20
1O
'()6-02
资助项目同家青年基金
(40906048)
作者简介
翁磊,另,硕
11
二,主要研究水矿
UJ
数值慎
J
以
.1
we
吨
@163.
com
1
1
判桌{,
'f
,~.、
n
早大学数型学院,
1
丰
j
京,
210044
翁磊
l
程国胜
l
卢长娜
l
。
引言
n
o
--
t
c
u
lu
o
I
t
n
y--
拦河筑坝,既可防洪又可用于灌溉、发电,对社会发展起积极作
用,但坝体一旦愤决,形成的洪水波会对下游地区造成严重灾害.洪
水波的速度影响着下游的预警和紧急疏散时间,决定受灾范围和程
度,因此,愤坝流数值模拟具有非常重要的研究意义.
、渍坝流常以浅水波方程组作为模型.愤坝流具有间断性,对浅水
波方程组离散时,要求数值方法具有稳定、高分辨率、能较好抑制振
荡等特性.近年来,用高分辨率数值格式建立一维溃坝流模型已有一
些研究.如
:Choi
等
lJ
基于迎风
TVD
,迎风
TVB
和
ENO
格式讨论了一
维溃坝问题并指出
TVB
格式精度较高,但
CPU
时间消耗较多
;Lin
等
2J
提出了一种二阶混合
TVD
有限体积格式,由于
TVD
格式本身固
有的限制,精度至多只能达到二阶,极值点附近降至一阶川;饱远林
等到基于
Boltzmann
方程,采用有限体积
KFVS
方法,成功模拟了下游
为干河道的一维溃坝流;
Mambretti
等川采用
M
c:
Corma
c:
k
-
Jameson
激波
捕捉法成功建立了一维带污染物两相流愤坝模型;文献
[6
J
应用
WE
NO
格式模拟了闸门瞬间开启后的水流,水位、流速与解析解拟合较
好;文献
[7J
建立了基于
WENO
格式的浅水间断流数值模型,得到了
有五阶精度的一维溃坝模拟算法.文献
[8
J
将一种基于
HWENO(Her
mit
Weighted
Essentially
Non-Oscillatorγ)
格式的限制器引人双由守恒
律方程的求解中,
HWENO
格式在重构边界点值时,选用
Hermit
插值
多项式,其优点是选用较少的模板就能达到较高的精度,同时减少计
算量,边界处理也变得简洁.本文将
HWENO
格式应用于一维愤坝流
进行数值模拟,所得结果与基于
WENO
格式的数值结果
7J
及精确
解
9
J
比较表明,本文建立的溃坝流模型,格式稳定,分辨率高,能更好
地处理类似愤坝这类浅水间断流问题.
1
浅水波方程组及时空离散
Shallow
water
wave
equations
and
spa
c:
e-time
dis
c:
retization
1.
1
有限体积法积分方程
在笛卡尔直角坐标系下,一维浅水波方程组守恒形式如下:
立
U+
立
F
=
O.
()t
ðx
、、,
F'
l
/,
E
飞、
文章编号:
1674-7070(2010
)05-0477-04
基于有限体积
HWENO
格式的一维溃坝流模拟
摘要
基于浅水波方程组建立一维渍坝
JA
模型,并给出数值模拟结果.其中,空间
离散采用
HWENO
(Hermit Weighted
Es-
盼
ntially
Non-Oscillat01Y)
格式,时问离散
采用四步
TVD
(Total
飞
V
怆咄
a
缸
nat
m
喀
g
引)
Rung
悖
e
←
ι-
1<
占
K
豆
Qω~l
凶川
u
川
1
川
tta
方法,模拟堤坝渍决时
洪
7
水扎
7
演突进过程.模拟结果表明:较采用
WENO
格式所得数值解更精确;同时,相
比
WENO
格式的相应算法,该算法解决
一维渍坝流问题能更有效地减弱振荡,
对间断具有更高的分辨率.
关键词
;责坝流;有限体积法;
HWENO
格式
中图分类号。
175.2
文献标志码
A
收稿日期
20
1O
'()6-02
资助项目同家青年基金
(40906048)
作者简介
翁磊,另,硕
11
二,主要研究水矿
UJ
数值慎
J
以
.1
we
吨
@163.
com
1
1
判桌{,
'f
,~.、
n
早大学数型学院,
1
丰
j
京,
210044
翁磊
l
程国胜
l
卢长娜
l
。
引言
n
o
--
t
c
u
lu
o
I
t
n
y--
拦河筑坝,既可防洪又可用于灌溉、发电,对社会发展起积极作
用,但坝体一旦愤决,形成的洪水波会对下游地区造成严重灾害.洪
水波的速度影响着下游的预警和紧急疏散时间,决定受灾范围和程
度,因此,愤坝流数值模拟具有非常重要的研究意义.
、渍坝流常以浅水波方程组作为模型.愤坝流具有间断性,对浅水
波方程组离散时,要求数值方法具有稳定、高分辨率、能较好抑制振
荡等特性.近年来,用高分辨率数值格式建立一维溃坝流模型已有一
些研究.如
:Choi
等
lJ
基于迎风
TVD
,迎风
TVB
和
ENO
格式讨论了一
维溃坝问题并指出
TVB
格式精度较高,但
CPU
时间消耗较多
;Lin
等
2J
提出了一种二阶混合
TVD
有限体积格式,由于
TVD
格式本身固
有的限制,精度至多只能达到二阶,极值点附近降至一阶川;饱远林
等到基于
Boltzmann
方程,采用有限体积
KFVS
方法,成功模拟了下游
为干河道的一维溃坝流;
Mambretti
等川采用
M
c:
Corma
c:
k
-
Jameson
激波
捕捉法成功建立了一维带污染物两相流愤坝模型;文献
[6
J
应用
WE
NO
格式模拟了闸门瞬间开启后的水流,水位、流速与解析解拟合较
好;文献
[7J
建立了基于
WENO
格式的浅水间断流数值模型,得到了
有五阶精度的一维溃坝模拟算法.文献
[8
J
将一种基于
HWENO(Her
mit
Weighted
Essentially
Non-Oscillatorγ)
格式的限制器引人双由守恒
律方程的求解中,
HWENO
格式在重构边界点值时,选用
Hermit
插值
多项式,其优点是选用较少的模板就能达到较高的精度,同时减少计
算量,边界处理也变得简洁.本文将
HWENO
格式应用于一维愤坝流
进行数值模拟,所得结果与基于
WENO
格式的数值结果
7J
及精确
解
9
J
比较表明,本文建立的溃坝流模型,格式稳定,分辨率高,能更好
地处理类似愤坝这类浅水间断流问题.
1
浅水波方程组及时空离散
Shallow
water
wave
equations
and
spa
c:
e-time
dis
c:
retization
1.
1
有限体积法积分方程
在笛卡尔直角坐标系下,一维浅水波方程组守恒形式如下:
立
U+
立
F
=
O.
()t
ðx
、、,
F'
l
/,
E
飞、