具有复杂计算域和地形的二维浅水流动数值模拟.pdfVIP

水 利 学 报 2005 年4 月 SHUILI XUEBAO 第36 卷 第4 期 文章编号: 0559-9350(2005)04-0439-06 具有复杂计算域和地形的二维浅水流动数值模拟 王志力，耿艳芬，金生 (大连理工大学 海岸和近海工程国家重点实验室，辽宁 大连 116024) 摘要：本文采用非结构化网格的有限体积方法，对具有复杂计算域和地形的二维河道进行数值模拟。采用Roe格式 的近似Riemann解计算界面通量，对地形变化源项采用特征分解平衡界面通量以保证格式的和谐性，对摩擦力源项 采用隐式或半隐式求解以增加格式的稳定性。通过与混合流算例理论解比较，验证了格式是和谐性的并具有良好 的间断捕捉能力与稳定性。应用此方法对地形复杂，存在多处险滩和深潭的某河段进行了实际模拟，数值计算结 果和模型试验结果吻合较好。 关键词：浅水方程；源项；Roe格式；有限体积 中图分类号：TV131.4 文献标识码：A 从20世纪50年代计算机广泛应用以来，浅水流动模拟广泛的应用到河道流场及洪水预报，溃坝决堤， 分滞蓄泄防洪规划，河道整治等方面[1～7] 。在进行河道，河口等实际的水流二维流场模拟时，由于地形、 计算域通常非常的复杂给数值模拟带来了很大的困难。为了克服这些困难，很多的计算水力学和计算流体 力学的学者们提出多种的解决方法，如采用贴体坐标[4]、非结构化网格[2，5，6]等，对于计算域变化问题提出 了动网格或变网格等方法。针对复杂计算域和地形，本文采用非结构化网格的有限体积法离散。空间离散 采用在空气动力学和计算浅水动力学中得到广泛应用的Godunov格式[6～13]，用Roe格式的近似Riemann解计算 界面通量。用Roe格式计算通量时，简单的方法处理源项会破坏格式的和谐性并影响稳定性，为避免源项 破坏格式和谐性很多学者提出了不同的方法处理源项，如平底模型，改变计算通量的水深，将底坡源项设 为0；Smolarkiewicz应用MPDATA方法处理源项[13]；Jenny[14]用特征分解的方法分解源项；此外还有 Glaister[15]，Bermudez[16]，Hubbart[17]以及Tomas[18]等也都提出各自的方法。本文采用特征分解平衡通量的 方法处理底坡源项，保证格式的和谐性。对摩擦源项采用隐式、半隐式离散，增加格式稳定性。 1 浅水方程 对于平面大范围的自由表面流动、水深尺度远小于平面尺度、无明显垂直环流、垂向流速小的浅水流 动，可以引入静水压力假设，并沿水深方向进行积分来简化方程，简化后的方程即为二维浅水方程。其守 恒形式为 U +E(U) +H(U) =S(U) (1) t x y ( ·) ，( ·) ，( ·) 分别为对时间，空间平面x，y方向的偏导数，式中 t x y 收稿日期：2004-04-07 作者简介：王志力(1977-)，男，河南商城人，博士生，主要从事计算水力学研究。E-mail:Wang-Z-L@ 1 水 利 学 报 2005 年4 月 SHUILI XUEBAO 第36 卷 第4 期   h hu  hv  0