Visualizzazione dei dati discontinui di Galerkin / elementi finiti

Vorrei visualizzare i risultati della simulazione, ottenuti usando l'approccio discontinuo di Galerkin (DG), all'interno di ParaView. Analogamente ai metodi a volume finito, il dominio problematico è diviso in celle a forma di cubo ("elementi"). A differenza dei metodi a volume finito, all'interno di ogni cella non esiste un solo valore per il vettore della soluzione , ma ogni cella contiene la soluzione u in più punti di integrazione di Gauss.$\mathbf{u}$$\mathbf{u}$

La mia domanda è se qualcuno ha esperienza con la visualizzazione di tali dati in modo efficiente con ParaView / VTK e quale approccio hai scelto per rappresentare i dati in VTK. Mi vengono in mente diversi modi, ma non so quale sia il più promettente:

(1) Usa voxel
Usa un voxel per ogni punto di integrazione.
Pro: tutti i plug-in che funzionano con i tipi di celle non strutturati VTK standard continueranno a funzionare senza cambiare nulla.
Contro: poiché i punti di integrazione non sono distribuiti uniformemente, potrebbe essere difficile trovare la posizione corretta dei vertici. Inoltre, la soluzione può essere definita due volte sulla superficie delle cellule, poiché il framework DG consente soluzioni discontinue. Inoltre, le informazioni gerarchiche (dominio diviso in elementi, ogni elemento contiene diversi punti) vengono perse.

(2) Usa poligonali
Usa un vertice per punto di integrazione.
Pro: più facile da implementare, facile da specificare più punti nella stessa posizione con soluzioni diverse.
Contro: la capacità di visualizzare i dati come "celle" viene persa, oltre agli stessi svantaggi di cui sopra.

(3) Usa schema di quadratura VTK
Utilizza il supporto integrato per schemi di quadratura.
Pro: implementazione piuttosto semplice, preserva tutte le relazioni e le proprietà della soluzione originale.
Contro: poiché si tratta di un tipo di cella completamente nuovo, molti (la maggior parte) dei plugin esistenti non funzioneranno più e probabilmente dovranno essere riscritti.

deal.II utilizza l'opzione (2) già da molto tempo con un ottimo successo. In altre parole, in 2d, ogni vertice appare 4 volte nel file di output, il che significa che ci sono più dati da scrivere ma i campi possono essere discontinui. L'unico altro inconveniente di cui sono a conoscenza è che alcuni programmi di visualizzazione non sono in grado di seguire le interfacce tra celle che non condividono vertici quando si eseguono streamline, anche se oggi sembrano essere migliorate.

Puoi anche dare un'occhiata a questa voce nelle nostre FAQ: 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_haf_too_m_fany_

Aggiornamento: 31 maggio 2018

Il link è stato spostato in: https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#the-graphical-output-files-dont-make-sense-to-me----they- sembrano da avere troppo molti gradi di libertà