CFD中文网

@coolhhh 不是在个人账户下写的，我自己的账户没动过，也能在本地跑。:xinlei:

@李东岳 在 各位老师，请问我是算出棋盘型压强场了吗？ 中说：

知识盲区了老铁。运动大体上没问题就问题不大。不过看起来有一些小bug没处理好，也不能说肯定是压力期盼分布

运动大体上差不多，但是肯定不完全对:xinlei:

@李东岳 在 VOF方法中如何准确计算自由表面曲率？ 中说：

这个方法在函数calculateK()里面，可直接调用。计算结果在3232网格上能得到光滑的结果，但不对。在6464网格上出现了数值振荡。

VOF结果在网格比较细的时候结果不对？这个结论很有挑战性

李老师，我用interIsoFoam算了一下，还是不对。在32$\times$32的网格上可以得到光滑的结果，但是和解析解不一样；在64$\times$64的网格上结果出现了振荡。interfoam也跑了一下，结果和interisofoam不一样，也不对。。。（这个算例的设置参照和解析解在主贴的参考文献[1]中）

算例在这里，大家有兴趣可以跑一下。要是算例设置有问题也可以指出来。另外，我是在v2006环境下跑的(这个链接可能需要注册坚果云，我没找到其他的上传方法:xinlei: )
capillaryWave

@cfdngu 收费也没办法，我之前找了一圈也没有找到比这个好用的

@李东岳 好的，谢谢李老师。我发现面上梯度可通过fvc::reconstruct这个函数重构出的体心梯度，我还是得研究一下这个函数。比如刚刚的重力源项，定义在面上的话，我还没掌握怎么把它放回控制方程里面:143:

单节点4核47秒，双节点每个2核72秒？

@bestucan 在 请问如何理解这一项？ 中说：

我猜是任意方向的加速度。

速度本身是矢量，做梯度就是个二阶张量了，前后乘以n和n的转置，应该是类似于对角化操作（基于某方向的对角化）。把各方向的加速度解耦出来，使之不在其他方向有分量。

这个“任意方向”就是右边那项分母上n的方向，咦，n就是法向哈:chouchou:

那应该就是了

但是要注意右边那项分子，速度有下标的，应该是n方向上的速度。
那这一项的意义应该是：“n方向上的速度，在n方向的加速度”

这也和左边对角化操作印证，嗯，就这么猜:xiezuoye:

理解了，这项应该是一个标量，即“沿自由表面法向的加速度”。

我水平有限，还是需要多研究。

非常感谢您的不吝赐教！

额 好像发现个问题，质量守恒方程应该是这个
7973d26e-efef-4af5-9e55-686a447c4996-image.png

其中𝑆(𝜌)与 𝑆(𝛾)形式一样，只不过自变量由相体积分数变成了密度，不知道这样对不对？

感觉有点弄乱了，求大佬指点，不胜感激:135:

@李东岳 在 interfoam控制方程的问题 中说：

@cfdngu http://dyfluid.com/interFoam.html 不是方程26么 这里面的推导跟nonNewtonianIcoFoam是一样的 那个更简单 你可以看看

好的，非常感谢东岳老师

@blookcfd
你好，感谢回复。

我发现我的问题可能是由于snappyHexMesh选定的区域与重叠区域重合导致的，现在已经解决了。

但是也有没解决的部分，我自己在原来overset网格的基础上加了一点东西，自己编译了另一个lib，这个lib与SHM还是有点冲突，运行的时候会报错，但实际上SHM执行了，网格也能用，所以我下一步好好看看我自己编的lib哪里有问题。

当用toposet + createPatch时，提示找不到stl文件，但是我的stl文件就在对应文件夹下放着，所以这个问题我还没找到原因。

我还在找寻替代方法，如果您知道其他方法能把stl放入计算域中，请不吝赐教，谢谢。

@bestucan 好的，非常感谢，我再研究一下~

@bestucan 非常感谢！

@cccrrryyy

行我再考虑下，谢谢

@Stan同学 非常感谢！

@CFDngu 而且forAll操作可以对空容器进行吧。。实在不行你可以给每个空容器加入第一个数-1，而网格编号不可能小于0，你循环的时候判断下网格编号是不是小于0就行了

@东岳

很难申请，不过我会努力的！

@cfdngu 好的，非常感谢您快速的解答，我再尝试尝试。

楼主做波浪的啊。我最近也读了一点文章，准备开始interFoam，想加一点浮冰啊，空气摩擦的影响进去，感觉可以交流一下。不过我看有些文章里面老是提SWAN，对这个你有了解吗

@cfdngu 楼主您好，我现在也有你之前的问题，算了一个和你类似的例子，波浪作用在矩形防波堤下三个自由度的响应，在沿波浪传播方向的那个自由度，无论是加密重叠区域网格，还是改进时间步和库郎数的关系，都没有好转，有了解贴体网格，但还没有采用做验证。以及您说的测点如何设置能在重叠网格上随物体一起运动。期待您的回答，祝您一切都好。

@bestucan 在 OpenFOAM在不同的系统下运行效率不一样？ 中说：

@CFDngu 我猜的，试试，总比再重装好点,,,,，(◔‿◔)

好的，十分感谢，我去研究一下