@wwzhao 在OpenFOAM-2.1的时候，程序里写的就是“CrankNicholson” ，所以controlDict里也应该写“CrankNicholson”，虽然这个拼写确实是不正确的。但是从OpenFOAM-2.3开始，这个名字改成了正确的拼写：“CrankNicolson”，相差一个字母。所以这就导致了2.1和2.3的算例会在这里不兼容。

@搬运工不好当

在对拉普拉斯项离散的时候，如果使用非正交网格会引入误差，因此需要非正交修正。所以nNonOrthCorr的作用就是消除非正交网格引入的误差。

你的arc定义错了，你觉得一个点能定义一条arc吗？

你得把错误内容贴出来，不然一句话看不出来。出错的原因太多了，而且不知道你修改了求解器，甚至case的哪部分...............

transportpProperties没有声明。

这几天放假没在学校，才看到回复。谢谢分享 。我研究一下。万分感谢

是的，可以忽略。

并没有用过pyfoam

@xpqiu 谢谢！

在liggghts/src下面找到这个fix_heat_gran_conduction.cpp文件，里面讲了颗粒和壁面之间的静态导热的计算，看看代码应该可以知道的。

用surfaceCheck检查STL文件的结果呢？

diagonal是对角矩阵求解器。

fvScalarMatrix rhoEqn ( fvm::ddt(rho) + fvc::div(phi) == fvOptions(rho) );

从这个方程里面可以看出，只有时间项是隐式离散，fvc和fvOptions都进入了矩阵右边。因此离散之后是一个对角阵。直接除过去就可以。并不需要迭代。因此迭代数为0.

这个原项是怎么弄得？fvOptions，之前我记得都是MRFSimpleFoam，SRFSimpleFoam，现在都整合到fvOptions里面了。原理是什么啊？

@mark 这种“提示文件找不到”的都是路径搞错了，我看你已经编译出来了，仔细查看查看路径。

把fvsolutuon里的tolerance由e-6改为e-8对结果影响大不。

像大家请教一个问题，算低速流体（10m/s以下）且求解域很小（1cm1cm1cm量级）当用可压求解器求解好像CFL不管设多小，算足够长时间都会发散，除开用不可压求解器以外还有什么其他解决方式么？

@李东岳 想明白了，多谢东岳大神！

是中国人，还没改国籍。

互联网的某个角落 or 提问。

@搬运工不好当 每个边界都要赋值，格式大概是
nonuniform List<scalar>
xxx
(
aaa
bbb
ccc
.....
);
xxx 要等于边界面元的数目，括号里的数字的个数要等于xxx。