@fangyuanaza :146: 非常感谢！

@fubianhanshu 好的我自己试一下，谢谢各位前辈 @星星星星晴 @fubianhanshu

例如下图所示的波浪-圆柱体耦合计算

波浪圆柱耦合.png

根据理论波浪力方程，圆柱所受到的波浪力（水平方向、垂直方向）为

理论波浪力.png

Fy表达式中的最后一项ρgV0为圆柱体的静水浮力（initial buoyancy in still water）；

请问：OF输出的Fy（垂直方向上圆柱受到的Force）里面包含ρgV0这部分吗？想要获得圆柱受到的垂直方向上的波浪力，是不是需要从OF获得的计算结果中手动减去ρgV0这部分？

@yfclark 谢谢！

十分感谢！@东岳 :146:

http://dyfluid.com/icoFoam.html 另外 这个例子你看过么 最后那个扰流

谢谢李老师一针见血的点评，新版本的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

@同学博 哈哈 好的