非正交曲线坐标下二维水流计算的 SIMPLEC

　2003年2月

水 利 学 报SH UI LI X UE BAO 第2期

收稿日期:2002201204

基金项目:国家杰出青年科学基金项目(50125924);高等学校博士学科点专项科研基金项目(2000014112);辽宁省自然科学基

金项目(2001101073)

作者简介:吴修广(1974-),男,山东阳谷人,博士生,从事环境水力学研究。文章编号:055929350(2003)022*******非正交曲线坐标下二维水流计算的SIMPLEC 算法

吴修广1,沈永明1,郑永红2,王平义

3(11大连理工大学海岸和近海工程国家重点实验室　辽宁大连　116023;21中国科学院广州能源研究所,广东广州　510070;

31重庆交通学院河海工程系　重庆　400074)摘要:本文采用Laplace 方程坐标变换方法生成正交曲线网格,并对浅水流动的控制方程进行坐标变换,方程离散时采用B 型交错网格。利用“水位扫描法”结合壁面函数法来处理移动边界,用SI MP LEC 算法解非正交曲线坐标下的k -ε双方程紊流模型,修正了由网格的非正交性引起的误差。通过对美国C olorado 洲Fall River 的资料进行流场验证,计算结果与实测资料基本符合,显示了本模型在不规则水域计算中的实用价值。

关键词:坐标变换;k -ε紊流模型;水位扫描法;壁面函数;SI MP LEC 算法

中图分类号:T V13114文献标识码:A

随着经济发展和社会进步,水利工程建设的步伐也在进一步加快,其中港航建设、大坝建设中的泥沙问题以及近来倍受世人关注的水污染问题已经成为制约水利发展的瓶颈问题,弄清河流、湖泊、海洋中水动力因素,是解决以上问题的重要基础。近年来,数学模型已逐步取代物理模型实验成为研究水流的重要手段,而浅水流动模型是处理大区域流场的一种非常有效的模型。它属于非线性方程组,在目前只能用数值方法求解,因此,有必要研究一种简单、高效的方法来求解浅水流动问题。自Patankar 和S palding [1]

发展了SI MP LE 算法以来,该方法被广泛应用于不可压缩流场的数值模拟,而且该方法还得到了进一步的发展,主要有SI MP LER 算法

[2]、SI MP LEC 算法[3]、SI MP LEX 算法[4]和SI M 2P LET 算法[5]等。这些模型均成功地应用于速度—压力耦合的流场计算,深度平均的浅水流动模型是在静压假定下导出的,一般流体模型中的速度—压力耦合也就转换成浅水流动模型中的速度—水深耦合[6]。

天然河流、海湾的边界曲折、地形复杂,采用坐标变换是解决问题的途径之一。目前多数N -S 方程的坐标变换中,流程全部采用逆变分量,这样就增加了方程的复杂程度。于是忽略掉方程中的非正交项,利用正交变换下的方程进行数值求解[7,8]。对于具有复杂边界的海湾及弯曲的河流,坐标变换中很难保证每个点都正交,特别是边界附近。水位变化是水力计算中难点之一,在目前的紊流数学模型中,多简单的利用“冻结法”,这样做将失去对边界出流动模拟的准确性。

本文研究中,采用正交曲线坐标变换生成数值网格,而数值计算中采用非正交曲线坐标下的k -ε双方程紊流模型,这样可以自动修正网格生成中的非正交项。流速除对流项中采用逆变分量,在其余各项中均采用原始分量,这样使得方程书写简单,有利于将各方程写成通用形式,编写的程序变得更规范。作者受Jian Y e 同位网格[9]

的启发,对普通交错网格做了修改,即采用B 型交错网格,使得u ,v ,k ,ε的计算布置在一个节点上,有利于节省计算程序代码,使程序书写更加规范。引入动边界扫描技术,结合紊流模型的壁面函数法,使壁面随着真实边界而变化。数值求解时,采用控制体积法离散方程,运用SI MP LEC 算法,使计算的流场更符合实际流场。

1　数值网格

本文对计算区域用Laplace方程实施坐标变换,生成正交的贴体网格,控制方程:

52x 5ξ2+52x

5η2=0

52y 5ξ2+52y

5η2=0

(1)

方程组(1)可采用S OR方法求解,然后用正交曲线网格边界正交化处理的三次样条插值边界滑动法[10],处理其边界处的网格,以提高其正交性。与其它方法相比,此方法虽在一定程度上提高了边界处网格的正交性,但难以保证完全正交,如果再把曲线坐标系下N-S方程中的非正交项去掉,由斜交网格引起的计算误差是不容忽视的。为了修正此项误差,本文采用非正交曲线坐标下的平面二维k-ε双方程紊流模型,模型能自动修正网格生成中的非正交项。

2　数学模型

笛卡儿坐标下深度平均的k-ε双方程紊流模型的通用微分方程:

5

5t(HρΦ)+5

5x(HρuΦ)+

5

5y(HρvΦ)=

5

5x(HΓΦ

5Φ

5x)+

5

5y(HΓΦ

5Φ

5y)+SΦ(2)

方程(2)转换到曲线坐标(ξ、η)下,仅在对流项中使用流速的逆变分量,而在其它项中使用原始变量,这样既简化了方程,又使所有方程仍可写为曲线坐标下的通用方程,模型的微分方程可写为如下通用形式:

5

5t(HρΦ)+1

J

5

5ξ(HρUΦ)+

1

J

5

5η(HρVΦ)=

1

J

5

5ξ(

HΓΦ

J

(αΦξ-βΦη)+

1

J

5

5η(

HΓΦ

J

(-βΦξ+γΦη))+SΦ(ξ,η)(3)

式中:Φ为所求问题的因变量;U和V分别为直角坐标下流速u和v的逆变分量,仅在对流项中出现;ΓΦ为扩散系数;SΦ为源项。当Φ表示某一特定量时,ΓΦ和SΦ对此特定量有特定的意义和表达式,这时方程(3)亦赋予特定的意义,模型的控制方程组如表1所示。

表1　模型控制方程组

方程ΦΓΦSΦ

连续100

ξ———动量uμ

eff

S u

η———动量νμ

eff

Sν

紊流动能k μ

eff

σ

k

H(G0k+G kν-ρε)

紊流动能耗率εμ

eff

σεH

ε

k

(Cε

1

G0k-Cε

2

ρε)+Gεν

S u=-1

J

ρgΗ(y

η

5z s

5ξ-yξ

5z s

5η)+

1

J

5

5ξ

Hμeff

J

y2η

5u

5ξ-yξyη

5u

5η

+1

J

5

5ξ

Hμeff

J

yξxη

5v

5η-xηyη

5v

5ξ+

1

J

5

5η

Hμeff

J

xξyη

5v

5ξ-xξyξ

5v

5η

+1

J

5

5η

Hμeff

J

y2ξ

5u

5η-yξyη

5u

5ξ-τbx

S V=-1

J

ρgΗ(x

ξ

5z s

5η-xη

5z s

5ξ)+

1

J

5

5ξ

Hμeff

J

x2η

5v

5ξ-xξxη

5v

5η