空腔回流区水沙特性的计算分析(董耀华)

热度130票 浏览67次 【共0条评论】【我要评论 时间:1999年4月01日 13:47

空腔回流区沙特性的计算分析*水利论文5B!k2Djc,qv|~

董耀华水利论文aU9W1CbK7j

长江科学院)水利论文.Q2l.k9W(^8Vy

摘 要建立了回流区水深平均平面二维水流泥沙数学模型,提出了平衡输沙模式下回流区泥沙淤积计算方法;采用长江科学院空腔回流泥沙淤积水槽试验数据,验证了回流区泥沙淤积计算方法的可行性和精度;通过数模计算分析,一方面加深了空腔回流区环流结构及含沙量分布的机理认识,另一方面研究了模型变率对回流区水沙特性的影响。

@3Y \-HC;i}&I0

关键词空腔回流区 二维水沙数模 环流结构 泥沙淤积 模型变率 紊流模型

Cqg*A0j;\}*r0

国家自然科学基金委员会和中国长江三峡工程开发总公司联合资助的自然科学基金重大项目(59493600-1-3)部分研究内容,长江科学院潘庆教高和余文畴教高审阅了本文。

w8IV6D F-qGV c9|0

1概述

c3E0e+Cqtt'FS0

  河道水流由纵向水流(主流)和次生流(副流)构成。当次生流与纵向水流结合在一起时,将形成螺旋流,弯道水流属于典型的螺旋流;当次生流基本上脱离纵向水流形成相对封闭的环流时,即产生了回流。回流流动形态上的封闭性是相对的,回流区内的水流与泥沙仍然与纵向主流区的水流和泥沙不断地进行着质量、动量及能量的交换,这种交换主要是通过回流区交界面(即紊动掺混区)进行的;回流区水流及泥沙的运动呈现出强烈的三维运动特性,如:竖轴环流(天然河道中最为常见的回流形态)不仅存在着平面上的环流运动,还存在着上下层水流和泥沙之间的强烈掺混;挟沙水流河道或水槽中的回流区既可能发生泥沙冲刷,也可能出现泥沙淤积,这主要取决于回流区环流强度、河床组成以及纵向水流含沙量等因素,天然河道中的回流区以泥沙淤积为主[3~5]水利论文%qk0K&KB W D*Y

  回流区环流结构和泥沙淤积特性的研究不仅具有理论意义,而且对于三峡工程泥沙问题的研究具有实际应用价值。三峡水库变动回水区河段呈现出宽窄相间、浅滩交错的河道形态,存在着许多港池回流区,如青岩子河段的沙湾,丝瓜碛河段的五羊背,铜锣峡河段的郭家沱、唐家沱和明月沱等;盲肠河段或渠段口门回流区的泥沙淤积也是三峡工程重要泥沙问题之一,已建的葛洲坝工程三江上下游引航道以及在建的三峡工程五级船闸引航道均存在着口门回流泥沙淤积及对策研究问题;两江汇合和分汊河道汇流区也会产生小及强度不等的回流区,它们对于河口拦门沙的形成及分汊河道的演变等具有重要作用;较为典型的回流区还包括汛期泄洪时坝前电厂两侧回流,丁坝、围滩及取排水口等工程局部区域回流等。水利论文Dt:oK$yO*wNn,T

  采用数学模型模拟回流区环流结构和计算回流区泥沙淤积是一种有效的方法之一。对于枢纽及河道泥沙实体模型,一方面模型设计以重力相似作为水流动力相似的首要条件,而让水流紊动粘性的相似产生偏离;另一方面部分模型采用了变态模型(特别是水库变动回水区模型和河口模型);加之,模型比尺一般较小,水流和泥沙测量等模型技术受到限制,因此,实体模型对于回流区的模拟难以做到相似和准确,均会出现不同程度的偏离,例如,三峡库尾河段变态模型的港池回流区泥沙淤积量和淤积高程普遍比原型偏高。数学模型则不受模型相似、变率、比尺及测量等方面的限制。水利论文:k/wXi.t+P!_r#V

  回流区泥沙淤积计算严格地讲应该采用三维水流泥沙数学模型,但对于天然河道中常见的“封闭”竖轴环流,当水深相对较小、上下层水体掺混相对较弱时,采用水深平均二维水流泥沙数学模型进行模拟计算是可行的。由于回流区紊动掺混较为强烈,水流运动控制方程(N~S方程)中的对流惯性项与紊动粘性项均是主要控制项,因此,数学模型应该采用紊流模型(如κ~ε方程)进行紊动粘性系数νt和紊动粘性项的精细模拟。

9eW]+n/W9sN&?0

2 二维水流泥沙数学模型水利论文"_3u'm#i4SN

2.1 水流模型

gY7G;Nuy8\9k0

  水流控制方程采用水深平均(沿水深积分)平面二维非恒定流浅水方程(守恒型)[2,5,6]

6{or&SK0

990206e1.gif (1318 bytes)水利论文J2Sup[Jy1dmQ"W

(1)水利论文g8g5KH ~YY+vnr

990206e2.gif (2267 bytes)水利论文'u F1[$?0yI2A}GJ

(2)

Tv`0j)P.j7x9F0

990206e3.gif (2248 bytes)水利论文&l/mP:gL-dQ

(3)

'Jxn~z"VA-R8z0

式中X、Y、t分别为平面二维坐标和时间坐标;u、ν分别为水深平均流速分量;z、h分别为水位和水深;n、νt分别为糙率和率动粘性系数。水利论文Z2kD[ yZ#A'P

  紊动粘性系数νt的精细模拟采用了水深平均平面二维κ-ε双方程模型(守恒型)[2,6]

6B/_7ex BT0

990206e4.gif (2165 bytes)水利论文1p1Be6o3i(us6A5q

(4)水利论文,JYC\r*Tn3U&xW

990206e5.gif (2212 bytes)

$\ [;npv){ I L I0

(5)水利论文^#LUu5opI2EIMd

990206e6.gif (1012 bytes)

k1d G|`6N4U0

(6)水利论文&e+zQ9j\ \1c

其中水利论文.N4N[z;`U&O;t@%w6D

990206e7.gif (1843 bytes)

#Y+Z&vC2`6e7f R,q0

(7)水利论文#F ]G K*A X#|&Ve

990206e8.gif (1231 bytes)

:j(Y'`Y*c }P0

(8)水利论文 [q:D0C;x#{Jq

式中κ、ε分别为紊动动能和耗散率;Ph为κ的产生率;Pκν、Pεν分别为附加源项;U*为摩阻流速;σκ、σε、C1、C2、Cμ、Cκ、Cε均为模型参数。

0K Z,LP/r0

对于恒定流,去掉方程中的时间惯性项;水流模型定解边界条件采用进口给定流速u、v、紊动动能κ、耗散率ε和紊动粘性系数νt的分布,出口按法向边界条件控制;模型进口各变量具体按下列经验公式给出

a+j U8aG4L CXm0

990206e9.gif (1106 bytes)水利论文yvV0CPUoV8vz

(9)水利论文!RJ3[7[,Y2[

990206e10.gif (1352 bytes)

G4}*l+?q:@0

(10)

+WO|"?&L_-}AM0

方程(1)~(3)数值求解基于SIMPLER算法[1,5,6];方程(4)~(6)数值求解采用有限控制容积法[1,6]水利论文:g;`!a*ku

  水流模型中主要参数的取值见表1

e m;a5yH5^"g-V0

表1 水流模型主要参数取值水利论文|f,Pd0i

Main parameters of flow model水利论文 pw,c6W]Z"d\7}


σκ水利论文Ki J8y |G2iQ

σε

0M:L3Vo l!s8VM#bd0

C1水利论文;hdOE9C8y

C2

S#D+T:Z3f`.foQ1u j0

Cμ水利论文6f?7qo*H X8| r

E水利论文)S1y|8UB%Q%R

Eκ水利论文l-H'K}n


1.0

j \5k4E*B]0

1.3

"rB Ztu2kO0

1.44水利论文)`UR!z*w8V;I

1.92水利论文!z'yy3?fW;B7J E2[

0.09

KmKyI3OG0

1/2~1/8水利论文.fdX$E@K

0.05~0.10

+k%M)`KeR'i)I0

2.2 泥沙模型水利论文L0D5NC{Gn

  含沙量分布的模拟采用了平衡输沙模式下水深平均形式的泥沙对流扩散方程[5,7]水利论文TZ7ph2_ \ @M F:JQ

990206e11.gif (1856 bytes)水利论文Z8]z.~Pn\@

(11)

,C:z+~7h+f!Na&l0

式中s为垂线平均含沙量;εs为泥沙紊动扩散系数(εs=κνt,h为经验系数)。水利论文R6G5N(a(A

  含沙量方程(11)边界条件仍采用进口给定分布,出口按法向边界条件控制;数值求解也采用了有限控制容积法[1,6]

TXH$j~H[ Ac,`2h0
  假定回流区外侧的纵向主流泥沙冲淤变化相对较小,采用平衡输沙模式可计算出整个回流区内的泥沙淤积量,计算公式为

990206e12.gif (1141 bytes)

O0RB8AF\Q0

(12)水利论文a"Ac s'P3ty

式中W为回流区泥沙淤积体积;Q为流量;Δt为计算时段;G1s为回流区上游进口断面的输沙率;G2s为回流区下游出口断面的输沙率;γ′为泥沙淤积物干容重。水利论文o0H2Bmm$fc+b

3 空腔回流区水沙特性计算分析水利论文"bUf`jIG

  利用长江科学院空腔回流泥沙淤积水槽试验资料,进行了同一水槽边界四种模型变率的水流流场、含沙量分布及回流区泥沙淤积计算。四组水槽的计算水沙条件、主要参数取值见表2。计算采用非均匀正交矩形网格,网格间距为5~10cm,计算节点为60×40(图1)。水利论文m {}7y2N+C'r"GQ r4V:d

990206t1.gif (3912 bytes)

`,p Qg[0

图1 空腔回流区二维水沙计算网格水利论文va,kLbp

Grid in 2-D model in cavity recirculation region

)A,M\-r E:iQ j0

990206t2.gif (4951 bytes)

m,A Ug5^ kh0

图2空腔回流区二维水沙计算流速矢量场

o @t^'P)B Q/q0k4E0

Vector field of calculated velocity

bekr7o c0

  图2和图3给出了第一和第四组水槽回流区流速矢量场和相对流速等值线分布,相对流速是指回流区二维计算节点水深平均流速占水槽进口断面平均流速的百分数;图4给出了四组水槽计算的无量纲化主槽流量沿程变化,无量纲化主槽流量为空腔主槽区流量与进口流量的比值。水利论文5d7m~:kC:vR9}+l(H

表2 空腔回流区二维水沙计算参数及成果水利论文w7e2v$f v

Calculated parameters and results of 2-D model

T]g3Hm#D0

计算组次水利论文eOd,?VX5NN

水利论文#C1n0@1i c,}

7bv-v U?w[ u0

}#z7ATe$[5Q$U_I `7l0

水利论文0QbN @B'B$_kt


水深h(cm)

VX p!\^"U6glw#o0

10

f6kQgWy0

25水利论文'u9g:ZS-dHH|Vt

50水利论文#x+WAH l @

70

-dJ1~)NlZ&u0

流量Q(l/s)

CD%loCR@;LYme0

6.17

oD*{;I!m h {*x0

24.38

a(QgYz*~j/z D0

68.96水利论文/G!W:rI+]6fft

114.23

2FX`!S,E']/O#h!MS%A0

含沙量S(kg/m3)水利论文 Lx3}!@ U)y t4g

7.13水利论文(Wd;]CFe9~"rQ2}

7.18

7XN-~8D ` d B(l0

7.18

jf\/?c+V0

7.43水利论文K5]$t!WMN D,{1i#T

泥沙中值粒径d50(mm)水利论文:d9`"gZ0l `&XT P

0.09水利论文@T:{s*Dd _"D)y!t3\ K*a

0.17

%Me7C CzM0

0.30水利论文B r;]My\#\r

0.35

.O(vWzM.PCSV4tI.BT0

放水历时t(min)

`q0ZV#i1f shjS0

360水利论文WL Gr_&w-zT

227水利论文0]({h8e pV

161

@_3BI.b2@6k#Ykm0

136

]s4tc+KIL0

进口流速U(cm/s)水利论文"R)MFq\5y

4.113水利论文S"v*Me0x7^:~B-Z#AE|

6.501水利论文v&]B,?nH!WO?

9.195水利论文r-S n4CJ

10.879

,M5RI0YQ0

模型变率e

,M9C X)V%E[,V]0

1.0

ng8mvt1V1G6XUMw0

2.5水利论文j(AA LX

5.0水利论文0^A,Zfe"RAB1M1C

7.0水利论文^D BD0j+R!X I(E

Froude数Fr

g5Nb*X"f(ZN\0

0.0415水利论文 EVK(Zuy

0.0415水利论文`Z+_#i @swM ~

0.0415

:v[(f?&mS0

0.0415水利论文:@Rz9D;Z4S:a2S2Y Y

Reynolds数Re

q9V4mpu0

4110水利论文r#x+|De

16250

vL{2F,o(Y2?%R0

45975

/b eDD e"D"Nc0

76150水利论文 \ ?N] A h h

水槽糙率n水利论文+hB J(XT(d

0.012

fT`.iOF I0

0.012水利论文E*dA4d(MT"O

0.012水利论文?N5\#V]$zo7v X){

0.012

Z#ktUfTn?Kj0

实测回流区淤积水利论文m8p4VNYP&W%t]

2.30水利论文:}8F{'@e-Az-EG$x

3.22

6UnG2nHo0

8.74水利论文 j)B.`9~4I Q

13.34水利论文'G$zZKJ*E I3PH@y

计算回流区淤积水利论文S MgS'f"}.o;H

3.66

E*~h!Ro3BJ#h hB:H0

6.27水利论文 z&?/H;c,_E \O

8.88水利论文L#K_j ~!_ES|I

10.64水利论文rk8?H.H^ c J


注: ① 无量纲化泥沙淤积:100×(泥沙淤积体体积/回流区水体体积)。水利论文7_Q"nmp*f3by(K

990206t3.gif (6487 bytes)水利论文aXE~T*K3\&`H

图3空腔回流区二维水沙计算无量纲化流速等值线水利论文k0o@F)inH

Calculated isopleth of velocity

#}}*GIn z%G^0

   表2给出了四组水槽回流区泥沙淤积验证结果;图5给出了第一和第四组水槽回流区含沙量等值线;图6给出了二维数模计算的空腔回流区泥沙淤积量与主槽流量亏损的关系,泥沙相对淤积量和主槽流量相对亏损分别指不同模型变率条件下无量纲化泥沙淤积和无量纲化主槽流量亏损与模型变率1时(第一组水槽) 对应值之比,无量纲化主槽流量亏损为进入空腔回流区最大流量占进口流量的百分比;图7给出了计算的空腔回流区泥沙相对淤积量与模型变率的关系。

?0XTeD{0

990206t4.gif (2872 bytes)水利论文CN6NvWJ3t

图4空腔回流区计算无量纲化主槽流量沿程变化

O.kh5~Z'y2D0Jc0

Calculated discharge along main channel

*m2J wp,e%YdM$p0

990206t5.gif (5641 bytes)水利论文 g.e.pt G#x*zm2p2q

图5空腔回流区二维水沙计算含沙量等值线水利论文 {0~K-k'@p jq

Calculated isopleth of sediment concentration

v*rc6zE.i0G0

990206t6.gif (2246 bytes)

'o5KB/p _XVne0

990206t7.gif (2231 bytes)水利论文R BMzCy{^"d_ V

图6空腔回流区泥沙淤积量与主槽流量亏损的关系

$QP |ensvo0

图7空腔回流区泥沙淤积量与模型变率的关系

P|H$v5wW0

Relation between sediment deposition and discharge loss in main channel水利论文8mQG{6iB(n!T }

Relation between deposition and distorted ratio of model水利论文m p7@3H2]y U

  上述计算成果不仅证明了回流区二维水流泥沙数学模型的可行性,验证了回流区泥沙淤积计算方法的精度,而且对于回流区环流结构、泥沙淤积以及模型变率对回流区水沙特性的影响等有如下几点认识:① 空腔回流属于竖轴环流,具有不对称性和不封闭性;回流区与主流区交界面(紊动掺混区)存在较大的流速梯度;水流从空腔下游侧壁进入回流区,以对流形式参与回流区环流运动;回流区环流运动受水流的紊动粘性项和对流惯性项的共同作用。② 回流区含沙量分布也呈现出不对称性;回流与主流交界面存在着较大的含沙量梯度,但其强度弱于流速梯度强度,表明回流区泥沙运动存在着滞后性;回流区泥沙淤积量与主槽流量亏损的关系表明:回流区环流特性决定了回流区泥沙特性及淤积量。③ 模型变率不同的空腔水槽,回流区环流形态及结构非常相似,但存在着环流中心及强度、进入空腔流量以及回流区水流紊动特性等方面的差异;模型变率不同的空腔水槽,回流区含沙量分布形态也非常相似,但回流区含沙量大小及泥沙淤积量差异较大;模型变率越大,回流区含沙量及水流挟沙力越小,泥沙淤积量越大;回流区泥沙淤积量与模型变率趋于线性关系,表明:变态模型均会出现回流区泥沙淤积量的不相似。

*dteL1\4D~0

4 结语

!P"m)FneuE%IO0

  1. 本文建立了回流区水深平均平面二维水流泥沙数学模型,提出了平衡输沙模式下回流区泥沙淤积计算方法;采用长江科学院空腔回流泥沙淤积水槽试验数据,验证了回流区泥沙淤积计算方法的可行性和精度;这一计算模式适用于水槽或天然河道最为常见的竖轴环流回流区环流结构及泥沙淤积计算。

2i7k:x+z\I/[F0

  2. 二维数模计算分析表明:空腔回流具有不对称性和不封闭性;回流区与主流区交界面存在着较大的流速和含沙量梯度,含沙量梯度强度弱于流速梯度强度;部分挟沙水流从空腔下游侧壁进入回流区,以对流形式参与回流区环流运动及泥沙淤积,回流区环流结构及泥沙运动受到水沙紊动粘性项和对流惯性项的共同作用;回流区环流特性决定了回流区泥沙特性及淤积量。

$j,F+|9?hW/r MA0

  3. 二维数模计算分析还表明:模型变率不同的空腔水槽,回流区环流结构及含沙量分布形态非常相似,但存在着环流中心及强度、进入空腔流量、回流区水流紊动特性以及回流区含沙量大小及泥沙淤积量等方面的差异;回流区泥沙淤积量与模型变率趋于线性关系,表明:变态模型均会出现回流区泥沙淤积量的不相似;因此,对于重力相似的变态模型,变率的选择应考虑回流区泥沙淤积相似的偏离,并应对回流区泥沙淤积量进行修正。

~D5e rqK~0

参考文献

}'g4sr.gD;N@o6K*b4{0

1 [美]S.V.帕坦卡著,张政译。传热与流体流动的数值计算。北京:科学出版社, 1984.

{.Z8Fp?%We9\0

2 [西德]W.罗迪著, 贺益英译。环境问题的紊流模型。北京: 水利电力出版社, 1987.

7fg:pdW7x0

3 张瑞瑾。论环流结构与河道演变的关系。武汉水利电力学院学报。1963(2):1-15.

/t9R vRC#C tAO0

4 刘青泉。盲肠河段口门掺混区的泥沙扩散。泥沙研究。1995(2):11-21.

}7luH-@)w$X0

5 董耀华。弯道水流的基本特性及数值模拟。长江科学院院报。1996(1):25-29.水利论文4@ui2\%Dy!^

6 董耀华。河道温度及浓度场的平面二维数模研究。长江科学院院报。1997(2):22-26.水利论文(HI n&g)Bjv

 

Xi1WqbG0
TAG: 沙特 董耀华 回流区 空腔
顶:13 踩:18
【已经有100人表态】
18票
极差
11票
很差
16票
较差
5票
稍差
12票
稍好
6票
较好
21票
很好
11票
极好
下一篇:黄河下游纵剖面自调整特性(尹国康)
上一篇:剖面二维非恒定悬移质泥沙扩散方程的数值方法(张耀新 吴卫民)
查看全部回复【已有0位网友发表了看法】