横流环境湍射流涡动力学特性数值模拟 numerical simulation of vorticity dynamics for turbulent jet in crossflow.pdfVIP

第21卷第3期 水科学进展 V01．21．No．3 INWATER 2010年5月 ADVANCES SCIENCE May，2010 横流环境湍射流涡动力学特性数值模拟 姜国强1，任秀文1，李炜2 (1．环境保护部华南环境科学研究所，广东广州510655；2．武汉大学水利水电学院，湖北武汉430072) 摘要：采用三维湍流模型及其混合有限分析解法进行横流中湍射流这一复杂三维流动的研究。在利用实验数据对 模型及其计算方法验证的基础上，对多种喷口形式和流速比工况下的流速场和涡量场进行了数值计算，模拟得到 了其旋涡结构发展演化特性。在射流初始阶段，横流在射流背流面形成绕流分离旋涡，其结构与射流喷口形式和流 速比有关。射流主轴沿流向布置的窄缝射流产生的旋涡最为特殊和复杂，在喷口侧面存在4个分离点，从而形成4 个旋涡。在射流远区，流动主要由反向旋转涡对所控制，并诱导出二次涡对。分析得出了射流喷口形式及流速比对 反向旋转涡对涡核位置和旋转强度的影响。 关键词：横流；湍射流；旋涡结构；水动力特性；数值模拟 中图分类号：TVl43 文献标志码：A 文章编号：100l-6791(2010)03-0307-08 横流流动环境下的湍射流广泛地存在于大量工程技术和科学领域，如环境科学领域中废水、废气的排 放；空气动力学领域中飞行器垂直起降气动性能的研究；热能动力领域中燃油发动机燃烧室内冷却气幕的引 入等。在上述3个工程技术和科学领域的流动都可概化为横流中的圆孔或有限长窄缝射流的旋涡运动，研究 复杂涡系对时均流动特性的影响，可揭示横流与射流两股流体之间复杂的相互作用，得出排放方式和环境条 件对湍射流发展影响的规律；对各种复杂涡系产生机理及演化方式的研究，将加深对湍射流物理机理的认 识。在试验研究方面，Yasuhiro…对不同流速比情况下的射流宏观特性进行了测量，发现基于射流流速场和 技术，试验研究了水槽底部排放的单孔和多孔射流，观测到射流出现分叉现象和马蹄涡结构。在数值研究方 面，SykesH。等运用一方程湍流模式，探讨了在高流速比下横流中射流的涡动力学特性，槐文信、李炜等”1 通过计算分析了横流中单圆孔紊动射流特性，发现射流出现了分叉现象；这些研究都表明旋涡结构，特别是 大尺度的旋涡结构控制着时均流动特性。本文在文献[6-9]的基础上，采用双方程湍流模型对三维流动基本 方程进行封闭，对不同喷口形式和流速比情况下的流速场和涡量场进行了较为详尽的计算，以期对不同喷口 形式和流速比情况下的涡动力学特性有一些规律性的认识。 1 数学模型 1．1 流动的基本控制方程 以圆孔射流孔径D或窄缝射流窄缝宽度曰为特征长度L，横流的来流平均速度砜为特征速度，可得出 量纲一化的控制方程，方程的具体数学表达式及模型常数取值见文献[9]。 1．2 定解条件 进口断面：“=Uo，秽=W=0，ko=0．001联，占o=0．001Vo／￡ 出口断面：塑ax=竺=面OW=-Ok=芸=0；对称平面：口=00xd名 Ox d菇 ，考=等=芸y=等=od， dy d d， 收稿日期：2009-06-02 基金项目：国家自然科学基金资助项目 万方数据 308 水科学进展 第2l卷 _l=3Zil孬：埘=0，idu=iov=警：警 az az dz dz 侧平面和下平面：采用壁函数法 exp(一1．7y2)／L 射流喷口(一L／2≤Y，L／2)：u=秽=0，埘=埘f，kj=o．06埘；exp(一1．7y2)，占，=0．06埘3i 对于不可压缩流动，有意义的是压力差而不是压力的绝对值，计算时选择横流入口断面某个节点，令其 压强为0，作为参照值，将其它节点处的压强作为相对于压强参照值的相对压强。 前述边界条件表达式中：u，移，W分别为石，Y，；3个方向上的速度分量；k为湍动