Autovalores do Laplaciano - Departamento de Matemática - UFMG
Autovalores do Laplaciano - Departamento de Matemática - UFMG Autovalores do Laplaciano - Departamento de Matemática - UFMG
Capítulo 6 Método de Elementos Finitos O método de elementos finitos é um outro método de discretização de equações diferenciais parciais baseado na reformulação variacional da equação. Por exemplo, como já vimos exemplos, encontrar a solução u de uma equação diferencial parcial dada é equivalente a resolver um problema de minimização F (u) = min F (v) v∈V onde V é um conjunto de funções admissíveis e F : V −→ R é um funcional. Em geral, a dimensão de V é infinita e portanto as funções em V não podem ser descritas por um número finito de parâmetros. Discretizar este problema através de elementos finitos é substituir o espaço de dimensão infinita V por um subespaço de dimensão finita Vh consistindo de funções simples (por exemplo, funções polinomiais). O problema discreto passa a ser encontrar o minimizador do funcional F sobre o subespaço Vh. Espera-se que este seja uma aproximação do minimizador de F sobre o espaço completo V , isto é, uma aproximação para a solução da equação diferencial parcial. 6.1 O Caso Unidimensional Nesta seção, desenvolveremos métodos de elementos finitos para resolver o problema de Dirichlet para a equação de Poisson em uma dimensão ′′ −u = f (x) em [0, 1] , (6.1) u (0) = u (1) = 0, onde f é uma função contínua. 6.1.1 Formulação Variacional Para obter uma formulação variacional deste problema, defina e V = v ∈ C 0 ([0, 1]) : v ′ é contínua por partes em [0, 1] e v (0) = v (1) = 0 F (v) = 1 2 1 0 (6.2) |v ′ (x)| 2 1 dx − f (x) v (x) dx = 0 1 2 v′ L2 − 〈f, v〉 L2 . (6.3) Veremos agora que uma solução para o problema de Dirichlet (6.1), que sabemos existir por integração simples, é solução tanto de um problema de minimização como de um problema variacional. 121
Rodney Josué Biezuner 122 6.1 Proposição. (Problema Variacional) Se u ∈ V é uma solução do problema (6.1), então u é a solução única do problema variacional Prova. Multiplicando a equação 〈u ′ , v ′ 〉 L 2 = 〈f, v〉 L 2 para todo v ∈ V. (6.4) −u ′′ (x) = f (x) por uma função teste v ∈ V e integrando sobre o intervalo (0, 1), obtemos Integrando por partes, temos Portanto, 1 0 − 1 0 u ′′ (x) v (x) dx = u ′′ (x) v (x) dx = u (x) v (x)| 1 0 − 1 1 u 0 ′ (x) v ′ 1 (x) dx = 0 0 f (x) v (x) dx. 1 u 0 ′ (x) v ′ 1 (x) dx = − u 0 ′ (x) v ′ (x) dx. f (x) v (x) dx. A unicidade de solução para o problema variacional (6.4) é facilmente determinada. Se u1, u2 satisfazem para todo v ∈ V , então 〈u ′ 1, v ′ 〉 L 2 = 〈f, v〉 L 2 , 〈u ′ 2, v ′ 〉 L 2 = 〈f, v〉 L 2 , 〈u ′ 1 − u ′ 2, v ′ 〉 L 2 = 0 para todo v ∈ V , em particular para v = u1 − u2, donde u ′ 1 − u ′ 2 L 2 = 0. Isso implica u1 − u2 = c para alguma constante c, e as condições de fronteira implicam que c = 0. 6.2 Proposição. (Problema de Minimização) u ∈ V é uma solução do problema variacional (6.4), se e somente se u satisfaz F (u) = min F (v) . (6.5) v∈V Prova. Suponha que u satisfaz (6.4). Dado v ∈ V , escreva w = u − v. Temos F (v) = F (u + w) = 1 2 u′ + w ′ L2 − 〈f, u + w〉 L2 = 1 2 u′ L2 + 〈u ′ , w ′ 〉 + 1 2 w′ L2 − 〈f, u〉 L2 − 〈f, w〉 L2 = 1 2 u′ L2 − 〈f, u〉 L2 + 〈u ′ , w ′ 〉 − 〈f, w〉 L2 + 1 2 w′ L2 = F (u) + 1 2 w′ L2 F (u) . Reciprocamente, suponha que u é um minimizador para o funcional F em V . Considere a função quadrática g : R −→ R definida por g (t) = F (u + tv) . Temos g (t) = 1 2 u′ L 2 + t 〈u ′ , v ′ 〉 + t2 2 v′ L 2 − 〈f, u〉 L 2 − t 〈f, v〉 L 2 = t2 2 v′ L 2 + t [ 〈u ′ , v ′ 〉 − 〈f, v〉 L 2] + F (u) .
- Page 71 and 72: Rodney Josué Biezuner 70 A norma l
- Page 73 and 74: Rodney Josué Biezuner 72 De fato,
- Page 75 and 76: Rodney Josué Biezuner 74 Além dis
- Page 77 and 78: Rodney Josué Biezuner 76 de uma ma
- Page 79 and 80: Rodney Josué Biezuner 78 donde n |
- Page 81 and 82: Rodney Josué Biezuner 80 3.10 Teor
- Page 83 and 84: Rodney Josué Biezuner 82 Da defini
- Page 85 and 86: Rodney Josué Biezuner 84 3.6.1 Esq
- Page 87 and 88: Capítulo 4 Métodos Iterativos par
- Page 89 and 90: Rodney Josué Biezuner 88 4.1.2 Mé
- Page 91 and 92: Rodney Josué Biezuner 90 ou seja,
- Page 93 and 94: Rodney Josué Biezuner 92 4.2.1 Con
- Page 95 and 96: Rodney Josué Biezuner 94 para t su
- Page 97 and 98: Rodney Josué Biezuner 96 4.7 Corol
- Page 99 and 100: Rodney Josué Biezuner 98 Por outro
- Page 101 and 102: Rodney Josué Biezuner 100 Suponha
- Page 103 and 104: Rodney Josué Biezuner 102 ou Se λ
- Page 105 and 106: Rodney Josué Biezuner 104 Para o c
- Page 107 and 108: Rodney Josué Biezuner 106 Reciproc
- Page 109 and 110: Rodney Josué Biezuner 108 Temos ω
- Page 111 and 112: Rodney Josué Biezuner 110 se ∆x
- Page 113 and 114: Rodney Josué Biezuner 112 A escolh
- Page 115 and 116: Rodney Josué Biezuner 114 a conver
- Page 117 and 118: Rodney Josué Biezuner 116 ou e k+1
- Page 119 and 120: Rodney Josué Biezuner 118 4.31 Cor
- Page 121: Capítulo 5 Métodos Multigrid 5.1
- Page 125 and 126: Rodney Josué Biezuner 124 é chama
- Page 127 and 128: Rodney Josué Biezuner 126 como vé
- Page 129 and 130: Rodney Josué Biezuner 128 No caso
- Page 131 and 132: Rodney Josué Biezuner 130 Prova. C
- Page 133 and 134: Rodney Josué Biezuner 132 onde é
- Page 135 and 136: Rodney Josué Biezuner 134 7.3 Coro
- Page 137 and 138: Rodney Josué Biezuner 136 Como seg
- Page 139 and 140: Rodney Josué Biezuner 138 Prova. O
- Page 141 and 142: Capítulo 8 Métodos Numéricos par
- Page 143 and 144: Rodney Josué Biezuner 142 o maior
- Page 145 and 146: Rodney Josué Biezuner 144 com q1
- Page 147 and 148: Rodney Josué Biezuner 146 e depois
- Page 149 and 150: Rodney Josué Biezuner 148 8.3.2 Im
- Page 151 and 152: Rodney Josué Biezuner 150 Pode-se
- Page 153 and 154: Rodney Josué Biezuner 152 Observe
- Page 155 and 156: Referências Bibliográficas [Asmar
- Page 157: Rodney Josué Biezuner 156 [Thomas2
Rodney Josué Biezuner 122<br />
6.1 Proposição. (Problema Variacional) Se u ∈ V é uma solução <strong>do</strong> problema (6.1), então u é a solução<br />
única <strong>do</strong> problema variacional<br />
Prova. Multiplican<strong>do</strong> a equação<br />
〈u ′ , v ′ 〉 L 2 = 〈f, v〉 L 2 para to<strong>do</strong> v ∈ V. (6.4)<br />
−u ′′ (x) = f (x)<br />
por uma função teste v ∈ V e integran<strong>do</strong> sobre o intervalo (0, 1), obtemos<br />
Integran<strong>do</strong> por partes, temos<br />
Portanto,<br />
1<br />
0<br />
−<br />
1<br />
0<br />
u ′′ (x) v (x) dx =<br />
u ′′ (x) v (x) dx = u (x) v (x)| 1<br />
0 −<br />
1<br />
1<br />
u<br />
0<br />
′ (x) v ′ 1<br />
(x) dx =<br />
0<br />
0<br />
f (x) v (x) dx.<br />
1<br />
u<br />
0<br />
′ (x) v ′ 1<br />
(x) dx = − u<br />
0<br />
′ (x) v ′ (x) dx.<br />
f (x) v (x) dx.<br />
A unicida<strong>de</strong> <strong>de</strong> solução para o problema variacional (6.4) é facilmente <strong>de</strong>terminada. Se u1, u2 satisfazem<br />
para to<strong>do</strong> v ∈ V , então<br />
〈u ′ 1, v ′ 〉 L 2 = 〈f, v〉 L 2 ,<br />
〈u ′ 2, v ′ 〉 L 2 = 〈f, v〉 L 2 ,<br />
〈u ′ 1 − u ′ 2, v ′ 〉 L 2 = 0<br />
para to<strong>do</strong> v ∈ V , em particular para v = u1 − u2, <strong>do</strong>n<strong>de</strong><br />
u ′ 1 − u ′ 2 L 2 = 0.<br />
Isso implica u1 − u2 = c para alguma constante c, e as condições <strong>de</strong> fronteira implicam que c = 0. <br />
6.2 Proposição. (Problema <strong>de</strong> Minimização) u ∈ V é uma solução <strong>do</strong> problema variacional (6.4), se e<br />
somente se u satisfaz<br />
F (u) = min F (v) . (6.5)<br />
v∈V<br />
Prova. Suponha que u satisfaz (6.4). Da<strong>do</strong> v ∈ V , escreva w = u − v. Temos<br />
F (v) = F (u + w) = 1<br />
2 u′ + w ′ L2 − 〈f, u + w〉 L2 = 1<br />
2 u′ L2 + 〈u ′ , w ′ 〉 + 1<br />
2 w′ L2 − 〈f, u〉 L2 − 〈f, w〉 L2 = 1<br />
2 u′ L2 − 〈f, u〉 L2 + 〈u ′ , w ′ 〉 − 〈f, w〉 L2 + 1<br />
2 w′ L2 = F (u) + 1<br />
2 w′ L2 F (u) .<br />
Reciprocamente, suponha que u é um minimiza<strong>do</strong>r para o funcional F em V . Consi<strong>de</strong>re a função quadrática<br />
g : R −→ R <strong>de</strong>finida por<br />
g (t) = F (u + tv) .<br />
Temos<br />
g (t) = 1<br />
2 u′ L 2 + t 〈u ′ , v ′ 〉 + t2<br />
2 v′ L 2 − 〈f, u〉 L 2 − t 〈f, v〉 L 2<br />
= t2<br />
2 v′ L 2 + t [ 〈u ′ , v ′ 〉 − 〈f, v〉 L 2] + F (u) .