邵仙水利枢纽全域流场的精细模拟.pdfVIP

刀日 Ⅶ 水～利～枢～纽一全。 域流～璺一的 ～精 ～细 邀拟一 王玲玲严忠民周春天罗龙洪 (河海大学水电学院，南京210024) 摘要 本文适过曲线同格生成技术在计算平面上生成了正交性良好的曲线网格．在该平 面用上采用交替方向隐式ADI法结台动边界处理技术，求解水深平均的圣一维南方程纽，井 运用该数学模型计算了苏北江都抽水站北侧的邵仙枢纽(包括邵仙闸、运盐闸及两闸之间的 拟新建套闸)上下游500m范围内多种工况的流场，计算得的流场形态和套闸分流比与物横 结果吻台良好．经分析比较，建议在套闸右堤堤头上游顺流向加建长约80m的隔水堤，可使 套闸引脘遵右堤堤头的横向流速分量明显减少．该数学模型还可用来计算河道、河El、大面 积海域的流场． 美键词 f圣一维南方程t 动边晃曼曲线阿倍 一入 图 1为苏北江都抽水站以北高水滴上的邵仙水利枢纽，枢纽现有邵拙闸和运盐闸各 一座，前者承担向北送水和通航的任务，后者则用于内河向长江排涝。因邵仙闸的过闸 流量较大．所以上游常发生沉船事故。为了协调航运和送水的矛盾，拟在邵仙阚和运盐 闸间再建一套闸满足航运之需。该套闸上游引航道即为邵仙闸上游引河，流遮较大，所 以需对套闸建成之后枢纽闸门区的水流流态进行深入绍致的研究。 ’‘。。●_～ 1曲线网格生成技术 由于图l所示物理区域形状复杂、极不规则。所以，在规则区域上推导得的水流运 动控制方程及其离散方法都不再适用。为能继续使用这些较为成熟的算法，就应当将物 理平面上的复杂区域转化到计算平面上，且形成规则区域，再将物理平面上的控制方程 转化到计算平面上，即以控制方程的复杂化换取计算区域的规则化、数值计算方法的简 单化。 计算网格的正交性程度对计算结果的精度影响极大，特别是对于差分类方法更是如 此，所以，正交曲线网格的生成技术已发展成为计算流体力学的一个重要分支。 利用流线和等势线正交的特性可知，在物理平面上求解一组泊松方程组 在区域内 Jf一+f一5尸(f，，，) (1) l口。+口∥；口(；．r／，J，，) ·117· -二=二‘：一=二。≮：一之沁／乒二二：二二一一 图l邵仙闸平面图 就可在计算平面上获得正交性良好的曲线网格。这里，x、y是物理平面上网格节点坐标： f、目是计算平面上网格节点坐标；P、Q是调节因子，其作用是调整物理平面上曲线网 格的形状及疏密程度。 求解方程(1)的边界条件是第一类Direhlet条件，即 J}=f(j，y) 在边界r上 Ir／=rl(x，y) 式(1)可改写成： 嘉(t菩)+专[t蒡]=。 Q’ 旦f上塑1+立f上丑1：0 (3) 其中 女： 、f毒：十毒；怛7j+吲qy=√砉蔷YVx；+； 因为物理平面上计算区域极不规则，在物理平面上求解上述方程组本身就很困难， 所以再将上述边值问题转化为计算平面上的边值问题： 锻 + +J + =O 筇 + +， + =O ，●●●●●●●●●●【 缈嚣 p∥聊 恤b $少 、¨、¨ (4) 式中 ，：一上掣些 口：一}掣 T d# 口=x；+y： ，=』；+Y； J=xfy一x_yf 求解方