河势贴体河道平面二维正交网格生成方法的研究及应用-长江科学院院报.pdfVIP

第1 8 卷第4 期 　江　科　学　院　院　报 V o .18 N o.4 2 0 0 1 年 8 月 Journal of Yangtze River Scientific Research Institute Aug.2 0 0 1 　　文章编号:1001-54 85 (200 1)04-0014-04 河势贴体河道平面二维 正交网格生成方法的研究及应用 董耀华 ( 江科学院河流研究所, 湖北 武汉　4300 10 ) 摘要:基于河势概念和Hermite 三次插值函数, 提出了河势贴体河道平面二维正交四边形网格的生成方法;与边界 层坐标系下水深平均流体力学控制方程和SIM P LER 算法联合使用, 建立了河道平面二维数学模型;进行了葛洲坝 枢纽至磨盘溪河段二维网格生成及水流泥沙数学模型的实例研究;讨论了该网格生成方法的主要特征。 关　键　词:正交网络;数学模型;河势;Hermite 函数;边界层;SIM P LER 算法 中图分类号:O242 　　　文献标识码:A 网格生成方法。采用这一网格生成方法及相应的河 道二维数学模型, 不仅可以进行天然河道大多数边 1　河势贴体网格 界条件下的水深平均平面二维水流、泥沙及温排和 河道平面二维数学模型网格生成方法研究中两 污排等计算分析及应用研究, 而且还为 河段河道 个关键问题是:(1)网格与河道拟合的贴体问题;(2) 二维及一、二维嵌套数学模型的研究和应用提供了 二维网格、控制方程和数值方法三者之间的匹配问 新思路。 题。 目前, 常用的河道二维正交网格生成方法是边 [1] 界拟合坐标系方法, 即河道T hompson 法 , 它主要 2 　网格生成方法 是通过物理平面(天然河道平面)与变换平面(数模 计算平面)之间Poisson 方程边值问题数值解实现二 [ 2] 张瑞瑾教授 认为:河道水流的流态(或河势) 维正交网格的生成。 具有很广的含义, 一切标志河道水流总体倾向的现 ξ+ξ=P (ξ, η) 象, 都被纳入这一概念之中。将河流动力学中基本 x x yy 　　(x , y )∈ D , (1) η +η = Q (ξ, η) 和核心的河势概念, 引入到河道二维网格生成方法 xx yy 式中:x , y 为物理平面网格坐标;ξ, η为变换平面 及其河道二维数学模型的研究之中, 具有更加坚实 网格坐标;D 为物理平面区域。 的理论依据和物理基础。 河道Thompson 法存在的主要问题有:(1)二维 河势贴体网格生成方法的总体思想为:放弃传 网络主要考虑与河道岸线的拟合, 多数情况下与河 统的网格生成方法中严格要求网格与河道岸线相拟 [3] 道的河势或主流线之间的拟合不理想;(2 )复杂洲、 合的思路, 采用Hermite 三次插值函数 , 生成河道 滩及岸线河道边界、非恒定流动岸和数值求解引起 沿程纵向与河势或主流线相拟合的河势拟合线(曲 的动边界等情况下, 网格与河道岸线之间的拟合同 线), 并使得河宽方向横向网格线(直线)与网格控制 样会出现实际偏离;(3