MIKE21 原理及快速入门.docx

1.1 MIKE 21简介MIKE 21是一个专业的工程软件包，用于模拟河流、湖泊、河口、海湾、海岸及海 洋的水流、波浪、泥沙及环境MIKE 21为工程应用、海岸管理及规划提供了完 备、有效的设计环境模型利用二阶精度的有限体积法对动态流的连续方程和动量 守衡方程求解最近几年发展起来的非结构网格模型，给用户提供了最大程度的灵活性，用户可以 在模拟区域自由选择计算单元网格的大小，比如，可以在工程区域/感兴趣区域使 用较小的网格单元，而在工程远区使用单位较大的网格单元，其结果是整体计算单 元较之传统的矩形网格大为减少，从而引起计算速度的提升MIKE 21水动力模拟由于各种作用力的作用而产生的水位及水流变化它包括了广 泛的水力现象，可用于任何忽略分层的二维自由表面流的模拟HD模块是MIKE 21软件包中的基本模块，为泥沙传输和环境水文学提供了水动力学的计算基础模型特点MIKE 21经过了超过25年的连续开发和改进，结合了世界上众多用户的使用经验DHI—直使用MIKE 21进行科研工作，同时不断地对模型进行一系列的开发和改 进MIKE21系统包括了如下四种主要建模网格： 单一网格：这是一种传统的矩形模型，是将研究区域划分成同一大小的矩形网格， 网格的大小（分辩率）由模拟区域大小及具体应用决定，网格越小计算精度越高， 但耗时越长。

嵌套网格：这也是一种矩形模型，只是在同一模型中可以有多种网格大小，在大网 格模型中可以嵌套小网格模型曲线网格：网格呈四边形或近似矩形，主要适用于蜿蜒河段的水动力学计算和河床 演变分析非结构网格：这是一种三角形/四边形混合网格，采用单元中心的有限体积解法 该网格能够很好地模拟弯道或水上结构物周围区域的流场所有的网格模块都支持先进的用户界面和高效的前后处理工具，包括： 网格地形生 成器、数据管理工具、结果浏览器等最近几年发展起来的非结构网格模型，给用户提供了最大程度的灵活性，用户可以 在模拟区域自由选择计算单元网格的大小，比如，可以在工程区域/感兴趣区域使用较小的网格单元，而在工程远区使用单位较大的网格单元，其结果使得整体计算 单元较之传统的矩形网格数量大为减少，从而引起计算速度的提升MIKE 21软件包中的水动力学模块（HD模块）可以模拟由于各种作用力的作用而 产生的水位及水流变化它包括了广泛的水力现象，可用于任何忽略分层的二维自 由表面流的模拟并且HD模块是MIKE 21软件包中的基本模块，它为泥沙传输和 环境模拟提供了水动力学的计算基础对于洪水计算而言，MIKE 21模型具有以下主要特点：MIKE 21 FM的方程离散方 法采用单元中心的显式有限体积法求解，保证了水量和动量在计算域内的守恒； 模型采用非结构三角形/四边形网格元，更有利于拟合复杂边界线，比如堤防； 利用干湿网格判断法处理潮滩移动边界，方便快捷。

支持并行计算，能够大幅度减 少计算时间尤其是对于几百公里长的河流，计算时间动辄几天，使用并行计算具 有非常大的时间优势控制方程对于水平尺度远大于垂直尺度的情况，水深、流速等水力参数沿垂直方向的变化较 之沿水平方向的变化要小得多，从而将三维流动的控制方程沿水深积分，并取水深 平均可得二维浅水控制方程组连续性方程：dh dhu dhv + + 二 hSdt dx dyX方向动量方程:dhu dhu2 dhvu2 dq h dp+ + 二 f vh — gh — 4 —d t d x d y d x p d xgh2 6p t t 1 + ―s^ — ―b^ ——2 p dx p p p0 0 0+ —(hT )+UhT)+ hu S dx xx dx xy0ds )+ 塔 dx丿Y方向动量方程:dhv dhuv dhv2 dn h dp=——f uh ——gh —— a —dy p dyds、yx + yydy dx 丿+ + dt dx dygh2 dP+T莎石 P P P0 0 0 0+ — ^hT)+2(hT Lhv S dx xy dy yyTsy —— by ——1 (ds方程中t为时间；x、y、z为右手Cartesian坐标系；"为水面相对于未扰动水面的高 度即通常所说的水位；h为静止水深；u、v、w分别为流速在x、y、z方向上的分 量；pa为当地大气压；P为水密度，P0为参考水密度；f = 2°sin申为科氏力参数（其中° = 0.729xlO-4s——1为地球自转角速率，e为地理纬度）；fv和和为地球自 转引起的加速度；Sxx、U、'yx、'yy为波辐射应力分量；Jx、Ty、［x、［y为水 平粘滞应力项，S为源汇项，（Us，I）源汇项水流流速。

边界条件开边界：n =n C)rrU = u (t)rv = v (t)rrn r、Ur、vr为开边界r上已知潮位、流速过程 闭边界：U = 0 或 v = 0 初始条件：n(x, y,t)| =n (x, y) / 0ug y, t = 0vlx, y, t = 0t=t00为计算初始时刻潮位空间分布函数数值方法模型采用非结构网格中心有限体积法进行求解，其优点为可拟合复杂地形，计算速度较快 对计算区域内滩地干湿过程，采用水位判别法处理，即当某点水深小于一浅水深 dry （如0.1米） 时，令该处流速为零，滩地干出；当该处水深大于£flood （如0.2米）时，参与计算，潮水上滩亚网格水工建筑物模拟水平尺度上的水工结构物通常较模型计算用的网格尺度小很多，因此结构物的影响 通常使用亚网格技术来处理模型中包含堰，涵洞，闸门，桥墩，涡轮机等五种不 同的结构物模拟闸门在模型范围中被定义为由点组成的断面（至少两个点是必需 的）,断面在数值模型的计算中定义为邻近两单元的一边闸门的运算是由一个控制 因子来主导如果控制因子是1, 闸门为开，如果控制因子是零，闸门为关。