@Shihang-Chen 哇！感谢，确实解决了！

有数据的PINN效果还行，但是完全的无数据的不太行

我不太懂这个。确认一下：你说的这个意思，是说PINN需要CFD算出多套参数的结果，然后这些多套参数和结果就是数据。比如先用CFD算100套进口条件，得出100套结果，把这100套进口参数+结果给到PINN，这样的话效果还行？

绝对到位！嘎嘎一嘎子！

另一种方法：仅添加到UEqn当中，压力方程不做处理，斥力的贡献依靠UEqn.H参与到后续的计算当中。结果就没有压力振荡了。

:146: :146: :146: :146:

Interesting！

好的老师，我算完给您回复。

谢谢岳哥，成了，牛牛牛

只要有浮力就会发散，为什么呢

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

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: