好的，谢谢李老师，我试一下

0时间步的U，是用湍流生成方法直接生成的，可以计算得到div(U)。然后就可以计算div(phi)： U -> phi = fvc::flux(U) -> div(phi)。测试结果是，这种方式计算的div(U)和div(phi)结果基本一样的，区别是小数点后五六位。此时div(U)=div(phi)=0

那这个算法很屌。不错。

@李东岳 在 CodedFixedValue边界设置问题 中说：

http://dyfluid.com/code.html 来这查一下

好的，谢谢东岳哥

@Calf-Z-DNS 在 用snappyHexMeshDict生成周期性边界问题 中说：

@东岳 谢谢东岳老师
我在cfd-online上面找到一种方法：
下面是原文：
“I encountered the same problem and I could solve it by defining rotational periodicity and then by declaring the corresponding vertices in ICEM as periodic. Then Right Click on Faces in the Blocking tree -> Periodic Faces -> ensure that really every periodic face is colored accordingly. createPatch went fine after this procedure.” ---by Fabian E.
然后我勾选了ICEM里面blocking里的faces的Periodic，再导入openfoam，发现就不会出现上述错误。
但我不太清楚为何这样会消除错误？

请问同学可以告知是怎么设置的吗，icem里面有face pericdic吗？我这里只设置了点对称

velocity.,txt
(
(0 (4 0 0))
(0.5 (4 0.0698 0))
(1 (4 0.0698 0))
(1.5 (4 0.0698 0))
(2 (4 0.0698 0))
(2.5 (4 0.0698 0))
(3 (4 0.0698 0))
(3.5 (4 0.0698 0))
(4 (4 0.1397 0))
(4.5 (4 0.2096 0))
(5 (4 0.2797 0))
(5.5 (4 0.3499 0))
(6 (4 0.4204 0))
(6.5 (4 0.4911 0))
(7 (4 0.5622 0))
(7.5 (4 0.6335 0))
(8 (4 0.7053 0))
(8.5 (4 0.7775 0))
(9 (4 0.8502 0))
(9.5 (4 0.9235 0))
(10 (4 0.9973 0))
(10.5 (4 1.0718 0))
(11 (4 0.9973 0))
(11.5 (4 0.9235 0))
(12 (4 0.8502 0))
(12.5 (4 0.7775 0))
(13 (4 0.7053 0))
(13.5 (4 0.6335 0))
(14 (4 0.5622 0))
(14.5 (4 0.4911 0))
(15 (4 0.4204 0))
(15.5 (4 0.3499 0))
(16 (4 0.2797 0))
(16.5 (4 0.2096 0))
(17 (4 0.1397 0))
(17.5 (4 0.0698 0))

)
入口设置
{
type uniformFixedValue;
uniformValue tableFile;
uniformValueCoeffs
{
file "velocity.txt";
outOfBounds repeat;
}
value $internalField;
}

@DY大世界 大佬，能给我转发一份吗？我最近也在提取重叠网格的位移，但是of自带的那个函数不知道为啥提取不出来，好像是有点问题，邮箱是1292153997@qq.com

@lcx 在 在coalChemistryFoam中加入MPPIC应力模型 中说：

在1基础上关闭MPPIC颗粒应力仅考虑颗粒所受曳力和重力得到较为合理的结果。颗粒没有出现过度堆积的情况。(gif见下楼)

结果如下
VEDIO3.gif

流场内的模型（stl文件）已经创建好，想在模型外面创建一个矩形的外流场，类似snappyHexMesh先利用blockMesh生成矩形网格。cfmesh有这个功能吗？

好的
我试试

@dxl 您好，在你的 https://www.cfd-china.com/topic/394/如果用dpmfoam求解稀相流会怎么样-误差大么/47 这篇文章中我也遇到了相同的情况可以讨论一下吗

@学流体的小明 感谢感谢:146: :146:

如题，在一个三维计算域中，均匀提取空间中各个位置的数据，除了使用probesDict文件添加监测点外，有没有其他方法实现？

@李东岳 谢谢李老师的回复，期待:xiexie:

多谢老师！

屌爆了，文章发出来之后我学习学习

谢谢啦

多谢老铁，$\frac{1}{\Delta t}$已更正为$\frac{\rho_\rP^{*}}{\Delta t}$

想复现溃坝一物理模型的数据，需要保持上游水头一直不变，请问各位大神这怎么操作

@麦迪文 你好，请问能参考下您的代码吗:xiexie:

@chengan-wang 怎么修正啊