谢谢李老师一针见血的点评，新版本的komegaSST模型是有一个共同的Base的，然后衍生出不同的komegaSST湍流模型，komega湍流模型也是这样的结构，我修改的时候这部分没有处理好。之前觉得就是给动量方程加一项，不懂C++应该也不难，现在发现C++知识严重不足，查了问题都看不懂，需要自己补补基础了。谢谢东岳老师的回复！

确实有很大改善，可能上方空气部分需要设置一个比较大的高度

@东岳 用interDyMFoam算潮流能水轮机，没有仔细看，计算速度大概有不到10%的提高

谢谢李老师的回复！因为雷诺应力不仅出现在动量方程里面，也出现在turbulent kinetic energy equation 与 turbulent frequency equation中，我的想法是如果我分别修改湍流模型的文件（kOmegaSSTBase）与动量方程的文件（linearEddyViscosity）就会要定义两次非线性项，所以我选择就在kOmegaSSTBase中更新turbulence->divDevRhoReff(U),在其中加入非线性项，参考的资料是OpenFoam里面的nonlineareEddyViscosity。由于这样修改需要在kOmegaSSTBase.C文件的constructors中加入eddyViscosity 也就是上图中蓝色后面的部分，我直接在原来的代码Turbulence后面加这一项会报错，想请问这怎么解决？非常感谢李老师耐心看我的问题，不知道我表述清楚没有？

解决了，把1换成1.0就好了，感谢 @麦迪文 师兄的帮助。

@东岳 是的VOF多相流，计算域三层，最上面是空气，红色是油，下面是水2020-10-22 21-10-44屏幕截图.png

@bestucan 你是行家呀，厉害厉害，终于知道原因了

如楼上所说，不指定就是六个自由度，重叠网格也不会出现网格变形的问题

晚上又测试了一下，在linux系统下是可以顺利的进行计算，谢谢东岳老师的帮助:xinxin:

@东岳 李老师，这样计算之后我发现计算结果不对，通过查找资料我发现这个边界条件可以用3b57a6e7-cef1-467b-a636-4f5d687a9af6-image.png 来代替，所以我用codedfixedValue编写了这个边界条件，现在的问题是我好奇codedfixedValue里面的value值必须要设置成$internalField吗？谢谢！

速度平均还是不会做，不过发现湍流强度对计算结果影响很大，大家算两相流的注意一下

比内能就是e啊，compressibleInterFoam顶层就有啊

@炎冰灬灬 检查网格看看呢，你也可以参考下这个帖子的设置
https://www.cfd-china.com/topic/664/cd计算不准的问题-les-re3900-三维圆柱绕流-pisofoam

@同学博 哈哈 好的

@bestucan 先说声谢谢，我去仔细看看

@东岳 用的komega湍流模型，有时在outlet处设置了intOut边界，value会设置成0，也尝试过zerogradient，没有value
，0文件夹的omega设置如下：

internalField uniform 53.6; boundaryField { INLET { type turbulentMixingLengthFrequencyInlet; mixingLength 8.52e-4; value uniform 53.6; } OUTLET { type zeroGradient; } BOTWALL { type omegaWallFunction; value $internalField; } UPWALL { type omegaWallFunction; value $internalField; } HEATWALL { type omegaWallFunction; value $internalField; } }

用pisoFoam计算了两百雷诺数的层流圆柱绕流问题。但是计算出来的升阻力系数出现了拍现象，有大神可以指点一下是什么问题吗。以下是升阻力系数图和算例设置求解器.png 离散.png 层流.png U.png P.png Cl.png Cd.png

@东岳 懂了~ 谢谢智慧的李老师~~

@东岳 刚性翼型用的pisoFoam，柔性翼型用的是fsiFoam。但是pisoFoam算出来的系数波动很大，而fsiFoam算出来系数波动就很小。不知道是为什么

@zhuyuxuan 您好，请问这里“已有的例子只提供了耦合后的icoFoam和simpleFoam求解器”例子是在哪里可以找到？