如何利用VTK高效写入固定网格上的时变场数据？

VTK和ParaView完全支持仅存储一次网格、后续仅写入场数据的需求，以下是几种可行的实现方案，按易用性和效率排序：

方案1：MultiBlock容器（.vtm）+ 独立网格/场文件 #

这是最直接且ParaView友好的方案，核心思路是用.vtm（VTK MultiBlock数据集文件）作为每个时间步的容器，静态网格只写入一次，场数据单独存储为小文件。

实现步骤： #

1. 写入静态网格：

   • 用vtkXMLUnstructuredGridWriter将固定网格（含完整几何拓扑，无场数据）写入static_grid.vtu。

   vtkNew<vtkXMLUnstructuredGridWriter> gridWriter;
   gridWriter->SetFileName("static_grid.vtu");
   gridWriter->SetInputData(yourStaticGrid);
   gridWriter->SetDataModeToBinary(); // 可选，压缩减少体积
   gridWriter->Write();
   

2. 写入单步场数据：

   • 每个时间步，创建仅包含当前场数据（点/单元数据）的vtkUnstructuredGrid，确保其点数量、单元数量和类型与静态网格完全一致，写入field_tX.vtu（X为时间步索引）。

   vtkNew<vtkUnstructuredGrid> stepGrid;
   stepGrid->ShallowCopy(yourStaticGrid); // 复用网格结构，不复制数据
   stepGrid->GetPointData()->AddArray(currentPressureArray); // 添加当前场数据
   // 可选：清空不需要的旧数据
   stepGrid->GetPointData()->RemoveArray("OldFieldName");
   
   vtkNew<vtkXMLUnstructuredGridWriter> fieldWriter;
   fieldWriter->SetFileName("field_t0.vtu");
   fieldWriter->SetInputData(stepGrid);
   fieldWriter->SetDataModeToBinary();
   fieldWriter->Write();
   

3. 创建时间步容器（.vtm）：

   • 每个时间步生成一个.vtm文件（如sim_t0.vtm），内容为引用静态网格和当前场数据的MultiBlock结构：

   <VTKFile type="vtkMultiBlockDataSet" version="1.0">
     <vtkMultiBlockDataSet>
       <Block name="StaticGrid">
         <DataSet file="static_grid.vtu"/>
       </Block>
       <Block name="TimeStepField">
         <DataSet file="field_t0.vtu"/>
       </Block>
     </vtkMultiBlockDataSet>
   </VTKFile>
   

4. ParaView可视化：

   • 创建.pvd文件管理时间序列（或直接用ParaView的"Load Time Series"选择所有.vtm文件），加载后ParaView会自动合并两个块的网格与场数据，显示动态效果。

方案2：手动生成仅含场数据的VTU文件 #

若想彻底避免场文件中重复存储网格，可手动构建仅包含场数据的VTU XML结构，无需写入几何拓扑节点。

实现步骤： #

1. 先写入完整静态网格到static_grid.vtu（同方案1步骤1）。
2. 每个时间步，手动生成仅含<PointData>/<CellData>的VTU文件，示例结构：

   <VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian">
     <UnstructuredGrid>
       <Piece NumberOfPoints="10000" NumberOfCells="8000"> <!-- 需与静态网格一致 -->
         <PointData>
           <DataArray type="Float64" Name="Velocity" NumberOfComponents="3" format="binary">
             <!-- 二进制场数据内容 -->
           </DataArray>
         </PointData>
       </Piece>
     </UnstructuredGrid>
   </VTKFile>
   

3. 在ParaView中，先加载static_grid.vtu，再用"Append Attributes"过滤器依次加载所有场数据文件，将过滤器输出设置为时间序列即可。

方案3：结构化网格专属优化（若适用） #

如果你的网格是结构化网格（如规则笛卡尔网格），可使用.vts格式：

• 结构化网格的拓扑由维度（如Dimensions="100 100 100"）定义，无需存储单元连接关系，因此首次写入.vts文件后，后续时间步的.vts仅需写入点数据，文件体积极小。
• 直接用vtkXMLStructuredGridWriter写入，ParaView加载时间序列时会自动复用拓扑信息。

辅助优化：压缩与二进制存储 #

无论哪种方案，都建议启用VTK的二进制存储和压缩，进一步减少文件体积：

writer->SetDataModeToBinary();
writer->SetCompressorTypeToZLib();
writer->SetCompressionLevel(6); // 压缩等级，1-9

内容的提问来源于stack exchange，提问作者Wesley Ranger