@东岳 谢谢东岳老师！:ok3:

@东岳 顶部#include已加

我找着我第一次发的方程改得，就是第二个block的代码，算出来的结果和MPPICInterfoam的基本一样。当然有可能我的粒子直径比较小，对计算结果影响不大？

我没明白您下面写的这个是什么意思。如果不看mrf的话， 你这个和interfoam相比就是减上了一项

fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)

这一项就是65293cb6-a527-49eb-bdc4-2355ec82c507-image.png
我不知道怎么推出的这一项？

https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/OpenFOAM/fields/GeometricFields/GeometricField/GeometricBoundaryField.C

这是边界场

@星星星星晴 谢谢！！

@红豆沙 太感谢了 啦 确实是这个命令，，半年没用OF 都给忘了 命令了 这错误犯的，，害 丢人列，，:mianmo:

@cccrrryyy 哦哦，好的。感谢大神的解答！

@pcxu 在 关于interFoam dambreak算例的一点疑问 中说：

defaultFieldValues
(
volScalarFieldValue alpha.water 0
);

看引用的部分
在上面setFieldsDict里，设置的默认场就是 alpha.water 0
所以water对应的就是0

你也可以把水设置成 alpha.water 1
这样水的参数对应的就都是1

@东岳 果然，linux版本的和windows版本的不一样。
linux在skip zero点掉以后，可以显示alpha.water的初始场了。
我之前用的windows版本的，0时刻就显示不出来。
至于paraFoam能不能编译出来，随缘吧

@多相流FOAM 热计算直接给定热流密度或者用对流换热算比较合适

@wangxiaowei hello，我可能理解错了你的问题，我先用imageJ软件从ct扫描切片中提取出土壤介质，生成了stl文件，然后用snappyHexMesh画出孔隙的网格

哦好像是定义没有弄明白，miu和nu是不一样的····

@东岳

我没有找到李老师说的“tutorial: decayIsoTurb”。不过在OF里找到了这个前处理程序7dbae584-0745-4a52-bc5f-ddfbc18a3618-image.png 正在看这篇文章，还看到了一篇中文的“一种改进的均匀各向同性湍流初始化方法”

我看的比较慢，有些基本概念还不太清楚。:136:

@D-Benjamin 具体一共都有哪些函数我也不是很清楚，具体函数和怎么用通常都是看看openfoam自带的边界或者求解器或者其它的应用有类似你想获得功能，就去翻源代码。gsum应该也可以对边界上的网格体积或者面进行求和，应该就能实现你说的功能。

edit -> find data
如果数据只在网格或者点，可以先用filter point data to cell 或者cell data to point筛一下。

感谢分享 :146: :146:

插值只对patch进行，和是否并行无关。并行对local patch进行插值。

@Izumi 这个没添加，我回去试试

@东岳 结果.tar
这个是发散的相分数，速度，压力图

个人理解，如果不改transport property的话，引用的nu应该要加上phase1的关键字，比如可能要输出的应该是nu.air。关键是看nu在你用的solver里注册成了什么名字，这个名字可能和输入的名字完全不一样（注：是注册的名字不是变量名）