sigma坐标系下的三维水流数值模拟探讨.pdf

doi：10．3969／j．issn．1001—9235．2010．02．00 2010年第2期 ·PEARLRIVER 人 民珠江 sigma坐标系下的三维水流数值模拟探讨 方神光，黄胜伟，崔丽琴 (珠江水利科学研究院，广东 广州 510611) 摘 要：建立了基于 sigma坐标系下的浅水三维水动力学数学模型，采用混合有限分析解法对其进行离散和求解。 分别对近海区域的表面风引起的回流、密度驱动回流和振荡潮波3种典型流态进行了数值模拟试验，并将计算结 果与理论分析解进行了对比验证。数值模拟试验准确地模拟了3种典型流态的流速和水位分布，证实了所建立的 数学模型和求解方法的可行性和有效性，为研究更为复杂的近海水动力学流态提供了一个有效的方法和工具。 关键词：风应力；密度梯度；振荡潮；三维数学模型；混合有限分析法 中图分类号：TV149．2 文献标识码：A 文章编号：1001．9235(2010)02-0003-03 河流、湖泊和近海区域的水流运动实际上是一种三维运 但由于自然界中各种影响因素的综合作用，导致流体运动状 动，但由于其水深较浅，与流场水平尺度相比，采用二维潮流 态非常复杂，同时在不同的工程研究问题中所关注的运动要 数学模型也能有效的模拟其运动流态。随着计算机容量、内 素不同，因此许多研究者以N—S模型为基础，建立了各种类 存和速度的迅速发展以及对潮流垂直结构的研究需要，20世 型的数学模型，典型的如雷诺应力模型和流函数涡量模型。 纪70年代出现 了一些三维动力学模型，如英 国海洋学家 此处建立的ro坐标系下的三维数值模型同样以该方程为基 Heaps 、Davies_2的谱模型和美国Leende~se等 的河 口和 础 ，采用静水压强分布和Boussinesq假定 ，且直角坐标系(， 海岸三维动力学模型。进入20世纪80年代后，美国Prince— y，z，t)与or坐标系( ，Y，Z，t)下的转换关系如下： ton大学的Blumberg和Mellorl4开发的ECOM和POM模型、 = ，Y=Y，= ( ，Y，ro，t)=o．H+田，t=t 德国Hamburg大学 Backhaus等 开发的HAMSOM陆架海 式中 日——当地水深H=h 叩；——水面水位；^——当地 洋动力模型等也得到了研究人员的广泛应用。三维潮流数 静水水深。ro=0时，z=钉表示水面；当ro=一1时，z 学模型中的数值模拟方法也是多种多样，按数值计算方法 = 一 h表示水下地形。 分，典型的如有限差分法、有限元法、有限体积法等；按使用 由此得到ro坐标系下的近海水域三维数学模型(略去 自 的坐标系分，有直角坐标系法和坐标系法 ，其中坐标变换法 变量上的撇号)： 以其在水深方向分层灵活且离散和求解相对其他方法容易 + + f1) 而获得了最为广泛的应用。针对这些数学模型和方法，李孟 Ot Ox Oy Oro 国和曹祖德 对海岸河 口数值模拟进行了详细的总结和评 詈+“+O。以uy+tO=击(【J)+专l(，号J)+妒1 述。 本文在流体力学原始方程N—S方程的基础上推导建立 ( )-gOx一pJ广~ft业Ox一鲁 】d+s(2 了坐标系下的近海区域三维潮流数学模型，并采用 Huang模 型 对坐标变换产生的扩散项余量进行 了简化处理。采用 詈+u+嚣+oqV=(【)+0＼NOd。vJ)+1 李炜 提出的混合有限分析法对数学模型进行了离散 ，应 用FORTRAN语言编写了计算求解程序，对封闭水体中的风 ( )-g一 [考一q,~Oy业OroJ1d+Sycs， 驱动回流、密度驱动回流和振荡潮波3种典型流态进行了数 值模拟计算，将计算结果与理论分析解进行了对 比验证 ，证 一 +十一 +十一