湍流理论和湍流模型(博士课程课件)

湍流理论和湍流模型西北工业大学2012年3月许和勇绕圆柱的理想流动：(a)无升力流动(b)有升力流动1(c)有升力流动2(d)有升力流动30Re44Re4040Re1903.5*105Re3*106103Re2*105绕圆柱的真实流动(P257)：Re=1.54Re=26Re=140粘性流体运动的两种流态------层流和湍流雷诺实验：1883年圆管内流动实验层流：管中水流稳定地沿轴向运动，流线之间层次分明、互不掺混，流体质点没有垂直于主流方向的横向运动；湍流：流体作复杂的、无规则的、随机的非定常运动，也称紊流;上临界流速：层流变湍流下临界流速：湍流变层流'ccVVdVdVRVcccc'''',dVdVRVcccc,)2320(cReRe)13800('ccReReRe)13800('cReRe流动为层流流动为不稳定的过渡状态流动为湍流工程上，将下临界雷诺数作为流态的判断依据。1.1湍流的不规则性湍流速度场是时间、空间坐标、实验次数的不规则函数)~,,(txuuii在不规则湍流中，流动的最小时间尺度和最小空间尺度都远远大于分子热运动的相应尺度，因此湍流运动产生质量和能量的输运远大于分子热运动产生的宏观输运，所以湍流场中质量和能量的平均扩散远大于层流扩散。随机变量的概率随机变量的概率密度TiNNuP/)(uuPup/)()(1.2湍流的统计随机变量：湍流速度变量u的实数集合，可表示为u(ω)事件集合：相同边界条件下不同初场演化出的所有流场状态系综：所有可能实现的事件集合举例：在相同边界条件下，N个真实初始条件产生N个实验流场（理论上N可以无穷大）是一个系综，其中某一次实验称为一个事件。概率的定义：规定全系综的测度为1，则随机变量u的概率P(x)定义为一切ux事件的测度M，又称为累积概率，可写作所有实验次数的实验次数出现xuxuMxP])~([)(概率密度的定义：如果累积概率P(x)是可微函数，则它的导数定义为概率密度，并用p(x)表示，即dxxdPxp)()(dxxpuuuPuu)()(2121累积概率可表示为联合概率：两个随机变量的累积概率所有实验次数的实验次数和同时出现和yvxuyvxuMyxP])~()~([),(联合概率密度：如果联合概率函数是可微的，则可定义联合概率密度函数为yxyxPyxp),(),(随机函数（或随机过程）的概率和概率密度：是随机变量中相应定义的推广，可以对每一时刻t给出u(ω,t)的概率所有实验次数的实验次数时刻出现在事件的测度时刻在xutxutMtxP)(),(xtxPyxp),(),(概率密度为所有实验次数的实验次数时刻出现和时刻出现在2221112121),;,(xutxutttxxP随机函数的联合概率湍流的统计量平均值：随机变量u依概率密度p(u)的加权积分称为u的期望值，在湍流中称为系综平均值TNiTiuNuduuupE1/)(全系综平均：NiiNuNuduuupEu1/lim)(随机函数或随机过程u(ω,t)的期望值或平均值是确定性变量t的函数，因为NtudutuuptuNiiN1)(lim),()(系综平均是确定性量，即，通过统计平均，不规则的信息已经全部消失，所以，系综平均可以看作一种低通过滤运算。脉动值：随机变量u和它的期望值或平均值之差是随机变量，称为涨落，在湍流中称为脉动脉动量的平均值等于零，因为0'uuuuEEEEu统计矩：随机变量u的n次幂的期望值或平均值称为随机变量u的n阶统计矩，在湍流中，称为n阶自相关量NudutupuuNininn1),(特征函数：概率密度的傅里叶变换称为随机变量的特征函数K(z)dxixzxpzK)exp()()(已知特征函数K(z)，通过傅里叶逆变换，可以求出概率密度dzixzzKxp)exp()(21)(自相关是用统计方法表示随机函数u(ω,t)在不同时刻之间的关系。确切地表示不同时刻的脉动的联系程度可以用自相关系数，定义为)',~(),~(')',;',(')',(tutududuttuupuuttRuu随机函数的自相关函数：随机函数u(ω,t)在时刻t和时刻t’的乘积的统计平均值，称为随机函数u(ω,t)的时间自相关函数，并用Ruu(t,t’)表示2/122/12)',~(),~(/)',(tututtRuu2/122/12)',~(),~(/)',~(),~(tutututu或者性质：),'()',(ttRttRuuuu0),(,0),(tRtRuuuu平稳过程：如果随机过程的自相关函数Ruu(t,τ)只和时间间隔τ有关，则称它为平稳过程，平稳过程有如下定理。各态遍历定理：设随机函数的涨落),~(),~(),~('tututu是平稳过程,即且有)(),~('),~('uuRtutudtRuu),(，则应有0)),~('1(lim02TTdttuT该定理表示平稳过程中随机变量的系综平均等于随机过程的时间平均，这一性质称为随机过程的各态遍历。意义：一次实验中u的时间序列几乎取尽了系综中所有可能出现的值。定常湍流：在时间历程上平稳过程的系综平均不仅可以用长时间平均来取代，而且平均值和时间无关，因此，可以把这种平稳过程简称为定常湍流。空间自相关如果随机变量和空间变量有关，则称它为空间上的随机过程，一般可以写作),~(xu，例如，圆管中湍流的脉动可写作),,,~('zru不同空间位置x1,x2上随机变量的自相关称为空间相关，空间自相关函数为),~('),~('),(2121xxxxuuRuu通常，令x2=x1+ξ，则),~('),~('),(),(11121ξxxξxxxuuRRuuuu如果x1=x2,或ξ=0,则空间自相关函数等于变量u’的2阶矩,即),~(')0,(121xxuRuu它又称为一点空间自相关。空间平稳过程的体积平均如果两点空间相关函数Ruu只和两点的相对位置有关，而和两点本身的空间位置无关，则称这种随机过程为空间平稳过程。即当),(),~('),~('121ξxxxuuRuu时，称),~('xu为空间平稳过程。类似于时间平稳过程中各态遍历定理，可证明0)),~('81(lim2321112233LLLLLLLduLLLixxuuu),~(),~('xx令，代入上式后得到)),~(81(lim112233321LLLLLLLduLLLuixx意义：空间平稳态中某一次实验在空间的分布值几乎遍历随机变量全系综的所有可能状态。空间平稳态的湍流称为均匀湍流。从动力学角度来看，完全均匀湍流必然是衰减的，但是有不少近似均匀湍流的例子，例如，风洞工作段的核心区。小结：一般情况下，湍流量的平均量是指系综平均NiiNuduuupu1/)(在定常湍流中，可以用长时间平均取代系综平均)1(lim0TTudtTu在均匀湍流中，可以用体积平均取代系综平均)),~(81(lim112233321LLLLLLLduLLLuixx时空自相关函数：例如，脉动速度的2阶时空自相关公式为),,~('),,~(),,,(1tututRuuξxxξx湍流的互相关函数：不同随机函数之间乘积的统计平均。湍流运动中流体速度ui、压强p、温度θ等都是随机函数。例如，两个速度分量u1,u2之间的2阶时空相关函数记作),(),(),,~(),,~(),,,(212121tututututRuuξxxξxxξx为简单起见，规定：（1）在以后系综平均表达式中，随机函数中表示系综事件的变量不再明确写出。（2）相关函数中的随机函数均指脉动函数，即平均值等于零的随机函数。1.3湍流脉动的谱1.3.1时间平稳态中的频谱定义：时间相关函数的傅里叶变换称为对应相关变量的频谱。2阶脉动速度的时间相关函数)()()(tutuRuu可变换到频率空间的脉动速度频谱diRSuuuu)exp()(21)(其逆变换为diSRuuuu)exp()()(dSuuu)(2当τ=0时,Suu(ω)表示湍动能在频带中的分布，它在所有频段上的积分等于湍动能的系综平均或时间平均值。时间相关函数与频谱是一一对应的，它们是统计量在时域和频域之间的转换。1.3.2均匀团流场中的波谱定义:空间相关函数的傅里叶变换称为对应相关变量的波数谱，简称波谱或谱。脉动速度的2阶相关函数)()()(ξxxξuuRuu的波谱为3213)exp()()2(1)(dddiRSuuuuξkξk其逆变换为当ξ=0时波谱Suu(k)表示脉动动能在波数段(k,k+dk)中的分布。332211eeekkkk是波数向量，ie是单位向量。321)exp()()(dkdkdkiSRuuuuξkkξ空间相关函数与波谱函数是一一对应的，它是统计量在物理空间和波数空间之间的变换。3212)(dkdkdkSuuuk小结：（1）频谱表示湍流脉动量在时间尺度上的分布频谱中高频成分表示快变的脉动（时间尺度小的脉动）低频成分表示慢变的脉动（时间尺度大的脉动）（2）波谱表示湍流脉动量在空间尺度上的分布波谱中高波数成分表示长度尺度小的湍流脉动低波数成分表示长度尺度大的湍流脉动总之,湍流脉动的谱(频谱和波谱)可以表示湍流脉动强度在各种尺度上的分布.1.5湍流脉动的测量原理湍流脉动的时间序列具有宽频带，测量仪器准确、响应特性好。测量点的脉动速度的时间序列测量方法：热丝风速计法、激光多普勒测速法脉动场的脉动速度的时间序列测量方法：统称为粒子图像测速法(PIV,paticleimagevelocimetry)数据采集的要求(1)测量精度：仪器精度+电子系统的高信噪比和宽频带的频率响应特性(2)采样频率：假如湍流脉动的最高频率为fh，则采样频率至少为2fh假如需要测量脉动量的n阶矩，则采样频率至少需要2nfh条件采样和统计方法思路：根据一定的准则检测湍流信号，当湍流信号满足条件准则时，开始记录一组或几组信号，然后对记录的数据进行统计分析。例如，湍流边界层外层，速度脉动并非始终具有很高的强度，而是间歇性地出现高强度脉动。最简单的条件采样时湍流间歇因子的测量。示性函数：)非湍流状态(01.0//',0)湍流状态(01.0//',1UuUuIUuUuIthth间歇因子：NiNI1/NiNiNtIIff11limNiNiNnIfIf11)1()1(lim湍流状态平均值和非湍流状态平均值：2.0Navier-Stokesequations推导取一质量为m的极小的运动流体单元为研究对象，对其运用牛顿第二定律:F=ma首先，分析其x方向的分量方程:Fx=max第二章湍流运动的基本方程zyxxpDtDuDtDuadxdydzmdxdydzzyxxpdxdydzzdxdzdyydydzdxxdydzdxxpppFzxyxxxxzxyxxxzxzxzxyxyxyxxxxxxxx,)(])[(])[(])[()]([充分小的运动流体单元x方向受力示意图DtDt称为随体导数或全导数，称为局部导数或当地导数。zyxzpDtDwzyxypDtDvzyxxpDtDuzzyzxzzyyyxyzxyxxx