CFD中文网

@石油cfder 理论上越小越好，但也不一定，可以参考文献的设置或者文献中的结果来选择

diagonal是对角矩阵求解器。

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

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

@allanzhong 看这个意思，应该和你说的差不多，理论上肯定为，但实际上从理论转换为数值算法的时候总会有些偏差，算法在具体实现的时候又会有一些偏差。如果是这样，这些偏差我认为是可以接受的。

@李东岳 好的老师:xiezuoye:

@bestucan 实在是太感谢老师了，刚刚去试了下，把单核变多核就能增加新相了。应该是因为服务器是用的多核计算，所以改参数也需要换成多核。

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

我这个是一个圆柱，分成1/26，成为周期对称模型。中间是空气，不能省略。

@xpqiu 谢谢！

这个文章没有讨论【一维压力波求解】

很多湍流模型的假定就是各向同性，这样模化起来、封闭起来比较简单。各向异性那种，就得上雷诺应力模型了，明显要复杂

@搬运工不好当 你先wclean一下试试:）

@李东岳 问题已经解决了，是我在录制脚本的时候进行了多余的不必要操作，再次感谢李东岳老师，您的建议大大的减少了我的工作量。

@东岳 :huahua: :huahua: :huahua: :huahua: :huahua: :huahua:

@闹市闲民 再给推荐一本粘性流体力学教材，个人亲测食用最佳:chitang: 7d7a70c8-f0cf-4676-b460-7866a9178ba8-image.png

@doubletrump 您好，我想向您请教，投影面积指的是？如果是一个圆柱，风垂直吹过来，指的是直径乘以高还是圆的面积？

曲线救国：建一个矢量场函数grad($$Velocity[0])，再去用这玩意儿的第2个分量..

@李东岳 谢谢东岳老师。是的，水气系统，但含有一些其它成分，所以液相粘度会增大一些。我能理解气含率增加，小幅增加比如0.25符合预期，但0.4/0.5有点过高

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

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

现在高版本只有Schnerr-Sauer模型和ZGB模型