二维水流数学模型mike21在潮流场数值模拟中的应用.docxVIP

二维水流数学模型mike21在潮流场数值模拟中的应用 目前，平面二维水流数学模型以垂直面上的平均流量因素为研究对象，平面水流场和详细测量的变化被广泛应用于建模中。20世纪90年代以来,国际上出现了不少成熟的平面二维水力学模型。丹麦水力学研究所开发的平面二维数学模型MIKE21,是其中应用较为广泛的商业模型之一。该模型可用于模拟河流、湖泊、河口、海湾、海岸及海洋的水流、波浪、泥沙及环境场,可为工程应用、海岸及规划提供完备、有效的设计条件和参数;该软件的高级图形用户界面与高效计算引擎的结合,使其在世界范围内成为很多专业河口海岸工程技术人员不可缺少的工具,并曾在丹麦、埃及、澳洲、泰国等许多国家和地区得到成功应用。如今,该软件已在国内的一些大型工程中得到广泛应用。如:长江口综合治理工程、杭州湾数值模拟、南水北调工程、重庆市城市排污评价、太湖富营养模型、香港新机场工程建设、台湾桃园工业港兴建工程等。 1 基于dhi的二维表面流模型 MIKE21是丹麦水力研究所(简称DHI)开发的系列水动力学软件之一,属于平面二维自由表面流模型。DHI采用MIKE21作为研究手段,不断在应用中发展和改进该软件。20多a来,MIKE21在世界范围内大量工程应用经验的基础上持续发展,其在平面二维自由表面数值模拟方面具有强大的功能。 1.1 评估功能的实现 1)用户界面友好,属于集成Windows图形界面。 2)具有强大的前、后处理功能。在前处理方面,具有根据地形资料进行网格划分功能;在后处理方面,具有流场动态演示及动画制作、计算断面流量、实测与计算过程验证、不同方案比较等分析功能。 3)热启动设置功能。用户因各种原因暂时中断了MIKE21模型,当再次开始计算时,只需将上次计算设置的热启动文件调入便可继续计算,极大地方便了计算时间有限制的用户。 4)干、湿节点和干、湿单元的设置功能,能较方便地进行滩地水流的模拟。 5)强大的卡片设置功能,可以进行多种控制性结构的设置,如桥墩、堰、闸、涵洞等。 6)多种类型水边界条件的定义功能,如流量、水位或流速等。 7)广泛应用于二维水力学现象的研究功能。可对潮汐、水流、风暴潮、传热、盐流、水质、波浪紊动、湖震、防波堤布置、船运、泥沙侵蚀和输移及沉积等进行研究,被推荐为河流、湖泊、河口和海岸水流的二维仿真模拟工具。 1.2 迭代参数的测量 1)数值参数:主要是方程组迭代求解时的有关参数,如迭代次数及迭代计算精度等。 2)物理参数:主要有床面阻力系数、动边界计算参数以及涡动粘性系数等。 2 维自由表面流模拟 MIKE21软件中的水动力学模块(HD模块)是其最核心的基础模块,可以模拟因各种作用力作用而产生的水位和水流变化及模拟任何忽略分层的二维自由表面流。由于其包括了广泛的水力现象,因此,该模块为泥沙传输和环境模拟提供了水动力学的计算基础。 2.1 深水流动质量及方程 对于水平尺度远大于垂直尺度的情况,由于水深、流速等水力参数沿垂直方向的变化比沿水平方向的变化要小的多,因此,将三维流动的控制方程沿水深积分,并取水深平均,可得到沿水深平均的二维浅水流动质量和动量守恒控制方程组。其连续性方程、X和Y方向动量方程,可分别表示为 式中H为水深,H=h+ζ,其中ζ、h分别为水位和水深;p、q分别为x、y方向上的流通通量;C为谢才系数;g为重力加速度;f为科氏力系数;ρ为水的密度;W、Wx、Wy为风速及在x、y方向上的分量;fw为风阻力系数;τxx、τxy、τyy为有效剪切力分量。 2.2 中心差分模式 采用隐式交替方向(ADI)技术对潮流模型质量和动量方程进行离散,所得的矩阵方程用追赶法求解,各微分项和重要系数均采用中心差分格式,防止离散过程中可能发生的质量和动量失真及能量失真,Taylor级数展开的截断误差可达到二阶至三阶精度。模型网格布置,见图1。 2.2.1 连续方程的离散格式 2.2.2 风阻力系数的计算 对动量方程逐项给出离散格式,此处仅给出x方向动量方程格式,y方向动量方程离散格式类似,不再列出。 1)时间偏导项,可表示为 2)重力项,可表示为 3)对流项,可分别表示为 其中:向下追赶时a=n+1,b=n;向上追赶时a=n,b=n+1,表达式可写成 4)底床阻力项,可表示为 6)风阻力项中各值已知,可用Smith和Banke方法计算风阻力系数,其计算式为 其中f0=0.000 63,W0=0,f1=0.002 6,W0=30;(16) 7)柯氏力项,可表示为 式中q*同式(12)。 2.3 有限差分法 2.3.1 时空上积分 模型采用ADI逐行法(Alternating Direction Implic)对质量及动量方程分别进行时空上的积分,每个方向及每个单独的网格线产生的方程矩阵用追赶法(Double Sweep)求解。各个