对短周期风洞中二维涡轮叶栅的非定常数值模拟.pdfVIP

中国工程热物理学会 热机气动热力学学术会议 对短周期风洞中二维涡轮芝￡婴的非定常数值模拟 邢秀清赵晓路 中科院工程热物理所，北京，100080 ● 摘要：为了对短周期风洞实验中，实验段建立起周期性稳定流动的时间进行估算，以 ● ■ 便在实验中把墨测量时同，保证测量数据的有效性和准确性，本文采用显式时间推进方 法求解二维Eu蛔方程对短周期风洞中二维涡轮叶栅的非定常流动进行了模拟．除预测 ● 建立起周期性稳定流动过程中的流场特征和达到周期性稳定流动所需的时间外，对叶片 的升力等冲击载荷在建立起周期性稳定流动过程中的变化也进行了估算，为三维涡轮部 件舜态实验方案的规划和实验设备的研制提供帮助。 关链词：短周期风涸，非定常流，敷值模拟 I引言 未来涡轮入口温度将达到甚至超过2000℃的高温，通过实验研究涡轮叶片在高温运 行下的热环境、冷却方法及效果，是提高涡轮设计质量的重要内容。短周期实验设备以 其易达到上述涡轮实际运行工况，同时具有能方便、独立地改变来流雷诺数、马赫数， 实验运行撮作简单，费用低廉等一系列优点，越来越受到航空发动机设计研究者们的青 璩．国外从七十年代起纷纷发展用于涡轮传热实验研究的短周期风洞和实验技术，其中 以荚田牛津大学和美国。邺PAN公司最为突出。他们发展的短周期实验设备，不仅有 二雏叶鼍实验风洞．也有用于全尺寸发动机透平部件的瞬春设备，为设计部门提供了大 量而有用的数据，促进了英、荑航空发动机的进步和发展【I肛J。 迄今为止，我国的短周期实验设备，只有二维叶栅实验风洞，需要建立一个用于涡 ■ 轮部件的瞬态实验台，为我国航空发动机的迅速发展提供帮助。由于短周期实验的实验 时间非常短，一般只有十几或几十毫秒时间，给测量带来很大困难，同时如此短的时间 内．实验段能否建立起周期性稳定的流动，叶片受到的冲击载荷有多大，以及由此带来 的轴承(特别是止推轴承)的受力情况如何，实验件强度、叶片强度和实验中测量时间 是否足够，是关键问题。 本文试图紧密结合实验研究和数值模拟，在实验方案规划、实验设备的研制过程中， 使用数值模拟进行定性乃至定量的估计和分析，以求尽量减少盲目性，达到缩短周期、 减少资金投入的目的。为此本文从对短周期风洞中涡轮叶播非定常流动的模拟开始，对 建立起周期性稳定流动过程中的流场特征和达到周期性稳定流动所需的时闻进行估算． 同时，借助非定常流数值模拟对叶片受到的冲击载荷进行模拟。最终实现对涡轮部件实 验中周期性稳定流动建立过程的数值模拟。为实验提供帮助。 2数值方法和边界条件 本文采用有限体积法，显式时间推进，四阶孙m辨-KⅡ曲法。为抑制数值振荡引入 了Jam∞∞=阶和四阶人工粘性项。计算过程每一点时间步长相等且恒等于全场的最小 时闻步长。对来流压力突升的两捧静止叶橱(导叶和动叶间无相对运动)的数值模拟求 解}“S方程，对来流压力突升的转、静干扰流动的数值模拟求解Ellh方程。 边界条件为：在所有固壁边界上，法向速度置为零．满足纯熟条件。在进出口边界． 采用胁s发展的二维非定常无反射边界条件。进口给定总温、总压、气流速度和气流 ● 角，出口给定背压13】。为简化计算并节省计算时间，本文假设转、静叶片数相等。 3算例及结果分析 ● ● 本文选用一工业燃气涡轮的末级叶中叶栅进行数值模拟【4I。计算中静叶叶橱中的网 格敷为周向25个节点，轴向49个节点，动叶叶橱中周向25个节点，轴向43个节点。 ● 采用H型网格。叶栅几何及工作参数见表一，时栅形状及计算网格见图l。