管道凸起激波湍流边界层干扰的数值模拟.pdfVIP

第25卷第i期 应 用 力 学 学 报 VoL25No．1 2008年3月 CHINESEJ(圯喂NALOFAPPLIEDM[Ea日渔NICS Mar．2008 文章编号：1000-4939(2008)01—0001-05 管道凸起激波／湍流边界层干扰的数值模拟。 任鑫 高歌 (北京航空航天大学100083北京) 摘要：采用基于部分平均的湍流方程组，模拟了管道凸起流动中两种不同强度的激波一湍流边界层 干扰现象，对流场的时均参数与实验值进行了比较，计算得到的壁面压力分布、摩阻系数分布、边界 层厚度和速度型与实验值比较吻合很好。结果表明基于部分平均的湍流方程组既能够准确模拟小 分离的激波／湍流边界层干扰流动又能够准确模拟大分离的激波／湍流边界层干扰流动，较好地预 测了典型的A激波结构。 关键词：可压缩流体；湍流；激波／湍流边界层干扰；基于部分平均的湍流方程组 中图分类号：0357．5 文献标识码：A 动算例进行模拟，进一步研究部分平均的湍流方程 1 引 言 组的性质以及探讨复杂湍流的模化的思路。 可压缩性带来的激波／湍流边界层干扰以及随 2控制方程 之而来的流动分离，对可压缩湍流的数值模拟提出 了严峻的考验。通常的两方程湍流模式难以准确模 造成湍流不封闭问题的罪魁祸首就是湍流中大 拟激波／湍流边界层干扰所引起的分离流、摩阻系 尺度旋涡的运动，大尺度旋涡运动导致了湍流运动 数、传热值[11。为了克服PANS方法碰到的困难，是一个强烈的非平衡过程；另外，流体的湍流运动仍 旧是连续介质力学的问题。既要体现大尺度旋涡运 Knight尝试采用大涡模拟技术来解决这一难题[2]， 尽管大涡模拟技术在二维高速流中可以应用，但是 动又要满足连续介质力学的要求，势必要求引入全 要想可靠地应用到三维高速流动情况下却十分困 新的平均方法对N-S方程平均后余下的脉动量再 难。近年来发展的非线性湍流模式，比如显示代数 进一步平均得到描述非平衡运动的控制方程。为 雷诺应力模式、非线性涡粘性模式，虽然在研究性的 此，文献[4]引入了部分平均方法将湍流运动分解为 计算上有了较好的结果，但是这些模式倚赖的经验 准平衡的平均运动和脉动的一阶统计平均运动，建 系数过多，限制了在工程应用的推广。由此可见，激 立平均运动和脉动一阶统计平均信息的描述方程， 波／湍流边界层干扰的数值模拟问题并没有较好的 耦合求解来描述湍流运动。如果采用Favre质量加 解决。 权平均方法，可以将部分平均湍流模式理论直接推 本文采用基于部分平均的湍流方程组[H]对标 广到可压缩流体，方程组中漂移流与不可压缩流体 准的管道凸起激波／湍流边界层干扰流动算例，即 的湍流模式基本相同。平均流增加了能量方程和状 Liu[5]的圆弧管道凸起激波／湍流边界层干扰流动算态方程。下面直接给出部分平均湍流方程组的平均 例和Deleryc6]的管道凸起激波／湍流边界层干扰流 流连续方程、动量方程、能力方程、漂移流连续方程、 -来稿日期：2006-06—15修回日期：2007-03-01 第一作者简介：任鑫，男，1979年生，北京航空航天大学能源与动力工程学院，博士生．E-mail：renxin960(园sohtl．eom 应用力学学报 第25卷 转换矩阵来求出