Memvisualisasikan data elemen Galerkin / elemen terputus

Saya ingin memvisualisasikan hasil simulasi, diperoleh dengan menggunakan pendekatan Disceruous Galerkin (DG), dalam ParaView. Demikian pula dengan metode volume hingga, domain masalah dibagi menjadi sel berbentuk kubus ("elemen"). Berbeda dengan metode volume hingga, dalam setiap sel tidak hanya ada satu nilai untuk vektor solusi , tetapi setiap sel berisi solusi u pada beberapa titik integrasi Gauss.$\mathbf{u}$$\mathbf{u}$

Pertanyaan saya adalah apakah ada orang yang memiliki pengalaman dalam memvisualisasikan data tersebut secara efisien dengan ParaView / VTK, dan pendekatan apa yang Anda pilih untuk mewakili data dalam VTK. Beberapa cara yang mungkin muncul di benak saya, tetapi saya tidak tahu mana yang paling menjanjikan:

(1) Gunakan voxels
Gunakan satu voxel untuk setiap titik integrasi.
Pro: Semua plugin yang bekerja dengan tipe sel VTK standar tidak berstruktur akan terus bekerja tanpa mengubah apa pun.
Con: Karena titik integrasi tidak terdistribusi secara merata, mungkin sulit untuk menemukan lokasi simpul yang benar. Juga, solusinya dapat didefinisikan dua kali pada permukaan sel, karena kerangka kerja DG memungkinkan solusi diskontinyu. Juga, informasi hierarkis (domain dibagi menjadi beberapa elemen, setiap elemen mengandung beberapa poin) hilang.

(2) Gunakan polyvertices
Gunakan satu titik per titik integrasi.
Pro: Mudah diimplementasikan, mudah untuk menentukan beberapa titik di lokasi yang sama dengan solusi yang berbeda.
Con: Kemampuan untuk memvisualisasikan data sebagai "sel" hilang, ditambah kerugian yang sama seperti di atas.

(3) Gunakan skema quadrature VTK
Gunakan dukungan bawaan untuk skema quadrature.
Pro: Penerapan yang lebih mudah, menjaga semua hubungan dan properti dari solusi asli.
Con: Karena ini adalah tipe sel yang sama sekali baru, banyak (sebagian besar) plugin yang ada tidak akan berfungsi lagi dan mungkin harus ditulis ulang.

deal.II menggunakan opsi (2) untuk waktu yang lama dengan kesuksesan yang sangat baik. Dengan kata lain, dalam 2d, setiap dhuwur muncul 4 kali dalam file output yang berarti ada lebih banyak data yang akan ditulis tetapi bidang dapat tidak kontinu. Satu-satunya kelemahan lain yang saya sadari adalah bahwa beberapa program visualisasi tidak dapat mengikuti antarmuka antara sel yang tidak berbagi titik ketika melakukan streamline, meskipun hari ini mereka tampaknya telah menjadi lebih baik.

Anda mungkin juga ingin melihat entri ini di FAQ kami: http://dealii.sourceforge.net/index.php/Deal.II_Questions_and_Answers#The_graphical_output_files_don.27t_make_sense_to_me_.E_80_4_they_se__dem_dem_dem_dem

Pembaruan: 31 Mei 2018

Tautan telah dipindahkan ke: https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#the-graphical-output-files-dont-make-sense-to-me----they- tampaknya memiliki terlalu banyak derajat kebebasan