轴对称紊流射流流场数值模拟.pdf

维普资讯 1995年 第 j9卷 石油大学学报 (自然科学版) Vo1．』 Ⅳ0．2 第 2期 JD口r¨ of,the r 坷 。 Pe~roteam．Chiua Ayr j995 轴对称紊流射流流场数值模拟 、 鼬石(油大学开发系，，山山东末营257062) 确 D 歹／ V 。摘要 选取标准的 ．亡二方程模型及优选的紊流模型系数，运用有限差分方法对淹没轴对称紊流 射流流场进行数值计算．计算结果表明，射流轴心速度是逐渐衰减的．其衰减规律和射流出口雷诺敷 无关 ·无园次等逮棱长度 (L／D)约为 {～5．和实验数据符台较好．对射瀛径向速度 P珊径 向截面 的 分 陋 射流基本殷内轴向速度分布的相似性进行了预测．L ； ， 主题词至多 堕璺墼堡墼 钧i靠 J {_牛 ) 大学博士研究生及石油开发系副教授、执事高压水射流数值模拟方面的研究． 0 引 言 前人在射流流动规律的研究领域 已做了大量的工作 ，着眼点大都基于理论分析和实验 研究．理论分析可给出一般的流动规律 ，但受边界条件的限制 ；实验可给出第一手 的数据 ． 但同样受测试手段和精度的限翩 ，并且所测数据有限．由于电子计算机应用的高速发展，数 值求解非线性方程组已成为可能．数值模拟的优点在于它可以正确给出复杂边界条件下的射 流流动规律，并能够预测出由于实验条件的限制所测试不到的流场信息．且花费低 本文选 取了 二方程模型以封 闭Jv—s方程组 ，运用Patankar方法对轴对称淹没射流流动进 行数 值求解 ．得出了一些有益结论，并由有关实验数据进行了验证． 1 控制方程及边界条件 1、1 控制方程 选取轴对称二维圆柱坐标下的微分方程式 ，为便于计算机程序的实现 ，控制方程可写为 如下的统一形式 + (pvr )一 (厂． )一 ( )一s。 (1) z d r ~／T d， 式中， 和 为轴 向和径 向速度分量； 除代表 和 外，还表示紊动能 ，紊流耗散 ；风为 相应方程的源项，式 (1)中厂．和 。的具体表达式可参阅文献 [53．标准 模型中的有关紊 流常敬取为 c．一O．∞，C=1．d4，0一192，=1、0，．=1．3． 标准的 胁 二方程紊流模型可以正确地预示多种紊流流动规律 ，但在描述轴对称射流时 却出现较大偏差．为此罗迪 (w．Rodi)给出了有关修正理论，并推荐 沛l艘 下式计算 ； 收稿 日期-199a一11一 ． 48 ． 维普资讯 l ．= 0．09-- 0．045 J =】．92--o．0667f … 去(等一I等1)r 此式反映出射流轴线速度 ”．减少时， ．和 ：的数值也相应减小． 在大量数值实验的基础上优选出 I为 1．505．该修正常数的合理性 已由实验数据予 以验 证 ． 1．2 计算网格的划分及边界条件 用一射流流入轴对称空腔来模拟轴对称淹没 射流 问题 ，当 d／D取得很小 时 (D为空腔 直 径)，计算将是合理的．计算区域划分为 30X30 个网格，轴 向长度采取渐扩式非均匀网格 ，长度 比 ／J 一1．05．计算 网格划分 示意 图见 图 I． 为求解椭圆型方程组 (1)，给定有关 的边界 条件 ： (1)进 口处 Ⅱ一 ， =0一k一5“(为紊 流强度)，