中科大FLUENT讲稿_第三章_湍流模型

24第三章，湍流模型第一节，前言湍流流动模型很多，但大致可以归纳为以下三类：第一类是湍流输运系数模型，是Boussinesq于1877年针对二维流动提出的，将速度脉动的二阶关联量表示成平均速度梯度与湍流粘性系数的乘积。即：2121xuuut3－1推广到三维问题，若用笛卡儿张量表示，即有：ijijjitjikxuxuuu323－2模型的任务就是给出计算湍流粘性系数t的方法。根据建立模型所需要的微分方程的数目，可以分为零方程模型（代数方程模型），单方程模型和双方程模型。第二类是抛弃了湍流输运系数的概念，直接建立湍流应力和其它二阶关联量的输运方程。第三类是大涡模拟。前两类是以湍流的统计结构为基础，对所有涡旋进行统计平均。大涡模拟把湍流分成大尺度湍流和小尺度湍流，通过求解三维经过修正的Navier-Stokes方程，得到大涡旋的运动特性，而对小涡旋运动还采用上述的模型。实际求解中，选用什么模型要根据具体问题的特点来决定。选择的一般原则是精度要高，应用简单，节省计算时间，同时也具有通用性。FLUENT提供的湍流模型包括：单方程（Spalart-Allmaras）模型、双方程模型（标准κ-ε模型、重整化群κ-ε模型、可实现(Realizable)κ-ε模型）及雷诺应力模型和大涡模拟。湍流模型种类示意图Zero-EquationModelsOne-EquationModelsSpalart-AllmarasTwo-EquationModelsStandardk-RNGk-Realizablek-Reynolds-StressModelLarge-EddySimulationDirectNumericalSimulation包含更多物理机理每次迭代计算量增加FLUENT提供的模型选择RANS-basedmodels25第二节，平均量输运方程雷诺平均就是把Navier-Stokes方程中的瞬时变量分解成平均量和脉动量两部分。对于速度，有：iiiuuu3－3其中，iu和iu分别是平均速度和脉动速度（i=1,2,3）类似地，对于压力等其它标量，我们也有：3－4其中，表示标量，如压力、能量、组分浓度等。把上面的表达式代入瞬时的连续与动量方程，并取平均（去掉平均速度iu上的横线），我们可以把连续与动量方程写成如下的笛卡儿坐标系下的张量形式：0)(iiuxt3－5jijllijijjijiiuuxxuxuxuxxpDtDu323－6上面两个方程称为雷诺平均的Navier-Stokes（RANS）方程。他们和瞬时Navier-Stokes方程有相同的形式，只是速度或其它求解变量变成了时间平均量。额外多出来的项jiuu是雷诺应力，表示湍流的影响。如果要求解该方程，必须模拟该项以封闭方程。如果密度是变化的流动过程如燃烧问题，我们可以用法夫雷（Favre）平均。这样才可以求解有密度变化的流动问题。法夫雷平均就是出了压力和密度本身以外，所有变量都用密度加权平均。变量的密度加权平均定义为：/~3－7符号～表示密度加权平均；对应于密度加权平均值的脉动值用表示，即有：~。很显然，这种脉动值的简单平均值不为零，但它的密度加权平均值等于零，即：0，0Boussinesq近似与雷诺应力输运模型为了封闭方程，必须对额外项雷诺应力jiuu进行模拟。一个通常的方法是应用Boussinesq假设，认为雷诺应力与平均速度梯度成正比，即：ijiitijjitjixukxuxuuu)(323－8Boussinesq假设被用于Spalart-Allmaras单方程模型和k双方程模型。Boussinesq近似的好处是与求解湍流粘性系数有关的计算时间比较少，例如在Spalart-Allmaras单方程模型中，只多求解一个表示湍流粘性的输运方程；在k双方程模型中，只需多求解湍动能k和耗散率ε两个方程，湍流粘性系数用湍动能k和耗散率ε的函数。Boussinesq假设的缺点是认为湍26流粘性系数t是各向同性标量，对一些复杂流动该条件并不是严格成立，所以具有其应用限制性。另外的方法是求解雷诺应力各分量的输运方程。这也需要额外再求解一个标量方程，通常是耗散率ε方程。这就意味着对于二维湍流流动问题，需要多求解4个输运方程，而三维湍流问题需要多求解7个方程，需要比较多的计算时间，对计算机内存也有更高要求。在许多问题中，Boussinesq近似方法可以得到比较好的结果，并不一定需要花费很多时间来求解雷诺应力各分量的输运方程。但是，如果湍流场各向异性很明显，如强旋流动以及应力驱动的二次流等流动中，求解雷诺应力分量输运方程无疑可以得到更好的结果。第三节，湍流模型3.3.1单方程（Spalart-Allmaras）模型Spalart-Allmaras模型的求解变量是~，表征出了近壁（粘性影响）区域以外的湍流运动粘性系数。~的输运方程为：YxCxxGDtDjbjj~~)~(1~2~3－9其中，G是湍流粘性产生项；Y是由于壁面阻挡与粘性阻尼引起的湍流粘性的减少；~和2bC是常数；ν是分子运动粘性系数。湍流粘性系数用如下公式计算：1~ft其中，1f是粘性阻尼函数，定义为：31331Cf，并且~。湍流粘性产生项，G用如下公式模拟：~~1SCGb3－10其中，222~~fdkSS，而1211ff。其中，1bC和k是常数，d是计算点到壁面的距离；Sijij2。ij定义为：jiijijxuxu213－11由于平均应变率对湍流产生也起到很大作用，FLUENT处理过程中，定义S为：),0min(ijijprodijSCS3－12其中，0.2prodC，ijijij，ijijijSSS2，平均应变率ijS定义为：jiijijxuxuS213－1327在涡量超过应变率的计算区域计算出来的涡旋粘性系数变小。这适合涡流靠近涡旋中心的区域，那里只有“单纯”的旋转，湍流受到抑止。包含应变张量的影响更能体现旋转对湍流的影响。忽略了平均应变，估计的涡旋粘性系数产生项偏高。湍流粘性系数减少项Y为：21~dfCYww3－14其中，6/1636631－15)(62rrCrgw3－1622~~dkSr3－17其中，1wC，2wC，3wC是常数，222~~fdkSS。在上式中，包括了平均应变率对S的影响，因而也影响用S~计算出来的r。上面的模型常数在FLUENT中默认值为：1335.01bC，622.02bC，3/2~，1.71C，~2211/)1(/bbwCkCC，3.02wC，0.23wC，41.0k。壁面条件在壁面，湍流运动粘性~设置为零。当计算网格足够细，可以计算层流底层时，壁面切应力用层流应力－应变关系求解，即：yuuu3－18如果网格粗错不能用来求解层流底层，则假设与壁面近邻的网格质心落在边界层的对数区，则根据壁面法则：yuEkuuln13－19其中，k=0.419，E=9.793。对流传热传质模型在FLUENT中，用雷诺相似湍流输运的概念来模拟热输运过程。给出的能量方程为：heffijjitpiiiSuxTtckxpEuxEt)(Pr)]([)(3－20式中，E是总能量，effij)(是偏应力张量，定义为：28ijiieffjiijeffeffijxuxuxu32)()(3－21其中，effij)(表示粘性加热，耦合求解。如果默认为分开求解，FLUENT不求解处effij)(。但是可以通过变化“粘性模型”面板上的湍流普朗特数（Prt），其默认值为0.85。湍流质量输运与热输运类似，默认的Schmidt数是0.7，该值同样也可以在“粘性模型”面板上调节。标量的壁面处理与动量壁面处理类似，分别选用合适的壁面法则。综上所述，Spalart-Allmaras模型是相对简单的单方程模型，只需求解湍流粘性的输运方程，并不需要求解当地剪切层厚度的长度尺度。该模型对于求解有壁面影响流动及有逆压力梯度的边界层问题有很好模拟效果，在透平机械湍流模拟方面也有较好结果。Spalart-Allmaras模型的初始形式属于对低雷诺数湍流模型，这必须很好解决边界层的粘性影响区求解问题。在FLUENT中，当网格不是很细时，采用壁面函数来解决这一问题。当网格比较粗糙时，网格不满足精确的湍流计算要求，用壁面函数也许是最好的解决方案。另外，该模型中的输运变量在近壁处的梯度要比k中的小，这使得该模型对网格粗糙带来数值误差不太敏感。但是，Spalart-Allmaras模型不能预测均匀各向同性湍流的耗散。并且，单方程模型没有考虑长度尺度的变化，这对一些流动尺度变换比较大的流动问题不太适合。比如，平板射流问题，从有壁面影响流动突然变化到自由剪切流，流场尺度变化明显。3.3.2标准k模型标准k模型需要求解湍动能及其耗散率方程。湍动能输运方程是通过精确的方程推导得到，但耗散率方程是通过物理推理，数学上模拟相似原形方程得到的。该模型假设流动为完全湍流，分子粘性的影响可以忽略。因此，标准k模型只适合完全湍流的流动过程模拟。标准k模型的湍动能k和耗散率ε方程为如下形式：MbkiktiYGGxkxDtDk3－22kCGCGkCxxDtDbkikti2231)(3－23在上述方程中，kG表示由于平均速度梯度引起的湍动能产生，bG是用于浮力影响引起的湍动能产生；MY可压速湍流脉动膨胀对总的耗散率的影响。湍流粘性系数2kCt。在FLUENT中，作为默认值常数，1C＝1.44，2C=1.92，09.0C，湍动能k与耗散率ε的湍流普朗特数分别为k＝1.0，＝1.3。可以通过调节“粘性模型”面板来调节这些常数值。293.3.3重整化群κ-ε模型重整化群κ-ε模型是对瞬时的Navier-Stokes方程用重整化群的数学方法推导出来的模型。模型中的常数与标准κ-ε模型不同，而且方程中也出现了新的函数或者项。其湍动能与耗散率方程与标准κ-ε模型有相似的形式：MbkieffkiYGGxkxDtDk3－24RkCGCGkCxxDtDbkieffi2231)(3－25kG表示由于平均速度梯度引起的湍动能产生，bG是用于浮力影响引起的湍动能产生；MY可压速湍流脉动膨胀对总的耗散率的影响，这些参数与标准κ-ε模型中相同。k和分别是湍动能k和耗散率的有效湍流普朗特数的倒数。湍流粘性系数计算公式为：~1~~72.132dCkd3－26其中，/~eff，100C对上面方程积分，可以精确得到有效雷诺数（涡旋尺度）对湍流输运的影响，这有助于处理低雷诺数和近壁流动问题的模拟。对于高雷诺数，上面方程可以给出：2kCt，0845.0C。这个结果非常有意思，和标准κ-ε模型的半经验推导给出的常数09.0C非常近似。在FLUENT中，如果是默认设置，用重整化群κ-ε