Método dos Elementos Finitos - Equação do Calor

Equação do Calor (1D)

A equação do calor em 1D descreve a evolução da temperatura $u(x,t)$ ao longo do tempo em uma barra unidimensional, sujeita a condições de contorno e a uma difusão térmica ($\alpha$). A equação diferencial que governa o comportamento do sistema é dada por:

Forma Forte

${\frac {\partial u(x,t)}{\partial t}}=\alpha {\frac {\partial ^{2}u(x,t)}{\partial x^{2}}}$

Condições Iniciais

$u(t=0,x)=sin(\pi x)$

Condições de Contorno

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle u(t,x=0)=0=u(t,x=L)}

Método dos Elementos Finitos (FEM)

O Método dos Elementos Finitos (FEM, do inglês Finite Element Method) é uma técnica amplamente utilizada para resolver equações diferenciais parciais, como a equação do calor. Para aplicar o FEM à equação do calor em 1D, a barra [0,1] é dividida em um número finito de elementos (intervalos) de comprimento Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Delta x} , e as variáveis contínuas são aproximadas por funções de base definidas localmente em cada elemento.

Forma Fraca

Apesar de o problema estar completo na forma forte, que é a mostrada acima, precisamos representá-lo na forma equivalente, chamada forma fraca. A forma fraca é preferível por diversos motivos, sendo um dos principais o fato de que ela não possui derivadas de segunda ordem, possibilitando o uso de funções lineares para fazer a aproximação, que na forma forte teriam essa derivada nula. Para obter a forma fraca a partir da forma forte, deve-se seguir os seguintes passos:

• Multiplica-se por uma função teste Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle v} do mesmo espaço de funções definido nos elementos.

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial u(x,t)}{\partial t}v = \alpha \frac{\partial^2 u(x,t)}{\partial x^2}v}

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \int_0^L \frac{\partial u(x,t)}{\partial t} v(x) \, dx = \alpha \int_0^L \frac{\partial^2 u(x,t)}{\partial x^2} v(x) \, dx}

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \int_0^1 \frac{\partial^2 u}{\partial x^2} v \, dx = \left[ \frac{\partial u}{\partial x} v \right]_0^1 - \int_0^1 \frac{\partial u}{\partial x} \frac{\partial v}{\partial x} \, dx}

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \int_0^1 \frac{\partial u}{\partial t} v \, dx + \int_0^1 \alpha \frac{\partial u}{\partial x} \frac{\partial v}{\partial x} \, dx = \left[ \alpha \frac{\partial u}{\partial x} v \right]_0^1}

O lado direito zera devido às condições de contorno

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \int_0^1 \frac{\partial u}{\partial t} v \, dx + \int_0^1 \alpha \frac{\partial u}{\partial x} \frac{\partial v}{\partial x} \, dx = 0}

Essa é a forma fraca, mas ainda precisamos discretizá-la tanto no espaço quanto no tempo para aplicar o FEM.