我又回来了，这次发现虽然写了用坐标去计算不同点的阻力系数，但是实际计算中发现还是均一的，以为是定义的变量被重新赋值了，重写了所有的中间变量发现也没办法解决，目前认为是坐标计算的不对或者openfoam porosityModel不支持非均一的阻力系数设置
如果想用场自定义阻力需要如何操作呢？直接改求解器加入源项，然后在0文件夹写初始场？

forAll(cells, i) { gamma_[zoneI][i] = Zero; scalar x = mesh_.C()[i].x(); scalar y = mesh_.C()[i].y(); scalar z = mesh_.C()[i].z(); scalar dir = y; if (xxx) { gamma_[zoneI][i].xx() = xxx; gamma_[zoneI][i].yy() = xxx; gamma_[zoneI][i].zz() = xxx; } D_[zoneI][i] = Zero; D_[zoneI][i].xx() += pow(gamma_[zoneI][i].xx(), 1/3) * dXYZ_.value().x(); D_[zoneI][i].yy() += pow(gamma_[zoneI][i].yy(), 1/3) * dXYZ_.value().y(); D_[zoneI][i].zz() += pow(gamma_[zoneI][i].zz(), 1/3) * dXYZ_.value().z(); D_[zoneI][i] = coordSys_.R().transformTensor(D_[zoneI][i]); // leading 0.5 is from 1/2*rho F_[zoneI][i] = Zero; F_[zoneI][i].xx() += 0.5 * fXYZ_.value().x(); F_[zoneI][i].yy() += 0.5 * fXYZ_.value().y(); F_[zoneI][i].zz() += 0.5 * fXYZ_.value().z(); F_[zoneI][i] = coordSys_.R().transformTensor(F_[zoneI][i]); }

@Ryo 前辈好！ 请问Gmesh能够对非结构网格划分block吗?

@李东岳 哈哈，好的，谢谢东岳老师

你是否比较Upwind格式下的残差曲线和linearUpwind 格式的差别呢？这两方法在OpenFOAM中的收敛性是不一样的。或者设置Porbes探针观察下数据采样点最后的波动情况。我觉的是收敛性影响的。有复杂地形或者建筑群的边界层计算时，OpenFOAM高阶阶方案的收敛性并不好，会影响计算精度@李东岳

谢谢李老师

@李东岳 好的哈，东岳老师，很期待，能不能请求一下也出个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: