湍流知识笔记.docVIP

湍流基础知识 0 引言 Reynolds在1883年在圆管流动中发现了自然界中两种不同的流动状态，第一种为流体运动比较规则，各层之间不会发生掺混，称为层流；第二种为流体运动呈现高度不规则状态，流体运动过程中各层之间发生掺混，称之为湍流。 在湍流流动中，物理量呈现高频的不规则运动，每个物理量都是随机函数，这种随机性主要具有两方面特点：1）在相同实验，或者外界条件相同的重复实验，空间中某点物理量随时间的变化关系不具有重复性；2）在相同试验，或者外界条件相同的重复实验，取出足够多样本进行统计平均，所得到的平均量与样本无关。 在实际问题中，与高频无规则而且无法充分的脉动相比，人们更关系湍流流动中可重复的平均量的变化。在实际应用中主要存在三种平均方法：1）样本平均：取出足够多样本进行平均；2）时间平均：在一次实验中，取物理量在某时间段随时间变化关系，并对其进行时间平均，上述时间段应该是远大于脉动时间尺度，而又远小于平均运动时间尺度的物理量，由于在湍流运动中，平均运动和脉动的时间尺度通常相差较大，因此该值在理论上存在；时间平均方法适用于定常流动情况，例如湍流边界层流动；3）空间平均：在一次实验中，取物理量在某空间范围的变化关系，并对其进行空间平均，上述区域应该是远大于湍流脉动的空间尺度，并且远小于平均运动的空间尺度；空间平均适用于均匀流动情况，如管流。 各态历经假设：假定在多次重复实验中出现的所有可能状态，在一次实验中（时间足够长或空间范围足够大）即可以相同概率出现，那么采用一次实验即可完成湍流统计平均量的研究，这样就大大减少了实验次数。 采用上述平均方法，那么湍流变量就可以分解为平均量与脉动量两部分，我们关系的是平均量的演化关系，而脉动量则需要更关系其平均值，实际上这种平均方法就可以知道，单一脉动量的平均值为0，不过脉动量之间的乘积的平均量就不为0，而且，这些值还会对平均量的运动产生影响，从而使得湍流运动与层流运动产生本质不同，那么这种不同到底是什么原因呢，雷诺通过将NS方程进行时间平均的方式进行了说明，并由此开始了湍流的研究。 雷诺平均方程 雷诺认为，湍流脉动只不过是时间尺度和空间尺度与通常所关注的尺度要小得多，其时间尺度和空间尺度还没有达到分子自由运动的程度。因此，湍流运动仍然满足NS方程，于是就可以采用将NS方程平均的方式进行湍流脉动量影响的研究。例如，对于不可压缩流动得分NS方程为， 采用时间平均方式进行变量平均： 代入控制方程，并对控制方程进行平均可得： 其中，。因此，进一步可得RANS方程为： 由此可见，与层流流动相比，湍流时均量仅在方程中增加了一项应力项，称之为雷诺应力。该项由脉动速度组成，表征了脉动速度对于时均量的作用效果，为应力量纲。在方程中引入雷诺应力以后，我们知道，方程中增加了三个方向的脉动速度，已经不再封闭，为了对方程进行求解，就必须对湍流应力进行模拟。湍流的数值模拟工作主要是对雷诺应力项的模拟，该模拟方法的不同，也对应了不同的湍流模拟方法，例如零方程（代数方程）模型，一方程模型，两方程模型，雷诺应力模拟等等。 雷诺应力的物理含义 根据上述推导过程，我们知道雷诺应力是脉动量对于时均量的影响，但是这种影响机制是什么，暂时还不得而知，其物理含义如何也不得而知。在平板湍流边界层流动中，我们可以对雷诺应力增加一些直观认识。 根据不可压缩湍流边界层实验结果，我们知道边界层内的流动速度从壁板开始逐渐增加，一直增加到边界层外层。湍流边界层可以分为内层和外层两部分：内层通常占据整个边界层的20%，其速度会恢复为来流速度的70%，内层分为粘性底层，该层粘性应力占据主导地位，y+在5~8以内，该层内；对数率层，该层湍流应力占据主导地位，y+在30~50以外，y+为1000或20%边界层厚度以内，该层速度分布为：。两者之间还存在过渡层，称之为缓冲层。外层也可以分为两个部分，尾迹率层，该层湍流应力已经开始减弱，流动开始偏离对数率，但仍然为全湍流状态，该层为20%到40%边界层厚度之间，速度分布为：；粘性顶层：该层显著特征是流动出现间歇性，非湍流的流动会间歇性地嵌入边界层内部区域，湍流与非湍流之间发生掺混，使湍流强度减弱，其边界是犬牙结构很不规则，该层位于40%到1边界层厚度之间，速度分布为：。 根据实验结果[吴子牛:空气动力学P404]，可以得到一些半经验公式。当雷诺数小于10^7时，整个不可压湍流边界层的速度分布可以近似描述为： 边界层厚度可以表示为： ， 壁面摩擦系数为： 根据边界层流动特征，均匀流的纵向梯度值远大于横向梯度值，横向速度远大于纵向速度。根据这个特点，对RANS方程进行化简，可得： 方程与层流运动相比增加了雷诺应力一项，方程中左边为对流项，对流项中表达了惯性力对流动的影响，根据连续性方程并对其进行量级分析，可知惯性力的