《FVCOM模型在潜堤附近水流运动模拟中的应用.》.pdf

第 35卷第 4期 水 道 港 口 V01．35 No．4 2014年 8 月 JournalofWaterwayandHarbor Aug．2014 FVCOM模型在潜堤附近水流运动模拟中的应用 吴毓儒 ，张庆河 ，严 冰 ，杨 华 (1．天津大学 水利工程仿真与安全国家重点实验室，天津 300072；2．交通运输部天津水运工程科学研究所 工程泥沙交通行业重点实验室，天津 300456) 摘 要：采用FVCOM中的导堤一丁坝模块 ，对上挑形淹没丁坝附近的水流进行数值模拟 ，并与已有物理 模型试验观测到的自由水面线和流速分布资料进行了比较。结果表明，FVCOM模型模拟结果与试验结 果有较好的一致性，模型能够合理描述潜堤顶部过水现象，反映潜堤附近三维环流特征。 关键词：淹没丁坝 ；窄堤算法 ；自由水面；三维流场；有限体积法 中图分类号：TVI41；O242．1 文献标识码：A 文章编号：l0O5—8443(20l4)04—0325一O6 潜堤 (本文中特指淹没式丁坝)是河流及海岸工程中常用的水工建筑物，可以起到调整水流方向、消减 波高、阻挡近底泥沙运动以减少航道淤积或减弱堤内水流促进泥沙落淤等作用。例如，黄骅港整治工程…、 长江 口深水航道整治工程均采用了一定形式的潜堤 。因此，掌握潜堤附近的水流与泥沙运动规律，对于实 现潜堤的建设 目的具有重要作用。 近年来，采用数学模型方法研究河流与海岸工程建筑物附近的流场和泥沙运动已经得到越来越广泛的 应用。在三维数学模型中，潜堤可以作为固定坡度的水下地形b，潜堤在水面变化时淹没或露出水面的情况 可通过干湿判断的动边界方法来处理。这种处理方法要求计算网格划分足够细，网格空间尺度应在几米的 量级 ，以便把潜堤顶部及边坡的地形变化描述出来。但对于实际河 口海岸问题，空间尺度动辄以几十乃至 上百千米计，上述潜堤处理方法对于三维水流、泥沙模拟的空间精度要求太高，计算量较大，在实际应用中 受到限制。 另外一种处理潜堤的方法是所谓的窄堤算法 ，这种算法忽略潜堤的宽度，将潜堤定义为网格的边界进 行计算，通过合理估算潜堤两侧计算单元之间的水流交换过程来实现潜堤的功能 。葛建忠等最近将这种 窄堤算法运用到FVCOM模型中 ，建立了基于有限体积法的非结构化网格导堤一丁坝计算模型，该模型能 够精确拟合不规则潜堤平面布置，并且能有效保证计算过程中潜堤两侧水体的质量守恒。但在文献 [5—6] 给出的流场验证算例中，基于窄堤算法的FVCOM模型只与出水堤附近流场进行了比较 ，尚未与潜堤附近的 流场进行比较。因此，本文拟利用文献中已有的上挑淹没丁坝附近水流试验数据 与FVCOM模型计算结 果进行比较，以进一步验证模型在描述潜堤周围水流运动方面的合理性。 1 模型介绍 1．1 FVCOM模型简介 FVCOM模型由美国马萨诸塞大学和伍兹霍尔海洋研究所联合开发，水平方向采用非结构化三角形网 格，垂直方向采用 坐标变换，数值方法采用有限体积法，其三维原始控制方程为 收稿 日期 ：2013—10—28；修回日期：2013—11-20 基金项 目：国家高技术发展研究863计划 (2012AA051709)；国家 自然科学基金 作者简介：吴毓儒 (1989一)，男，福建省泉州人 ，硕士研究生，主要从事海洋数值模拟研究。 Biography：WUYu一111(1989-)，male，masterstudent． 2014年8月 吴毓儒 ，等 FVCOM模型在潜堤附近水流运动模拟中的应用 327 心，其他标量定义在网格节点上。定义在形心上的M和 在动量控制体MCE中计算获得，定义在节点上的标 量则在质量控制体TCE中计算获得。MCE与TCE在非结构化网格中的定义如图i所示 ，其中，图1-a代表 一 般单元的MCE与TCE，图1一b代表潜堤单元的MCE与TCE。 在潜堤计算中，设 为潜堤的绝对高度 ，经过潜堤两侧的TCE可以分成左右两个部分分别计算通量，即 Qz鲁=一(加)一 )HfD)dy—fT(o)l