如何使用VTK和Python在时间步中生成字段的快照

Question

我正在使用OpenFoam安装（v1812）的cavity示例。执行它（icoFoam）后，得到以下结果：

$ ls
0  0.1  0.2  0.3  0.4  0.5  constant  system

我想知道是否有可能通过脚本对时间步进行后处理，以生成具有给定字段值（例如压力）的像素矩阵。

除了每个时间步长的数组（可以用numpy处理）之外，我还想显示这样的热图：

Answer 1

您可以找到使用cavity的here示例vtkplotter的后处理示例。

Answer 2

我认为您也可以通过函数对象来完成此操作。通过在每个时间步骤写入和提取所需的数据，您可以拥有任意数量的数据。通过使用pyFoam或任何其他python代码，您可以将其用作数组。首先尝试搜索functionObjects。