python进行OpenFOAM流场后处理

[田畔的风]

@尚善若水 没能理解你的第二个问题

[尚善若水]

@田畔的风 比如我画出的了OH云图，但是我同时想在其上面绘制当量混合分数Zmix=constant这条线。

[田畔的风]

@尚善若水 画单纯的等值线可以参考 https://matplotlib.org/3.1.1/gallery/images_contours_and_fields/tricontour_smooth_delaunay.html#sphx-glr-gallery-images-contours-and-fields-tricontour-smooth-delaunay-py

[尚善若水]

@田畔的风 好的，谢谢。

[李东岳]

[fangyuanaza]

@田畔的风 这个流场图真的很漂亮。请教一下老师，我的网格是O型+H型，Fortran代码算的，能导入paraView。想请教一下，这种情况如何导入到python画图呢？如果直接用contourf 或者tricontourf都会导致网格变形。比如在paraview中为：

在python中直接
ax1.tricontourf(x,y, u,linewidth=0.25, levels=levels,cmap='coolwarm',vmin=0, vmax=1),显示出来为：

造成的原因可能是因为您说的Matplotlib只能处理三角形的网格。想请问在不清楚网格具体信息的情况下，能不能从paraview导出成三角网格？尝试过您代码中的Slice，但由于本身是2D，不能切片。请问有什么建议么？

[田畔的风]

我不太清楚你表述的网格扭曲是什么意思，如果只是纵横比不正确的话，可以使用set_aspect进行调整。比如我上面代码段中的ax.set_aspect('equal')，表示x轴和y轴始终保持1:1的比例。

[fangyuanaza]

@田畔的风 感谢您的回复～

不是纵横比的问题，调整纵横比之后：

而是因为计算域边界是曲线，当调用tricontourf时，在一个矩形域内进行了插值，导致边界变形，成了直线。同理，用contourf时，首先需要对mesh进行插值，也会有同样的问题。

这个问题已经解决了，对于计算域的两个block, 分开处理。先reshapre到对应的block网格，然后用contourf来画。可以得到物理量的云图。

[yuhxFoamer]

@田畔的风
感谢老师分享非常实用的代码，本人也是初次使用pvpython处理数据，请问老师一个简单的问题，当我运行您的程序时，出现了以下的错误，我想这个错误应该是非常常见的，但是目前我尝试的方法有：1 重新install pathlib(用pipe 和 conda都试过) 2. 更新python的版本，这两个方法都没有将此错误排除，求问老师有没有其他方法？

Traceback (most recent call last):
  File "./readData.py", line 3, in <module>
    from pathlib import Path
ImportError: No module named pathlib

[田畔的风]

@yuhxfoamer

1. 首先pathlib这个模块在最终的代码中并没有被使用到，所以你可以将这个模块和对应的代码移除，这并不会影响代码运行结果。
2. 然后谈谈你这个报错。readData.py这个文件是用pvpython执行的，这个python解释器是paraview自带的，对应的第三方库也是在安装paraview时配置好的，所以错误来源于你的paraview附带的pvpython解释器（可能是安装问题，可能是系统环境问题，也可能是软件版本问题，请自行排查），和系统python解释器无关。你使用pip和conda只是在安装你系统python解释器的第三方库。
3. 尽管最终版本这个模块被弃用，不过这里再补充一下：pathlib在本代码中的用途是创建空文件夹（pathlib.Path），在上面的讨论中也提到了，因此可以使用os.mkdir替代。

[知行]

@落花风 老师您好,请问下您是怎么安装的 paraview包呢,我在pypi里面没有发现这个包唉.非常感谢您.

[落花风]

@知行 是的，好像并没有这个包。我是在paraview绿色版里找到python的相关包，直接复制过去的

[知行]

非常感谢

[尚善若水]

@田畔的风 尝试在wsl下的Ubuntu中使用代码，可以顺利写出，如下：

但是meshio读取时出现问题。如下：

读取代码

请问是什么原因呢？

[田畔的风]

@尚善若水 是不是有多边形网格？先尝试打开切片三角化吧

[尚善若水]

@田畔的风 是用blockMesh生成的三维的结构化网格，不过射流出口做了类似处理：

设置 Triangulatetheslice=True 可以正常读取了，但是这里的代码应该怎么修改呢？直接注释并取消相应参数不行。

[尚善若水]

@田畔的风 再次尝试之后（将三角化代码全部注释掉），可以输出了，但显示还是不正确

triang = tri.Triangulation(mesh.points[:, 0], mesh.points[:, 1])
tpc = ax.tricontourf(triang, postData, cmap='jet', levels=16)

正确结果

[田畔的风]

@尚善若水 很显然你并没能成功构建matplotlib所定义的三角网格单元，可以看看我在13楼中针对cell变量的定义。直接用三角网格数据绘图的话，此处你直接设定cell = mesh.cells[0].data应该就行（我现在没电脑，无法验证）。

不正确设定cell和node的对应关系的话，tri.Triangulation会根据你输入的点集进行自动Delaunay三角化，造成你所描述的绘制错误。

[尚善若水]

@田畔的风 是的，感谢大佬，果然可以了。得抓紧学python了，批量后处理急需。