Openfoam使用SST k-omega湍流模型相关问题。

[Joker]

我在做射流模拟的时候，使用了SST k-omega湍流模型，控制了y+～1，但模拟出来的边界层和自由剪切层与实验相比厚了两倍，没有使用壁面函数，这种情况下该从哪个方面着手去解决呢？湍流边界条件omega是根据文献中建议的值去给定的。

[EricLiu]

@Joker 可以讨论一下吗，我在做液膜射流问题，也打算换k-omega做，但是初始值边界条件的确定遇到了点问题

[Joker]

@EricLiu 什么问题？

[EricLiu]

@Joker 目前我在用hybridCentralFoam这个求解器算超音速下液膜演化问题，就是分析超音速来流对下壁面流出的液膜的影响，这个求解器本身默认的湍流模型是laminar，下面是算出来的液相体积分数，红色就是液膜，很平滑，没有体现出超音速气流对液膜的剪切破碎作用。

所以我打算把laminar换成看k-omega模型，算出来的液膜是下面这样的
感觉液相突然间少了许多，就像凭空消失了一样，而且温度场也跟层流下的差异很大
我是按照公式大概估算了入口的k,nut,omega这些值，所以不知道出了什么问题

[Joker]

@EricLiu 我觉得层流和湍流差异很大很正常吧，而且你要模拟液膜被切碎我感觉已经出来了，有没有文献对照一下？看看流速啥的，还有你的湍流边界条件具体咋给的？

[EricLiu]

@Joker 现在的问题就是文献对不上，只是有类似的样子，但是具体的破碎特点、温度、速度分布都有差异，这个问题我自己再研究研究。您说的湍流边界条件是这样的：
入口的k,omega都是用公式算的fixed value，nut 给的是calculated条件
壁面都用的是壁面函数
出口的k,omega是zero gradient，nut 是calculated。
目前我只根据湍流模型修改了0文件，fvScheme和fvSolution我没有改动，因为这块还是不太懂

[Joker]

@EricLiu 你的这些设置没有问题，我觉得是网格的问题，你从层流过渡到湍流，网格要更细一些，你算y+了吗？

[EricLiu]

@Joker 我按照超音速来流算的y+，特征长度取的是计算域的x方向长度的10%，跟目前网格的第一层厚度基本上一致，

[Joker]

@EricLiu 你用openFoam自带的那个函数去监测y+，别用你自己算的那个，不准。我个人觉得哈

[EricLiu]

@Joker 好的，感谢大佬建议，我去试一下这个