干湿边界的斜对角笛卡尔方法在平面二维水沙数学模型中的应用.pdfVIP

泥 沙 研 究 年 月 第 期 干湿边界的斜对角笛卡尔方法在平面二维水沙数学模型中的应用 刘斌 方红卫 段杰辉 清华大学水沙科学教育部重点实验 北京 中国水利水电科学研究院 北京 摘要 复杂边界 尤其是发生干湿变化的动边界 一直是二维水沙数学模型计算的难点 也是控制计算收敛和 精度的关键 本文主要介绍在二维浅水方程和泥沙数学模型中采用斜对角笛卡尔方法结合干湿边界点的变 化来处理复杂边界问题 并将这一方法应用到北京首都国际机场扩建工程跑道高程设计和华能丹东电厂海 水蓄水库泥沙淤积的模型计算中 计算表明 该方法拟合精细 处理简单 计算高效 与实测结果吻合较好 关键词 干湿边界 斜对角笛卡尔方法 二维水沙数学模型 中图分类号 文献标识码 文章编号 前言 浅水运动方程是目前平面二维水流计算的主要工具 如计算河口海岸 漫滩 分洪区和湖泊水库等 处的水流 在处理有以上问题的各类实际工程的计算中 复杂边界是计算的主要难点之一 如果水位还 随水流过程发生变化 边界将不仅复杂而且可变 对这种发生干湿变化的动边界的模拟要求就更高 因 此如何准确有效地对复杂边界和干湿边界进行处理是浅水运动方程应用的关键 对于复杂边界的处理 已经有很长的研究历史和大量的研究成果 目前常用的方法主要有贴体坐标 法 非结构网格法 多重网格法以及直角坐标系统下的锯齿法 插值法 混合法等 贴体坐标法由 提出 该方法的主要优点是用曲线坐标逼近复杂边界 在边界模拟上具有较高的精度 但是该方法在较复杂的边界上不易形成正交贴体网格 如果是非正交的曲线坐标系 就会产生大量的交 叉导数项 给求解带来很大的困难 非结构网格法使用三角形和四边形的计算单元 使得在复杂区域形 成计算网格更加灵活和方便 但是往往要建立一套标明网格节点之间关系的复杂系统 因此在编程计算 处理时比较麻烦 特别是针对自由面流动和干湿边界变化的问题时 这种困难将显得更加突出 多重网 格法则由于各网格系统之间的相对独立性 难以保证整个系统的守恒性 传统的笛卡尔坐标方法被应 用的历史最长 由于其简单 易形成计算网格 正交等特点 在工程界应用得最为普遍 但在边界复杂的 地方容易形成锯齿状边界 导致边界的拟合精度不高 因此该方法在大尺度计算中应用较多 对于干湿边界问题的模拟 许多数学模型也做出了探索与尝试 早期的模型将水陆交界处设置成 固壁边界 完全忽略水位升降时发生干湿变化的区域 也有在初始化时将全计算域都赋上一薄层水 而 这种方法往往会导致错误的流动 进一步采用有限差分法和有限元法的数学模型 使用固定的计算网 格监控计算节点的干湿性来跟踪干湿边界 另有一些模型使用了随干湿边界移动的计算网格 等 尽管这种方法看起来很美观 计算却很费力 而且对复杂地形处的水流模拟得也不够充分 因为边界每动一次 计算网格都必须重新生成 而且经常在水位升高时要添加计算节点 退水时又要删 收稿日期 基金项目 国家自然科学基金委杰出青年基金资助 作者简介 刘斌 女 黑龙江人 清华大学水利水电工程系硕士研究生 主要从事水力学及河流动力学数值模拟方面的 研究 减 才能保证网格不扭曲变形 和 对此法做了改进 将始终处于湿状态的网格固 定 对发生干湿变化的网格点再进行加减 然而用有限元法离散浅水方程只能保证整个计算域的质量 守恒 却不能确保每个单元内或每个节点上的守恒性 一些研究者还将有限体积法和固定网格追踪格 式结合起来模拟发生干湿变化的浅水区域 魏文礼和金忠青 汤立群 等人采 用 冻结法 动边界技术进行具有复杂边界的二维模拟 较好地解决了边界随水位的变动问题 本文的研究是在处理复杂边界斜对角笛卡尔方法 黄炳彬等 的基础上 进一步集成由于水位变化而使计算边界发生变化的干湿边界处理技术 这两种技术的结合不 仅保留了锯齿法正交 易形成计算网格 计算效率高等特点 同时又能有效地