output - 如何从 FiPy 输出 .pvd 或 xdmf 文件

Question

在示例文档中，有以下部分：

FiPy 不会绘制或输出任何内容，除非您告诉它：

 if __name__ == "__main__":
     viewer = Viewer(vars=(phi,), datamin=0., datamax=1.)

我了解当前配置会导致使用 matplotlib 或 Mayavi 查看器打开查看器。但是，我希望能够导出 .pvd 或 .xdmf 文件以合并模拟。

谢谢你的帮助！

Answer 1

FiPy 目前没有这样的能力，虽然我已经对 XDMF 和多时间步 Gmsh MSH 文件进行了一些实验，但需要找时间来恢复它。您可以使用VTKViewer类保存单个数据快照，但不能保存时间序列。

vw = fp.VTKCellViewer(vars=(phi, psi))
vw.plot(filename="myFile.vtk")

如果有特别想看的，请更新issue #132。