可视化不连续的Galerkin /有限元数据

我想在ParaView中可视化使用不连续Galerkin（DG）方法获得的仿真结果。与有限体积方法类似，问题域被划分为立方体形的单元（“元素”）。与有限体积方法相反，每个像元中不仅有一个解向量，而且每个像元在多个高斯积分点处都包含解。$\mathbf{u}$$\mathbf{u}$

我的问题是，是否有人有使用ParaView / VTK有效地可视化此类数据的经验，以及您选择了哪种方法来表示VTK中的数据。我想到了几种可能的方法，但我不知道哪种方法最有前途：

（1）使用体
素每个集成点使用一个体素。
Pro：所有与标准VTK非结构化单元类型一起使用的插件将继续工作，而无需进行任何更改。
缺点：由于积分点分布不均匀，因此可能难以找到顶点的正确位置。同样，由于DG框架允许使用不连续的溶液，因此可以在细胞表面上定义两次溶液。同样，层次信息（域划分为元素，每个元素包含几个点）也会丢失。

（2）使用多
顶点每个积分点使用一个顶点。
专业：最容易实现，易于使用不同的解决方案在同一位置指定多个点。
缺点：能力来可视化数据为“细胞”丢失，加上相同的缺点如上。

（3）使用VTK正交方案
使用对正交方案的内置支持。
优点：相当简单的实现，保留了原始解决方案的所有关系和属性。
缺点：由于这是一种全新的单元格类型，因此许多（大多数）现有插件将不再起作用，可能必须重写。

deal.II已长期使用选项（2）并取得了很好的成功。换句话说，在2d中，每个顶点在输出文件中出现4次，这意味着要写入的数据更多，但字段可以不连续。我知道的唯一另一个缺点是，某些可视化程序在进行流线处理时无法遵循不共享顶点的单元格之间的接口，尽管今天它们似乎已经变得更好。

您可能还想在我们的常见问题解答中查看此条目：http : //dealii.sourceforge.net/index.php/Deal.II_Questions_and_Answers#The_graphical_output_files_don.27t_make_sense_to_me_.E2.80.94_they_seem_to_have_too_many_degrees_ofs_of。

更新：2018年5月31日

链接已移至：https : //github.com/dealii/dealii/wiki/Frequently-Asked-Questions#the-graphical-output-files-dont-make-sense-to-me----they-似乎有太多的自由度