Métodos Matemáticos da Fı́sica 1
Equação do Calor - Parte 1
A. Resolva os problemas do calor abaixo. Identifique os tipos de condições de contorno e
interprete fisicamente o que elas significam.


ut = 2uxx , 0 < x < 1
ut = uxx , 0 < x < 2






 u(0, t) = u(1, t) = 0, 0 < t < ∞
u(0, t) = 0, 0 < t < ∞
.
3.
1.
.
1
1
u(2, t) = 4, 0 < t < ∞



u(x, 0) = sin(4πx) + sin(6πx),



3
5
u(x, 0) = x2 , 0 ≤ x ≤ 2

0≤x≤1

ut = 3uyy , 0 < x < L




u(0,
t) = 1, 0 < t < ∞


u
=
u
,
0
<
x
<
π
4.
.
u(L,
t) = 1, 0< t <∞
 t
xx


ux (0, t) = ux (π, t) = 0, 0 < t < ∞ .
2.

 u(x, 0) = 2 sin 3πy , 0 ≤ x ≤ L


u(x, 0) = x, 0 ≤ x ≤ π
L
B. Analise as afirmações abaixo, referentes a problemas simples associados à equação do calor.
Formule o problema correspondente a cada uma delas e resolva-os. A idéia é perceber a relação
entre fenômenos facilmente observáveis e os termos da solução da equação do calor.
1. Uma barra isolada com temperatura inicial constante permancerá a temperatura constante
para t > 0.
2. Uma barra inicialmente a temperatura constante T0 , e com temperaturas idênticas T0
especificadas nas fronteiras permanecerá com temperatura constante T0 para t > 0.
C. Vamos analisar um pouco mais em detalhe a solução do problema de Neumann para a equação
do calor do tipo:

ut = κuxx , 0 < x < L



ux (0, t) = α, 0 < t < ∞
,
u (L, t) = β, 0 < t < ∞


 x
u(x, 0) = f (x), 0 ≤ x ≤ L
com α, β ∈ IR. Considere as questões abaixo.
1
1. Suponha, inicialmente, α = β = 0. Mostre que lim u(x, t) =
t→∞
L
de f (x) no domı́nio [0, L].
Z
L
f (x)dx, que é a média
0
2. Suponha, agora, que α = β = Q 6= 0. Explique fisicamente o que está acontecendo neste
problema (atenção na interpretação!). Fazendo u(x, t) = w(x, t) + U (x) e, escolhendo
adequadamente w(x, t) e U (x), encontre os problemas que devem ser resolvidos para estas
novas incógnitas.
3. Resolva o problema para U (x). Note que sua escolha anterior não especifica unicamente
U (x). Neste ponto, não há nada que você possa fazer para determinar unicamente esta
função. Isto faz sentido?
Métodos Matemáticos da Fı́sica 1 - 02/14 - Prof. Yuri Dumaresq Sobral - MAT/UnB
4. Suponha, agora, que f (x) = A, com A constante. Resolva o problema para w(x, t) e
encontre, finalmente, u(x, t). Seu resultado é único, isto é, há algum problema em não
podermos determinar unicamente U (x)? Determine lim u(x, t) e comente seu resultado.
t→∞
5. Suponha que, agora, tenhamos α 6= β e que α, β 6= 0. Mostre que não é possı́vel resolver
este problema com o procedimento sugerido no item 2. Por quê? Aprenderemos a resolver
este problema mais para frente.
D. (Condições de contorno distintas em cada extremidade) Considere que uma barra
unidimensional de comprimento π e difusividade térmica κ = 1 é dividida em duas partes de
tamanhos iguais. Uma parte é resfriada a uma temperatura de T = 0 e a outra é aquecida a
uma temperatura T = 1. No instante t = 0, as duas partes da barra são colocados uma ao
lado da outra, em contato, formando novamente a barra de comprimento π. A extremidade da
parte que estava resfriada (x = 0) é mantida a T = 0, enquanto a outra extremidade (x = π) foi
isolada. Formule este problema matematicamente e encontre a solução T (x, t). Note que este
é um problema com condições de contorno de tipos diferentes em cada uma das extremidades.
Desta forma, fique atento na hora de resolver o problema de autovalores-autofunções! Use o
computador para traçar a solução obtida e analise seu resultado.
E. (Condições de contorno de Robin) Considere o problema da equação do calor abaixo,
em que umas das condições de contorno é uma condição do tipo mista (Robin):

ut = uxx , 0 < x < 1



u(0, t) = 0, 0 < t < ∞
u(1,
t) + ux (1, t) = 0, 0 < t < ∞



u(x, 0) = x, 0 < x < 1
1. Aplique o método de separação de variáveis u(x, t) = F (x)G(t) e formule o problema de
valor de contorno para a função F (x). Mostre, então, que os autovalores σ são aqueles
que satisfazem a equação transcendental tan σ = −σ. Trace um gráfico das funções tan σ
e −σ e identifique quantos autovalores existem. As três primeiras raı́zes da equação transcendental acima são σ1 ≈ 2.02, σ2 ≈ 4.91 e σ2 ≈ 7.98. Note que os autovalores σ não são
mais múltiplos inteiros de π e, portanto, as autofunções associadas a cada um deles não
desenvolvem um número inteiro de perı́odos no intervalo [0, 1]. (Se quiser, verifique isto
num computador)
2. Obtenha que o problema de valor inicial requer a função u(x, 0) seja escrita em uma série
de Fourier de senos. Usando a ortogonalidade das funções sin(σn x), mostre que:
Z
1
x sin(σn x)dx
0
bn = Z
=
1
sin2 (σn x)dx
2σn
sin(σn ) − σn cos(σn )
.
σn − sin(σn ) cos(σn )
σn2
0
3. Calcule os três primeiros bn e mostre que
u(x, t) ≈ 0.24e−4t sin(2.02x) + 0.22e−24t sin(4.91x) + 0.03e−63t sin(7.98x) + · · ·
Métodos Matemáticos da Fı́sica 1 - 02/14 - Prof. Yuri Dumaresq Sobral - MAT/UnB
Download

Métodos Matemáticos da F´ısica 1