湍流理论学习

湍流理论学习1．层流和湍流粘性流体的运动存在着两种完全不同的流动状态：层流状态和湍流状态。雷诺首先于1883年通过做圆管内流动实验观察到层流与湍流现象。当圆管中流动速度较小时，管中的流线之间层次分明，互不掺混，这样的流动称为层流。当流速增大后，流体作复杂、无规律、随机的不定常运动，称为湍流。流动状态与雷诺数eR、下临界雷诺数ecR和上临界雷诺数ecR有关。当eecRR时，流动为层流；当eceecRRR时，流动为不稳定过渡状态；当eecRR时，流动为湍流。湍流是在连续介质范畴内流体的不规则运动，它有别于物质分子的不规则运动。具体来说，在极不规则的湍流中，流动的最小时间尺度和最小空间尺度都远远大于分子热运动的相应尺度。因此湍流运动产生的质量和能量的输运将远远大于分子热运动产生的宏观输运。2．湍流的平均化、雷诺粘性应力经典的湍流理论认为，湍流是一种完全不规则的随机运动，湍流场中的物理量在时间和空间上呈随机分布，不同的瞬时有不同的值，关注某个瞬时的值是没有意义的。因此，雷诺首创用统计平均方法来描述湍流的随机运动，即对各瞬时量进行平均得到有意义的平均值。从N-S方程出发，利用平均化运算的法则推导平均物理量满足的方程组。只考虑不可压缩流体情形，假设体力可以忽略，此时，N-S方程具有下列形式1110uuuupuvwutxyzxvvvvpuvwvtxyzy(1)运用(1)式中的连续性方程，运动方程可改写为2221110uuuvuwputxyzxvuvvvwpvtxyzywuwvwwpwtxyzzuvwxyz(2)对方程组(2)中各式两边进行平均化运算，并利用平均化运算法则得到____22____22____221110uuuvuwuuvuwputxyzxyzxvuvvvwuvvvwpvtxyzxyzywuwvwwuwvwwpwtxyzxyzzuvwxyz(3)考虑到方程组(3)的第四式，方程组(3)中的头三个方程可改写成另一种形式，把脉动项移到右边，得到____2____2uuvuwuuuupuvwutxyzxxyzvuvvwvvvvpuvwvtxyzyxyz____20wuwvwxyzuvwxyz(4)将上式和应力形式的运动方程对比ddivdtVP其中P是应力张量，则有2pPISP其中I是单位张量。S是平均运动变形速度张量，P是对称的二阶张量，它的表达式如下___________2___________2___________2xxxyxzyxyyyzzxzyzzuuvuwuvvvwuwvwwP(5)由(5)式得到在湍流运动中除了平均运动的粘性应力外，还多了一项由于脉动所引起的应力，称为雷诺粘性应力。平均流动元所满足的方程(4)称为雷诺方程。上述雷诺方程是不封闭的，方程个数只有4个，而未知函数有10个，即V、p及六个雷诺粘性应力分量。为了使方程组封闭，必须在雷诺应力即平均速度之间建立补充关系式。3．湍流模型所谓湍流模型理论，就是根据理论和经验，对雷诺平均运动方程和雷诺应力项建立表达式或方程，然后对雷诺应力方程的某些项提出尽可能合理的模型和假设，使方程组封闭求解的理论。3.1零方程模型所谓零方程模型，就是在雷诺方程和连续方程以外，不需要另外再加任何方程式来使方程组封闭。即雷诺粘性应力能直接用某些物理量和物理常数表达出来，没有引进高阶微分方程，又称代数模型。(1)Boussinesq涡粘性系数法一般的方法利用Boussinesq假设把雷诺压力和平均速度梯度联系起来：_____23jiiijttijjiiuuuuukxxx(6)(2)Cebeci-Smith模型该模型通常写成涡粘性系数形式，在内层有22()1exp()()mDyUyly(7)式中Dl是衰减长度系数，当壁面不存在质量引射时为常数。在外层有*0.0168mU(8)式中*为边界层位移厚度，U是外部自由流速度，为间隙因子。3.2一方程模型(1)TheSpalart-Allmaras模型Spalart-Allmaras模型的变量中v是湍流动粘滞率除了近壁区域，方程是：221()()()()ivbvvivjjjvvvvuGvCYStxxxx(9)这里vG是湍流粘度生成的，vY是被湍流粘度消去，发生在近壁区域。对于解决动力漩涡粘性，S-A模型是相对简单的方程，不必要去计算和剪应力层厚度相关的长度尺度。(2)简化的二方程模型由于方程难模化和精度低，可以在保留k方程的基础上舍弃方程，将表示为3/2kl(9)式中，l可视为混合长度，必须根据具体流场情况另外给出。c二方程模型(1)标准k模型标准k模型是个半经验公式，主要是基于湍流动能和扩散率。k方程是个精确方程，方程是个由经验公式导出的方程。湍流动能方程k方程为()()()tikbMkijkjkkkuGGYStxxx(10)扩散方程方程为2132()()()()tikbijjuCGCGCStxxxkk(11)方程中kG表示由层流速度梯度而产生的湍流动能，bG是由浮力产生的湍流动能，MY由于在可压缩湍流中，过渡的扩散产生的波动。(2)标准k模型标准k模型是一种经验模型，是基于湍流能量方程和扩散速率方程。()()ikkkkijjkkkuGYStxxx(12)()()iijjuGYStxxx(13)在方程中，kG是由层流速度梯度而产生的湍流动能。G是由方程产生的，k和表明了k和的扩散率。kY和Y由于扩散产生的湍流。4．雷诺应力模式雷诺应力模式简称为RSM，由于雷诺平均运动方程所包含的雷诺应力为一未知项，所以最直接的办法是建立雷诺应力_____ijvv的方程。(1)雷诺应力方程和湍流动能方程得到_____ijvv的方程步骤为：把N-S方程化成()ijDvvDt的形式，将瞬时速度表示成平均速度与脉动速度之和，即iiivvv、jjjvvv，将其代入()ijDvvDt形式方程，减去()ijDvvDt平均运动的形式_____()ijDvvDt，最终得到雷诺应力_____ijvv方程为___________________________________________________()ijijijkkjijijkikijkkkjiikjkkkDvvvvvvvDttxvpvvvpvvvxxvvvvvvxx湍流扩散项分子扩散项_____________________________2jjiikkjivvvvpxxxx耗散项压力变形项产生项(14)在模式理论中要用到湍动能___2/2ikv，在式(14)中令i=j得到___________________________11222jjiijkiijikijijijjDkkkvDttxvpvvvvsvvsssx(15)这就是湍动能方程。(2)经模化后的雷诺应力方程和湍动能方程在雷诺应力方程(14)和湍动能方程(15)中存在未知项，如_______ijkvvv、_____ivp等，因此，需要建立这些未知量的表达式，这一过程称为湍流模式建立的过程。最终经模化后的雷诺应力方程和湍动能方程为：_______________2_____12232233ijijijkijijlllijijijijkDvvvvvvkCDtxxxCvvkCkPPP(16)2kklllDkkkkCDtxxxP(17)其中_____________________,jiiijikjkkilkklvvvvvvvvvxxxPP(18)式中经验常数120.09~0.11,1.5~2.2,0.4~0.5kCCC。以上方程中出项的湍流耗散率仍是未知项，建立方程，进行模化得到__22_____12iilllllvDkCCvvCDtxxxkxk(19)式中经验常数120.07~0.09,1.41~1.45,1.9~1.92CCC。联立封闭微分方程组，这就是雷诺应力封闭模式。