圆射流粘性稳定性的谱元法计算.pdfVIP

中国工程热物理学会 学科：传热传质学 学术会议论文 编号：083221 ** ** 圆射流粘性稳定性的圆射流粘性稳定性的谱元法谱元法计算计算 圆射流粘性稳定性的圆射流粘性稳定性的谱元法谱元法计算计算 谢明亮, 周怀春, 张引弟 华中科技大学煤燃烧国家重点实验室, 武汉, 430074 Tel：0278505， Email：mlxie@mail.hust.edu.cn 摘要摘要：本文研究了圆射流线性稳定性问题，数值方法采用谱元法，通过与前人计算结果的比较，验证 摘要摘要 了本文所采用的计算方法切实可行，同时也给出了一阶螺旋状扰动的拟序结构，为进一步研究圆射流 场的流动结构控制方法奠定了基础。 关键词关键词：射流；流动稳定性；谱方法；有限元方法 关键词关键词 0 前言 射流是流体运动的一种重要形式，在航空航天、动力、化工、自动控制、农业、环 境工程等方面有广泛的应用。因此对射流研究并进行控制具有理论价值和实际意义。圆 射流是指流体从圆形出口射出流入周围另一流体域内的流体运动。流体射出时，只要射 出速度与周围流体速度不同，就会与周围流体之间形成速度不连续的间断面。间断面一 般不可避免地受到干扰，于是间断面要失去稳定而产生旋涡。早在上世纪 60 年代，人 们就开始了研究圆射流场拟序结构的产生与演化机理[1-3] ，现有研究发现，圆射流被大 尺度拟序结构所控制，这些结构起初是剪切层中的不稳定波，后来波的振幅增大，形成 拟序的涡结构，当涡结构往下游发展时，临近涡结构发生合并，这些过程决定了射流场 的发展。由于对射流的控制在工程中有非常重要的作用，如提高燃烧效率、加速化学反 应过程、控制火眼扩散、改进喷射效果等，目前对射流稳定性问题的研究正方兴未艾[4-6] 。 由于圆射流的轴对称性，描述圆射流的控制方程通常采用圆柱坐标系。这种处理往 往带来数学上的奇异性，在数值计算中，一般需要采用附加的物理约束，如对称性边界 条件，正则化约束等[7-9] 。由于有限差分方法在处理边界条件上的灵活方便性，因而在 圆柱坐标系 （或极坐标系）中得到了广泛的应用[4-6] 。此外，射流还是一种无界问题， 这都给谱方法在射流流动稳定性问题的准确计算增加困难[10] 。 本文充分利用谱方法[11]具有速度快、占用机器内存小，精度高等优点，结合Petrov Galerkin 有限元方法在处理边界条件上的优势，计算圆射流稳定性方程，并通过与高精 度紧致差分方法的结果进行比较，证明本文方法的切实可行。 1 控制方程 圆射流线性稳定性方程[12] ： *本文得到国家自然科学基金 和“111”引智计划 （B06019 ）资助 1  2 1 n2 +1 2 2  −ikcu =−Dp + D u + u − u −k u − inu −ikU u (1)