Autovalores do Laplaciano - Departamento de Matemática - UFMG

Autovalores do Laplaciano - Departamento de Matemática - UFMG Autovalores do Laplaciano - Departamento de Matemática - UFMG

09.05.2013 Views

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) .

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) .

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!