龙口垂向二维流场数值模拟.pdfVIP

1 第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 龙口垂向二维流场数值模拟龙口垂向二维流场数值模拟∗ 龙口垂向二维流场数值模拟龙口垂向二维流场数值模拟 赵庚润 1 吴 卫 1 刘 桦 1 （1. 上海交通大学船舶海洋与建筑工程学院，上海，200240 ） 摘要摘要: 以长江口青草沙水库主龙口施工过程为背景，在恒定来流及出水口入水口水位差给 摘要摘要 定的情况下，应用基于 RANS 方程和 VOF 方法的 CFD 技术对龙口垂向二维水流流态进行了 数值模拟。计算结果揭示了龙口下游回流区范围、最大流速出现位置，以及龙口上下游和龙 口处自由液面的变化过程。比较了在不同龙口底槛高程的情况下龙口局部流场的变化特征。 关键词关键词: 垂直二维；数值模拟；龙口；最大流速 关键词关键词 1 引言 在围垦工程中，堵口截流工程最为艰巨。龙口水力计算的准确性，包括龙口的布置与防 护，堵口水力计算及水力条件分析，堵口的顺序与方法，将对堵口施工起到关键性的指导意 义。对施工强度安排、经济投入、进度要求都有决定性的作用。 在龙口的合拢过程中，其水流为束窄河道水流，随着戗堤进占加剧，使过水断面在水平 方向和垂直方向产生收缩，产生局部能量损失，同时，过水断面减小，流速加大，部分势能 转变为动能，形成水面跌落（贝努利水丘）出现宽顶堰流[1] 。 Chan、Huang 和 Hsieh[2]在粘 [3] 性模型下对入射波在潜坝的作用下的流动特性进行了验证。迟耀瑜等 介绍了龙口下游水力条 [4] 件对立堵截流影响的物理模型试验工作。季永兴等 对浦东国际机场围海大堤龙口进行了物理 模型试验和平面二维数值计算。 本研究利用 FLUENT 软件，结合上海青草沙水库龙口实际地形情况，对水库龙口垂向二 维流场进行了数值模拟。采用了κ − ε 湍流模型，水气两相流的 VOF 模型（Volume of Fluid ） 进行自由液面的追踪[5] 。 重点讨论不同龙口底槛高程的情况下龙口的局部流场特征。 2 数学模型 ∗ 国家自然科学基金项目资助. - 986 - 1 2 第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 数学模型的控制方程为： ∂ρ u ∂ρ u u ∂p ∂ ∂u ∂u 2 ∂u ∂ i + i j = − + [µ ( i + j − δ k )]+ (−ρ u u ) + ρ f + F (1) ∂t ∂x ∂x ∂x ∂x ∂x 3 ij ∂x ∂x i j i i j i j j i k j ∂ρ + ∂ (ρ ui ) = 0 (2) ∂t ∂x i 其中，对三维流动