可压缩湍流边界层的基本方程-中国力学学会

1、Copyright2014版权所有中国力学学会地址：北京市北四环西路15号邮政编码:100190Address:No.15BeisihuanxiRoad,Beijing1001901）第八届全国流体力学学术会议CSTAM2014-S010372014年9月18-21日甘肃兰州可压缩湍流边界层平均场的多层结构对称性吴斌*，毕卫涛*，2），张又升+，佘振苏*（北京大学工学院力学与工程科学系，湍流与复杂系统国家重点实验室，北京100871）+（北京应用物理和计算数学研究所，北京100094）摘要应用余振苏等近期创建的壁湍流结构系综理论（StructuralEnsembleDynamics或SED）研

2、究了零压力梯度超音速平板湍流边界层平均场的多层结构特性。通过分析不同马赫数、雷诺数和壁面温度的可压缩湍流边界层（CTBL）的直接数值模拟数据，证明了在CTBL中，速度和温度混合长存在SED理论所预言的多层结构对称性。基于该多层结构对称性给出了速度和温度混合长的解析表达式，获得了对CTBL平均速度和平均温度剖面的误差小于1%的理论预测。测量了速度和温度混合长的多层结构参数，分析了其随马赫数、雷诺数和壁面温度的变化。结果表明，多层结构参数能够有效刻画CTBL的马赫数效应及壁面温度效应，具有显著的物理意义。关键词边界层，可压缩湍流，多层结构，对称性，平均场引言可压缩湍流边界层是航空航天与能源动力工程

3、中广泛存在的基本流动。构建可压缩湍流边界层的基础理论，对发展湍流工程计算模式、提高气动设计水平具有重要意义。然而，与不可压缩湍流边界层相比，可压缩湍流边界层有着更多的复杂性，使其基础理论的研究非常困难。比如对于最简单的所谓规范湍流边界层（零压力梯度光滑平板湍流边界层），可压缩流动在雷诺数效应之上又增加了马赫数和壁面温度等效应，几种效应（或涡、声和熵模式1）的复杂非线性相互作用，增加了理论分析的难度。已有的理论成果主要是在不可压缩湍流边界层的理论中，以“恰当”的方式引入了气体平均物性参数变化的影响，即采用了Morkovin假设的思想1，2。这些理论存在若干不足，比如通常具有唯像性和经验性的特征，

4、定量上精确度有限、可应用的参数范围不足够明确，也难以对理论进行修正和拓展，等等。最近，余振苏等提出了湍流的结构系统理论（StructuralEnsembleDynamics或SED），为构建复杂湍流的基础理论提供了新的思路3，4。SED理论认为湍流具有内在多尺度、多自由度的（第一类）复杂性和外部多状态、多环境的（第二类）复杂性，是力学系统中典型的复杂系统，应采用复杂系统学的认识论和方法论来研究5，6。简言之，SED理论构建了一元二面的湍流本体论模型和多维多层次的湍流研究方法论7。SED理论认为，湍流的平均和脉动场构成了湍流的“静”“动”二面；其一元为自组织性，体现为湍流场具有的对称性。在SED

5、理论中，对称性被认为是湍流的基本原理，也是湍流不封闭问题中所缺失的定解原理。这种对称性由一个称为“序函数”的变量来表达。序函数是平衡态统计物理学中序参量概念在非平衡系统中的推广，是最能反映脉动场对平均场状态影响的宏观变量。在湍流中，序函数是联系脉动场（湍流未封闭项）与平均场的桥梁，尤其是反映了不同脉动结构系综的标度行为的变化6。因此，识别序函数、确定其对称性成为湍流研究的关键。在这方面，SED理论提出了多维多层次的湍流研究方法论，以及湍流的序函数分析方法，为开展一系列具有工程意义的复杂湍流的研究提供了框架和解决方案。SED理论在规范壁湍流（不可压缩的圆管、槽道和零压力梯度平板边界层）的研究中取

6、得了很大成功6。余振苏等发现在规范壁湍流中混合长等变量具备序函数的特征。通过创新性的应用李群理论，他们构造了三类拉伸群1）基金资助项目：国家自然科学基金（11372008,11221062）2）Email：.cn不变解形式，从而获得了混合长的全场多层结构形式的解析解，由此得到了平均速度剖面、表面摩阻系数等的理论预测，与直接数值模拟（DNS）和实验结果高度一致（精度达99%）。一个标志性的成果是发现了0.45的普适的卡门常数，回答了学界关于幂次律和对数律，以及关于卡门常数取值的争议。SED理论也针对规范壁湍流之外的其它复杂湍流开展了有成效的研究，包括R-B热对流，以及有粗糙度效应8、压力梯度效应

7、和可压缩效应的壁湍流9,10等。同时，SED理论在湍流脉动场的研究、湍流模式理论研究等方面也取得了初步而重要的成果4。鉴于SED理论在不可压缩的零压力梯度平板湍流边界层的研究上取得了极大成功，我们展望其会对可压缩湍流边界层的基础理论研究产生推动。该推动体现在两个层次上，其一是仍然基于Morkovin假设的框架，但由于不可压缩的湍流边界层有了精确的平均场理论，可压缩湍流边界层的理论有望得到提升。在这方面，我们已经开展了一些工作，获得了关于马赫数效应和壁面温度效应的更高精度的理论刻画9，10。其二是发展可压缩湍流边界层的SED理论，即辨识决定可压缩湍流边界层的各个物理维度的序函数，确认其多层结构对

8、称性，研究多层结构参数的演化规律，并且探索背后的物理机制和原理。这个可压缩湍流边界层的SED理论有希望最终从定量上和理论深度方面对经典的Morkovin假设理论形成突破，回答既有理论未能明确的一系列关键科学问题。本研究即依据这个思路开展了初步的探索，通过对零压力梯度可压缩平板湍流边界层的DNS数据开展序函数分析，分别研究了速度和温度混合长的多层结构对称性。1直接数值模拟简介本研究基于对可压缩平板湍流边界层的DNS数据的分析。DNS采用了中科院力学所李新亮研究员等研制的CFD软件OpenCFD-SC。OpenCFD-SC是一个开源的以科学计算为目的的全N-S方程高精度有限差分求解器。利用Open

9、CFD-SC，我们先后在上海和天津的超算中心完成了大规模的并行计算。计算中，N-S方程的时间推进采用了3阶TVD-RK格式；无粘对流项的求解采用了Steger-Warming分裂法和7阶JS-WENO格式；粘性项的计算采用了8阶中心差分格式。本文计算的是空间发展的可压缩湍流边界层11。来流为层流边界层，经吹吸扰动后边界层发生转捩，逐步发展为充分发展的湍流边界层，之后进入出口缓冲区。壁面为无滑移等温边界条件，远场为零压力梯度无反射边界条件，出口为无反射边界条件，展向为周期性边界条件。在流向的充分发展湍流区和法向的近壁区，都对网格进行了加密。流动控制参数、计算区域和网格参数如表1所示。表1DNS算

10、例的基本参数算例MaRetT/TT/TM2.252.25500-1000w81.90wr1.00M4.5T04.5325-6004.390.95M4.5T14.5600-8602.500.54M4.5T24.51950-23501.000.22M6.06.0275-6006.980.93算例LxXLyXLzNxXNyXNz8+X8+X8+M2.2512.23X0.5X0.1754000X85X256xywz13.23X0.9X6.15M4.5T012.46X0.56X0.194500X120X3005.25X0.43X2.74M4.5T114.33X0.56X0.194500X120X3007

11、.97X0.89X5.64M4.5T210.31X0.56X0.0664600X150X30020.54X1.15X7.9M6.013.57X0.6X0.205000X150X3205.63X0.45X2.812可压缩湍流边界层的基本方程-uv+du+(8)零压力梯度的可压缩平板湍流边界层的平均流向动量方程和平均能量方程可以写作1：dUdUdCpu+p=-puv+口二+h.o.t.dxdydyduTOC o 1-5 h zdy丿(1)dHdHddudTpu+pv=(-pHv+_u+k)dxdydydydy+h.o.t.(2)其中上标-表示雷诺平均；上标和”分别表示Favre平均和相应的脉动量：

12、X=pX/-，u2X=XX；h为总焓：H=CT+-；其p2它量的定义遵循惯例。从壁面(y=0)起沿壁面法向将(1)和积分,有：与(8)相比,可压缩湍流边界层的平均密度和平均粘性系数是温度的函数,沿壁面法向有剧烈变化,因此,动量方程(6)需要与能量方程(7)联立求解。这是可压缩湍流边界层的复杂之处。图1画出了同一雷诺数下,不同马赫数的近绝热壁可压缩湍流边界层的无量纲总应力的剖面。总应力的马赫数效应非常微弱,这个结果与MorkovinScaling相符2。SED理论给出了总应力的一个近似模型(亏损标度律模型)：TOC o 1-5 h z1-(料)1-2二I匕)1-2ReOt99由图1可见,模型(9

13、)与DNS符合良好。(9)中的指数1-2存在一定的雷诺数效应,需要进一步研究。0.00.11101001000y+TOC o 1-5 h z一_du_-puv+_=_dy=-Jy(pu+pv)dy+Jh.o.t.(3)0dxdy一_dudTpHv+_u+k=qdydy=-q-Jy(+)dy+Jh.o.t.(4)w0dxdy其中厂为总应力,q为总热流，q=kwdy丿w为壁面平均热流。由于H沁CT+uu，将P代入(4)，可将(4)改写为：dT-CpTv+k+uT=-qpdy(3)和(5)以壁面参数无量纲化后有：-+_du+_dy+idT(5)(6)+tegdubnoitcnufygren图1总应力

14、剖面，蓝线为SED模型。151050-5-10-15Ma=2.25kdTdy.-C+q+updu/dy-u0.11y+qupdu/dy-UMa=4.5ikdTdy-Cqupdu/dy-uMa=6.0+kdTdy*-C101001000-pTv+k+(y-1)Mu+t+Prdy+T=Bq+q.其中，M=u/訂RT为摩擦马赫数，B=q/(CpuT)为壁面传热因子。qwpwTw不可压缩湍流边界层的动量方程的无量纲形式为：(7)图2能量方程(5)中各项的剖面，其中总应力分解为粘性应力和雷诺应力。图中DNS为同一雷诺数(Ret=550),不同马赫数近绝热壁可压缩湍流边界层。能量方程(5)的各项如图2和图

15、3所示。图2为同一雷诺数,不同马赫数的近绝热壁的情形。在近绝热等温壁的条件下,总热流近似为基金资助项目：国家自然科学基金(11372008,11221062)Email：.cn零；在边界层主流区存在一定的马赫数效应。图3是不同壁温的情形，显示壁面存在热流时，能量方程的各项均发生了显著变化。我们也可以针对总热流建立一个与(9)类似的SED经验模型(图3)：q+=1-(二)1(10)ReTT/T=4.39,Re=550对于可压缩湍流边界层，Prandtl混合长如何恰当的定义？是否仍然满足多层结构对称性？这些都是未知的问题。此外，由于气体物性参数的变化，用于封闭雷诺应力的Prandtl混合长不足以求

16、解平均场方程(6)和(7)，还需要针对湍流热流引入相应的序函数。本节将针对这些具体问题进行初步的探讨。3.1速度和温度混合长+tegdubnoitcnufygren3020100-10-20kdTdySEDmodel口c-Cqudu/dy-uT/T=2.50,Re=550-kdTdy-Cqudu/dy-UT/T=1.OO,Re=1750在可压缩湍流边界层中，Prandtl混合长这里称速度混合长)可以写作：1V|l-pwv+m.p+du+/By+(12)+kdTdy+-C+q+udu/dy0.1110100100010000y+图3能量方程(5)中各项的剖面，其中总应力分解为粘性应力和雷诺应力。

17、图中DNS为Ma=4.5，不同壁温的可压缩湍流边界层。3序函数分析不可压缩湍流边界层的SED理论确定了Prandtl混合长是一个序函数，满足如下的三层结构标度律6：1+=b()Y1-(1+()讣2/nimy+y+subsub(1-rm)-(1+(y+buffM)2办3/叮m(1-r)(11)针对零压力梯度绝热壁的可压缩湍流边界层，我们最近还提出了马赫数不变的混合长9：1+=云+1+=-puV+MIm1Bu+(13)由于因子、p+n+t0.261，1+仅是对1+的微MIm弱修正。将(12)代入(6)有：-f-r*du+du+p+/+2()2+g+=V+mBy+By+引入温度混合长1：-pTv+B

18、u+BT+p+-By+By+(14)(15)温度混合长1+和速度混合长1+之间存在关系：Um1+2m=其中l+=m阪+/切为Prandtl混合长；b为近1U+puv+BT+/By+pTv+Bu+/By+=Pr(16)壁系数：G=Ky+J/y+，K为卡门常数；subbuffy+和y+分别为粘性底层和缓冲层的厚度；subbuff丫、丫、丫为标度指数；n、n为跃迁强度；12312r二1-y/8为外区坐标；5为边界层外区尺MM度(外区的拉伸对称中心)；m=4。由(11)、(8)和SED的总应力亏损模型可以预测不可压缩湍流边界层的平均速度剖面。该理论预测得到了美国斯坦福、法国里尔、瑞典KTH和澳大利亚墨

19、尔本等一系列国际知名研究机构的实验数据的支持，理论与实验的误差在2%以内4，6。其中P：为湍流普朗特数。将(15)代入(7)有：2Bu+BT+1BT+p+1+2+k+(V-1)M2u+r+UBy+By+PrBy+T二-Bq+(17)q由于气体物性参数p+、和k+是温度T+的确定性函数，在已知速度混合长和温度混合长的情况下，(14)、(17)和(9)、(10)构成了一组封闭的平均场方程，能够提供平均速度和平均温度分布的理论预测。3.2混合长的多层结构对称性我们仍以不可压缩湍流边界层中Prandtl混合长的三层结构模型来对速度和温度混合长进行模拟。图4给出了不同马赫数近绝热壁可压缩湍流边界层的速度

20、混合长与其SED三层结构模型的比较。模型参数见表2。模型参数的确定过程将另文阐述。图4显示，速度混合长的SED三层结构模型能够在全流场与DNS数据较好地相符（误差小于4%），表明该速度混合长满足SED理论所预言的多层结构对称性。此外，模型给出了一个新的边界层厚度5mo5m是边界层的一个外区尺度，在模型中为外区拉伸对称中心的位置。在有限雷诺数下，5m明显小于传统定义的边界层厚度.9。图5给出了不同马赫数近绝热壁可压缩湍流边界层的温度混合长与其SED三层结构模型的比较。模型参数见表3o图4不同马赫数近绝热壁可压缩湍流边界层的速度混合长及其多层结构模型。图中算例的Ret=550。表2速度混合长的SE

21、D多层结构模型参数参数M2.25M4.5T0M6.0O7Ysub+101010Ybuff+514443Yi丫丫3-1-1-1ni444n2444久/%0.690.690.69m444y+图5不同马赫数近绝热壁可压缩湍流边界层的温度混合长及其多层结构模型。图中算例的Ret=550o表3温度混合长的SED多层结构模型参数参数M2.25M4.5T0M6.0o6Ysub+101010Ybuff+494141YiYY3-0.9-0.9-0.9ni444n24440.690.690.69

22、m444图6马赫数2.25的近绝热壁可压缩湍流边界层的湍流普朗特数及其SED模型。图中算例的Ret=550o-0.041101001000y+(a)温度混合长与速度混合长的模型参数仅有微小的区别，区别主要反映在标度指数y3上。以往的研究表明可压缩湍流边界层的湍流普朗特数约在0.7-0.9之间，明显偏离110。鉴于湍流普朗特数是速度混合长与温度混合长之比的平方(公式(16)，温度混合长与速度混合长的SED三层结构模型也提供了湍流普朗特数的多层结构理论模型。图6给出了马赫数2.25边界层的湍流普朗特数的DNS与SED理论模型的比较，二者非常一致。数据分析表明，湍流普朗特数偏离1与温度混合长的标度指

23、数Y偏离-1有直接关系，其机理有待进一步研究。(b)图7(a)马赫数2.25的近绝热壁可压缩湍流边界层的平均速度剖面及其SED理论预测结果；(b)理论与DNS的相对误差。图中算例的Ret=550。由速度和温度混合长的SED三层结构模型，以及动量与能量方程(14)和(17)、总应力与总热流的亏损模型(9)和(10)，可以得到可压缩湍流边界层的平均速度剖面和平均温度剖面的理论预测结果。图7是仅由速度混合长的SED理论模型预测的平均速度剖面(气体的物性参数采用了DNS值)与DNS结果的比较，二者非常一致，误差在全流场都小于1%(图7b)。该结果表明可压缩湍流边界层的混合长序函数的多层结构理论模型对平

24、均场有极为准确的预测能力。4结论本文尝试将最新发展的规范壁湍流的结构系综理论拓展到零压力梯度可压缩平板湍流边界层。初步的研究表明速度和温度混合长均存在SED理论发现的多层结构对称性，是能够封闭平均场方程的序函数。速度和温度混合长的SED多层结构模型能够提供平均场的精确理论预测，其多层结构参数反映了可压缩湍流边界层的流动状态的变化。下一步的研究内容包括进一步扩展可压缩湍流边界层DNS的参数范围(提高雷诺数和马赫数)，观测序函数的多层结构模型参数随流动控制参数的演化，进而构建完整的可压缩湍流边界层的结构系综理论。致谢：在数值计算上受到李新亮研究员的指导，在此表示感谢。同时受益于与陈曦博士的讨论。也

25、感谢上海和天津超算中心提供了计算资源。感谢来自国家自然科学基金的资助。参考文献SmitsAJ,DussaugeJP.Turbulentshearlayersinsupersonicflow.NewYork:SpringerVerlag,2006,119138MorkovinMV.Effectsofcompressibilityonturbulentflows.MechaniquedelaTurbulence,GordonandBreach,367380Zhen-SuShe,NingHu,YouWu.Structuralensembledynamicsbasedclosuremodelforwa

26、ll-boundedturbulentflow.ActaMechanicaSinica.2009,25(5):731-736佘振苏，飞行器气动力学与光学设计中的关键湍流问题,973项目验收报告，国家科技报告服务系统，20135张志雄.基于结构系综的湍流复杂系统研究.Ph.D.thesis,北京大学,20096陈曦.壁湍流的多层结构理论.Ph.D.thesis,北京大学，20127佘振苏.复杂系统学新框架.中国科学出版社,2012Zhen-SuShe,YouWu,XiChen,FazleHussain.Amultistatedescriptionofroughnesseffectsinturbu

27、lentpipeflow.NewJournalofPhysics,14093054,2012.You-ShengZhang,Wei-TaoBi,FazleHussain,Xin-LiangLi,andZhen-SuShe,Mach-Number-InvariantMean-VelocityProfileofCompressibleTurbulentBoundaryLayers,PhysicalReviewLetters,109,054502,2012You-ShengZhang,Wei-TaoBi,FazleHussain,Xin-LiangLi,andZhen-SuShe.Agenerali

28、zedReynoldsanalogytheoryforthecompressiblewall-boundedturbulence,JournalofFluidMechanics,2014,Vol.739,pp.392-420GaoHui,FuDe-Xun,MaYan-Wen,LiXin-Liang.Directnumericalsimulationofsupersonicturbulentboundarylayerflow.ChinesePhysicsLetters.2005,22(7):1709MULTILAYERSCALINGOFMEANVELOCITYANDTHERMALFIELDSOF

29、COMPRESSIBLETURBULENTBOUNDARYLAYERSWUBin1BIWeitao1ZHANGYousheng2SHEZhen-Su1(1StateKeyLaboratoryforTurbulenceandComplexSystems,DepartmentofMechanicsandEngineeringScience,CollegeofEngineering,PekingUniversity,Beijing100871,China)(2InstituteofAppliedPhysicsandComputationalMathematics,Beijing100094,China)AbstractRecently,asymmetry-basedstructuralensembledynamics(SED)theorywasproposedbySheetal.forcanonicalwall-bounde