Візуалізація переривчастих даних Галеркіна / кінцевих елементів


15

Я хотів би візуалізувати результати моделювання, отримані з використанням перерваного підходу Галеркіна (DG), в межах ParaView. Аналогічно методам обмеженого обсягу, проблемна область поділяється на клітиноподібні клітини ("елементи"). На відміну від методів обмеженого обсягу, всередині кожної комірки існує не одне значення для вектора рішення , але кожна комірка містить рішення у кількох точках інтеграції Гаусса.уу

Моє запитання - чи має хто-небудь досвід ефективної візуалізації таких даних за допомогою ParaView / VTK, і який підхід ви обрали для представлення даних у VTK. Мені спадають на думку кілька можливих способів, але я не знаю, який із них є найбільш перспективним:

(1) Використання вокселів
Використовуйте по одному вокселю для кожної точки інтеграції.
Про: Усі плагіни, які працюють зі стандартними неструктурованими типами комірок VTK, продовжуватимуть працювати, нічого не змінюючи.
Con: Оскільки точки інтеграції розподілені не рівномірно, може бути важко знайти правильне розташування вершин. Крім того, розчин може бути визначений двічі на клітинних поверхнях, оскільки рамки DG дозволяють розриви рішень. Також втрачається ієрархічна інформація (домен, розділений на елементи, кожен елемент містить кілька точок).

(2) Використання полівертиків
Використовуйте одну вершину на точку інтеграції.
Pro: найпростіший у здійсненні, легко вказати кілька точок в одному місці за допомогою різних рішень.
Con: Можливість візуалізації даних як "комірок" втрачена, плюс ті ж недоліки, що і вище.

(3) Використання квадратурної схеми VTK
Використовуйте вбудовану підтримку квадратурних схем.
Про: Досить проста реалізація, зберігає всі відносини та властивості оригінального рішення.
Con: Оскільки це абсолютно новий тип комірок, багато (більшість) існуючих плагінів більше не працюватимуть і, ймовірно, доведеться переписати.


Чи міг би хтось із правами на редагування додати (нові) теги paraview, vtk, перерваний-galerkin? Дякую.
Michael Schlottke-Lakemper

Гаразд, додав теги.
Арон Ахмадія

Відповіді:


6

deal.II довго використовує варіант (2) з дуже хорошим успіхом. Іншими словами, у 2d кожна вершина відображається 4 рази у вихідному файлі, що означає, що потрібно записати більше даних, але поля можуть бути розривними. Єдиний інший недолік, про який я знаю, - це те, що деякі програми візуалізації не можуть слідувати інтерфейсам між клітинками, які не поділяють вершини під час потокових потоків, хоча сьогодні вони, здається, стали кращими.

Ви також можете ознайомитись із цим записом у наших поширених запитаннях: 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_mde

Оновлення: 31 травня 2018 року

Посилання було перенесено на: https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#the-graphical-output-files-dont-make-sense-to-me----they- здається, занадто багато градусів свободи


Дякую, Вольфганг, саме таку відповідь я шукав.
Michael Schlottke-Lakemper

Нова адреса запису "FAQ": code.google.com/p/dealii/wiki/…
Вольфганг Бангерт
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.