Me gustaría visualizar los resultados de la simulación, obtenidos usando el enfoque discontinuo de Galerkin (DG), dentro de ParaView. De manera similar a los métodos de volumen finito, el dominio del problema se divide en celdas en forma de cubo ("elementos"). A diferencia de los métodos de volumen finito, dentro de cada celda no hay un solo valor para el vector de solución , sino que cada celda contiene la solución en múltiples puntos de integración de Gauss.
Mi pregunta es si alguien tiene experiencia visualizando dichos datos de manera eficiente con ParaView / VTK, y qué enfoque eligió para representar los datos en VTK. Se me ocurren varias formas posibles, pero no sé cuál es la más prometedora:
(1) Use voxels
Use un voxel para cada punto de integración.
Pro: todos los complementos que funcionan con los tipos de células no estructuradas VTK estándar continuarán funcionando sin cambiar nada.
Con: Ya que los puntos de integración no están distribuidos de manera uniforme, puede ser difícil encontrar la ubicación correcta de los vértices. Además, la solución se puede definir dos veces en las superficies celulares, ya que el marco DG permite soluciones discontinuas. Además, la información jerárquica (dominio dividido en elementos, cada elemento contiene varios puntos) se pierde.
(2) Use polyvertices
Use un vértice por punto de integración.
Pro: más fácil de implementar, fácil de especificar múltiples puntos en la misma ubicación con diferentes soluciones.
Con: La capacidad de visualizar datos como "celdas" se pierde, más las mismas desventajas que las anteriores.
(3) Use el esquema de cuadratura VTK
Use el soporte incorporado para esquemas de cuadratura.
Pro: implementación bastante sencilla, conserva todas las relaciones y propiedades de la solución original.
Con: Dado que este es un tipo de celda completamente nuevo, muchos (la mayoría) de los complementos existentes ya no funcionarán y probablemente tendrán que reescribirse.