tahts difficult to answer, I would suggest to get to know the algorithm first, if you want to know the settings, the knowledge should come to the first place. So, get some books to have a look at the algorithm, afterwards the settings shouldn't be a problem.

@东岳 是的，粘度调高也可以出现抛物线，感谢 感谢。我把网格尺度变小了，也能观察到。

:143:

我的上述计算中没有用到估算的方法，那论文里的是不是因为边界层网格采用的估算的行为，所以才会导致friction Reynolds numbe因人而异，出现微小的差别呢！

@东岳 李博，我准备试试，之前写过fluent多面体的网格转换程序，不过这个发不了论文:mihu:

@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: