Я хотел бы визуализировать результаты моделирования, полученные с использованием прерывистого подхода Галеркина (DG), в ParaView. Как и в методах конечных объемов, проблемная область делится на кубические ячейки («элементы»). В отличие от методов конечного объема, в каждой ячейке имеется не только одно значение для вектора решения , но каждая ячейка содержит решение u в нескольких точках интегрирования Гаусса.
У меня вопрос, есть ли у кого-нибудь опыт эффективной визуализации таких данных с помощью ParaView / VTK, и какой подход вы выбрали для представления данных в VTK. На ум приходит несколько возможных путей, но я не знаю, какой из них наиболее перспективный:
(1) Используйте воксели.
Используйте один воксель для каждой точки интеграции.
Pro: Все плагины, которые работают со стандартными неструктурированными типами ячеек VTK, будут продолжать работать без каких-либо изменений.
Против: Так как точки интегрирования не распределены равномерно, может быть трудно найти правильное расположение вершин. Кроме того, решение может быть определено дважды на поверхностях ячейки, поскольку структура DG допускает прерывистые решения. Также иерархическая информация (домен разделен на элементы, каждый элемент содержит несколько точек) теряется.
(2) Используйте поливертики.
Используйте одну вершину на точку интегрирования.
Pro: Легко реализовать, легко указать несколько точек в одном месте с различными решениями.
Против: Возможность визуализировать данные как «ячейки» теряется, плюс те же недостатки, что и выше.
(3) Использовать квадратурную схему VTK
Используйте встроенную поддержку квадратурных схем.
Pro: Довольно простая реализация, сохраняет все отношения и свойства исходного решения.
Против: Поскольку это абсолютно новый тип ячеек, многие (большинство) существующих плагинов больше не будут работать и, вероятно, их придется переписать.