高含沙水流运动特性综述(刘兆存,徐永年)

热度116票 浏览61次 【共0条评论】【我要评论 时间:2001年7月01日 09:49

含沙水流运动特性综述

"TqJl3FE$KUJ0

刘兆存徐永年水利论文kKT8z!~I
(中国水利水电科学研究泥沙研究所)水利论文QR,buL6V I

摘要:本文综述了高含沙水流的阻力特征、流速分布、紊动强度分布、能量分布等流动过程中所具有的特性,并和清水、低含沙水流作了对比分析。系统论述高、低含沙水流所共同具有的规律。水利论文;?YhIT(l&Ux.guP

关键词:高含沙水流;特性;低含沙水流基金项目国家自然科学基金资助项目(59890200)。

*hN_3P6\!f6A bW/T9cx0

作者简介:刘兆存(1970-),男,中国水科院博士。
p2VQ5_2]#t2Q#W0收稿日期:1999-03-01

*b#r&n_R3fZ%o0

1 引言水利论文FD y?W A3D3i

  和清水水流、低含沙水流运动特性相比,高含沙水流运动有其自身的规律[1~4]。本文作为低含沙水流运动特性的继续[5],重点探讨高含沙水流运动时的阻力特性、流速分布、脉动强度分布、能量分布等其运动过程中所固有的规律,并和低含沙水流、清水水流的相应特性做了对比分析。

cQ%V(d)B\O*TT0

2 流速分布

b1n n S(N0

  文献[1]研究了明渠高含沙水流的阻力规律,对于二元明渠实测表明,当含沙量较高时,典型的泥流层流流速分布如图1所示(u0为最大流速)。

7S~t(@1[1v+c)X8\0

  在高含沙水流中,实测证明[1]紊流边界层虽是多层结构,但不论是对过渡区或阻力平方区,实测证明充分发展的外流区均可用对数分布近似表达,其方程为(d50=0.042, d10=0.01,单位:mm,含沙量变化范围为0~650kg/m3)水利论文w.@vGr

水利论文 p^:{1Gl+{ E L

(a)层流流速分布

U W vdEMZ0
(b)有微动层的层流流速分布

图1 流速分布
A y/\ G'q7n;l0
Fig.1 Velocity distribution

q-p5P!u%Cxb nt0

u/u*=1/κlny/Δ+A水利论文k)^yG+E@0~j;]s

(1)

S$i Q:J&h_0rS'x J0
式中 u为距床面为y处的点流速,κ为卡门常数,Δ为床面突起高度。在低含沙水流中,类似Image781.gif (1112 bytes),令
Image782.gif (1187 bytes)

E1=u2*Image783.gif (860 bytes)/k2(1/y1-1/y2)水利论文%aGh!dJ"]J `!G3sg

(2)水利论文UU]?g!c2Bg9bR'ST

Image783.gif (860 bytes)为η的垂线平均值,文献[2]的研究表明,既可适用于高含沙宾汉流体,又可适用于清水和高含沙牛顿流体的恒定二维均匀明渠紊流的流速分布公式为

umax-u/水利论文5yR4s} q Kq4_K

Image784.gif (1049 bytes)=1/klnah/y

(3)

8T3g]/Q4z!M#I5F0

其中,a=y0/h即为相对非流核厚度。当水流为有流核的高含沙水流时

kx6O,mD7b;|0y0

umax=q/κ+a/κ水利论文2ku3m LP6O(rk0y

Image784.gif (1049 bytes)

(4)水利论文-yX6G0s ` LgN r

当水流为清水或无流核的高含沙水流时,τB=0,a=1,此时

a;^;_5x gM_U.{P0

umax=q/h+u*水利论文2Q.N)t1~%P c0F:i

(5)水利论文v7{\xS \U ]2At

其中q为单宽流量,卡门常数为含沙量的函数(κ的最小值对应于S=200kg/m3左右),关系如图2示。

K8PJ5t@ S)\7N-q UE0

  文献[3],[4]研究后认为,在高含沙水流中流速分布的统一表达式为(这里以均质、光滑边壁的宾汉流动为例)水利论文5\f2h Du

Image599.gif (854 bytes)/u*=Image785.gif (853 bytes)l(Image786.gif (854 bytes)/u*)lt(Image786.gif (854 bytes)/u*)水利论文4[9Z9vK^$O-m

(6)

V%YVB~ {0

9ZMD%m8Ex0

图2 k~S关系水利论文2fj~cu2u
Fig.2 k~S relation
水利论文 Jwv{-xOC%ko

其中Image786.gif (854 bytes)/u*示无量纲平均流速分布(u*为摩阻流速),Image785.gif (853 bytes)l、γt分别为层流、紊流发生率,(Image786.gif (854 bytes)/u*)l、(Image786.gif (854 bytes)/u*)t分别为层流、紊流时的流速分布,具体表达式和清水时类似,可参看[3],[4]。

  由此表达式知式(1)、(3)均为其特例。分析知,在紊流充分发展的情形下,宾汉流体和牛顿流体的流速分布表现出类似的分布规律,具有[5]述的类似特性。据文献[6]的研究,泥浆粗颗粒两相管流的流速分布为水利论文fD xfCO%Y ]&_

Image787.gif (3049 bytes)水利论文7Sc@P"W(d

(7)

其中系数Image42.gif (855 bytes)=0.8,卡门常数κ为κ=κ(Svf,Svc,u*),Δ为近底层流层厚度,ym为垂线最大流速点处的y值,Svf、Svc分别为二相流中细、粗参加了个远程的Microsoft认证辅导,考过了我就做soho颗粒泥沙浓度。顺便指出,在稀相流中,颗粒速度比液相的要小[7]。水利论文g"`.l we Hn

  在清水水流和低含沙水流中,层流时的流速分布研究得相对较为成熟,而紊流时的流速分布特性可参看[5]。在高含沙水流中,层流区流速分布和低含沙水流有明显的不同,而在紊流区,如上所述,二者流速分布有类似的特性。造成这种现象的原因是,高含沙水流的本构方程为:τ=τ0+ηdu/dy,τ0,η的实验结果表明,虽然高含沙水流具有触变特性,由τ、η的自身物理特性和随着所处的运动状态的演化而演化的特性,共同决定了高含沙水流流动时的流速分布。和清水及低含沙水流情形时相同,高含沙水流运动特性之一的流速分布是空间点坐标和时间的函数,在从层流向紊流过渡时,高含沙水流的转捩特性应和清水水流相同。在有细颗粒的高含沙水流运动中,在絮网的形成和破坏的动态平衡和随运动状态的演化而演化的过程中,一旦完成转捩形成充分发展的紊流状态时,流动的流体内部的流体质点间有和清水水流类似的流动关系和制约形式。换言之,流体质点运动形式随流体运动状态(不同的雷诺数)的演化有类似的演化机制。从能量观点看,流速分布是流体运动过程中能量分布形式的体现,由于动量传递属于一种扩散现象,流速的连续性本质上是依靠内界面的内摩擦阻力来维持的,从动力学的观点上看,流速分布一方面具有对演化历史的记忆功能,另一方面又具有反映演化状态的不可逆的耗散性。高含沙水流的阻力平方区的对数流速分布其实质上在于高含沙水流和清水水流的演化机制的相似性,这个问题下面还要讨论。

2~3\;w`)_5Ks0

3 紊动特性水利论文(L-d0{1c d)dq

  文献[8]、[9]对高含沙水流的水流结构及运动机理进行了深入的研究。研究认为(试验用沙取自于黄河花园口放淤区,中值粒径d50=0.021~0.027mm,小于0.01mm的细沙约占全部沙重的7%~17%,试验水流的含沙量S=129~666kg/m3),在非流核区,高含沙水流的相对纵向紊动强度σx/u*沿垂线的分布规律基本上与清水相似,令σx、ux分别为当地沿流向的紊动值(Image788.gif (1007 bytes))和流速值,σx/ux表示相对紊动强度,则无论对于清水水流或高含沙水流,其σx/ux数值变化于0.021~0.08之间,与含沙量和雷诺数并无明显的关系,即单位时均流速所对应的бx基本上是同一数量级的,令ηα表示相对高度,则其脉动强度分布如图3所示(a 组含沙量529<S<666kg/m3,b组含沙量90<S<472kg/m3,c组为清水)。

w6_s\RTw0

  文献[10]、[11]研究了宾汉体泥浆湍流的结构特征,其中Rem=ρu(4H)/η(1+τBH/2ηu, 其中ρ为浆体密度,u为断面平均流速,H为泥深。对脉动流速的量测结果表明,Rem<2000时为层流,2000<Rem<8000时为过渡湍流,Rem>8 000时为充分发展的湍流,若定义y+

q an0aT$z\&k c:j0

y+=ρyu*/η(1+τBH/2ηu)水利论文-w*TI&Q*kb[ {:y

(8)

试验发现泥浆湍流的近壁层流层范围为y+=0~1,与牛顿体湍流的粘性底层不同,泥浆湍流的近壁层流层内完全保持层流状态。若定义间隙系数φ为

,Wo2S]5F_?0

φ=流速发生脉动的时间/统计测量的总时间(9)水利论文 U7Um&U WC [c.NW3g

则依据时均流速u/u*、脉动流速Image789.gif (929 bytes)(u*为摩阻流速)泥浆过渡区湍流沿垂直床面方向可以分成五层,即(1)近壁层流区(2)间隙湍流区(3)强脉动层;(4)扩散层;(5)表面层。

:Ms6A)\;E P#r0

,c$G*zA'~6D Bc0

图3 脉动强度分布图
,d })d]2uzn0Fig.3 The distribution of fluctuating intensity

@ d)hUBU.QN0

  在近壁层流层(0<y+≤1)中,完全没有脉动,这一点和牛顿流体的近壁层性质不同。水利论文Zj U wWPk{'wN

  y+=1~4为间隙湍流层,该层相当于牛顿流体的粘性底层,并具有类似的猝发过程,该层紊动强度较大,对整个湍流的形成有重要意义。水利论文-T z q&y-EI7I

  强紊动层,从y+=5一直到流核区,或者流速梯度较小的流区的下方。该层中φ≈1,紊动旋涡在该层进一步发育和分解,各种频率的脉动都在此出现,脉动强度也较大。

|-_!o T~W1K2k8u0

  在流核区或时均流速梯度很小的中上层,流体剪切变形率很小,不会产生紊动,该层测到的流速脉动都是由下层扩散而来的紊动旋涡所引起的,因而称之为扩散层。扩散层脉动强度较小,甚至为0。与强脉动层和间隙湍流层不同,虽然二者的φ值较小,但发生湍流时,流速脉动较大,有较大的脉动强度,实测结果表明,充分湍流区泥浆流速分布和清水没有区别,过渡湍流区上部仍有流核,下部流速分布,也是对数分布,清水加推移质和泥浆加推移质的湍流时均流速分布仍为对数分布。水利论文+R[*RUD!C

  令G(n)示R(t)=Image790.gif (1107 bytes)的频谱比密度,实验结果知,泥浆湍流的紊流能量集中于低频脉动,推移质可以吸收低频脉动,并转化为高频脉动,因而低频脉动能量比例大为降低,中高频脉动能量比例增加2~4倍,推移质运动还使泥浆湍流中含能和耗能旋涡尺度均减小,而在清水中,推移质运动减小平均流速,增大脉动流强度,使湍流脉动微观上随机性增大,客观上相关性加强,充分湍流区清水和泥浆的脉动频谱分布如图4示。显然,G(n)满足∫0G(n)dn=1。实测结果表明:泥浆湍流(宾汉体)中紊动产生于间隙湍流层上部和强紊动层下部,间隙湍流层和扩散层的紊动有间歇,造成分布曲线的双峰。泥浆浓度增大,自相关系数增大,低频脉动能量增加。

  应用紊流随机理论研究均质宾汉流体的脉动结构(公式较繁,具体结果参看[3],[4])发现,充分紊动时其结构和清水时类似,表明二者具有类似的特性。

]\&sm0dV*w/^ C0

7@[-i E"a2p]0

图4 充分湍流区清水和泥浆的脉动频谱分布
5~6u9_,N g5i&Oh0
Fig.4 Fluctuating spectral density distribution水利论文RNrG+p
of developed turbulent range of clear-water and mud flow
水利论文'o-h4G$yVT`p

4 阻力关系水利论文 y%I3y ?;l&A#].D

  文献[12]对管道高含沙水流的阻力特性的研究表明(泥沙的中值粒径D50分别为0.031,0.042mm,含沙量变化范围为0~630kg/m3),无论是在管道或明渠中,在由层流向紊流转变的过渡区内增加含沙量时,阻力损失相对减小,称之为过渡区减阻,其原因是由于增加含沙量使水流粘性增大,抑制了水流紊动,从而减小了阻力。即使如此,实测结果清楚显示了与同流速的清水相比含沙水流是增阻而不是减阻,即使在整个过渡区,甚至包括临界点在内,也不小于清水的阻力损失。[13]从高含沙运动的物理图象出发,概念清晰地指出,高含沙水流中的紊流的阻力方程为

IK G:moe/X0

τ=τBT+ηdu/dy+τTm

@c9^cAHVw0
(10)

其中τTm为结构流中紊动应力,即

;no(GR E,ic0

τTm=-ρmImage791.gif (940 bytes)

4B3G\Pl6jayHY1o0
(11)

u′m认为是由内聚力作用所引起的作用力,式(10)垂线平均后为水利论文VU9YCai4nf9o

Image614.gif (845 bytes)=1/2Image792.gif (872 bytes)hJmBT(1-y0/2h+η(umax)m/h+1/h∫0h-y0τTmdy水利论文0jROk-L_6E$z'\

(12)

式中 h为水深,y0为流核厚度,umax为水面最大流速,比较相同水深条件下,清浑水的平均切应力可得水利论文^*w&g9jA@_F N!C

Jm/J0BT(1-y0/2h)h+η(umax)m+∫0h-y0τTmdy/μ(umax)0+∫h0τT0dy(13)

式中 脚标0表示清水,在相同条件下

{2o t&Sc[0

λm0=Jm/J0水利论文6t|%y3z-~

(14)

在层流区时,相同条件下

:v Q.}&mUx:Ff0

λm0=[2(ue-η)(2-α)/η(3-α2)+(1-α)2/1-1.5α+0.5α3]vm/v0>1

iex.zTX)Nc&`B0
(15)

式中水利论文+@xW4R N9sr

a=τBww=rmhJm

@/u^sAJ7nVn6G0
(16)

在阻力平方区时

g:C%M#a2ye:D7Q Q0

λm0=η(umax)m+∫0hτTmdy/μ(umax)0+∫h0τT0dyImage793.gif (861 bytes)0

D$U0ZHd^P)pDL0
(17)

分析知上式中λm0<1。如令ET表示紊动损耗的能量,kt表示紊动损失占总能量损耗的比例数,则有

;E9t.Vhdp0

ET=v/h20h-y0τTmdy=kTrmvJm水利论文Q7D0^? k)P#t!\fe

(18)

kT=1/2-η(umax)mBT(1-1/2y0/h)h/τwh水利论文`Z1swdO

(19)

以阻力公式代入知水利论文;x#r;c&q E(Q.Yk

kT=1/2-[8λm/Reu·η/ue·(umax)m/v+τBTw(1-1/2y0/h)]

Y x$rOx*D0
(20)

ET由三部分组成,由于粘性耗散为热能EH,由于克服内聚能部分EB,抗拒重力作用的悬浮功ES,则水利论文/r7Gl)Q(S{ {

ET=EH+EB+ES

\&a.M HH0
(21)

  减阻和增速是一个问题的两个方面。高含沙水流的阻力和清水相比,在层流区时阻力损失较同条件下的清水为大,这似乎已成共识。也有人认为紊流粗糙区时二者阻力近似相同,而对过渡区,似乎未取得一致意见,这可能与比较标准不同有关[6]。不同浓度条件下的高含沙水流运动情形下的减阻,和泥沙的浓度及电化学、物理特性密切相关,从紊动来看,和猝发的周期、流条的空间尺寸有关。虽然不同浓度条件下的高含沙水流运动减阻问题在管道输运和治黄中得到了足够的重视,但关于减阻机理目前还不太清楚(一般情形下得到一定的含沙浓度和同条件下的其它含沙浓度相比具有减阻作用,或含沙浓度和阻力的定量关系),起码和高分子减阻情况相比是这样,因而其研究尚需进一步深入。水利论文CEZ&d2h

  需要指出的是,关于阻力问题,杨文海、赵文林,戴继岚、万兆惠等都进行了研究,窦国仁对此进行了分析并得到了一些理论结果[3,4],但关于阻力问题是一个远没得到很好解决的问题,因而上述关于阻力问题的讨论仅是相当初步的。

6er3F!rKb oW l!N?0

5 能量特性

5P${7s,e(b"vE0

  两相流中的能量关系一直是人们感兴趣的研究方向,继维里坎诺夫、张瑞瑾、吴德一先后做了有益的研究后,陈立、刘兴年,王光谦、倪晋仁张红武、舒安平等都对此作了有益的探索(这里介绍他们研究的共性部分,结合[5],以期对能量关系有整体性的认识)。文献[14]从固液两相流理论出发,推导了二维两相挟沙水流紊动能量时均平衡方程式

6sM'k |m [k0

Image794.gif (1894 bytes)水利论文/bA ZN[ D*m

(22)
式中 ue为混合相有效粘度,ue=um+uT,其中um为混合相动力粘性系数,uT为紊动粘性系数,σB为混合相普朗特(prandtl)数,Bt=1/2Image795.gif (918 bytes)为单位质量混合相紊动动能。

  总流方程式水利论文Wi4T-\ Qv4d

Image796.gif (1677 bytes)

da#z)a$B*w'R0
(23)

其中τ为水流总切应力,τBT为宾汉切应力,ux,uy为x,y向脉动流速,以es表示悬移质运动效率系数。水利论文f \ U5I7|@,\

定义es=Es/E1=泥沙悬浮功/单位时间紊动能=

9iCfoU6R5b8O0
Image797.gif (1235 bytes)  (24)

用kt表示紊动能转化率,定义

&}|WzSH%W9u0

kt=E1/E=单位时间紊动能/单位时间总能量=水利论文7d;?x'z7[YE v

Image798.gif (1411 bytes) 

(25)水利论文q"^-A {uZu)E

易知es实为理查逊数。推导知水利论文9_!xR.T }Zi9QF

  效率系数es的理论表达式为水利论文%^H[x o[4q)Q0v

es=φ1/κ2κt(fm/8)1/2水利论文j8^4S4`Z7dw u.F

(26)
  而

φ=1.64[1-cucd(1+c21)3/23s水利论文9Y:J/K.U,S'{f]C-x

(27)

φ=f1r)f2(κ,fmmsm,u2/gκ,w/u)

3U|7X2o1w]O Z8F(T2Y0
(28)

其中 fm为浑水的阻力系数,cu,cd,c1为系数,ηs为无因次高度,κ为卡门常数,w为沉速,f1,f2为形式函数表达式。由实测资料知es、sv关系为如下图5所示。水利论文&\5~n}&@9up N

  经过推导知kt表达形式为

ZL1@;|s(F"\k0

kt=1-1/τ0h[τBT(2h-y0)+2umkdv水利论文%o:g5FZ5_#o

(29)

  其中y0BT0h为高含沙水流流核厚度,h为水深。令umax=kdu,kd、u分别表示流速分布形状系数和垂线平均流速。一般情况下kt<1,只有在充分紊动的挟沙水流中可近似取kt≈1。

L%B7t8tn'~J[0

  在清水和低含沙水流中曾论述,卡门参数是空间点坐标和时间的函数,是流体运动过程中流体质点所处的运动状态的一个系统整体参数,和雷诺应力有关。[15]指出,在清水中加入单一粒径的泥沙(d=0.16mm),粗略地讲,当含沙浓度S>50kg/m3时,浑水的紊动强度小于相应清水的紊动强度,否则,则大于。由于流体运动过程中平均流能量和紊动能是两种不同性质的能流,图2实际上反映了加入泥沙后泥沙对水流结构的影响。水利论文0aq6L/t5BPK

  在[5]中,我们详细计算了流体运动过程中的能量分配关系,诸如粘性耗散能,紊动能产生和紊动能等等并和实验结果作了对比分析,由式(24)和图5容易估计泥沙悬浮功的大小。这样,对实验室水槽中的悬移质运动和紊动能转化等各种能量间的相互关系有了一个清晰的了解。水利论文R*s p$x(\[dG

5Y@ev!G!z0

图5 实测效率系数es随含沙浓度sv的变化
H%q{z&Z8n,C0
Fig.5 Variation of measured effective coefficient with sediment concentration

-GPT6t%C%p ?0

  其实,出现在阻力特性和这里的紊动能转化率是同一个系数,它们都是理查逊数的函数在不同情形下的具体表达形式。对于含有悬移质的明渠中的能量关系,钱宁早就指出,由于悬浮泥沙的能量(悬浮功)取自于紊动动能,而其最终来源是由有效势能转化而来的,一般情况下,悬浮功不过占有效势能的百分之四、五,或者更小;即使在高含沙条件下,一般也不超过百分之十。紊动产生能的传递方向是自下向上,而有效势能的传递方向是自上向下。结合第三部分对紊动特性和脉动频谱的讨论对流动过程中的能量关系就会有更清楚的认识。6 讨论和结论水利论文Mx-BEB%["wr

  在低含沙水流中曾指出如下事实,悬浮泥沙的尺寸增加,将使底壁水流更容易猝发,流动雷诺数增加将抑制猝发;泥沙的纵向流速比水流纵向流速大得愈多,愈不容易导致流体猝发的形成,但当泥沙的纵向流速比水流的纵向流速为小时(天然情形中常常这样),较易导致猝发的产生;泥沙尺寸在一定范围内愈大,将使猝发次数增加,使湍流度和雷诺应力增加,泥沙含量增加,使猝发频率和流条平均空间尺寸不变。观察表明,似乎是近壁区的低速流条的上升和破裂控制着泥沙的传输,泥沙大部分在低速流条区积聚,后被流体上举带离床面,开始运动的。猝发历时随壁粗糙度的增加而减少。脉动是由猝发引起的,从直观上讲,相同条件下,猝发频率越高,则脉动也越强。对于单个颗粒在流体中的运动特性,其具体细节尚不完全清楚。结合实验结果可知,推移质运动使其惯性平衡区的特征尺度变小(这里借用流体运动的惯性平衡区),因为含能旋涡主要集中于大尺度旋涡,耗能旋涡主要集中于小尺度旋涡,这和流体粘性减小的效果是一致的(仅在所论特性方面)。出现这种现象的原因可能和平均流和脉动流的关系有关,传统认为,平均能不直接转化为紊动能,而是先蜕化为紊动产生能,再影响脉动。

'\zjW^Mz0

  牛顿流体的本构方程为

MJ1l2U'|"P vH'M pF0

Image799.gif (1558 bytes)

2f;FyYp#a}0
(30)

其中 p为动水压强,μ为动力粘性系数,τij为应力(这里符号和材料力学中符号相一致),若为不可压时,有水利论文Dhl5zh f&? _;rf&b

τij=-pδij+水利论文{)a [Dg:V;h a

Image800.gif (1235 bytes)

     (31)水利论文Ctb F_:^7m7q

应用牛顿第二定律,在连续介质假定成立时,有水利论文.AN*ZZ&MI1HBt

Dui/dt=fi+1/ρ水利论文 \2F9eK+^(b

Image801.gif (967 bytes)    (32)

其中fi为质量力。

2D9xVN*J\0

将(31)代入(32)式知有N-S方程成立水利论文&WvJ"rUlp

Image802.gif (1211 bytes)

4}'u#aS0K-U&ka u0
(33)

对于宾汉流体,依据实验成果假定应力关系符合下式水利论文(O"zn3fN-L6N

τij=-pδij+

q ~ro+ZR7V Q0
Image800.gif (1235 bytes)ij0

(34)

{-E?c+G&?/w1I0

将(34)代入(32)式知有下述方程成立水利论文1m+h[ rxY2I'O

Image803.gif (1909 bytes)水利论文~U|N^

(35)

  比较(35)、(33)两式知,对于空间某点而言,若η为常数,τij0为常数,则(35)、(33)两式相同,此时在相同边界条件下,高含沙和泥浆的宾汉体表现出和清水相同的运动规律。若η=η(x,y,z,t),τij0ij0(x,y,z,t)时,即流变系数是宾汉体运动状态的函数(这在实际情形中常会发生),此时高含沙和泥浆的宾汉体表现出和清水相似的运动规律(其异同如前述,此不赘述)。

E X| ^i0

  由于方程(31)成立的基础是连续介质假定成立,由上述知,决定清水和高含沙水流运动特性异同的基础是二者之间本构关系的异同。在某些情形下,流变参数是随流体运动状态的演化而变化,使高含沙水流表现出和清水水流相似但又具有其自身的运动特性。水利论文&\ _azDfDh9j e[

  关于脉动产生的本质,即紊流的发生问题,是一个没有很好解决的问题,根据近年来湍流研究的成果结合实验结果,笔者认为,清水加推移质、泥浆加推移质中紊流发生的本质和清水、泥浆流中紊流发生的本质相同,差别仅在于流体特性及推移质对流体运动的不同影响造成的。其实,如果撇开具体细节,就紊动发生的本质,或者脉动在平均流中产生的机制,二者有统一性的一面。因而实验结果显示了清水、泥浆、或二者加推移质后所表现的特性在本质有共性一面的基础上,表现出特别的和本身物理特性相符的各自的紊动特性来。经过分析,笔者认为在低含沙水流中的平均流和脉动流的能量关系的图象原则上也适用于高含沙水流。水利论文c w7q7s,sfr6b

  如何从实验出发,探讨泥沙状况和水流结构间的本质关系及内在联系,有待今后工作的继续深入。

dW7D,~k5^D0

  本文系统分析、比较、论述了高含沙水流和低含沙水流、清水水流时相应的流速分布特性、紊动特性、能量间分配关系等流动现象的内部规律及其异同,并指出了高含沙水流所固有的特性,为进一步研究高含沙水流的其它特性奠定了基础。水利论文S^a-HP_

致谢:衷心感谢万兆惠、匡尚富教授对本工作的指导,衷心感谢王力、李文武等全室对本工作的支持。水利论文x)q&E1Z*vxIq+c~-P

参 考 文 献水利论文Ny9pR f&b/Ak

[1] 中国科学A, 1982.6.水利论文2bUg'^6x| |

[2] 王明甫,段文中等。高含沙水流的水流结构及运动机理。中国科学A,1987.5.[3] 黄亦芬,窦国仁。均质高浓度紊流的时均结构和阻力系数。水利水运科学研究。1989.

j"m1nbq0

[4] 黄亦芬,窦国仁。均质高浓度流从层流向紊流的过渡及紊动强度。水利水运科学研究。1989.

-^5gzkN.e#xzX0

[5] 刘兆存,徐永年。低含沙水流运动特性综述。泥沙研究,1999,(2).

G,Z qd`0

[6] 宋天成,万兆惠,钱宁。细颗粒含量对粗颗粒两项高含沙水流流动特性的影响。水利学报。 1986,(4).

D8A4uzq0

[7] 王光谦,倪晋仁。固液两相流流速分布特性的试验研究。水利学报,1992.11.

QW0Z*E!we"@$X0

[8] 王明甫,詹义正。高含沙水流紊动特性的试验研究。第二次河流泥沙国际学术讨论会论文集。1983,南京。

o!H{!y8_o0

[9] 王明甫。高含沙水流及泥石流。水利电力出版社,北京,1995.

1D!q!B?;RO)UR)f0

[10] 王兆印,张新玉。推移质运动对清水和宾汉体泥浆湍流结构的影响。泥沙研究, 1991,(2).

x2Spg t0WS0

[11] 王兆印,任裕民,王兴奎。宾汉体泥浆湍流的结构特征。水利学报, 1992,(12).水利论文;dR3P)nW

[12] 王兆印,任裕民,王兴奎。宾汉体泥浆湍流的统计特征及脉动频谱分布。水利学报,1993,(4).

{7?5Vw;v*a,uA0

[13] 任增海。管道高含沙水流的阻力特性。泥沙研究, 1985,(3).水利论文&n kYxb'q'q;O

[14] 乐培九。高含沙紊流的阻力及其输沙特性。泥沙研究,1986,(2).水利论文TR7Eod_%k|

[15] 舒安平。高含沙水流挟沙能力及输沙机理的研究。清华大学博士论文,1994.

N` JXhB2b%U0

[16] 梁在潮。粒子相的分类及对紊流场的作用。泥沙研究。1996,(3).

;X[H[)eC(i0

[17] 段文忠,刘建军。过渡区光滑明渠高含沙水流阻力的实验研究。泥沙研究,1991,(1).

d UN6h\d(^bC0

[18]曹如轩,陈诗基,卢文新,任晓枫。高含沙异重流阻力规律的研究。第二次河流泥沙国际学术讨论会论文集。南京。水利电力出版社。1982.

L k9O8d;e u['uBY0

[19] 钱宁主编。高含沙水流运动。北京。清华大学出版社。1989.水利论文` Tm Kc~7`

[20] 费祥俊,朱程清。高含沙水流运动中的宾汉切应力。泥沙研究,1991,(4).

x\+{ozA+E0
TAG: 含沙水流 刘兆存 特性 徐永年 运动
顶:9 踩:19
【已经有88人表态】
12票
极差
8票
很差
11票
较差
8票
稍差
15票
稍好
9票
较好
11票
很好
14票
极好
下一篇:风洞中挟沙气流水平集沙实验研究(倪晋仁, 李振山)
上一篇:长江中游界牌河段过渡段控制工程数值分析研究(陈晓云,高凯春)
查看全部回复【已有0位网友发表了看法】

广告投放

广告投放