流动环境中平面负浮力排放的数值模拟.pdfVIP

水 利 学 报 2003 年07 月 SHUILI XUEBAO 第7 期 文章编号：0559-9350 (2003) 07-0058-06 流动环境中平面负浮力排放的数值模拟 1 1 杨中华 ，槐文信 (1.武汉大学 水利水电学院，湖北 武汉 430072) 摘要：建立流动环境中平面负浮力倾斜射流的二维k- ε湍流模型，采用D M Shahrahani 和J D Ditmars(1976)的试 验资料进行检验，并且对平面负浮力倾斜射流流动特性进行了数值预报。给出了射流的分区及分区的界限，同时亦 给出了涡心断面的物理量的分布及典型断面上湍动能平衡图。 关键词：负浮力射流；流动环境；平面倾斜射流；数值模拟；内部结构 中图分类号：O358 文献标识码：A 化肥厂产生的含石膏废水及海水淡化厂的高盐水通过扩散器的深水排放等流动可以看作是负浮力射 流。负浮力射流的研究对于废水排江扩散器的设计和环境影响评价具有一定意义。DMShahrabani 和 ［1］ JDDitmars 对平面负浮力射流作过实验研究，研究的主要内容是中心线上的最大上升高度和撞击点到喷 口的距离以及这两个位置的稀释度等等，并且得到静止环境中，45°倾角排放的水平稀释范围(撞击点到 喷口的距离)最大。槐文信等［２］用二阶k- ε湍流模型来模拟静止环境中倾斜负浮力射流，并且采用混合有 限分析法进行数值计算。流动环境中45°倾角排放的水平稀释范围是否最大还有待进一步研究。本文将在 文献［２］的基础上对流动环境中平面45°倾斜负浮力射流进行数值模拟，来对比两种环境中的流动特性。 1 数学模型 如图 1 所示，一股密度为ρ 的射流从宽度为 B 的狭缝中， 0 以初始速度u 沿倾角45°倾斜地流入密度为ρ的均匀环境中， 0 a ρ＞ρ，上游有恒定来流，流速为 u 。取 x 轴沿水平方向，y 0 a a 轴铅直向上。 设流动是定常的，由布辛涅斯克型基本方程并用k- ε湍流模 型进行封闭，可得：连续性方程 图1 坐标示意 ∂u ∂v + 0 (1) ∂x ∂y 收稿日期：2002-04-01 基金项目：本文得到国家自然科学基金的资助 作者简介：杨中华(1977-)，男，山东成武人，博士研究生，从事环境水力学、计算流体力学方面的研究。 58 水 利 学 报 2003 年07 月 SHUILI XUEBAO 第7 期 x 向动量方程 ∂u ∂v p u   u v  1 ∂ ∂  ∂  ∂ ∂ ∂ u