Fluent壁面函数的选取依据
责任编辑:Cherry     时间:2024-05-31     来源:转载于:CAE交流之家
责任编辑:Cherry
时间:2024-05-31  来源:转载于:CAE交流之家
分类: 技术分享
浏览量: 471

1. Fluent壁面函数

  前面介绍了壁面函数的由来及相关的理论,这里我们介绍Fluent中壁面函数的选取依据。牢记:使用壁面函数的前提是y+>15

  Fluent在两种湍流模型中需要选择壁面函数分别是k-e模型和Reynolds Stress雷诺应力模型,其他的湍流模型不必考虑壁面函数的问题,同时也不必考虑y+问题,我们后面会详细说明。

640.png

Fluent提供了四种壁面函数以供选择,分别是:

  Standard Wall Functions 标准壁面函数

  Scalable Wall Functions  扩展壁面函数

  Non-Equilibrium Wall Functions 非平衡壁面函数

  User-Defined Wall Functions  自定义壁面函数

2. 标准壁面函数

2.1 Standard Wall Functions

  标准壁面函数是由Launder and Spalding提出的,广泛应用于工业流体流动,是Fluent默认的壁面函数。但是我们计算时尽量不要使用这种壁面函数。

标准壁面函数使用典型的对数律:

640.png

其中

640.png

  式中,κ= 0.4187为卡门常数;E= 9.793为经验常数;Up为紧邻壁面网格中心速度;kp为紧邻壁面网格中心湍动能;μ为动力粘度。U*和y*是Fluent中的无量纲物理量,等同于u+和y+。详情可查看四十五、壁面函数理论及y+的确定

  标准壁面函数对y+要求非常严格,y+必须大于15,如果低于这个值,求解结果准确性会变得很差。其他的壁面函数对于y+的要求有所宽松,但尽量还是保证y+>15。

  上述公式是对边界层内速度的近似,对于能量方程、组分方程和湍流方程同样有近似规律,这里不做介绍。

2.2 使用限制

  标准壁面函数基于壁面恒剪切应力和局部平衡假设,因此当近壁流动受到很大的压力梯度的影响时(边界层分离),即流动处于非平衡状态时,预测结果可能会不准确。

  其实就是因为标准壁面函数忽略了粘性底层,当粘性底层存在大压力梯度时,计算结果肯定是有问题。

3. 扩展的壁面函数

3.1 Scalable Wall Functions

  顾名思义,Scalable Wall Functions在标准壁面函数的基础之上进行了扩展。当y*<11时,标准壁面函数是无法使用的,而Scalable Wall Functions可以正常使用。

  Scalable Wall Functions对y*进行了一定的限制。

640.png

  式中y*limit=11.25

  如果y*>11.25,就取它自身的值。如果y*<11.25,就直接令y*=11.25。也就是说,y*>11.25时,Scalable Wall Functions和standard wall function功能相同。但是如果y*<11.25,那么小于11.25的网格结果相同,且都等于11.25时的规律。

  举个例子,如果我们有三层网格的y*<11.25,那么这三层网格的y*直接按照等于11.25计算的,计算结果肯定是相同的。

3.2 使用限制

  在y*问题上,Scalable wall functions比标准壁面函数应用范围要广,但是和标准壁面函数相同,当近壁流动受到很大的压力梯度的影响时(边界层分离),当流动处于非平衡状态时,预测结果可能会不准确。

4. 非平衡的壁面函数

4.1 Non-Equilibrium Wall Functions

  由于Standard Wall Functions和Scalable wall function对于壁面压力梯度较大时都不适用,因此需要提出一种新的方式来解决这个问题。Fluent提供了Non-Equilibrium Wall Functions。

  Non-Equilibrium Wall Functions基于两层假设来计算壁面剪切应力τw、湍动能k和湍动能耗散率e。而和压力梯度相关性不大的物理量如能量方程、组分方程等则和标准壁面函数保持一致。

640.png

  yv为粘性底层的厚度,y为网格到壁面的距离。y处于粘性底层和粘性底层之外时,分别使用不同的公式来描述流动。

4.2 使用限制

  通过这种方式,Non-Equilibrium Wall Functions能够弥补标准壁面函数的缺陷,适用于分离、撞击等复杂流动

640.png


5. 标准壁面函数 VS扩展的壁面函数 

  为了对比上述的内容,我们使用一个案例加以说明,案例的源文件在公众号文章chapter45中。

5.1 网格情况

  模型为二维平板,长2.5m,宽0.1m,进口流速为6m/s,物性参数保持默认。通过计算可知y+>15,第一层网格高度取2.1E-03m;y+=1时,第一层网格高度取8E-5m。分别对上下两个壁面进行不同的y+网格划分如下图

640.png

  以下使用标准壁面函数进行计算。

5.2 两壁面Y+

  上壁面y+基本等于1,很小;下壁面y+整体都大于15。

640.png


5.3 速度云图

  直观上看,上下壁面的速度云图并不相同,主要还是因为标准壁面函数对y+非常敏感。通过前面的分析,我们应该知道下壁面的速度分布更加合理

640.png

  通过x=1截面上的曲线图也能看出,两壁面附近速度有所区别。

640.png

5.4 扩展的壁面函数对比

  标准壁面函数和扩展的壁面函数明显有所不同。两种壁面函数在y+>15处速度曲线几乎重合,而在y+很小时差距变大。Scalable wall functions在y+较小时相对更准确一些。

640.png

6. 壁面函数的使用限制

  尽管基于标准壁面函数做了很多改进如Scalable Wall Functions和Non-Equilibrium Wall Functions,但是壁面函数仍然存在一些问题。

壁面函数基于对数律,要么忽略粘性底层,要么对粘性底层进行修正,对于粘性底层的求解仍然不够精确,因此对于以下问题,壁面函数并不适用:

  很低的雷诺数流动,如毛细现象

  壁面相变问题,如壁面沸腾现象

  大压力梯度导致的边界层分离现象

  依靠体积力驱动的流动,如自然对流,浮力等

  对于3D模型,边界层歪斜度较大也不适用壁面函数

640.jpg

  既然壁面函数存在一些适用不了的工况,那么我们就想研究这样工况应该怎么办呢??

  还记得四十五、壁面函数理论及y+的确定文章,对于边界层细节捕捉问题,其实是有两种处理方法的。第一种就是我们刚刚介绍的壁面函数的方式,第二种是我们刚开始就想到的加密网格的方式。

  本来为了减少网格数量,我们想到使用壁面函数。现在壁面函数无论如何满足不了需求了,我们就只能回归老本行,通过加密网格的方式来捕获细节。

  Fluent提供了两种方式用来专门捕获壁面处细节Enhanced Wall Treatment和Menter-Lechner。这部分我们下篇文章再详细讲解。

7. 壁面函数总结

  1) 壁面函数只会出现在k-e模型和Reynolds Stress雷诺应力模型

  2) Standard Wall Functions:适用于高雷诺数流动,要求y+>15

  3) Scalable Wall Functions:也适用于高雷诺数流动,但对于y+要求比较宽松,但尽量满足y+>15。不要用Standard Wall Functions,而尽量选择Scalable Wall Functions

  4) Non-Equilibrium Wall Functions:适用撞击、分离等问题,y+<15也可以使用。

  5) 如果不想考虑那么多,就直接使用Scalable Wall Functions

  虽然进行了壁面函数的推荐,但实际上对于k-e模型和Reynolds Stress雷诺应力模型,Fluent推荐不要使用壁面函数,而使用近壁面处理。限于篇幅,下篇文章详细讲解。



来源:转载于:CAE交流之家

回复:

Copyright © 2021 .长沙麦涛网络科技有限公司 All rights reserved. 湘ICP备20015126号-2
联系我们