可压缩平板湍流边界层拟序结构直接数值模拟.docVIP

可压缩平板湍流边界层拟序结构直接数值模拟 　　摘要： 文中直接数值模拟了来流马赫数（M∞）为2.25的空间发展的可压缩平板湍流边界层。研究了可压缩平板湍流边界层内的拟序结构的形成过程。指出在粘性底层出现低速条带是触发拟序结构的第一个信息。条带升起，开始发生振动，然后破裂，并伴随着产生强烈的湍流脉动，这便形成了湍流的猝发过程。因此壁湍流拟序结构的重要意义在于它是生成湍流的重要机制。条带的破裂还伴随着上抛和下扫现象，这是壁湍能生成和耗散的关键。 　　Abstract： A Direct numerical simulation of a spatially evolving compressible flat-plate turbulent boundary layer with incoming Mach number （M∞） of 2.25 was performed in this paper. The coherent structure forming process is studied. The low-speed strip appearing in the viscous sublayer is the first information to trig coherent structure. Strip rises， begin to vibrate， and rupture， and accompanied by strong turbulence?? which will form the turbulent bursting process. Therefore the wall turbulent coherent structure is significant in that it is the important mechanism to generate turbulence. Strip rupture accompanied by up and down sweep phenomenon， it is the key of wall turbulent energy generation and dissipation. 　　关键词： 可压缩；平板湍流边界层；直接数值模拟；拟序结构 　　Key words： compressible；flat-plate boundary layer；direct numerical simulation；coherent structure 　　中图分类号：O35 文献标识码：A 文章编号：1006-4311（2013）10-0003-02 　　0 引言 　　对湍流拟序结构的研究，使得人们在揭示湍流的物理本质方面寻觅到了一个新的有效的方向，从而得以有能力预测和控制湍流。在航空航天、海洋、化工等工程领域存在着大量的壁湍流问题，其中平板湍流边界层是一种典型的壁湍流。目前，该方面的主要研究成果可参见文献[1-9]。可以了解到，其中研究较多的是不可压缩流动问题。比如文献[3]中利用大涡模拟方法研究了三维槽道完全发展湍流近壁区流动的拟序结构。相对而言，人们对可压缩壁湍流近壁区流动的拟序结构的研究成果还比较少。文献[9]对来流马赫数0.7，来流Reynolds数250000/Inch，锥角为20°的尖锥边界层的整个空间转捩过程进行了直接数值模拟，显示了近壁湍流的高、低速条带结构。 　　本文采用高精度差分方法求解三维可压缩Navier-Stokes方程，直接数值模拟了空间发展的可压缩平板湍流边界层，主要分析了拟序结构的形成过程，利用脉动速度分布显示了近壁条带结构的特征。 　　1 控制方程和数值方法 　　直接数值模拟可压缩平板湍流边界层流动的数值方法是从三维可压缩N-S方程出发，采用五阶精度的迎风紧致格式离散对流项，六阶精度的对称紧致格式离散粘性项，三步三阶精度的R-K方法离散时间导数项[10]。 　　本文所模拟的流动问题、流动参数、网格划分、边界条件的具体处理方法等可以参见文献[11]。 　　2 拟序结构的分析 　　在文献[11]中，通过对流向涡的分析，以及一些湍流统计平均结果与分析，可以知道在本文的计算条件下，平板湍流边界层已经得到了充分地发展，所采用的计算域和网格密度是满足流场的要求的。因此，我们可以对可压缩平板湍流进行更深入的研究。在以下讨论中，x、y和z分别代表流向、平板法向和展向的坐标方向。 　　近壁速度的条带结构是壁湍流最重要的拟序结构之一。图1呈现了脉动速度u在充分发展湍流区（8.5≤x≤9.6）不同y+处的（x，z）平面上的等值线分布，图中黑色的条纹对应u0。这里，u0表示当地流体速度大于平均速度，并称为高速区。 　　从