浅水环境中含污染物横向射流三维浓度分布特性及研究.pdfVIP

浅水环境中含污染物横向射流三维 浓度分布特性研究。 樊靖郁 张燕 王道增 (上海大学，上海市应用数学和力学研究所，上海200072) 摘要本文采用RNG湍流模型对浅水环境中含污染物横向射流近区的三维流场和浓度 场进行了详细的数值计算，并与相应的实验结果进行了比较。数值计算结果表明在近区冲击 区内射流冲击底壁后形成的壁射流与环境横流的相互作用具有强烈的三维性，在滞止点上游 存在一高浓度聚集区，其纵向和横向范围受到环境水深和流速比的影响。滞止点下游壁射流 区的浓度衰减随环境水深的增大而增大，随流速比的增大而减小。当流速比相对较小时，过渡 区内在靠近对称面以及底层壁射流与环境横流的横向边界附近出现横向高浓度聚集区，其横 向位置和范围随下游距离的增加逐渐沿横向扩展。 关键词浅水，横向射流，冲击射流，浓度分布 1 引言 横向射流(transversejet)在实际工程问题中有许多重要应用，如环境工程中的河流或近 海排污，化学工程中的混合容器和燃烧工程中的燃料喷射等。在某些实际工程问题中，射流的 起始密度大于环境流体的密度，如环境工程中化肥厂产生的含石膏废水、海水淡化厂的高盐水 的排放以及航道疏浚工程、水底结构工程和大规模填海工程中的高浓度泥浆的倾泻和排放 等[】]。由于射流受到与重力方向一致的负浮力作用，为充分利用环境水体的稀释能力有时需 要采用表面排放的方式。对于垂直向下射人浅水环境中的含污染物动量射流，在射流初始动 量和负浮力的作用下，射流的发展和演化会受到底壁的阻碍形成冲击射流。底壁限制下的含 污染物横向射流涉及到射流、横流和底壁的相互作用，射流边界的剪切层和底壁边界层的存在 使得近区范围内的流动性质和污染物浓度分布具有复杂的三维性，特别对于冲击区和过渡区 内的浓度分布特征在以往的研究中还未见详细报道口]。本文采用RNG湍流模型对浅水环境 中含污染物横向射流近区的三维流场和浓度场进行了详细的数值计算，并分析了射流对称面 以及近区范围内不同区域横断面的污染物浓度分布特性。 · 上海市重点学科资助项目 控制方程和数值计算方法 控制方程采用三维形式的雷诺平均N—S(RANS)方程和湍流扩散方程： 等+蠡(瓦，=。 P-。万3ui+玩等)一一差+甭；+者(p罄一乃丽) (2) (3) d￡ adC￡+u，亳一矗(矾蠢一面)+厨’瓦一瓦‘um石一。“，)十M 其中“一”表示雷诺平均值，“；(i=1，2，3)分别为流向、横向和垂向的流速分量，lD为密度，P 为压力，c为浓度，g为重力加速度，D，为分子扩散系数，M为源项。 本文数值计算采用RNGk-e湍流模型，控制方程采用有限体积法进行离散，对流项采用 相同口]。射流出口直径D=0．005m，U和砜分别为射流出口和横流初始流速，相应的流速 比R(定义为R=U，／u。)分别为7．79和15．6。射流出口初始密度B为1170kg／m3，初始密 度Froude数Fr=】0．21。 计算区域采用与实验完全一致的几何条件，长L为lm，宽2B为0．5m，水深H为 个不同区域分别布置非均匀网格，其中射流出口共布置128个网格，邻近射流出口及近壁区域 进行局部加密。综合考虑计算精度和效率，垂向网格数依赖于水深而略有差异。当相对水深 x40；当相对水深较高(H／D=9)时，总的网格 较低(H／D一5)时，总的网格数为180×60 数为180×60×50。 计算区域内边界条件包括射流进口，横流进口及下游出口，对称面，自由表面和固壁。其 中射流和横流进口边界条件分别采用预先计算充分长的管流和明渠流的平均流速和紊动参 数，简单地给定均匀边界条件会给冲击射流近区的计算结果带来较大误差。本文对固壁边界 采用壁面函数法来进行处理，同时对自由表面进行了简化，没有考虑射流入水时对自由表面变 形的影响，采用自由滑移边界条件。下游出口边界采用～般的均匀出流零梯度边界条件。 3结果和分析 3．1对称面浓度分布特征 图1