压力梯度作用下近壁湍流相干结构相互作用.pdfVIP

574 2003空气动力学前沿研究论文集 压力梯度作用下 近壁湍流相干结构的相互作用 张强陆利蓬 (北京航空航天大学，北京100083) 摘要近壁区的相干结构是湍流边界层中最活跃的大尺度运动，它是湍流发生和维 持的重要基础。工程实际中压力梯度特别是逆压梯度造成的边界层分离常给流体机 械带来负面影响。而相干结构在时间上和空间上都是多频率出现的，并且它们之间 有很强的相互作用。因此研究压力梯度作用下两个相干结构的相互作用非常必要。 本文采用直接数值模拟的方法研究压力梯度作用下两个相干结构之间相互作用的演 化问题。结果表明逆压梯度增强了相干结构之间的作用。 关键词 相干结构 湍流边界层压力梯度直接数值模拟 中圈分类号文献标识码 1引言 相干结构理论模型的建立，不仅对于了解湍流发生的机理有重要意义，而且对湍流的计算 也有帮助。在工程实际中，特别是航空发动机内部的叶栅通道的流动中，边界层中的湍流经常 处于非平衡状态，流动变得更加复杂。特别是逆压梯度造成的边界层分离会造成叶轮机效率 的降低，严重的造成失速，导致喘振等严重后果，所以研究逆压梯度下的流动具有工程实际价 值。实际流动中相干结构时间上和空间上都是多频率出现的，而且它们之间有很强的相互作 用。因此有必要研究逆压梯度作用下相近的相干结构的相互影响问题。 2理论分析及数值方法 在近壁区由于黏性作用很强，湍流的小尺度脉动很弱，可以将瞬时速度场和压力场分解为 基本流项和相干结构项： H=U+Ⅱ。，P=P+P。 (1) 其中，U，P表示基本流项；H…P表示相干结构项。将它们代入N—S方程可得出相干结构项 满足的扰动方程： (2a) 酱+(【，·V)‰+(¨v)v=去V2峨一vpc-(¨V)虬 V·H，=0 (213) 基本流【，(_(z，Y)，i(z，Y)，0)由解N—S方程直至稳定而得，其初始场为经法向压缩的层流 剖面一部分，边界条件为在，+=100处与湍流平均速度光滑连接。其合理性见文献[1，2，3J。 压力梯度作用下近壁湍流相干结构的相互作用 575 y，z，t)的Fourier谱展开式可表示如下形式： 爿1 妒(。，y，：，f)：∑‰(。，y，￡)e—zm＆m=一了N，…，可N一1 (3) 一 m=譬 其中，N是截断误差的阶数；卢为展向波数；‰(z，y，t)为P(z，y，z，f)的谱展开系数。对非 线性对流项进行Fourier谱展开为 掣- 2难 (4) (虬·v)坼=∑Fm[(u。·V)‰]e m=譬 个二维方程组： (5a) (5b) V。·U。=0 N N． m2一百，‘一，—i l V。={瓦3，￡，一i邢}，V。2=毒+尊一m2萨 对流项的空间离散采用傅德薰、马延文(1996)的5阶迎风紧致格式计算。压力梯度项中通过 ／+v。．H二=o得到V乙“=里乎令6=m2