明天给你看一下 老铁

感觉就是不同的研究者对这个的规定不一样，论文里面推导过程正确就行。
当然如果按照张量的知识，向量并积计算的话，我觉得应该是
$$
\nabla {\bf{U}} = \frac{{\partial {U_j}}}{{\partial {x_i}}}
$$

谢谢李老师

感谢分享

@liujm 在 编译overWaveDyMFoam失败请求援助 中说：

@chenboyao 在 编译overWaveDyMFoam失败请求援助 中说：

@liujm 你对照当前的wave2foam中waveFoam求解器和interFoam求解器，然后你看看修改了哪儿，你就把修改的地方复制到overInterDyMFoam中，在把修改好的文件复制到wave2Foam存放求解器的地方，重新编译一下wave2foam就行

@chenboyao

Hello, 我说下我的编译方法，避免出错，您看看哪里有问题：

把multiphase/interFoam/overInterDyMFoam文件夹复制到solvers1912_PLUS目录下（与原有waveFoam同级），并将overInterDyMFoam改名为overWaveDyMFoam。

在overWaveDyMFoam.C的基础上对照waveFoam.C，填补缺少的代码。

做到第二步时我有两点疑问：

您提供的文件中没有完全将waveFoam.C中的代码补充到overWaveDyMFoam.C中，例如部分头文件以及如下等代码。请问略去的代码是不需要吗？ #include "readDyMControls.H" phi = mesh.Sf() & Uf(); MRF.update(); 两个文件夹中的部分同名文件也有较大区别，例如correctPhi.H等，这些该如何处理？对于waveFoam中有的文件，而overWaveDyMFoam没有的文件，是否应该复制到后者处？

我这个代码对应的OpenFOAM版本是2206，不知道1912版本适不适配

记录一下，由于fe40太老了，在ubuntu2204下面编译以下问题，给我卡住了1个小时..

sigFpe.C文件要改malloc scotch要用6.0.4，不能用6.0.0 scotch要自己编译 fe41默认lib创建so连接，要不然还是找不到

已更新：http://dyfluid.com/download.html

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

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

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

@李东岳 谢谢李老师回复，问题已经解决了:high:

@五好青年 大佬，这个问题解决了吗:135:

@lizhisongsjtu 我不太关心涡这个东西，有几个疑问，邀请李老师解惑，给我开开眼：

不同的涡判据，比如liutex、lambda、Q之类的，都是用来判别涡的，实际使用中，是哪个看起来更流弊就用哪个么？ :134:

如何判定那个涡识别方法更流弊呢？:134:

@队长别开枪 太流弊了大佬，真羡慕啊

@李东岳 太好了！学习中！:142: :142: :142: :ok3:

好的，谢谢李老师

@yuky 文章里看到过，https://doi.org/10.1016/j.csite.2022.102560 ，记忆中以前的老版本好像有人开发过，可以搜一下

@zhaoym 老铁，后来这个做的怎么样了，有进展么

interPlicFoam求解器已经并入扩展包geometricVofExt，仓库地址：https://github.com/daidezhi/geometricVofExt

@秋山结月 在 200万网格并行算力测试（OpenFOAM版本） 中说：

@tidedrinker 我之前也觉得可能是因为这个，但是事实上多插了两根内存条并没有改善现状。插满内存条的经费我现在也没有，只有先等等看以后能不能重新测个数据了。

4根内存条也不至于这么慢，插4根应该是每个三缓1根，我记得这个型号一共有4块三缓，8个核心共用一块。插4根内存，8核也应该在135s左右。没改善估计是位置插错了，内存通道插满，32核成绩应该在70s左右。

并行计算时候，每个进程掌握的网格就固定那些。
举例b点只存在于processor x1，不在processor x2当中，processor x2自然无法获取b点的信息。
要想processor x2获取b点信息，可以用reduce函数，也就是进程之间通讯一下。
首先初始化intU为零向量，然后赋值，再reduce。

point b (1,1,1); label cellID = mesh.findCell(b); vector intU(0,0,0); intU = U[cellID]; reduce(U,sumOp<vector>());

具体代码你再看一下，我这个是思路，代码可能哪里不对。