mike21软件在数值模拟中的应用.docxVIP

mike21软件在数值模拟中的应用 20世纪末，全国各流域委员会和大学联合研究了平面三维数学模型。当时，二维平面数学模型中使用的网格是正矩阵网络。适用于大范围和宏观区域的洪水模拟。由于小范围和宏观区域的洪水模拟无法正确概括，因此无法正确总结边界。如在弯曲河道上修建一座桥梁, 正交矩形网格对圆柱形桥墩只能用一个矩形网格概化, 对弯曲的河道也无法准确概化出河道边界。传统的概化方法是缩窄桥墩位置的河道过水断面或者将圆形桥墩用加密的单个网格近似模拟, 但上述概化方法只是把圆形桥墩当做正交的方形桥墩来模拟, 相同墩厚的方形桥墩的阻水影响比圆形桥墩大得多, 此种概化方法与圆形桥墩局部水流的实际情况有所出入, 无法准确模拟圆形桥墩的绕流特点。 随着城市发展步伐的加快, 涉河工程越来越多, 圆形桥墩以其阻水作用小、抗冲性强等优点被越来越多的使用, 城市的现代化发展也对圆形桥墩阻水的计算手段和模拟精度提出了更高的要求。 丹麦水力研究所开发的平面二维数学模型MIKE21是一款功能较为强大、应用较为广泛的商业模型, 它采用非正交曲线网格, 忽略了垂向水流加速度, 以垂向平均的水流因素为研究对象, 模拟计算海洋、湖泊、河道、蓄滞洪区的流场、流速、水位的变化。该模型曾经在丹麦、埃及、澳洲、泰国等国家及中国的香港、台湾地区得到成功应用。目前其在中国诸多大型工程中得到广泛应用, 如:长江口综合治理工程、南水北调工程、天津子牙新河口和独流减河口防洪影响论证、杭州湾数值模拟、重庆市城市排污评价、太湖富营养模型、香港新机场工程建设、台湾桃园工业港兴建工程等。 1 生长二维水力学模型的功能 MIKE21是丹麦水力研究所 (简称DHI) 开发的系列水动力学软件之一, 属于平面二维自由表面流模型。DHI采用MIKE21作为研究手段, 在应用中发展、改进和推广该软件。20多年来, MIKE21在世界范围内大量工程应用经验的基础上持续发展起来, 在平面二维自由表面流数值模拟方面具有强大的功能。 (1) 用户界面友好, 属于集成的Windows图形界面。 (2) 具有强大的前、后处理功能。在前处理方面, 能根据地形资料进行非正交矩形和三角形网格剖分, 将散点数据插值至网格中, 处理迅速;在后处理方面, 能够自动生成流场、流速、水位等值线、动态演示和动画生成, 并可自动处理不同时刻的最高、最低和平均的水位和流速等值线。 (3) 当用户因各种原因需暂时中断MIKE21运行时, 可将格式为dfsu的结果文件作为初始条件赋值给模型, 截取边界条件后, 继续进行计算。 (4) 能进行干、湿节点和干、湿单元的设置, 用户可自定义干湿节点转换的精度, 能较方便地进行平面水流的模拟。 (5) 可以进行多种控制性结构的设置, 如桥墩、堰、闸、涵洞等。 (6) 可以定义多种类型的水边界条件, 如流量、水位或流速等。 (7) 可广泛地应用于二维水力学现象的研究, 如潮汐、水流、风暴潮、传热、盐流、水质、波浪紊动、湖震、防浪堤布置、船运、泥沙侵蚀、输移和沉积等, 被推荐为河流、湖泊、河口和海岸水流的二维仿真模拟工具。 该模型也存在一定的局限性: (1) 三角网格剖分时, 无法进行过渡剖分, 如从河流入海口到远海实现网格剖分从小到大渐变是无法实现的。 (2) 模型是在许多水力条件假设的情况下进行的, 如垂向的水流加速度忽略不计, 不适用于急流和三维水流模拟。 2 基本方程和解算方法 2.1 合自由表面模拟技术 对于桥墩阻水引起的壅水影响, 采用较为成熟的二维紊流雷诺平均应力模型方程结合自由表面模拟技术进行。模型主要由以下方程组成:描述水流运动的方程、描述紊动应力模型方程、描述自由表面方程。 在二维模型水流模拟计算中, 水流阻力主要是通过河床的糙度 (床面的阻力高度) 、水流强度和断面形态综合考虑进行计算。 2.2 数值方法选取 本次利用MIKE21 FM 非结构网格模型进行模拟计算。非结构网格模型中采用的数值方法是单元中心的有限体积法。控制方程离散时, 结果变量U、V位于单元中心, 跨边界通量垂直于单元边。有限体积法中法向通量通过在沿外法向建立单元水力模型并求解一维黎曼问题而得到。 3 模拟段 虚拟一条长500 m, 上口宽100 m、下口宽70 m的梯形断面河道, 边坡1∶3, 下游河底高程0 m, 上游河底高程0.1 m, 纵坡1/5 000, 河道内共布设横纵各4排圆柱形桥墩, 墩厚2 m, 跨度20 m, 前后排桥墩墩距8.5 m。 模拟时段为2011年1月1日0∶00至2011年1月1日5∶00, 时间步长为30 s, 模拟步数为600。模型的构建过程如下。 3.1 计算方法的稳定性 计算区域的网格剖分是河道水流模拟的基础性工作。网格剖分的型式决定了采用的计算方法和模型的概化准确度, 网