@Zhy2022

有没有什么进展或者文章见刊啊？

试试reconstructParMesh之后再reconstructPar？
如果直接处理的话可以运行 touch case.foam && paraview case.foam或者paraFoam -builtin，但是paraview没法直接处理拉格朗日场，可以试试用tecplot

@leilei operation pressure是0的话，只要有压力下降就会有负数

我又回来了，这次发现虽然写了用坐标去计算不同点的阻力系数，但是实际计算中发现还是均一的，以为是定义的变量被重新赋值了，重写了所有的中间变量发现也没办法解决，目前认为是坐标计算的不对或者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]);
        }

收到了+1，谢谢李老师

@TINGHAIK 下载了还没用过，感觉算气动应该挺不错的了，看里边的案例做的还挺有意思的，不过界面上应该还可以下点功夫，不涉及工程计算用现有的工具可以了。

@Samuel-Tu 最近几天看了看salome，感觉照着官方的tut操作一遍，使用上应该没啥问题了，除了操作有点繁琐，其他感觉用起来还挺顺手的，我先学学，学会了可能在b站录个系列教程。
salome的第一个tut就可以创建只有一层的网格，楼主有需要可以自己去看看。

@xpqiu u 好的，感谢！

@Samuel-Tu 好滴，gmsh和salome哪个好一些呢？