河势贴体河道平面二维正交网格生成方法的研究及应用.docVIP

河势贴体河道平面二维正交网格生成方法的研究及应用 目录 TOC \o 1-9 \h \z \u 目录 1 正文 2 文1：河势贴体河道平面二维正交网格生成方法的研究及应用 2 1 　河势贴体网格 2 2 　网格生成方法 3 3 　河道二维数学模型 5 二维水流泥沙控制方程 5 4 　应用实例 6 5 　讨　论 7 文2：河道平面二维水沙数学模型的有限元方法 8 1 基本方程 9 (3)  9 (5)  10 (7)  10 2 　有限元方程 10 （10） 11 (11) 11 (12) 11 (13) 11 (14) 11 (15) 11 (16) 11 (17)  11 3　数值解法 12 3 　数模计算有关问题处理 13 三角形网格划分 13 5 　计算实例 14 6　结论 16 参考文摘引言： 17 原创性声明（模板） 18 文章致谢（模板） 18 正文 河势贴体河道平面二维正交网格生成方法的研究及应用 文1：河势贴体河道平面二维正交网格生成方法的研究及应用 1 　河势贴体网格 河道平面二维数学模型网格生成方法研究中两个关键问题是： (1) 网格与河道拟合的贴体问题； (2) 二维网格、控制方程和数值方法三者之间的匹配问题。目前，常用的河道二维正交网格生成方法是边界拟合坐标系方法，即河道 Thompson 法 [1] ，它主要是通过物理平面 ( 天然河道平面 ) 与变换平面 ( 数模 计算 平面 ) 之间 Poisson 方程边值问题数值解实现二维正交网格的生成。 式中： x ， y 为物理平面网格坐标；ξ，η为变换平面网格坐标； D 为物理平面区域。 河道 Thompson 法存在的主要问题有： (1) 二维网络主要考虑与河道岸线的拟合，多数情况下与河道的河势或主流线之间的拟合不理想； (2) 复杂洲、滩及岸线河道边界、非恒定流动岸和数值求解引起的动边界等情况下，网格与河道岸线之间的拟合同样会出现实际偏离； (3) 控制方程变换引起的数值求解困难和数值诱发耗散问题； (4) 对于常见的宽、窄相间的河道平面形态，二维网格不均匀间距可能导致的数值计算精度问题； (5) 与河道一维断面数学模型的嵌套和联合存在接口困难等等。 为了尽可能地避免河道 Thompson 法的上述问题和困难，本文提出了河势贴体河道平面二维正交网格生成方法。采用这一网格生成方法及相应的河道二维数学模型，不仅可以进行天然河道大多数边界条件下的水深平均平面二维水流、泥沙及温排和污排等计算分析及应用研究，而且还为长河段河道二维及一、二维嵌套数学模型的研究和应用提供了新思路。 2 　网格生成方法 张瑞瑾教授 [2] 认为：河道水流的流态 ( 或河势 ) 具有很广的含义，一切标志河道水流总体倾向的现象，都被纳入这一概念之中。将河流动力学中基本和核心的河势概念，引入到河道二维网格生成方法及其河道二维数学模型的研究之中，具有更加坚实的理论依据和物理基础。 河势贴体网格生成方法的总体思想为：放弃传统的网格生成方法中严格要求网格与河道岸线相拟合的思路，采用 Hermite 三次插值函数 [3] ，生成河道沿程纵向与河势或主流线相拟合的河势拟合线 ( 曲线 ) ，并使得河宽方向横向网格线 ( 直线 ) 与网格控制断面相吻合，从而构造出平面二维正交四边形网格。由此生成的二维网格，一方面避开了河道复杂和变动的洲、滩及河岸岸线，另一方面体现了控制河道水流运动的河势概念。 Hermite 插值函数 Hermite 插值函数不仅要求插值节点上的函数值相等，而且还要求节点处一阶甚至高阶导数相等。两个插值节点情况下的 Hermite 三次插值函数可表述如下 [3] ： 需要满足的节点条件为 函数表达式为 系数表达式为 式中： H 3 (x) ， H 3 ′ (x) 为 Hermite 三次插值函数及其导数；α，β为插值系数； x ， y ， m 分别为插值节点的坐标、函数和一阶导数。 网格生成步骤 河势贴体河道平面二维正交四边形网格生成方法包括如下三个步骤。 确定网格控制断面和节点 选取研究河段的进出口断面、河段内水位 / 水文站点或测流断面、河势控制断面以及需要重点研究河段的控制断面等作为网格控制断面；在所选取的网格控制断面上确定网格控制节点，这些节点可以任意选取在控制断面的左右岸、深泓点、主