关于OpenFOAM中K-Omega SST模型使用壁面函数时的Y+问题
-
最近在画网格的时候碰到了一个问题,K-Omega SST模型使用壁面函数时要求Y+在哪个范围内?我看到一些资料说是30~300,然后也有说上限为200的。我比较好奇的是,Y+下限是30吗?因为我的网格画出来Y+有一部分小于30的,是否会影响模拟结果?
另外,对于水轮机、风机这类速度沿展弦方向变化很大的情况,边界层厚度是否应该根据速度变化而变化,还是保持一个定值?对于不同区域边界层厚度对应Y+不同(但都在范围内)是否可以?
还有一个问题,OpenFOAM自带的工具好像没办法计算二相流的Y+,就是那个postProcessing -func yPlus命令,是否有什么其他工具可以办到? -
最近我也在用SST k-w模型计算翼型周围的流场,开始用y+大于30,但是我有一部分网格y+会小于30,算出来升力并不能收敛。。而是处于周期性振荡情况。看到researchGate 上面说openfoam里面的sst k-w模型有两种模式?https://www.researchgate.net/post/Should_the_RANS_model_kOmegaSST_be_fully_resolved_at_low_Re-number_Is_it_possible_to_use_it_with_a_wall_function
一种,是低雷诺数模型模式,y+要求在0~5,另一种,是高雷诺数模型模式,y+在30~300;用的边界条件不一样。
但是,原本的sst k-w模型应该是一种低雷诺数模型,应当对边界层求解,高雷诺数模式有人试过吗?是否靠谱?
二相流的y+不能计算吗。。我算的空化,直接在controlDict里设置yplus就可以输出呀。。 -
在controlDict里设置yplus设置?意思是只能在计算中输出?
而且,OpenFOAM的K-Omega SST模型壁面函数好像与其他软件有所不同,可以自行判断使用哪种模式? -
我的算例是在计算中输出的,没什么问题,你可以试试。。
第二个问题也是我想问的,OpenFOAM的sst k-w模型是否会自动识别用哪种模式?感觉此处需要@一些大神 ^_^ -
@lllwonderliquid 想问同样的问题,用sst k-w是否还需要指定用低Re壁面函数?
-
@lllwonderliquid 感谢大神的回复,我这就去试下
-
@lllwonderliquid 你好,请教下,使用壁面函数不就是为k方程壁面边界条件为kqRwallFuction, omega边界条件为omegaWallFuction, nut边界条件为nutkWallFunction吗,这里又说“我后来的做法是将网格做到y+<1,然后不采用壁面函数”。“网格做到y+>30,然后指定壁面函数也是可以的”这里又该如何指定壁面函数呢,望赐教!!!
-
@lllwonderliquid 可以请教下如何确保区域里近壁面的y+都小于1吗,谢谢!
-
@lllwonderliquid 方便po一下controlDict中的yPlus吗,我试了但是没有成功,谢谢!
-
我后来的做法是将网格做到y+<1,然后不采用壁面函数,k方程壁面边界条件为kqRwallFuction, omega边界条件为omegaWallFuction, nut边界条件为nutkWallFunction
另外,我觉得将网格做到y+>30,然后指定壁面函数也是可以的同样有疑问,指定了k w 和nut 的边界条件不就是采用了壁面函数吗?
还是说壁面函数有什么另外的指定方法?
(还是我认知错误) -
@hoversoar 在 关于OpenFOAM中K-Omega SST模型使用壁面函数时的Y+问题 中说:
我后来的做法是将网格做到y+<1,然后不采用壁面函数,k方程壁面边界条件为kqRwallFuction, omega边界条件为omegaWallFuction, nut边界条件为nutkWallFunction
另外,我觉得将网格做到y+>30,然后指定壁面函数也是可以的同样有疑问,指定了k w 和nut 的边界条件不就是采用了壁面函数吗?
还是说壁面函数有什么另外的指定方法?
(还是我认知错误)不知道理解的对不对,当y+<1的时候,可以不采用壁面模型。实际上此时计算使用的是近避免模型,当Y+在1~5区间的情况下,此时是速度是线性规律,30往上是对数率。一般不推荐y+取值在5~30之间。不确定这个答案是否正确,仅仅是我的一些推断。
-
关于OpenFOAM中在采用K-Omega SST模型时,使用/不使用壁面函数的Y+及边界设置问题,可以参考这个链接链接文本,解释的很详细。