时域谱方法用于内流非定常流动计算研究.pdfVIP

中国Ii程热物理学会 热机气动热力学 学术会议论文 -t 编号：062103 时域谱方法用于内流非定常流动计算 苏欣荣，袁新 清华人学热科学‘j动力T．程教育部重点实验室，北京，100084 摘要：奉文使用时域潞方法对内流气动声学基准算例和．二维喷管算例进行丫计算，初步计算表明相 对千常规的1r定常时间推进方法。对于刷期性1F定常流动，时域谱方法口l以缸州『问方向}：达到很扇的 精度l川时减少计算耗时。从而町以应用fl-：nl’轮机械椎定常流动计算中并减小计算规模． 关键词：时域潜方}左，周期性非定常流动，商精度格式 1．引言 叶轮机械内部的流动在本质上是1卜定常的，受客观条件限制，对1F定常流动现象的 较为细致的研究，包括实验研究和数值计算在最近才蓬勃发展起来。由丁．相邻叶排之间 的相对周向运动，叶轮机械内部的大部分北定常流动都可以看成是周期性的。现有的非 定常计算方法通常在时间方向上使用二阶Euler后差格式【1】等时间推进方法。当使用时 间推进方法计算周期性1卜定常流动时，由丁．给定的计算初值一般井不严格满足1仁定常控 制方程，冈此计算过程要经过一段过渡态的计算才能使流场中各点的流动参数呈现周期 性变化，计算规模越人，需要计算的过渡过程时间越长，从而计算耗时越多。 在计算叶轮机械内部的非定常流动时，对丁．周期性1仁定常流动应当充分利用流动的 周删性。早期的基丁线性化的处理方法[2】把1卜定常成分看作是周期性的小扰动，从而在 线性化控制方程的基础上求解1F定常流动。He在此基础上考虑了1E线性效应【3】，发展 harmonic 出了Nonlinear 均中对雷诺应力项的处理，并把该方法Jl=Ij丁考察多级压气机的定常和非定常特性。 Jameson等【4，5】把谱方法应用在时间方向上，而对空间项则采用有限体积方法处理。 在该方法中，只要求所求解的流动是周期性的，并不需要做小扰动假设或者线性假设。 为了考察时域谱方法在内流流动计算中的可行性，本文使用时域谱方法对一个气动声学 基准算例和二维喷管算例进行了计算，初步的计算结果表明对于周期性非定常流动时域 谱方法可以达到较高的精度同时人人减少计算时间。 2．数值方法 计算对象分别为一维变截面喷管和二维喷管，目Ii{『的进行的计算都是无粘计算。一 维喷管取白rICASE／LaRC计算气动卢学专题[6】，二维喷管为Ⅱ音速对称喷管。为了简 便，在卜面的讨论过程中用(2．1)式米分别代替二维Euler方程和准一维Euler方程。 a，‘ 豢+厂(“)=0 ·(2．1) ar 2．1． 时间格式 住时间方向上分别使心了两种格式来逼近叫’问偏导数项：时域谱方法和二阶Euler 肝蓐方法。谱方注的基本原理不再叙述，参考文献【7，8】。把谱方法席川剑时问方向上， 774 鉴于流动在时间方向上的周期性，取基函数和权函数均为三角函数。臼变量为 U=“(f，．．．)且在时间方向上具有周期性，则zl可以用有限项级数和的形式近似表示为： “Ⅳ(f)=∑反e拗 (2．2) 七 借助丁．基函数和权函数的正交性，可以得到控制方程在频域中的表达式： ikco五l+以=0，尼=一N／2，…N／2 (2．3) 往上式中，缈为基频频率。在中，玩和Z是通过F0urier变换计算得到的：在实际中， f的形式可能很复杂，此时使用离散Fourier变换同样可以得到相同的形式。对于角频 借助丁级数求和，“(o)可以表示为： 、． “，=