旋转盘腔中紊流的数值模拟.pdfVIP

中国航空学会推进系统气动热力学1}业委员会 cS，LA9扣019 江苏(苏鲁皖)航空航天学会动力裂鲞：专业委员会 中国人民解放军空军工程学院 ’99学术会议 旋转盘腔中紊流的数值模拟 吉洪湖张靖周梁波 (南京航空航天大学动力工程系，南京210016) 摘要本文对一个旋转盘腔的流场进行了数值计算研究．谊盘腔由一个旋转盘，一 个静止盘厦静止的外围屏组戍．在盘腔的核m区采用标准k-e模型，在近壁区分别采 s模型，并且在第二个模型中考虑了固体旋转引起的耗散率降低的修正．计算结果表 明第二个模型比单方程模型能更准确地再现实验数据． 关键词紊流模型，旋转盘胜 l引言 先进的发动机中有各种各样的涡轮盘腔结构，其中一种常见的结构是由转盘和静盘 形成的转静盘腔。一般有高压冷却气流穿过这些盘腔，对其进行冷却并起防止高温燃气进 入涡轮盘腔的封严作用。深入地丁解涡轮盘腔中流场的基本特性，对改进冷却结构设计， 减少对冷却气量的需求，从而提高发动机的性能是十分有意义的。 本文对旋转盘腔流场进行了数值分析研究。所计算的旋转盘腔的几何结构与文【l】实 验研究中的结构相同，如图I所示。图中横坐标与纵坐标分别为用转盘半径b无因次化的 为静止的环形围屏。为了清晰起见，图中的网格没有完全画出来，而是每五条线中画出一 条线。在计算中采用有限容积法对雷诺平均的流动控制方程进行离散，用SIMPLE算法对 其求解。计算结果与实验数据进行了比较。为了节省篇幅，本文仅对所采用的紊流模型进 行简要的描述．其他细节不再给予讨论。 ro‘罐骧枷s 景o¨睦要塞篓匡霉雾蓁弱啪蛔‘ 。蘑茎綦蓦嘻喜争摹薹蓦碧眦or 图I计算域 2紊流模型 本文采用常用的紊流模拟方法，把流场区域分为完全紊流区和低雷诺数近壁区两部 分。在完全紊流区采用高雷诺数模型．在近壁区，由于壁面的存在极大地影响了该区域的 紊流运动。使其有效紊流粘度变小，因此分子粘度的作用变得显著起来。在紊流模拟时须 采用低雷诺数模型。这两类模型各有许多种，在本文中采用的高雷诺数是标准的k．￡模型， 吉洪湖张靖周粱波 旋转盘腔中紊流的数值模拟 模型。 2．I高雷诺数k．e紊流模型 在高雷诺数区域中采用k-￡模型。根据有效粘度近似紊流应力可由下式模拟 丽；b“盟尸iIo敏u，,+百OUj】 (1) 的以下微分方程求解 掣=毒陋期卦矿一胆 ∞ 掣专陋期讣伽妻P-c62P譬 ④ 式中，Fm为普朗特一施密特数，下标中代表k或s；紊流动能生成率|P定义为： P=-—u,u—jOU，／gx．。它反映了通过丈尺度的紊流涡和平均应变场之间相互作用引起的从平 均过程到脉动过程的能量传递。方程中的常数 表1方程中的常数 c“t％，a。，c∥c。2由表1给出。 d 2．2低雷诺数紊流模型 卜 k 盯P CPI ％l 2．2．1单方程模型 0．09 1 1 22 1．44 1．92} 单方程模型的紊流应力和紊流动能的方 程与标准I一￡模型中的类似，分别是方程(1) 和(2)。耗散率e由给定的近壁区长度尺寸，通过代数方程}=^7／，f求出。式中，长度尺 离。紊流粘度根据胁=pcJ．uki计算，此处‘；2．55Y0一exp(-0．016y‘m· 2．2．2Launder-Shar