有关Mike21模拟计算的文章.pdf

丹麦 MIKE21 模型在桥渡壅水计算中的应用研究 袁雄燕 徐德龙 （长江水利委员会水文局，湖北 武汉，430010） 摘要：丹麦 DHI 公司开发的 MIKE21 数学模型属于平面二维表面流数学模型，可广泛地 应用于潮汐、水流，风暴潮等二维水力学现象的研究；本文将 MIKE21 应用于阿深北高速公 路兴建巴河大桥前后水位、流场变化情况的二维数值模拟中，以确定 MIKE21 模型在桥渡壅 水计算中的适用性；经检验，计算成果令人满意，其精度可以满足预测要求。 关键词：MIKE21；平面二维；数学模型；桥渡壅水；巴河 中图分类号：TV147 文献标识码：A 平面二维数学模型以垂线平均的水流因素作为研究对象，模拟计算河段平面流场及河床 细部的变化情况，在短河段短时期的河床细部数值模拟计算中应用广泛，如桥渡壅水问题的 研究。20 世纪 90 年代以来，国际上出现了不少成熟的平面二维商业水力学模型，丹麦水力 研究所开发的平面二维数学模型MIKE21，是其中应用较为广泛的一款商业模型，曾经在丹麦、 埃及、澳洲、泰国及中国香港、台湾等国家和地区得到成功的应用。目前该软件在国内的应 用发展很快，并在一些大型工程中广泛应用：长江口综合治理工程；杭州湾数值模拟；南水 北调工程、重庆市城市排污评价、上海太湖富营养模型；香港新机场工程建设；台湾桃园工 业港兴建工程 [1]等。 1 MIKE21 模型简介[2][3] 丹麦水力研究所（Danish Hydraulic Institute，简称 DHI 公司）是丹麦一所私营研究 和技术咨询机构，成立于 1964 年，MIKE21 是该公司开发的系列水动力学软件（DHI Software） 之一，属于平面二维自由表面流模型。丹麦水力研究所不断采用 MIKE21 作为研究手段，在 应用中发展和改进该软件，MIKE21在 20 多年持续发展和世界范围内大量工程应用经验的基 础上发展起来，在平面二维自由表面流数值模拟方面具有强大的功能： （1）用户界面友好，属于集成的 Windows 图形界面。 （2）具有强大的前、后处理功能。在前处理方面，能根据地形资料进行计算网格的划 分；在后处理方面具有强大的分析功能，如流场动态演示及动画制作、计算断面流量、实测 与计算过程的验证、不同方案的比较等； （2）可以进行热启动，当用户因各种原因需暂时中断 MIKE21 模型时，只要在上次计算 时设置了热启动文件，再次开始计算时将热启动文件调入便可继续计算，极大地方便了计算 时间有限制的用户； （3）能进行干、湿节点和干、湿单元的设置，能较方便地进行滩地水流的模拟； （4）具有功能强大的卡片设置功能，可以进行多种控制性结构的设置，如桥墩、堰、 闸、涵洞等； （5）可以定义多种类型的水边界条件，如流量、水位或流速等； （6）可广泛地应用于二维水力学现象的研究，潮汐、水流，风暴潮，传热、盐流，水 质，波浪紊动，湖震，防浪堤布置，船运，泥沙侵蚀、输移和沉积等，被推荐为河流、湖泊、 河口和海岸水流的二维仿真模拟工具。 作者简介：袁雄燕，女，长江水利委员会水文局，助理工程师。 - 1 - 该模型也存在一定的局限性： （1）MIKE21矩形网格计算模块采用矩形网格有限差分法，在海岸或防波堤等不规则边 界時，边界常被处理成齿状，计算出的结果不尽理想。但 MIKE21 的其他模块（如曲线网格 MIKE21C 模块基于贴体坐标系的有限差分法）可以解决该问题，或者可以通过细化矩形网格 来改善该问题。 （2）MIKE21矩形网格计算模块难以进行小尺寸局部水工建筑物的绕流模拟。 （3）MIKE21 模式是在许多水力条件假设的情况下进行的，如垂向的水流加速度忽略不 计，属大范围平面二维数学模型，不适用于近区三维问题或急流等流动。 MIKE21 用户界面如图 1 所示，其计算参数包括两类，一类是数值参数，主要是方程组 迭代求解时的有关参数，如迭代次数及迭代计算精度；第二类是物理参数，主要有床面阻力 系数，动边界计算参数以及涡动粘性系数等。 2 基本方程及离散格式[4] 2.1 基本方程 水流连续方程