数值切割单元法在河流数值模型中的应用前景展望 Prospects of the sharpen immersed boundary method in the River Fluvial Modeling.pdfVIP

水 利 学 报 2013年2月 SHUILI XUEBAO 第44卷 第2期 文章编号：0559-9350（2013）02-0173-10 数值切割单元法在河流数值模型中的应用前景展望 黑鹏飞 ，方红卫 ，陈稚聪 ，夏建新1 2 2 1 （1. 中央民族大学 生命与环境科学学院，北京 100081； 2. 清华大学 水沙科学与水利水电工程国家重点实验室，北京 100084） 摘要：天然地表水边界复杂，且处于不断的冲淤变化中。数值切割单元法可任意追踪拟合复杂动边界，且网格生 成简单，在计算流体中具有重要的研究和应用价值。本文对数值切割单元法常用的直接插值法、虚源反射法和直 接源项法的基本原理、适用条件及关键难点解决方面等进行了系统总结。结合地表水流的动力特性，给出数值切 割法在复杂边界地表水模拟中的具体应用方法，提出在二维模拟中，可直接采用基于有限体积法的直接插值法进 行数值切割，而在三维模拟中，联合采用sigma-网格及直接源项法可同时实现模型的准确性及高效性。最后，通 过对比目前数值切割法在其它领域典型算例，提出数值切割单元法在崩岸、河工建筑和植物绕流模拟、高精度紊 流模拟等地表水模型难题方面的应用前景。 关键词：数值切割单元法；淹没边界法；河流动力；数值模型；复杂边界 中图分类号：TV143 文献标识码：A 1 研究背景 河流数值模拟中，边界拟合直接决定着河岸边界的剪切 应力和动压计算的准确性，从而决定着河岸冲刷及河流水沙 计算的准确性。但传统网格在河流边界拟合方面都存在难以 突破的缺陷。曲线贴体网格生成过程较为复杂，且很难实现 对复杂边界区域的拟合（如图1所示的长江南京段三江口窝崩 形态）。非结构网格虽可灵活逼近边界，但离散所得的代数方 程无法采用线/块迭代解法，从而降低求解的效率，更为重要 的是，动边界计算中，非结构网格无法在计算过程中自动追 踪拟合边界。 因此复杂动边界追踪拟合成为河流数值模拟的主要难题 ［1］ 图1 三江口窝崩发生前后河道形态 之一。但是，1972年以来，复杂动边界拟合一直是国际计算 ［2］ 流体力学的热点，最为典型的是IBM （immersedboundarymethod）非贴体网格方法（据不完全统计， SCI收录的IBM相关文章已达591篇，其中Journalofcomputationalphysics 近170篇）。IBM方法主要包 括经典IBM、IIM（immersedinterfacemethod）、直接源项法（directforcemethod）和数值切割网格法。其 中前3种方法的主要思想是将边界的作用转化为等效力源项进行实现，适用于流固耦合、颗粒流以及 刚体绕流的计算。而数值切割单元法不仅可实现2阶以上精度，且边界力无耗散，而且适用于高Re 流动的计算中，因此尤其适用于河床演变的研究，但是该方法尚未在河流动力模拟中得到重视。 长江三峡工程蓄水以来，水库下游河床冲刷强度加剧，由蓄水前的“冲槽淤滩”转变为“滩槽均 收稿日期：2012-03-17 基金项目：国家自然科学基金；中央民族大学“985工程”项目（MUC98507-08）；中央民族大学自主科研项目（1112KYQN30） 作者简介：黑鹏飞（1979-），陕西人，讲师，主要从事水力学及河流动力学研究。E-mail：heipf06@ — 173 — ［1］ 冲”，局部崩岸时有发生，较为典型的是长江南京河段的三江口窝崩 。2008年11月18日，长江南