Taxas de convergência de elementos finitos para problemas mistos


8

Eu codifiquei um problema do Stokes Flow usando elementos finitos e estou no processo de verificar se ele funciona. Só não tenho certeza de qual taxa de convergência devo esperar enquanto refino a malha globalmente.

Conheço problemas escalares usando funções de base linear que eu esperaria convergência de ordem ( é o tamanho do elemento), e usando funções de base quadrática eu esperaria convergência de ordem na norma e uma potência a menos em o seminário . O problema que estou tendo agora é que, ao codificar o fluxo de Stokes, usei o elemento Taylor-Hood, que usa lineares para a pressão e quadrático para os componentes de velocidade. É tão simples quanto as velocidades convergentes em e a pressão na ordem ?h2hh3L2H1h3h2

Publiquei isso no mathoverflow primeiro e me disseram que poderia ser mais adequado para este fórum.

Respostas:


7

A aproximação de Taylor-Hood do fluxo de Stokes é um método misto de elementos finitos, para o qual as estimativas de erro geralmente têm a forma que é o exato solução e é a aproximação. Para o fluxo de Stokes, e (com valor médio zero). Para o elemento - -Taylor-Hood, consiste em polinômios quadráticos por partes contínuos e

(1)uuhV+pphMC(infwhVhuwhV+infqhMhpqhM),
(u,p)V×M(uh,ph)Vh×MhV=H01(Ω)dM=L2(Ω)P2P1VhMhde polinômios lineares contínuos por partes, para os quais os dois termos do lado direito podem ser delimitados por erros de aproximação quadrática usando argumentos padrão (por exemplo, regras de transformação e lema de Bramble-Hilbert): (regra geral "número de derivadas à esquerda" potências de número de derivadas a direita"). Inserir isso em gera (supondo que a solução exata seja de fato bastante regular).
infwhVhuwhH1Ch2uH3infqhMhpqhL2Ch2pH2
+h =(1)
uuhH1+pphL2Ch2(uH3+pH2).

Como o problema de Stokes contínuo e o discreto são um sistema linear acoplado triangular superior da forma a estimativa de erro (que explora que a diferença das soluções satisfaz um sistema semelhante para e e usa a de no kernel de ) para depende em geral do erro em .

Au+Bp=fBu=0
AhBhABuuhpph

No entanto, se você olhar para a prova, há um furo de loop: se o espaço nulo de estiver contido no espaço nulo de , o termo do acoplamento será eliminado e você obterá uma estimativa de erro envolvendo apenas : A partir daí, você pode aplicar o truque de Aubin-Nitsche padrão (se a equação adjacente estiver bem colocada, que é o caso se o domínio for regular o suficiente - um polígono convexo em 2D ou tiver um limite que possa ser parametrizado por uma função diferenciável de Lipschitz) para obter uma taxa de convergência para o erro de uma ordem maior: BhBu

uuhH1Ch2uH3
ΩL2
uuhL2Ch3uH3

Você pode encontrar esses resultados em Ern, Guermond: Teoria e Prática de Elementos Finitos , Springer, 2004 . (As estimativas de erro são coletadas no Teorema 4.26, enquanto a regularidade necessária para está definida no Lema 4.17; infelizmente, as provas estão espalhadas pelo livro e acho que não é verificado em nenhum lugar explicitamente.)ΩkerBhkerB


Parece que você trabalhou muito na compilação desta resposta. Eu realmente aprecio isso. Você pode esclarecer onde obteve a estimativa: , ou isso é apenas dado? uuhH1+pphL2Ch2(uH3+pH2)
Lukas Bystricky

1
Não tem problema, feliz se for de alguma ajuda. A estimativa vem da inserção dos erros de interpolação padrão na primeira estimativa; Eu adicionei essa etapa.
Christian Clason

Experimentalmente (usando um domínio de caixa que não pode ser parametrizado por uma função diferenciável), estou obtendo convergência em e convergência em para todo o sistema. h2L2hH1
Lukas Bystricky 16/07/2014

1
Hmm, então eu não tenho muita ideia no momento. Essa é realmente mais uma pergunta nova, que as pessoas mais familiarizadas com os problemas de Stokes estão melhor equipadas para responder. Por que você não faz uma nova pergunta (use o link abaixo), onde descreve exatamente qual problema está solucionando (lado direito, condições de contorno, geometria, solução exata), o que espera ver e o que está obtendo ? (Se o código deal.II é pequeno o suficiente, você pode postar isso também, alguns dos principais desenvolvedores são contribuintes freqüentes aqui.)
Christian Clason

1
Acontece que eu estava olhando para os componentes incorretos na minha solução. Quando olho para os componentes corretos, recebo em e em para as velocidades e 1 potência menor para a pressão. h3L2h2H1
Lukas Bystricky
Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.