Visualização de dados descontínuos de Galerkin / elementos finitos

15

Gostaria de visualizar os resultados da simulação, obtidos usando a abordagem descontínua de Galerkin (DG), no ParaView. Similarmente aos métodos de volume finito, o domínio do problema é dividido em células em forma de cubo ("elementos"). Ao contrário dos métodos de volume finito, dentro de cada célula não existe apenas um valor para o vetor de solução , mas cada célula contém a solução u em vários pontos de integração do Gauss.vocêvocê

Minha pergunta é se alguém tem experiência em visualizar esses dados com eficiência com o ParaView / VTK e qual abordagem você escolheu para representar os dados no VTK. Várias maneiras possíveis me vêm à mente, mas não sei qual é a mais promissora:

(1) Use voxels
Use um voxel para cada ponto de integração.
Pro: Todos os plugins que funcionam com os tipos de células não estruturados VTK padrão continuarão funcionando sem alterar nada.
Contras: Como os pontos de integração não são distribuídos uniformemente, pode ser difícil encontrar o local correto dos vértices. Além disso, a solução pode ser definida duas vezes na superfície da célula, pois a estrutura DG permite soluções descontínuas. Além disso, as informações hierárquicas (domínio dividido em elementos, cada elemento contém vários pontos) são perdidas.

(2) Use polvertices
Use um vértice por ponto de integração.
Pro: Mais fácil de implementar, fácil de especificar vários pontos no mesmo local com soluções diferentes.
Con: A capacidade de visualizar dados como "células" é perdida, além das mesmas desvantagens acima.

(3) Usar esquema de quadratura VTK
Use o suporte interno para esquemas de quadratura.
Pro: implementação bastante direta, preserva todas as relações e propriedades da solução original.
Contras: Como esse é um tipo de célula completamente novo, muitos (a maioria) dos plugins existentes não funcionarão mais e provavelmente terão que ser reescritos.

Michael Schlottke-Lakemper
fonte
Alguém com direitos de edição pode adicionar as (novas) tags paraview, vtk, descontínua-galerkin? Obrigado.
Michael Schlottke-Lakemper
Ok, adicionei as tags.
Aron Ahmadia

Respostas:

6

deal.II usa a opção (2) há muito tempo, com muito bom sucesso. Em outras palavras, em 2d, todo vértice aparece 4 vezes no arquivo de saída, o que significa que há mais dados a serem gravados, mas os campos podem ser descontínuos. A única outra desvantagem que conheço é que alguns programas de visualização não podem seguir interfaces entre células que não compartilham vértices ao executar linhas de fluxo, embora hoje elas pareçam ter se tornado melhores.

Você também pode consultar esta entrada em nossas Perguntas frequentes: 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_reo_manydegrees.24

Atualização: 31 de maio de 2018

O link foi movido para: https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#the-graphical-output-files-dont-make-sense-to-me----they- parece ter muitos graus de liberdade

Wolfgang Bangerth
fonte
Obrigado, Wolfgang, este é exatamente o tipo de resposta que eu estava procurando.
Michael Schlottke-Lakemper
Novo endereço da entrada da FAQ: code.google.com/p/dealii/wiki/…
Wolfgang Bangerth