平面二维非恒定流数学模型的研究及其在三峡上游引航道中的应用.pdfVIP

平面二维非恒定流数学模型研究 及在三峡上游引航道中的应用 黄国兵王才欢 郑邦民 (长江科学院水工研究所) (武汉水利电力大学) 提要本文采用沿水深平均的平面二维数学模型，分析计算了三峡水利枢纽工程因船闸充水而在上游引 航道中引起的非恒定流运动。成果表明，该数学模型的结果与物理模型试验规律一致．在精确性、收敛速 度、稳定性方面有其优越性，可供设计部门作方寨比较时参考。 关键词三峡工程 非恒定流 数值模拟 1 问题的提出 三峡工程通航建筑物采用双线连续五级船闸和垂直升船机。当船闸充水时，在上游引航 道内会产生非恒定流波动．从而影响船队在引航道内的航行和停靠．而航运部门对上游引航 道水力学指标有严格规定，要求： 船闸闸前水位最大波幅≤ojom波峰减波谷．以下同)； 升船机前水位最大波幅≤0．30m： 航道内横向流速≤O15n以； 航道口门区最大纵向流速≤20n以，最大横向流速≤O．3m／s。 为此，设计部门根据工程实际及国内数家辩研单位的建议，提出了三类不同布置方案； “小包”方案、“大包”方案及“全包”方案。所谓“小包”、“大包”及“全包”主要针对 航道右侧长隔流堤的位置而言。“小包”方案右隔流堤只把船闸上游航道包围在内：“大包” 方案则将船闸航道和升船机航道包围在一起共用一条引航道；“全包”方案则将升船机右侧 临时船闸改建的冲沙闸也包围在一起，使航道水域面积更大。作者对上述几种方案进行了数 值模拟计算，提出了研究成果报告。本文为研究成果的一部分，主要对“大包”方案不考 虑大江泄流时由船闸充水而引起的非恒定流水流运动特性进行分析论证{】]。 “大包”方案船闸与升船机共用一条引航道，引航道底部开挖高程130m，泥沙淤积控 为l：2。引航道口门处渠底宽度为220m，全航道底宽均大于220m。另外，在船闸右侧设 660m短隔流堤。船闸取水口设在短隔流堤外(简称侧向取水)，见图l。 2数学模型【2】【3】 2．1坐标变换与网格划分 采用曲线坐标I’(‘，q)代替笛卡尔坐标xJ(xJ)以适应上游引航道的不规则形状 并生成适体坐标的计算网络，控制方面即是解曲线坐标的Laplace方程： 图l “大包方案”计算区域图 V2p。去旁慨”势 ㈩ 式中矿为度量张量的逆变分量，即妒=(岛)一，而g=峨l为度量张量gi的行列式值。 在二维情况下 其中 铲艮舅 gn=『g∞g一：]土 L92lgIIjg 晶t=乓+兵 g一2：gzl：-‘+雎一g∞=毫+Z 盯于止父凹线坐称gl=-ga20 g。=一1￡] g“=[i；Z。]吉 如果函数p是几何坐标№，则对正交曲线坐标有 嚣。参+后专c瓜“，罄％g g哦：喝、‰+矗毛瓯}瓜、x： 一矗毫幅|正％增 (2l gnyg强ty xq+瓜熹t佤{瓯)y： 一压嘉(√i，√云撕=90 (3) 当‘，自变量为计算网格坐标，就可阻用(2)、(3)式求出正交曲线网格上的K值。№ 为适体坐标，P，Q为贴体坐标疏密控制函数。 2．2基本方程 由于“大包”方案上游航道底高程控制在139m，当上游水位为145m时，水深只有6m， 而整个计算水域水面宽大于250m，宽深比大于40。因此，采用了沿水深平均的平面二维数 学模型进行计算，控制方案州为： 詈+嗉(㈣+导(肋)；o