明渠非恒定流传播特性及流速分布研究.
第11章_明渠非恒定流

Q vwA vwB
其中, B ( B1 B2 ) / 2
22
二、动量方程
A1 P1
v1-vw A2 v2-vw P2
l
A2 (v2 vw )[(v2 vw ) (v1 vw )] F P P2 1
A2 (v2 vw )(v2 v1 ) F P P2 (11.50) 1
不可压缩流体,ρ
=const
Q=vA
(10.20)
A Q 0 t s
(11.3)
(3) Q / s 0
A / t 0 Q 沿程不变,恒定流
10
11.2 明渠非恒定渐变流的基本方程式
一、连续性方程
( A) ( vA ) 0 t s
不可压缩流体,ρ
=const
Q=vA
(10.20)
A Q 0 t s
(11.3)
A ( Av ) 0 t s
式(11-4)是明渠非 A v A A v 0 (11.4) 恒定流连续性方程 t s s 的另一种表达式
矩形断面明明渠,A=bh
h v h h v 0 t s s
(2) Q / s 0
上 a b t2 ds t1
A / t 0 Z 随 t 下降,落水波
下
a b 下
上
这说明如果流进的 流量少,流出的流 量多,微分区间内 水位将随时间而下 降,明渠中会产生 落水波。
9
11.2 明渠非恒定渐变流的基本方程式
一、连续性方程
( A) ( vA ) 0 t s
11.1 明渠非恒定流的特性及波的分类
一、特性
1、水力要素如u、Q、A、Z或h等都是时间t和 位置s的函数,它是非恒定的非均匀流动。
明渠恒定非均匀流渠道流量实时监测方法研究与探讨

பைடு நூலகம்
EXPLoRATI oN AND STUDY oN THE REAL —TI E oNTI M M oR ETHoD F M Ol S TEADY NoN —UNⅡ U FLoW N oI I oPEN CHANNEL
2 .Wa r osrac ueuo e—kuds c,hn og ogig 5 20,hn ; t ne nyB ra f eC v H o ir tS adn ,D n y 7 0 C ia ) t i n2
Absr c : d a lc rn i e o o c lu ain, c a a trsis o ta y u i r a d n n— u i r o t a t Hy ru i sp cplsf r f w a c lto i l h r ce t f se d nf m n o i c o nf m f w, o l c re tp o lms o o mo io n n o e h n e r e c b d.F rt o v nin e o o c l u ain. u r n r b e fr f w n tr g i p n c a n lwe e d s r e l i i o he c n e e c ff w ac l t l o
a d f r n e e u t n w su e o mo i r t e w trlv l c a g s o ta y n n—u i i o p n c a n l i ee c q ai a s d t nt h a e e e h n e f e d o o o s n f T f w i o e h n e o nl n i s a fdf rn il q ai n .T i a e p l d mo e n t c n lge omo i r t et o w tr e e aU s o n t d o i e e t u t s h sp p ra p i d r h oo ist n t h ae v l l e f e f ae o e e o w l V
明渠水流特点

当 A = 一定,要求:Q → Qmax
梯形水力最佳断面条件
i、n、A及m = Const
(b, h) min
d 0 dh d2 0 2 dh
A A(b, h) (b mh )h 2 b 2 h 1 m
A 2 mh 2h 1 m h
由均匀流条件:P1=P2 ;v1 =v2 ,则
G sin F f or i sin Ff G 0
B P1 v1
C α G τ0 Ff
2
v2 P2 D
(Z1>Z2)
A
v1 P2 v 2 Z1 Z2 hw 2g r 2g
P1
2
必要条件
恒定流
流量沿程不变(无分叉和汇流情况)
h
水面宽度 过水面积
B=b+2mh
A (b mh)h
b 2h 1 m 2
R A (b mh)h b 2h 1 m 2
湿
周
水力半径
1 2 / 3 1 / 2 1 A5 / 3 1 / 2 Q VA AC Ri A头线互相平行
J
h v i Jp
v
h
J Jp i
θ
底坡、水面坡度、总水头线互相平行
二、 明渠均匀流产生条件
充分必要条件
•力学条件
渠壁摩擦阻力与水重力在流动方向的分力始终平衡(大 小相等,方向相反) •能量条件 水面下降值始终等于水头损失。
由动态平衡,得
P1 G sin P2 F f 0
明渠均匀流输水能力的影响因素,决定于过水断面的大小和形状 (i随地形而定,n决定于壁面材料),当i、n和A一定时,使通过 流量最大的断面形状,称为水力最优断面。
非恒定流垂线流速分布规律的初探

式 中 : 为水 面 比降 。 . ,
令
() 2
将 上 述有 关 项 代 入 式 ( —1 即 得 : 1 )
+ ) +警 _ P P
式 () 非恒定流剖 面流微方程 。 3即
() 3
2 次 生 流运 动 方程
2 1 方 程 的 建 立 .
在恒定均 匀流条件 下 , 即
天 然河 速 分 布 常 常 偏 离 对 数 分 布 规 律 。涨 潮 时 底 层 特 大 其 流 速 大 , 层 流 速 小 ; 潮 时 表 层 流 速 加 大 , 层 流 速 减 少 。在 有 径 流 加 盟 时 , 表 落 底 由于 盐 水 楔 活 动 , 种 现 象 不 这 难 理 解 。在 无 盐 水 楔 或 其 他 类 型 异 重 流 的情 况 下 也 有 此 类 现 象 。这 里 试 图从 水 流 的 非 恒 定 和非 均 匀 性 所 引 起 的水 流 结 构 变 化 去 寻 求 解 答 。
垂 线 流 速 分 布 一 般 为 B型 和 C 型 。 B型 分 布 在 主 流 区 , 布 梯 度 大 于 A型 ; 分 C型 分 布 与 A 型 相 反 , 布 分 梯 度 一 般 为负 。 关键 词 : 恒 定 流 ; 定 非 均匀 流 ; 速分 布 非 恒 流 中 图 分 类 号 : V 3 . T 13 2 文献 标识 码 : A 文 章 编 号 :0 5—84 (( )2—0 5 10 432  ̄2 0 0 4—0 6
水力学专题 明渠非恒定流

z z z z zw zw udz vdz u w v w w u b v b w 0 z z x b y b y y x z z w x z zb
结合边界条件:
z w (hu ) (hv ) 0 t x y
三、水波的分类 明渠非恒定流是一种具有自由水面的波动 (一)按主要作用力分 重力波(gravitational wave) :主要恢复力为重力――洪水波 表面张力波(capillary wave) :风波成长初期的涟波 压力波(compression wave) : (少见――水下爆破引起等) (二)按质点运动方式分 振动波:波运动时,其水质点沿封闭的轨迹运动――海洋波浪 运行波:波运动时,水质点不断行进――洪水波
u
z zw
u
v
w
z zw udz u w z x b x
dw
zb x
z zb
z zw vdz v w z y b y
v
z zw
zb y
z zb
zw
zb
w dz 0 z
w w w dx dy dz x y z
这里,我们只给出了两个方向的动量方程,z 方向的动量方程省去了。事实 上,由于垂向上的加速度和重力加速度相比可以忽略,垂向的速度 w 很小,所以
w 亦可以忽略。因此,垂向上的重力加速度和压力梯度平衡:
p g z
直接积分得: p g ( z w z )
一般地,河道中的水流均为紊流,所以要研究河流中的平均运动,就要用 雷诺时均方程,雷诺时均方程和 NS 方程相比,只是增加了紊动切应力项。 由于脉动流速在各方向上具有相同的数量级,所以雷诺应力在各方向上也 具有相同的数量级。但只是他们的梯度出现在方程中,而雷诺应力在垂向变化的 距离比纵向和横向上的小得多,故垂向上的梯度变化是最重要的。 鉴于此,描述自由表面流动的雷诺方程可简化为:
水力学系统讲义第八章-明渠流动

流量Q=25.6m3/s,过水断面宽5.1m,水深3.08m,问渠底坡
度应为多少?并校核渠道流速是否满足通航要求(通航允
许流速[v] ≤1.8m/s)
解:
将Q AC
Ri K
i写成i
Q2 K2
Q2 A2C 2R
R A 5.1 3.08 1.395m
5.1 2 3.08
b 0.83h 1.64m
(2) 2 b / h
A (b mh)h (2h h)h 3h2
b 2h 1 m2 2h 2h 112 4.828h
R A 0.62h
Q A R2/3i1/2 1.542h8/3 n
h 1.55m
超高
m
3.2m
b
解:按均匀流计算,当超高为0.5m时,渠中水深 h=3.2-0.5=2.7m,此时断面要素为:
A (b mh)h (34 1.5 2.7) 2.7 102.74m2
b 2h 1 m2 34 2 2.7 11.52 43.( 1 m2 -m)=2( 112 -1)=0.83
h
A (b mh)h (0.83h h)h 1.83h2
又水力最优 R h 2
Q AC Ri A R i 2/3 1/2 0.815h8/3 n
h ( 5 )3/8 1.98m 0.815
1h
h
其中 为宽深比
m
b
b 2mh 2h 1 m2
梯形水力最优断面的水力半径: R A (b mh)h (b mh)h h
b 2h 1 m2 b b 2mh 2
工程流体力学中的流速分布与流动特性

工程流体力学中的流速分布与流动特性工程流体力学是研究流体在各种工程应用中的力学性能和行为的学科。
在工程流体力学中,流速分布和流动特性是非常重要的概念和参数。
流速分布描述了流体在横截面上的速度分布情况,而流动特性则涵盖了流体的速度、压力、阻力、粘性等特性。
首先,我们来讨论流速分布在工程流体力学中的重要性。
流速分布是指流体在流动过程中的速度分布情况。
在工程应用中,了解流速分布可以帮助我们预测和控制流体的运动方式和行为。
例如,在石油工业中,了解油管内流体的流速分布可以帮助我们优化管道布局和提高油井的产量。
在水力学中,了解水流的流速分布可以帮助我们设计和优化水电站的水轮机,提高能量转换效率。
接下来,我们来探讨流动特性在工程流体力学中的重要性。
流动特性是指流体在流动过程中的速度、压力、阻力、粘性等特性。
了解流动特性可以帮助我们分析和预测流体的运动规律和行为。
在风力发电中,了解风的流动特性可以帮助我们选择合适的风力发电机组和优化风力发电场的布局。
在航空航天工程中,了解气流的流动特性可以帮助我们设计和优化飞机的机翼和气动外形,提高飞机的操纵性和燃油效率。
在工程流体力学中,还存在一些常见的流动特性。
例如,当流体通过管道或孔隙时,会发生收缩和加速的现象,称为孔隙流动。
在孔隙流动中,流速分布呈现逐渐加速的特点,其原理可以用伯努利方程来解释。
此外,还存在旋转流、透明流、层流和湍流等不同的流动特性。
这些特性的不同对于工程应用具有重要意义。
除了流速分布和流动特性,工程流体力学还研究流体的压力变化和阻力特性。
压力变化是指流体在流动过程中由于速度变化引起的压力变化。
在管道输送中,可通过流速分布和流动特性的研究,分析流体在管道中的压力变化情况,从而优化管道的设计和运行参数。
而阻力特性则是流体在流动过程中与物体碰撞并产生的阻力。
在流体力学中,研究流体的阻力特性可以帮助我们计算流体对物体的压力以及物体所受到的力。
总结来说,流速分布和流动特性是工程流体力学中的重要参数和概念,对于工程应用具有重要意义。
明渠非恒定流的特征线解法

a ueya do e lwst eb n c tl n v f o h a k.M o e v ri t eisa c eiu ,i m a a s h a k c atd.Th ril r o e h n tn ei sr s t f s o y c u t b n o se e e ea t e c
文章编号 :00—9 3 2 0 0 —0 2 10 83(0 6)1 0 0—0 3
明渠 非 恒 定 流 的特 征 线 解 法
万 晖
( 河海大学 水利水 电工程学 院, 南京 209 ) 10 8
摘
要: 长距 离输水 明渠 中, 常遇到 闸门启 闭过程 中, 中水面波动 比较剧烈 , 渠 导致漫堤 的情况 , 重的可能会 引起 渠堤失稳滑坡 。 严 利用特征线法对描述 明渠 中水流运动状态的圣维南方程组进行求 解 , 编制程 序求解 了简单输水 模型 中, 闸门启 闭情况下
0 引
言
为一个个短 的渠段 , 水流通过节制闸, 从上游渠段 向
下 游 渠段输 送 。每个 渠 段 类 似 于 一 个 小水 库 , 距 长
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1 明渠非恒定流传播特性明渠非恒定流传播过程中,存在波的变形,不但沿时间存在横向变形,而且纵向波高上也在变化。
通常来说,周期性非恒定流在传播时,随着传播距离的增加,其上升段往往越来越陡,下降段越来越缓,甚至出现波的破碎现象,但总的周期保持不变;而在波的高度上,随着传播距离增加存在坦化现象,波幅越来越小,波峰变矮而波谷变高,非恒定流有逐渐均匀化的趋势(如图1所示)。
非恒定流的横纵变形,导致非恒定流在传播过程中沿程水深和流量变幅及波形不一致,使得非恒定流的传播表现出与恒定流不同的性质。
从诸多工程问题来说,波速如何确定,波峰和波谷的高度如何计算等是大家较为关心的,但是,1.1 试验条件1.1.1 非恒定流过程的概化非恒定流试验都采用正弦函数表示的周期性非恒定流过程。
对于周期为T 的非恒定流,给定如下形式的非恒定流过程:)22()(210ππ-+=T t Sin x Q Q Q b (1)20 t(s)Q ,H1 明渠非恒定流传播过程波形变化图示)22()(210ππ-+=T t Sin x A A A b (2)式中:0Q 为基流;b Q 为非恒定波的波幅,即流量的变幅;t 为时间;0A 为基流的过水面积,b A 为面积的变幅。
1.1.2 水槽底坡的确定为了简化问题寻求规律,同时也方便同均匀流对比,本文试验均在正坡条件下进行,坡度为3‰、5‰,恒定流时在自由出流状态下能形成均匀流,非恒定流状态下当周期较大则趋近于均匀流状态。
1.1.3 试验段的选择水槽的4#~7#水尺之间,在恒定流时能形成均匀流,而在非恒定流状态下, 4#~7#水尺之间最大最小水深相差很小,最大水深连线和最小水深连线基本水平。
因而,非恒定流状态下,4#~7#水尺形成一种特殊的、相对较为稳定的非恒定流,其波高基本不变,4#~7#水尺之间即为本文非恒定流的试验段。
x(m)H (m )图2 比降3‰时非恒定流波的沿程分布 (T=20s,Q=15~40L/s)1.1.4 水流条件从前面的分析也可知,在3‰、5‰坡度下3~40L/s 之间的流量时,恒定流在28m 水槽4~7#水尺之间能形成均匀流,非恒定流下则能形成相对稳定的非恒定流状态。
本文在这两种坡度下和5~40L/s 的流量变幅之间,根据不同的流量、比降和周期组合,共进行了40组次的28m 水槽非恒定流试验,给定流量过程均为正弦函数分布,各组次的试验条件分别见表1。
表1 28m 大水槽非恒定流试验水力参数表x(m)H (m )图3 比降5‰时非恒定流水深的沿程分布(T=20s,Q=15~40L/s)1.2 非恒定流传播速度采用4#和7#两把水尺之间的距离除以各自波峰达到的时间差来计算非恒定流的传播速度:)#4()#7(#4#70peak peak t t x x dt dxC --==(3)分别代入max min H H 、利用上式计算波速,将其与实测结果在图5.10中进行对比发现,(3)式计算的结果与实测结果差别较大,之间还相差了某一特征流速。
本文的非恒定流过程的非恒定流传播不仅是微波传播,而且是整个水体向前运动,因而传播速度不能单独用微波速度来表示,需要另外加上一个水流运动速度,以表示如下:000U gh C += (4)分别代入max min H H 、和max min Q Q 、,通过(4)式对40组次的实测试验数据进行分析,结果也统计于图4。
从对比结果来看,0h 和0U 都与非恒定流传播时的最大流量对应的水深和流速更为接近:maxmaxmax 0BH Q gH C +=(5) 式中,max Q 为非恒定流过程的流量峰值,max H 为最大水深,B 为水槽宽度。
图4统计对比了实测的波速与利用式(5)计算的波速,计算值与实测值二者符合程度较高。
1.3 水深外包线的确定非恒定流传播过程中,水深变化过程即:)22()(21)()(0ππ-+=T t Sin x H x H x H b (6) 式中,)(0x H 为水深基值;)(x H b 为非恒定流水深的变幅。
由于本文的非恒定流试验段选择在4~7#水尺之间,即水槽的12~22m 之间,在该区段内波峰波谷的连线沿程基本不变,因而可以认为在某一级非恒定流下试验段内max H 、min H 、0H 、b H 相同:)(21min max 0H H H +=)(min max H H H b -=研究非恒定流波高的变化,需要研究0H 及b H 的变化,这涉及到试验段内最大、最小水深的计算,即非恒定流波峰、波谷的水面外包线。
对于最大流量为max Q 、最小流量为min Q 和周期为T 的非恒定流,定义如下参数:max H 表示非恒定流波峰的水深;min H 表示非恒定流波谷的水深;max Q H 表示max Q 的均匀流水深;min Q H 表示min Q 的均匀流水深;H 表示)(21max min Q Q +时的均匀流水深。
非恒定流在传播过程中,其最大水深(波峰)和最小水深(波谷)与非恒定流本身的参数有关,即和非恒定流的流量最大值、最小值以及非恒定流的周期有关,光滑水槽内的明渠非恒定流试验同样也表明了这一点。
图5~6分别统计了实测的3‰、5‰坡度时不同组次非恒定流在试验段内波峰、波谷的水深变化情况。
试验结果表明,在确定了非恒定流的最大、最小流量max Q 、min Q 以后,非恒定流波峰、波谷的水深随周期T 的变化呈如下特点:(1)非恒定流下的波峰水深不会大于均匀流状态下最大流量对应的水深,而是随周期在max Q H ~H 之间变化;波谷水深不会小于均匀流状态下最小流量对应的水深, 而是随周期在min Q H ~H 之间变化。
(2)当周期减小时,非恒定流波峰水深max H 和波谷水深min H 逐渐靠近平均流量)(21max min Q Q 的均匀流水深H ;(3)当周期变大时,波峰水深max H 和波谷水深min H 逐渐分散,波峰更高,靠近max Q 的均匀流水深max Q H ,波谷则更低,靠近min Q 的均匀流水深min Q H 。
测数据显示非恒定流传播时水深外包线的规律性较强,所以可以构造关于),,(min max Q Q T 的函数,确定非恒定流的波峰、波谷水深。
图5 非恒定流波峰、波谷变化规律统计(J=0.003,Q=15~25L/s )根据对多组试验结果的分析,可以构造如下型式的函数来分别描述波峰和波谷的水深:波峰 )90(2)(max max ππT arctg H H H H Q -+= (7)波谷 )90(2)(min minππT arctg H H H H Q --= (8)从图7~8中可以看出,用式(5.8)和式(5.9)描述的波峰、波谷趋势线能很好的反映非恒定流传播过程中波峰、波谷的变化规律。
从图5.19和5.20波峰波谷实测值与计算值的对比图中可以看出,二者测符合程度较高。
图7 波峰实测值与公式(7)计算值的比较图6非恒定流波峰、波谷变化规律统计(J=0.003,Q=5~40L/s )1.4 出口流量流量的外包线非恒定流传播过程中,随着传播距离的增加,能量沿程衰减,流量的波幅也逐渐衰减,波峰流量越来越小,而波谷流量则越来越大。
要知道上游发生的非恒定流传播一定距离以后的情形,必须要知道基流0Q 和波幅b Q 如何变化。
本文在28m 水槽内进行的试验研究无法测得水槽中部的流量过程,因而主要针对出口流量的变化过程进行研究。
图9统计的是典型组次水槽进出口流量的变化对比情况,从40组试验结果比较来看,非恒定流传播过程中出口流量的变化有如下特点:(1)非恒定流传播时流量的外包线变化与水深外包线变化规律基本相似,波峰流量、波谷流量在进口的波峰波、谷流量之间变化。
(2)波峰流量随周期增加而增大并逐步靠近进口的波峰流量max Q ,随周期减小而减小,并向平均流量靠近。
(3)波谷流量随周期增加而减小并逐步靠近进口的波谷流量min Q ,随周期减小而增大,并向平均流量靠近。
(4)进出口最大最小流量的平均值不变,即非恒定流传播中,基流沿程不变,存在这样的关系式:)(21)(21min max min max (出口)(出口)=Q Q Q Q ++。
从图9和10统计的进出口平均流量对比可知,出口实测的平均流量与进口平均流量几乎完全一致。
对多组试验结果分析,可以构造类似于水深的函数来分别描述出口的波峰流量和波谷流量:波峰 )90(2)(max max ππT arctg Q Q Q Q -+=(出口) (5.9) 波谷 )90(2)(min min ππT arctg Q Q Q Q --=(出口) (5.10) 出口波峰、波谷的流量计算值与实测值分别统计于图11中,从图示对比的结果可以看出(5.9)、(5.10)两式的计算结果和实测值的符合较好。
(a) Q=15~25L/s,J=0.003(b) Q=30~40L/s,J=0.005 图9 进出口波峰波谷流量变化对比图图10 水槽进出口的平均流量对比1.5 非恒定流的变形特征周期为T 的非恒定流在传播过程中,随着沿程距离的增加,非恒定流的波形在时间方向上会发生变形,一般波的上涨段变短、变陡,而下降段变长、变缓。
但是每一个波形的总时间仍然不变,维持在一个周期T 。
定义一个周期内某一距离x 处非恒定流上涨段持续时间为)(x t u ,下降段持续时间为)(x t d ,则有:T x t x t d u =+)()( (11)定义T x t x T u /)()(=α为波形随时间的变形率,可知5.0≤αT,而且,随着传图12 10s 周期非恒定流沿程波形的变化播距离的增加,αT 会越来越小。
由40组试验数据的相关分析可知,非恒定流的变形率与传播距离之间存在如下的关系:x k B T x t u 0/)(-= (5.12)在非恒定流传播过程中,)(x t u 会越来越小,而在初始位置,即非恒定流开始的地方,)0()0(===x t x t d u ,5.0=B 。
从40组试验结果看,当不考虑测量误差时,该数字也基本在0.5附近。
所以,(5.12)式即为:x k T x t u 05.0/)(-= (5.13) 式中的斜率0k 待定,其值与非恒定流参数有关。
Q=15~40L/s,T=10s,J=0.0050.10.20.30.40.510203040x(m)Q =15~25L/s,T=5s,J=0.0030.10.20.30.40.50.6102030x(m)图13 非恒定流周期变形与传播距离x 的简单相关从40组次的试验数据分析,0k 值是表征非恒定流参数的max Q 、min Q 和周期T 的函数关系式。
图14分别是典型试验组次的0k 值与max min max /)(Q Q Q -和T 的相关性分析,0k 值与这二者的成很好的函数关系,因而可以利用这两个参数构造0k 的表达式。