rans方程和附面层方程耦合求解转捩位置的方法 transition prediction for 2d flows based on coupling rans equation and boundary layer equation.pdfVIP

第39卷第1期 航空计算技术 V01．39No．1 2009年1月 Aeronautical Jan．2009 ComputingTechnique RANS方程和附面层方程耦合求解转捩位置的方法 邓磊1，乔志德1，熊俊涛2 (1．西北工业大学翼型、叶栅空气动力学国防科技重点实验室，陕西西安710072； 2．加利福尼亚大学欧文分校机械与航空与太空工程系，CA92697—3975) Navier—Stokes)方程和附面层方 摘要：研究了基于线性稳定性分析的RANS(Reynolds—Averaged 程耦合求解转捩位置的方法。RANS方程计算得附面层外边界的速度分布作为附面层方程计算的 输入，求解附面层方程得到稳定性分析需要的附面层参数。由于附面层的解在分离点处的奇异性。 用交互式附面层方程求解方法。用基于线性稳定性分析的e“方法计算转捩位置，并考虑了Toll． 位置收敛。RANS求解中使用转捩过渡区模型，避免了点转捩模型引入的数值扰动。通过对NA． CA0012翼型的计算并和实验结果及Xfoil计算的结果进行比较，吻合较好，有较好的工程实用价值。 关键词：Navier—Stokes方程；附面层方程；e“方法；转捩预测 中图分类号：V211．3 文献标识码：A 文章编号：1671．654X(2009)01—0027．04 引言 的解，在分离点垂直方向的速度分量无穷大，这称为附 面层方程在分离点的奇异性问题。这时需要给定附面 层流一湍流流动中准确的转捩位置的预测有非常 层的动量厚度分布，将附面层外边界的速度或者压力 重要的实际意义。在RANS方程求解中，一般使用全 作为解的一部分进行计算，这种方法称为逆问题求解 湍流计算或者使用固定转捩计算。转捩位置的变化。 方法。这样既解决了附面层方程在分离点的奇异性， 虽然对翼型表面的压力分布影响不大，但是对翼型表 又满足了分离点下游“回流”的存在。另外，已经证明 面的摩擦系数有较大的影响。因此转捩位置的准确预 对有分离气泡的流动，用逆问题计算的方法求得解和 测对亚音速翼型阻力计算有重要意义。 N—S方程的解吻合很好。 根据翼型的外形、攻角及雷诺数()，附面层或 RANS方程求解中，一般使用“点转捩”模型，即在 者在最小压力点之后很快进入失稳转捩或者先分离而 转捩点设置湍流模态开关。这样由于转捩点附近的位 后进入转捩。从层流到湍流的转捩问题，一般有几个 移厚度明显减小，引入了一个当地的扰动，这样导致一 方法：比如基于抛物化稳定性方程(PSE)的方法，或者 个对上游流动有明显影响的有粘／无粘迭代区域。对 基于非定常Naver—Stokes方程的方法(直接数值模拟