@bestucan 好的，我试一下，谢谢前辈

@bestucan 好的谢谢。感觉你说的挺有道理的。之前尾部网格是很扁的矩形（因为存在边界层），结果微小的幅度旋转都会使网格拉扯出现负体积，所以我才逐步减小尾部网格的横纵比，让网格看起来不会那么扁，成功实现了8°范围的俯仰，然后为了做更大的俯仰角度，继续减小网格横纵比，试图可以让他多转一些角度，当网格的横纵比减小到现在的情况时，并没有出现网格畸形，也出现了发散，就不知道是为什么了。

@星星星星晴 可以留个联系方式

@cccrrryyy 非常感谢~

@东岳 好的，谢谢李老师。今天早上破案了，OpenFOAM7成功编译了。

@东岳 感谢李老师！！是我边界条件设的不合理加上网格密度取值的问题。已经更正。谢谢李老师耐心解答！！不可压的alpha就是groMetricOneField，我想太简单了，以为稳流模型把空气自动过滤掉了。

@东岳 好的，谢谢老师～～那主要还是从模拟结果上分析最为有效吧～～

@东岳 好的东岳老师，谢谢东岳老师，实在不好意思，上次断网没发出去刚发现。

@东岳 对的，我也正在琢磨在cfd-online和stackoverflow上推一下的，此插件支持双语版本，根据用户的vscode界面语言选择使用中文还是英文显示

@东岳 谢谢东岳老师！: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画出孔隙的网格