河势贴体河道平面二维正交网格生成方式-研究及应用

1、河势贴体河道平面二维正交网格生成方式-研究及应用            摘要：基于河势概念和Hermite三次插值函数，提出了河势贴体河道平面二维正交四边形网格的生成方法；与边界层坐标系下水深平均流体力学控制方程和SIMPLER算法联合使用，建立了河道平面二维数学模型；进行了葛洲坝枢纽至磨盘溪河段二维网格生成及水流泥沙数学模型的实例研究；讨论了该网格生成方法的主要特征。 关键词：正交网络 数学模型 河势 Hermite函数 边界层 SIMPLER算法  1河势贴体网

2、格河道平面二维数学模型网格生成方法研究中两个关键问题是：(1)网格与河道拟合的贴体问题；(2)二维网格、控制方程和数值方法三者之间的匹配问题。目前，常用的河道二维正交网格生成方法是边界拟合坐标系方法，即河道Thompson法1，它主要是通过物理平面(天然河道平面)与变换平面(数模计算平面)之间Poisson方程边值问题数值解实现二维正交网格的生成。函数表达式为式中：H3(x)，H3(x)为Hermite三次插值函数及其导数；，为插值系数；x，y，m分别为插值节点的坐标、函数和一阶导数。2.2网格生成步骤河势贴体河道平面二维正交四边形网格生成方法包括如下三个步骤。确定网格控制断面和节点选取研究河

3、段的进出口断面、河段内水位/水文站点或测流断面、河势控制断面以及需要重点研究河段的控制断面等作为网格控制断面；在所选取的网格控制断面上确定网格控制节点，这些节点可以任意选取在控制断面的左右岸、深泓点、主流点、中心点等处，所选择的网格控制节点即为数学上Hermite三次插值函数的计算节点。生成河势拟合曲线利用上述Hermite三次插值函数，可以生成一条既通过网格控制节点，又垂直于网格控制断面的河道纵向网格控制曲线。通过多次调整网格控制断面和节点，使得所生成的网格控制曲线与研究河段的河势或主流线相拟合，将最终生成的河道纵向网格控制曲线确定为河势拟合曲线。这一步是河势贴体网格生成方法的关键和核心。构

4、造平面二维网格选取和调整纵向和横向网格间距，构造由平行于河势拟合曲线的曲线簇和垂直于河势拟合曲线的直线簇(包括网格控制断面)的河势贴体平面二维正交四边形网格。    3河道二维数学模型河势贴体河道平面二维正交四边形网格生成方法与边界层坐标系下水深平均流体力学控制方程4以及合适的数值方法(如SIMPLER算法5)之间联合使用，可以建立河道二维数学模型。本文给出了河道二维水流泥沙数学模型的控制方程和数值方法。3.1二维水流泥沙控制方程采用边界层坐标系下简化的河道二维浅水控制方程和泥沙对流扩散方程4，6图1葛洲坝枢纽至磨盘溪河段河势图(1997年9月河道地形)

5、图3葛洲坝枢纽至磨盘溪河段二维水沙数模计算的流速等值线 二维水流控制方程的数值离散和迭代求解基于SIMPLER5算法，泥沙对流扩散方程的数值离散和迭代求解基于有限控制容积法5，离散方程迭代求解方法具体包括以三对角追赶法(TDNA法)为核心的逐行法和高斯塞德尔点迭代法，并配合块修正和欠松驰修正技术等。计算过程中采用了动边界模拟技术，具体处理措施包括：每次迭代根据二维网格节点的计算水深值，均要判断和区分水域和陆域节点；对于陆域节点采用边界隔墙法5处理，并让陆域节点始终保持一个较小的富余水深等。4应用实例河势贴体网格生成的主要目的是为河道平面二维数学模型的研究提供二维计算网格及离散节点。作

6、为应用实例，本文给出了葛洲坝枢纽至磨盘溪近坝河段(图1)二维网格生成的具体过程及二维水流泥沙数学模型研究的部分成果。该项研究的主要目的是分析论证葛洲坝下游壅水工程措施(如胭脂坝左汊布设潜坝)的壅水效果及其对防洪和航运的影响，而采用壅水工程措施的主要目的是为了解决三峡工程运用初期枯水位时葛洲坝下游引航道通航水深不足问题。依据上述河势贴体网格生成方法及步骤，首先选取0(坝轴线    2007-04-23        )，Y34(庙咀)，Y37(宜昌)，Y39，Y41(宝塔河)，

7、Y44，Y46，Y49和Y50(磨盘溪)共9个河段河势控制或水文/水位断面作为网格控制断面，选取9个断面的中点作为网格控制节点；然后，由Hermite三次插值函数生成了本河段的河势拟合曲线A(图1)；最后，通过确定纵向x，横向y的网格间距，构成本河段河道平面二维计算网格(图2)。最终生成的二维计算网格节点数为137×41，河势方向(x方向)网格间距为100200m，断面方向(y方向)网格间距为50m。采用生成的河势贴体河道平面二维正交网格，利用上述的河道二维水流泥沙数学模型，即可进行葛洲坝枢纽至磨盘溪近坝河段的二维水流泥沙数模计算分析研究。计算程序为自编“HELIU11”程序，有关研

8、究内容及成果请详见长江科学院“九五”三峡工程泥沙问题研究子题分报告：“葛洲坝枢纽下游近坝段整治二维水流泥沙数学模型研究”。本文仅给出了宜昌流量13500m3/s时该河段二维数模计算的流速等值线(图3)和1993年11月1日至1997年8月31日冲淤验证计算中的冲刷部位等值线(图4)。5讨论本文提出的河势贴体河道平面二维正交网格生成方法具有如下几个主要特点：(1)引入河流动力学的河势概念，凸现了河道二维数学模型研究中的主要矛盾，避开了河道岸线拟合这一复杂但相对次要的矛盾；(2)生成的二维正交网格与边界层坐标系下河道水深平均流体力学控制方程及相应合适的数值方法(如SIMPLER算法)三者之间匹配较好；(3)网格保留了河道一维断面形式