由于我在用SPH方法仿真时用的是FORTRAN语言，并且没有找到直接输出vtk文件的代码，因此偷懒通过MATLAB转换一下数据。

用到的Matlab子程序可通过一下链接找到。

Matlab处理数据导出Paraview可读的vtk文件（一）

假设我每个粒子输出的格式按照下面输出的

x(i), y(i), z(i), vx(i), vy(i), vz(i), p(i), rho(i), phase(i)

含义为每个粒子的坐标，速度，压强，密度，相。

1vtkwrite(‘output.vtk’,'unstructured_grid',x(:),y(:),z(:),'vectors','velocity',vx(:),vy(:),vz(:),'scalars', 'pressure',p(:), 
2         'scalars', 'rho',rho(:),'scalars', 'phase',phase(:));

解释：

1，输出文件名

2，定义非结构网格，还有其他选项

3，后面跟每个点的坐标

4，‘vectors’定义向量，后面‘velocity’是自定义的名字，比如速度，然后后面跟速度的三个分量

5，‘scalars’定义数值，后面引号中自定义名字，比如压强和密度等

6，向量和数值最后在paraview中都有对应的显示方式。

将生成的vtk文件拖进paraview就可以在Representation中选择Point Gaussian。

 在Coloring中选择要显示的数据

二、显示向量

 1、选择菜单中filter->Alphabetical->Glyph

或者点下面

/2、设置你要显示的向量