44第四章,湍流流动的近壁处理壁面对湍流有明显影响。在很靠近壁面的地方,粘性阻尼减少了切向速度脉动,壁面也阻止了法向的速度脉动。离开壁面稍微远点的地方,由于平均速度梯度的增加,湍动能产生迅速变大,因而湍流增强。因此近壁的处理明显影响数值模拟的结果,因为壁面是涡量和湍流的主要来源。实验研究表明,近壁区域可以分为三层,最近壁面的地方被称为粘性底层,流动是层流状态,分子粘性对于动量、热量和质量输运起到决定作用。外区域成为完全湍流层,湍流起决定作用。在完全湍流与层流底层之间底区域为混合区域(Blendingregion),该区域内分子粘性与湍流都起着相当的作用。近壁区域划分见图4-1。图4-1,边界层结构第一节,壁面函数与近壁模型近壁处理方法有两类:第一类是不求解层流底层和混合区,采用半经验公式(壁面函数)来求解层流底层与完全湍流之间的区域。采用壁面函数的方法可以避免改进模型就可以直接模拟壁面存在对湍流的影响。第二类是改进湍流模型,粘性影响的近壁区域,包括层流底层都可以求解。对于多数高雷诺数流动问题,采用壁面函数的方法可以节约计算资源。这是因为在近壁区域,求解的变量变化梯度较大,改进模型的方法计算量比较大。由于可以减少计算量并具有一定的精度,壁面函数得到了比较多的应用。对于许多的工程实际流动问题,采用壁面函数处理近壁区域是很好的选择。45如果我们研究的问题是低雷诺数的流动问题,那么采用壁面函数方法处理近壁区域就不合适了,而且壁面函数处理的前提假设条件也不满足。这就需要一个合适的模型,可以一直求解到壁面。FLUENT提供了壁面函数和近壁模型两种方法,以便供用户根据自己的计算问题选择。4.1.1壁面函数FLUENT提供的壁面函数包括:1,标准壁面函数;2,非平衡壁面函数两类。标准壁面函数是采用LaunderandSpalding[L93]的近壁处理方法。该方法在很多工程实际流动中有较好的模拟效果。4.1.1.1标准壁面函数根据平均速度壁面法则,有:**1ln()UEyk4-1其中,1/41/2*/ppwUCkU,1/41/2*ppCkyy,并且k=0.42,是VonKarman常数;E=9.81,是实验常数;pU是P点的流体平均速度;pk是P点的湍动能;py是P点到壁面的距离;是流体的动力粘性系数。通常,在*30~60y区域,平均速度满足对数率分布。在FLUENT程序中,这一条件改变为*11.225y。当网格出来*11.225y的区域时候,FLUENT中采用层流应力应变关系,即:**Uy。这里需要指出的是FLUENT中采用针对平均速度和温度的壁面法则中,采用了*y,而不是y(/uy)。对于平衡湍流边界层流动问题,这两个量几乎相等。根据雷诺相似,我们可以根据平均速度的对数分布,同样给出平均温度的类似分布。FLUENT提供的平均温度壁面法则有两种:1,导热占据主要地位的热导子层的线性率分布;2,湍流影响超过导热影响的湍流区域的对数分布。温度边界层中的热导子层厚度与动量边界层中的层流底层厚度通常都不相同,并且随流体介质种类变化而变化。例如,高普朗特数流体(油)的热导子层厚度比其粘性底层厚度小很多;对于低普朗特数的流体(液态金属)相反,热导子层厚度比粘性底层厚度大很多。1/41/2*()wppPTTcCkTq4-246=1/41/2*2*1/41/2221PrPr21Prln()1PrPrPr2pptptptcCkyUqEyPkCkUUq****()()TTyyyy4-3其中P的计算采用下列公式[L93]1/41/2Pr/41sin(/4)PrPrttAPrPk4-4其中,fk是流体导热系数;是流体密度;pc是流体定压比热;q壁面热流;pT近邻壁面控制体温度;wT壁面温度;Pr/pfck为分子普朗特数;Prt是湍流普朗特数,壁面取0.85;A=26,是VanDries常数;k=0.42,是VonKarman常数;E=9.793,是壁面函数常数;cU是**Tyy时的平均速度大小。Fluent中,当选择了流体介质后,就可以根据流体介质的物理性质,计算出分子普朗特数,热导子区厚度*Ty,存储备用。在求解的时候,根据*y与已经存储的*Ty之间大小关系,判断是采用线性法则还是对数法则来计算壁面温度WT或热流率q。在采用雷诺应力模型或k双方程模型时,包括壁面近邻的控制体的湍动能都要计算,其边界条件为湍动能在壁面法向方向上梯度为零。湍动能产生项kG及耗散率是湍动能输运方程的源项组成部分,根据局部平衡假设来计算。根据这一假设,与壁面毗邻的控制体种湍动能及其耗散率是相同的。则湍动能产生率为:pp4-5耗散率不需要求解输运方程,直接用如下公式计算:pppkykC2/34/34-6以上所介绍的标准壁面函数是FLUENT程序的默认设置。标准壁面函数包含了定常剪切和局部平衡假设条件,如果壁面有很强的压力梯度,并且很强的非平衡性,则我们可以选择非平衡壁面函数方法。4.1.1.2非平衡壁面函数47在非平衡壁面函数方法中,平均温度的壁面法则与标准壁面函数中相同。而对数分布的平均速度对压力梯度更加敏感:ykCEkkCUw2/14/12/14/1ln1/~4-7式中,22/1*2/1*ln21~ykkyyyykkydxdpUU4-8y是物理粘性底层厚度,用下式计算:2/14/1*pkCyy4-9其中,225.11*y。非平衡壁面函数在计算近壁控制体湍动能时采用了双层的概念,并且需要求解湍动能k。假定与壁面毗邻的控制体积是由粘性底层和完全湍流构成,则湍流量由如下公式得到:wt0yyyyppkkyyk2)(yyyyyCkykl2/322yyyy4-10式中,4/3kCCl,y是有量纲的粘性底层厚度,2/14/1*pkCyy。利用上面的公式,近壁控制体里面的控制体平均湍动能产生率及其耗散率就可以计算出来。这里我们可以看出,非平衡壁面函数抛弃了标准壁面函数中的局部平衡假设,从而可以考虑非平衡的影响。标准壁面函数对于高雷诺数流动问题,有壁面作用的流动过程等有较好的计算结果;非平48衡壁面函数则把壁面函数方法推广到有压力梯度和非平衡的流动过程中。但是,如果流动情况偏离了壁面函数的理想条件,则壁面函数就不合适了。如:高粘度流体流过狭窄的通道,壁面由渗透的流动,大压力梯度并导致边界层分离的流动,由强体积力的流动,近壁区域三维性很强的流动问题。如果要成功解决上述问题,必须采用改进模型的方法来模拟近壁流动。FLUENT提供了双层区模型(Two-LayerZonalModel)。4.1.1.3双层区模型在双层区模型中,认为近壁流动只分两个区域,即粘性影响的区域和完全湍流,用基于到壁面距离y的雷诺数yRe来区分两个区域。ykyRe4-11其中,y是计算网格到壁面的垂直距离;FLUENT中,y是到最近壁面的距离:wrrrywwmin4-12式中,r是点在流场中的位置矢量;wr是在边界上的位置矢量;w是所有壁面边界的集合;这样,我们可以去处理流场里有复杂边界的问题。而且,这样定义y跟网格的形状没有关系,对非结构网格也同样适合。在完全湍流区域(200Rey),采用雷诺应力模型或者k模型;在粘性影响区域(200Rey),采用Wolfstein[L181]的单方程模型。动量和湍动能输运方程跟前面介绍的没有区别,但计算湍流粘性系数的方法不同。这里湍流粘性系数计算公式如下:utlkC耗散率计算49lk2/34-13上面的长度尺度根据参考文献[L29]的方法计算:]1[ReuyAlueycl4-14]1[ReAlyeycl4-15如果所有的计算区域都在粘性影响的区域以内(200Rey),耗散率的输运方程并不需要求解,而是用上面的代数方程来就得。上面长度尺度计算过程中的模型常数采用ChenandPater[L29]的结果。4/3kCcl,70uA,lcA2表4-1,几种壁面处理方法比较优点缺点标准壁面函数应用比较多,计算量较小,有较好精度适合高雷诺数流动,对低雷诺数流动问题,有压力梯度,强体积力及强三维性问题不适合非平衡壁面函数考虑了压力梯度,可以计算分离,重附及撞击问题对低雷诺数流动问题,有较强压力梯度,强体积力及强三维性问题不适合双层区模型不依赖壁面法则,对于复杂流动,特别是低雷诺数流动很适合要求网格密,因而要求计算机处理时间长,内存大。第二节,湍流计算中近壁处理对网格的要求一个成功的湍流计算离不开好的网格。在许多的湍流中,空间的有效粘性系数不同,是平均动量和其它标量输运的主要决定因素。因此,如果需要有足够的精度,这就需要保证湍流量要比较精确求解。由于湍流与平均流动有较强的相互作用,因此求解湍流问题比求解层流时候更依赖网格。你可以用后处理面板去画出y,*y和yRe的值来检查网格是否满足自己的计算要求。需要指出的是计算出来的y,*y和yRe并不是只与几何参数有关的固定量,它们也和最后的收敛解解有关系。所以,如果你把网格加密一倍(到壁面的距离减少一半),计算得到的y并不一定是加密前计算的y的一半。对于近壁网格而言,不同的近壁处理对网格要求也不同。下面对常见的几种近壁处理的网格要求做个说明。采用壁面函数时候的近壁网格:50第一网格到壁面距离要在对数区内。通常计算的距离为y(/yu或*y。如果网格在对数区内,y和*y的值差不多大小。我们知道,对数区的y30~60。FLUENT在y12.225时候采用层流(线性)准则,因此网格不必要太密,因为壁面函数在粘性底层更本不起作用。对数区与完全湍流的交界点随压力梯度和雷诺数变化。如果雷诺数增加,该点远离壁面。但在边界层里,必须有几个网格点。壁面函数处理时网格划分采用双层模型时近壁网格要求当采用双层模型时,网格衡量参数是y,并非*y。最理想的网格划分是需要第一网格在y=1位置。如果稍微大点,比如y=4~5,只要位于粘性底层内,都是可以接收的。理想的网格划分需要在粘性影响的区域内(200Rey)至少有十个网格,以便可以计算粘性区域内的平均速度和湍流量。采用双层区模型时网格划分采用Spalart-Allmaras模型时的近壁网格要求该模型属于低雷诺数模型。这就要求网格能满足求解粘性影响区域内的流动,引入了阻尼函数,用以削弱粘性底层的湍流粘性影响。因此,理想的近壁网格要求和采用双层模型时候的网格要求一致。51采用大涡模拟的近壁网格要求对于大涡模拟,壁面条件采用了壁面法则,因此对近壁网格划分没有太多限制。但是,如果要得到比较好的结果,最好网格要细,最近网格距离壁面在y=1的量级上。第三节,用FLUENT求解湍流问题设定1,击活粘性模型面板上的湍流模型(Spalart-Allmaras,k-epsilon,ReynoldsStressorLargeEddySimulation),如果选择k-epsilon模型,将需要继续选择采用标准k模型、重整化群k模型或可实现k模型中的一种。2,如果流动问题中包含壁面,选择k或者雷诺应力模型,在粘性模型面板上选择近壁处理方式。近壁处理方式包括:标准壁面函数;非平衡壁面函数和双层区模型。3,Option选项设置4,变量的边界条件设置5,全场变量赋初始值任意选项设置如果选择Spalart-Allmaras模型,会出现如下选项:.Vorticity-basedproduction.Strain/vorticity-basedproduction.Viscousheating(alwaysa