谢谢李老师

@李东岳 好的哈，东岳老师，很期待，能不能请求一下也出个petsc安装教程，感觉装petsc坑更多一些

@李东岳 至少论文里面是这样说的。The strain rate tensor S is treated implicitly in the modified flow solver.

@小考拉 这个不错。之前用230版本的gamg一直报错，有了 cacheAgglomeration这些关键词，终于跑通了

其实也不难，你用grep大概能猜出个所以然来。
grep -rn "validate"
其他都是空函数只有一个函数有点意思

template<class BasicTurbulenceModel> void Foam::ReynoldsStress<BasicTurbulencemModel>::validate() { correctNut(); }

然后我又grep 了correctNut()
常用的是kOmega 我看了下

template<class BasicTurbulenceModel> void kOmega<BasicTurbulenceModel>::correctNut() { this->nut_ = k_/omega_; correctBC; BasicTurbulenceModel::correctNut(): }

我虽然没有一行行debug， 大概猜出来就是湍流初始化。最开始给定k omega 初始值计算出第一步的nut 初始值给k omega方程，即previous time step，欢迎指正。

遇到过，忘记怎么解决的了，换个paraview低版本？

@史浩 按照老师说的用alpha识别相界面，不管我设置的那个横截面是否把网格切分，比如我设置z轴0.0002m一下为液体区域，当z轴方向50个网格时，初始计算他的液相含量是0.2；当我z轴方向51个网格时，初始计算他的液相含量变成了0.196078，导致0文件中alpha非0即1，因此并没有识别出相界面。

@chengan-wang 对

2004里面g9已经有了，如果装2012，换回g9就可以，任意切换使用

那可能是你的组分里根本就没有定义N2，你把inertSpecie N2换成inertSpecie H2O试试。

@bestucan 我去试试，感谢大佬::140:

@李东岳 东岳老师，我在使用multiphaseEulerFoam（OF8）的时候出现了这种情况，有没有可能是局部固相分数太小造成的这种发散？

謝謝老哥，後來用的歪門邪道搞定的，你説的這個方法不錯 @bestucan

@李东岳 在 关于cht的网格划分问题 中说：

writeZones

万分感谢李老师的解答，转化格式（写入命令语句）为
fluentMeshToFoam xxx.msh -writeZones，
本人OpenFoam小白，望老师多多包涵:chouchou:

可以啊 这个跟几何没关系啊

@bestucan 万分感谢！

@李东岳 对的，大概有十几公里，那谢谢东岳老师，我再使劲努努力钻研一下。个人觉得根据文献植入算法太难了，虽然能看懂离散和求解过程，但是想在openfoam里面实现觉得难度挺大的，当时老师您带我们植入的单相可压缩的，目前这是两相可压缩，我努力努力，再次感谢东岳老师的持续关心:xinxin:

问题已解决，勾选paraView中的“Use VYkPolyhedron”选项后交错就会消失

@bestucan 我看了目前用到这个 -I. 的地方，包括reactingParcelFoam跟sprayFoam，都有子文件夹，比如sprayFoam下的egineFoam，应该跟compile子求解器有关

@李东岳 东岳老师是想诱惑我再报一次您的辅导班:141:

@cresendo 我也是看了亿遍~~~~~