Existem duas abordagens gerais para representar soluções no método descontínuo de galerkin: nodal e modal. Modal : As soluções são representadas por somas de coeficientes modais multiplicados por um conjunto de polinômios, por exemplo, onde geralmente é polinômio ortogonal , por exemplo,...
Gostaria de saber se existe uma maneira rápida de calcular a distância euclidiana de dois vetores na oitava. Parece que não há função especial para isso, então devo usar apenas a fórmula com
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...
Estou aprendendo a teoria por trás dos métodos DG-FEM usando o livro Hesthaven / Warburton e estou um pouco confuso sobre o papel do 'fluxo numérico'. Peço desculpas se essa é uma pergunta básica, mas procurei e não encontrei uma resposta satisfatória. Considere a equação da onda escalar linear:...
No artigo http://www.sciencedirect.com/science/article/pii/S0045782509003521 , uma equação local elemento-HDG é descrita na página 584 equação (4), com uma das equações da seguinte forma: - ( uh, ∇ q)K= - ⟨ u^h⋅ n , q- q¯⟩∂K−(uh,∇q)K=−⟨u^h⋅n,q−q¯⟩∂K-(u_h,\nabla q)_K = -\left\langle\hat{u}_h \cdot...
Sempre ouvi dizer que a fácil paralelização era uma das vantagens dos métodos DG, mas não vejo realmente por que nenhuma dessas razões também se aplica ao Galerkin
I ter implementado um esquema ADER-Galerkin Descontínuo para a resolução de sistemas lineares de leis de conservação do tipo de e observou-se que a condição de CFL é muito restritiva. Na bibliografia, um limite superior para o intervalo de tempo Δ t ≤ h∂tvocê+ A ∂xvocê+ B ∂yvocê=...
Em breve implementarei o método descontínuo de Galerkin nodal, e depois disso, antes de conhecer as matrizes básicas de indexação que precisarei computar, dados dados de malha e polinomiais. O problema que encontrei no código anterior foram erros sutis que cometi ao computar coisas como indexação...