Visualisation des données discontinues de Galerkin / éléments finis


15

Je voudrais visualiser les résultats de simulation, obtenus en utilisant l'approche Galerkin discontinue (DG), dans ParaView. De façon similaire aux méthodes à volume fini, le domaine problématique est divisé en cellules en forme de cube ("éléments"). Contrairement aux méthodes de volumes finis, dans chaque cellule, il n'y a pas qu'une seule valeur pour le vecteur solution , mais chaque cellule contient la solution à plusieurs points d'intégration de Gauss.uu

Ma question est de savoir si quelqu'un a de l'expérience dans la visualisation efficace de ces données avec ParaView / VTK et quelle approche vous avez choisie pour représenter les données dans VTK. Plusieurs pistes me viennent à l'esprit, mais je ne sais pas laquelle est la plus prometteuse:

(1) Utiliser des voxels
Utilisez un voxel pour chaque point d'intégration.
Pro: Tous les plugins qui fonctionnent avec les types de cellules VTK standard non structurés continueront de fonctionner sans rien changer.
Con: Étant donné que les points d'intégration ne sont pas distribués uniformément, il peut être difficile de trouver l'emplacement correct des sommets. De plus, la solution peut être définie deux fois sur les surfaces des cellules, car le cadre DG permet des solutions discontinues. De plus, les informations hiérarchiques (domaine divisé en éléments, chaque élément contient plusieurs points) sont perdues.

(2) Utiliser des polyvertices
Utiliser un sommet par point d'intégration.
Pro: plus facile à mettre en œuvre, facile à spécifier plusieurs points au même endroit avec différentes solutions.
Inconvénients: la capacité de visualiser les données sous forme de "cellules" est perdue, plus les mêmes inconvénients que ci-dessus.

(3) Utiliser le schéma de quadrature VTK
Utilisez le support intégré pour les schémas de quadrature.
Pro: Implémentation plutôt simple, préserve toutes les relations et propriétés de la solution d'origine.
Con: Comme il s'agit d'un type de cellule complètement nouveau, de nombreux (la plupart) des plugins existants ne fonctionneront plus et devront probablement être réécrits.


Quelqu'un avec des droits d'édition peut-il ajouter les (nouvelles) balises paraview, vtk, discontinuous-galerkin? Je vous remercie.
Michael Schlottke-Lakemper

D'accord, a ajouté les balises.
Aron Ahmadia

Réponses:


6

deal.II utilise déjà depuis longtemps l'option (2) avec un très bon succès. En d'autres termes, dans 2d, chaque sommet apparaît 4 fois dans le fichier de sortie, ce qui signifie qu'il y a plus de données à écrire mais que les champs peuvent être discontinus. Le seul autre inconvénient que je connaisse est que certains programmes de visualisation ne peuvent pas suivre les interfaces entre les cellules qui ne partagent pas les sommets lors de la rationalisation, bien qu'aujourd'hui ils semblent être devenus meilleurs.

Vous pouvez également consulter cette entrée dans notre 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_have_too_many_frey

Mise à jour: 31 mai 2018

Le lien a été déplacé vers: https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#the-graphical-output-files-dont-make-sense-to-me----they- semblent avoir trop de degrés de liberté


Merci, Wolfgang, c'est exactement le genre de réponse que je cherchais.
Michael Schlottke-Lakemper

Nouvelle adresse de l'entrée FAQ: code.google.com/p/dealii/wiki/…
Wolfgang Bangerth
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.