首页 理论教育模型验证:深槽水流和泥沙计算分析

模型验证:深槽水流和泥沙计算分析

【摘要】:图5-112中性物质在弯道中图 ,图 和图 连续排放的平均浓度验证图5-112中性物质在弯道中图 ,图 和图 连续排放的平均浓度验证图5-113试验深槽以及垂向计算网格示意图图5-113试验深槽以及垂向计算网格示意图图5-114计算深槽水流结果与Van Rijn 试验及Alfrink 计算结果比较5.5.5.3深槽水流、 泥沙计算验证图5-114计算深槽水流结果与Van Rijn 试验及Alfrink 计算结果比较5.5.5.3深槽水流、 泥沙计算验证Van Rijn 在进行了多组不同尺寸和流量的水流跨过深槽的水槽实验。

5.5.5.1 弯道水流计算验证

图5-110 实验弯道示意图 (单位:m)

为了研究弯道螺旋流对物质的混合过程,Chang(1971)对弯道上水流速度和中性物质浓度分布进行了测量,实验水槽(如图5-110)由两个90°弯道组成,弯道半径8.53m,弯道之间由4.27m 直水槽过渡,弯道进出口分别由长为2.13m 的直段过渡,弯道横断面为矩形,宽2.34m。进口平均流速为0.366m/s,平均水深为0.115m,水面比降S=0.00035。流速测量断面位置如图5-110中虚线所示。模型计算网格为354×24×22个,网格大小为0.1m 左右,计算时间步长0.1s。数值试验表明迭代计算1000步接近稳定。

模型验证了弯道13个测量断面上的流速,限于篇幅只列出第二弯道尾部 (PI/2)横断面的验证结果(图5-111)。可见,模型计算纵向流速和横向流速与实测值吻合较好,尤其是在弯道中部附近的流场计算精度相对较高。在弯道凹岸靠近自由面附近区域横向流速计算值和实测值偏离较大,主要是在该区域存在一个与主二次环流旋转方向相反的次生二次环流,而紊流模型不能够模拟次生二次环流。此外,河道边壁网格疏密和壁函数取值可能对边壁流场也有一定的影响。

图5-111 PI/2断面沿主流和垂直主流方向流速沿水深分布

y—断面横向坐标;z—与河底距离;D—水深;B—宽度;

5.5.5.2 弯道输运和泥沙计算验证

Chang (1971)试验在第一个弯道进口的左、中、右分别设置中性物质排放点。试验在两弯道过渡段和第二个弯段上13个断面测量中性物质浓度。试验用中性物质密度与水相同且在输运过程中不发生降解,因此中性物质输运方程的源项取为0。对物质在第一个弯道进口中部、右和左岸排放,模型计算横断面三维浓度分布垂向平均值与实测资料对比如图5-112所示。对中部排放,计算值和实测值吻合较好,但在弯道凹岸有一定偏离,这可能是弯道螺旋流误差引起。两侧连续排放,采用S1断面上垂向平均浓度作为进口条件,得到垂向物质平均浓度计算值和对应实测值吻合较好。

图5-112 中性物质在弯道中图 (a),图 (b)和图 (c)连续排放的平均浓度验证

图5-113 试验深槽以及垂向计算网格示意图

图5-114 计算深槽水流结果与Van Rijn (1987)试验及Alfrink (1983)计算结果比较

5.5.5.3 深槽水流、 泥沙计算验证

Van Rijn (1987)在进行了多组不同尺寸和流量的水流跨过深槽的水槽实验。深槽几何尺寸见图5-113。试验进口平均流速为0.395m/s,水深为0.204m,床面泥沙平均高度为0.006m,有效的粗糙高度为0.02m。深槽水流测量位置如图5-113中标注。模型的计算网格为(180×12×22),中间纵向立面计算网格分布如图5-113所示。试验测量了各个断面的纵向流速,紊动动能和切应力分布。同时,Alfrink (1983)对这一试验也进行了模拟计算,结果如图5-114 所示。同时,根据三维模型计算河床摩阻流速:U*=

摩阻流速与实测结果比较见图5-115。可见,本计算和试验流场资料以及他人的计算结果都比较吻合。对于泥沙计算的验证,仍选用Van Rijn (1987)深槽系列实验。深槽几何尺寸同前面,河床表面铺设的泥沙特征参数为:D50=0.16mm,D90=0.2mm。床沙厚为0.2m,进口悬沙输沙率为0.03 kg/sm,推移质输沙率为0.01kg/sm。

图5-115 模型计算深槽底部摩阻流速和实测资料等的验证比较

图5-116 模型计算泥沙浓度分布和Van Rijn (1987)实测资料的比较

泥沙浓度计算与实测的比较 (图5-116)可见,含沙量分布在深槽前后边坡和实测值吻合良好,而在深槽中部悬沙浓度计算值与实测值相差较大,误差达到50%左右。水流经过深槽时,由于断面增大、流速减小,河床发生淤积,不断的淤积沙洲随着过程的发展向下游传递;在深槽恢复正常水深后,由于流速减速加大,深槽后部的边坡发生一定冲刷。图5-117给出了在T=0,2,4,6,8,10,12,15h时计算地形的分布情况和T=15h Van Rijn 实测的地形资料及Van Rijn的立面二维模型计算的地形的对比。模型能模拟出深槽地形随时间的变化过程和实测地形比较吻合。

图5-117 模型计算河槽冲淤过程和Van Rijn (1987)实测和计算结果的比较