@Marcllory 对，氢气才两个质子的大小。空气中氧气、氮气、二氧化碳二三十个质子中字。当空气是个窟窿眼很大的网，氢气在里面穿梭。

如果想确定这种简化误差有多大，可以关注这两个数
Screenshot_20201106_221545.jpg
Screenshot_20201107_095249.jpg

忽略的就是氢气的密度（动量运输）和氢气浓度差造成的扩散（扩散速率）。如果这两个占的比重跟对流比起来很小那就可以忽略。

以上都是我猜的，我没做过这个:chigua2:

有偿请教，问一下有关拉格朗日颗粒喷射的问题，我往喷射模型（manualInjection模型）动态的传递喷射点坐标，Positions文件坐标每一个时间步都会发生变化，怎么实现时刻读取新的喷射点坐标并往出喷粒子呢？

@刘雄国 赞，太感谢了。我赶紧去补一下vector<>:146:

这个问题感觉和我前几天提的问题一样，我已经解决了，原因就是低粘度液体在半径大的管道中就会出现这种速度分布的平台状。
解决办法有两个：提高流体的粘度，或者减小管道的半径尺寸。
你看一下两个模型的尺寸设置是不是一样的或者流体粘度设置的是不是一样。

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就好了，感谢 @麦迪文 师兄的帮助。