含有紊动黏性项浅水方程的数值求解 Numerical solution to the shallow water equation with diffusion motion.pdfVIP

第3期 水利水运工程学报 No．3 AND 2010年9月 HYDRO-SCIENCEENGINEER矾G Sep．2010 含有紊动黏性项浅水方程的数值求解 何杰1”，辛文杰12 (1．南京水利科学研究院，江苏南京210029；2．水文水资源与水利工程科学国家重点实验室，江苏南京 210029) 摘要：紊动黏性项对水流方程数值求解的稳定性至关重要．采用基于无结构网格的有限体积法对方程进行离 散，对流项的界面通量采用Roe格式的近似Riemann解计算，紊动黏性项采用单元交界面的平均值估算，底坡源 项则采用特征分解法处理．应用该模式模拟了圆柱绕流，模拟结果与实际情形较相符．对港珠澳大桥的海中人 工岛水域的潮流运动进行了有效模拟，结果表明该模式能够较好地模拟出水建筑物附近水域复杂的水流结构． 关 键 词：紊动黏性项；浅水方程；有限体积法；无结构网格；人工岛 中图分类号：TVl31．4 文献标志码：A 计算浅水动力学在数值离散方法上先后经历了有限差分法、有限元法和有限体积法等．有限体积法是一 种结合有限单元法改进的有限差分法，通过对每一个控制体积积分可以严格满足物理量守恒．在基于无结构 网格的有限体积法方面，浅水方程的数值求解研究经历了从一维、二维到三维的过程，其中单元问的界面通 量计算成为浅水方程求解的重点．采用数值方法模拟实际水流运动时，尤其是河口海岸地区大范围的水体运 动，往往忽略紊动黏性项在浅水方程中的作用．这主要是因为在模拟大水体运动过程中，计算网格较粗，浅水 方程中的对流项占主要作用，有无紊动黏性项模拟计算结果的差别很小．但对于小范围的水体，如模拟出水 建筑物(尺寸以米计)周围的水体运动时，紊动黏性项在水流方程中所起的作用将会突显出来．因此，在模拟 尺寸较小建筑物周围水体的运动时，需考虑有紊动黏性项时浅水方程的数值求解． 基于无结构网格离散计算区域、采用有限体积法求解浅水方程的过程中，数值通量的计算非常重要…． 的源项处理也很关键，若简单处理则可产生多变地形下静水虚假流动现象M-，较常见的做法是参照对流项 特征向量的分解形式对源项进行分裂”o．紊动黏性项中含有偏微分二次项，对计算单元积分后将面积分转 化为沿单元边界的线积分后仍含有一次微分项，这是紊动黏性项在浅水方程求解过程中不易处理的主要原 因之一．本文采用单元交界面的平均值计算该界面紊动黏性项的数值通量，对于微分项则采用类似有限元的 方法处理，从而实现带有紊动黏性项的浅水方程的数值求解．耿艳芬旧。等采用不考虑紊动黏性项的浅水方 程对桥墩等出水建筑物进行了水流模拟．本文将采用类似的计算方法结合考虑浅水方程的紊动黏性项计算 模拟出水建筑物附近的水流形态． 1 控制方程 在平面尺度远大于垂直尺度情况下，可用平面二维浅水方程描述计算水域的水流运动．考虑到水流紊动 因素的影响，水流控制方程可表示为： iou+vE：5+VE6 、(1) af 收稿日期：2009—10—09 基金项目：水利部科技创新项目(SCX2003—03) 万方数据 水 利 水 运 工 程 学 报 2010年9月 ‘2) 簟≥型： F=[d“2≥2／2]，c=[d缈：乏：／2] ^I !d=矗+，7 式中：全水深d=h+，7，h为水平面以下水深，田为水面波动 (见图1)；“，移分别表示石，Y方向的流速． 图1水深相关量 Related factorsabout