Capı́tulo 4

Equações diferenciais parciais

4.1. Introdução

O problema de equações diferenciais parciais aparece quando em determinadas aplicações aparecem

modelos matemáticos em que existem duas ou mais variáveis independentes. Nesse tipo de aplicações, o
modelo matemático é representado por equações diferenciais parciais.
Para resolver uma equação diferencial apresentamos o método de separação de variáveis. A idéia
central desse método consiste em substituir a equação diferencial por um conjunto de equações diferenciais
ordinárias que são resolvidas sujeitas a condições iniciais ou de contorno. Portanto, para desenvolver o método
de separação de variáveis devemos revisar as propriedades básicas de problemas de valores de contorno para
equações diferenciais ordinárias, assim como das propriedades de séries de Fourier em senos e/ou cosenos.

4.2. Problemas de valores de contorno para fronteiras com dois pontos

Quando estudamos as equações diferenciais ordinárias da forma:

00 0
y + p(t) y + q(t) y = g(t) (4.1)

normalmente informamos duas condições iniciais em um mesmo ponto, geralmente em t = to , da

seguinte forma:

0 0
y(to ) = yo y (to ) = yo (4.2)

Entretanto, existem casos práticos de problemas em que o valor da variável dependente y ou de sua
derivada é especificado em dois pontos diferentes. Esse tipo de informação é chamada de condições de
contorno que são tratadas de forma diferente às condições iniciais que especificam os valores de y e de
0
y no mesmo ponto.
Uma equação diferencial junto com as condições de contorno adequadamente especificadas forma um
problema de valores de contorno com dois pontos como a seguinte:

00 0
y + p(t) y + q(t) y = g(t) (4.3)

1
junto com as condições de contorno:

y(α) = yo y(β) = y1 (4.4)

A solução do problema (4.3)-(4.4) assume a forma y = φ(x) que satisfaz (4.3) no intervalo α < x < β e que
tem valores especificados yo e y1 nos extremos do intervalo. Em geral encontramos primeiro a solução geral
da equação diferencial e depois usamos as condições de contorno para encontrar os valores das constantes
arbitrárias.
Observações: As seguintes observações são muito importantes:

1. Se g(t) = 0 para todo t e se os valores de yo e y1 também são nulos então o sistema (4.3)-(4.4)
é homogêneo. Em caso contrário, o problema é não homogêneo.

2. Os problemas de valor inicial (4.1)-(4.2) e de valores de contorno com dois pontos (4.3)-(4.4) apre-
sentam soluções que diferem em aspectos importantes. Assim, o problema de equação diferencial de
valor inicial geralmente tem uma única solução. Por outro lado, as equações diferenciais de valores
de contorno com dois pontos podem ter uma única solução, não ter solução ou ter uma infinidade de
soluções. Neste aspecto, os problemas de equações diferenciais lineares com valores de contorno são
parecidas com as equações algébricas lineares.

3. Relembrando os sistemas lineares consideremos o seguinte sistema algébrico linear:

Ax=b (4.5)

onde A é uma matriz de dimensão n × n, b é um vetor de dimensão n × 1 e x é um vetor de dimensão

n × 1. Nesse contexto temos o seguinte:

Se A tem inversa então o sistema (4.5) tem solução única.

Se A é singular (não tem inversa) então o sistema (4.5) não tem solução ou tem infinitas soluções.

4. Seja o sistema algébrico linear homogêneo:

Ax=0 (4.6)

obtido de (4.5) com b = 0. Nesse contexto temos o seguinte:

O sistema (4.6) sempre tem a solução trivial x = 0.

Se A tem inversa então x = 0 é a única solução.
Se A é singular (não tem inversa) então o sistema (4.6) tem uma infinidade de soluções não nulas.

5. Solução de uma equação diferencial de segunda ordem com coeficientes constantes homogênea:
Um dos tópicos principais analisados em equações diferenciais ordinárias é a solução geral de uma
equação diferencial de segunda ordem com coeficientes constantes homogênea. Essa equação tem a
seguinte estrutura:

00 0
ay +by +cy =0

em que a, b e c são parâmetros conhecidos. A solução geral desse sistema assume a seguinte forma:

2
 y(x) = c1 y1 (x) + c2 y2 (x)

onde y1 (x) e y2 (x) são soluções da equação diferencial de forma que o wronskiano dessas funções seja
diferente de zero. Neste contexto, um wronskiano diferente de zero significa que y1 (x) e y2 (x) devem
cumprir a seguinte relação:

0 0
y1 (x) y2 (x) − y1 (x) y2 (x) 6= 0

Uma estratégia para encontrar y1 (x) e y2 (x) é identificar uma forma geral da solução e, a partir dessa
estrutura, encontrar duas soluções que cumpram com a exigência de que o wronskiano seja diferente de
zero (soluções linearmente independentes). Pode-se verificar que uma solução para a equação diferencial
assume a seguinte forma genérica:

y(x) = erx

cujas derivadas assumem a seguinte forma:

0 00
y (x) = r erx y (x) = r2 erx

Substituı́ndo essas relações na equação diferencial temos o seguinte:

h i
a r2 erx + b r erx + c erx = erx a r2 + b r + c = 0 =⇒

a r2 + b r + c = 0

A relação anterior é chamada de equação caracterı́stica da equação diferencial. A solução desse sistema
assume a seguinte forma:
√ √
−b + b2 − 4ac −b − b2 − 4ac
r1 = e r2 =
2a 2a

Neste caso existem três tipos de soluções:

√
Quando: b2 − 4ac > 0:
Nesse caso as duas soluções que podem ser obtidas, isto é, y1 (x) = er1 x e y2 (x) = er2 x cumprem
com a exigência do wronskiano e, portanto, a solução geral assume a seguinte forma:

y(x) = c1 er1 x + c2 er2 x

√
Quando: b2 − 4ac = 0:
Nesse caso as duas soluções são iguais. Em outras palavras temos apenas uma solução. Seja
b
α = 2a . Pode-se mostrar que se a primeira solução é y1 (x) = e−αx então a segunda solução pode
ser y2 (x) = x e−αx . Assim, a solução geral assume a seguinte forma:

y(x) = c1 e−α x + c2 x e−α x

3
 √
Quando: b2 − 4ac < 0:
√
Nesse caso as duas soluções são chamadas de complexas conjugadas. Seja β = 4ac − b2 . Assim,
os valores de r1 e r2 assumem a seguinte forma:

r1 = −α + i β r2 = −α − i β

Nesse contexto, pode-se mostrar que a solução geral assume a seguinte forma:

y(x) = e−α x [c1 Sen(βx) + c2 Cos(βx)]

Exemplo 1: Equação diferencial com solução única:

Resolver a equação diferencial com valores de contorno em dois pontos:

00
y + 2y = 0 y(0) = 1 y(π) = 0

A equação caracterı́stica da equação diferencial assume a seguinte forma:

√
r2 + 2 = 0 =⇒ r2 = −2 =⇒ r = ± 2i

Portanto, as soluções assumem a seguinte forma:

√ √
y1 = Cos( 2x) y2 = Sen( 2x)

Assim, a solução geral assume a seguinte forma:

√ √
y = c1 Cos( 2x) + c2 Sen( 2x)

Aplicando as condições de contorno temos o seguinte:

1 = c1
√ √
0 = c1 Cos( 2π) + c2 Sen( 2π)

√ √
=⇒ c1 = 1 =⇒ Cos( 2π) + c2 Sen( 2π) = 0

√
=⇒ c2 = −Cotag( 2π) = −0,276

Dessa forma, a solução final assume a seguinte forma:

√ √ √
y = Cos( 2x) − Cotag( 2π) Sen( 2x)

Deve-se observar que a equação diferencial do exemplo 1 tem uma única solução.
Exemplo 2: A equação diferencial não tem solução única:
Resolver a equação diferencial com valores de contorno em dois pontos:

4
 00
y +y =0 y(0) = 1 y(π) = a

A equação caracterı́stica da equação diferencial assume a seguinte forma:

r2 + 1 = 0 =⇒ r2 = −1 =⇒ r = ±i

Portanto, as soluções assumem a seguinte forma:

y1 = Cos(x) y2 = Sen(x)

Assim, a solução geral assume a seguinte forma:

y = c1 Cos(x) + c2 Sen(x)

Aplicando as condições de contorno temos o seguinte:

1 = c1 =⇒ c1 = 1
a = c1 Cos(π) + c2 Sen(π)

=⇒ c1 Cos(π) + c2 Sen(π) = a =⇒ c1 Cos(π) = a

=⇒ c1 = −a

Dessa forma, a solução da equação diferencial assume a seguinte forma:

c1 = 1 e c1 = −a c2 qualquer valor

Para a solução encontrada temos o seguinte:

Se a 6= −1 o sistema não tem solução.

Se a = −1 o sistema aceita infinitas soluções com c1 = 1 e c2 assumindo qualquer valor e, nesse

contexto, a solução geral assume a seguinte forma:

y = Cos(x) + c2 Sen(x)

Exemplo 3: A equação diferencial tem solução única:

Resolver a equação diferencial com valores de contorno em dois pontos:

00
y + 2y = 0 y(0) = 0 y(π) = 0

O problema anterior é uma equação diferencial homogênea da seguinte forma:

5
 00 0
y + p(x) y + q(x) y = 0 (4.7)

com condições de contorno:

y(α) = 0 y(β) = 0 (4.8)

que tem a solução trivial y = 0 (solução trivial nula) mas pode ter soluções não triviais que pode ser de
interesse encontrar.
A forma geral da solução é a mesma do exemplo 1 e, portanto, a solução geral assume a seguinte forma:
√ √
y = c1 Cos( 2x) + c2 Sen( 2x)

Aplicando as condições de contorno temos o seguinte:

0 = c1 =⇒ c1 = 0
√ √ √
0 = c1 Cos( 2π) + c2 Sen( 2π) =⇒ c2 Sen( 2π) = 0 =⇒ c2 = 0
Portanto, y = 0 é a única solução da equação diferencial.
Exemplo 4: A equação diferencial tem infinitas soluções:
Resolver a equação diferencial com valores de contorno em dois pontos:

00
y +y =0 y(0) = 0 y(π) = 0

A equação diferencial é a mesma do exemplo 2 e, portanto, a solução geral assume a seguinte forma:

y = c1 Cos x + c2 Sen x

Aplicando as condições iniciais temos o seguinte:

0 = c1 =⇒ c1 = 0
0 = c1 Cos(π) + c2 Sen(π) =⇒ c2 Sen(π) = 0 =⇒ c2 qualquer valor
Portanto, a solução assume a seguinte forma:

y = c2 Sen x

que mostra que a equação diferencial tem infinitas soluções.

Observações: As seguintes observações são muito importantes:

1. Os exemplos apresentados mostram uma relação entre problemas de equações diferenciais com valores
de contorno com dois pontos homogêneos e não homogêneos com o que acontece com os sistemas
algébricos lineares homogêneos e não homogêneos. Se um problema de equação diferencial com valores
de contorno não homogêneo tem uma única solução e o problema homogêneo correspondente tem
apenas a solução trivial. Além disso, se um problema de equação diferencial não homogênea não tem
solução ou tem infinitas soluções e o problema homogêneo correspondente tem (infinitas) soluções não
triviais.

6
 2. Sabemos que o sistema algébrico linear:

Ax=λx (4.9)

tem a solução trivial x = 0 para todo valor de λ. Adicionalmente, para determinados valores de λ,
chamados de autovalores, existem soluções não triviais. Um comportamento semelhante acontece nas
equações diferenciais.

4.3. Autovalores e autofunções em equações diferenciais

Consideremos a equação diferencial:

00
y + λy = 0 (4.10)

com as seguintes condições de contorno:

y(0) = 0 y(π) = 0 (4.11)

Nesse contexto, os valores de λ para os quais existem soluções não triviais do problema (4.10)-(4.11) são
chamados de autovalores e as soluções não triviais correspondentes são chamados de autofunções.
Vamos encontrar os autovalores e as autofunções do problema (4.10)-(4.11). Para isso, consideramos três
casos: (a) λ > 0, (b) λ = 0 e (c) λ < 0. Deve-se observar que em cada um desses casos a solução assume
uma estrutura diferente.

1. Quando λ > 0:
Neste caso mudamos a notação para λ = µ2 para facilitar os cálculos. Assim, temos o seguinte:

00
y + µ2 y = 0 (4.12)

A equação caracterı́stica dessa equação diferencial assume a seguinte forma:

r2 + µ2 = 0 =⇒ r = ±i µ

Portanto, a solução geral assume a seguinte forma:

y = c1 Cos(µx) + c2 Sen(µx) (4.13)

Aplicando as condições de contorno temos o seguinte:
0 = c1
0 = c2 Sen(µπ)
As soluções não triviais acontecem com c2 6= 0 =⇒ Sen(µπ) = 0. Portanto, os valores permitidos para
µ são os seguintes: µ = 1, 2, 3, . . .. Como λ = µ2 então os valores de λ que representam soluções da
equação diferencial são os seguintes:

7
 λ1 = 1
λ2 = 4
λ3 = 9 (4.14)
..
.
λ n = n2

Assim, os autovalores da equação diferencial assumem a seguinte forma:

λn = n2 n = 1, 2, 3, . . . (4.15)

As autofunções são obtidas de (4.13) e assumem a seguinte forma:

yn = c2 Sen(n x)

Assumindo c2 = 1, as autofunções da equação diferencial assumem a seguinte forma:

y1 = Sen x
y2 = Sen 2x
y3 = Sen 3x
..
.
yn = Sen nx

Usando notação compacta, as autofunções da equação diferencial assumem a seguinte forma:

yn (x) = Sen(n x) (4.16)

2. Quando λ < 0:
Fazemos λ = −µ2 < 0 =⇒ µ > 0. Assim, em (4.10) a equação caracterı́stica assume a seguinte forma:

r2 − µ2 = 0 =⇒ r2 = µ2 =⇒ r = ±µ

Nesse contexto, a solução geral da equação diferencial assume a seguinte forma:

y = k1 eµx + k2 e−µx

1 1
Fazendo: k1 = 2 (c1 + c2 ) e k2 = 2 (−c1 + c2 ), a relação anterior assume a seguinte forma:
· ¸ · ¸
1 ¡ µx ¢ 1 ¡ µx ¢
y = c1 e − e−µx + c2 e + e−µx
2 2

A relação anterior assume a seguinte forma:

y = c1 Senh(µ x) + c2 Cosh(µ x) (4.17)

8
 Aplicando as condições de contorno na solução geral mostrada temos o seguinte:
0 = c2 =⇒ c2 = 0
0 = c1 Senh(µπ) =⇒ c1 = 0 porque Senh(µπ) 6= 0
Das relações anteriores podemos concluir que y = 0 é a única solução e não existem soluções não
triviais. Portanto, quando λ < 0 não existem soluções não triviais.
3. Quando λ = 0:
Neste caso a equação caracterı́stica assume a seguinte forma:

r2 = 0 =⇒ r = ±0

que mostra que as soluções são reais e iguais. Nesse caso a solução geral assume a seguinte forma:

y = c1 e0 + c2 x e0 =⇒ y = c1 + c2 x (4.18)

Aplicando as condições de contorno na solução geral temos o seguinte:

0 = c1 =⇒ c1 = 0
0 = c1 + π c2 =⇒ c2 = 0 =⇒ y = 0
Os resultados encontrados mostram que a equação diferencial tem apenas a solução trivial (y = 0).
Portanto, quando λ = 0 não existem soluções não triviais.

A análise anterior da equação diferencial representada pelo sistema (4.10)-(4.11) nos permite concluir
que essa equação diferencial apresenta soluções não triviais apenas quando λ > 0 e nesse caso o sistema
(4.10)-(4.11) que tem a seguinte estrutura:

00
y + λy = 0 y(0) = 0 y(π) = 0 (4.19)

tem uma sequência infinita de autovalores positivos da forma:

λn = n2 n = 1, 2, 3, . . . (4.20)

e as autofunções correspondentes tem a seguinte estrutura:

yn = k Sen(n x) (4.21)

sendo k uma constante.

Exemplo 5: Um caso de aplicação do tópico anterior:
Considere o sistema (4.10)-(4.11) em que π = L. Para esse caso particular encontre os autovalores e as
autofunções dessa equação diferencial.
Neste caso temos a seguinte equação diferencial com as condições de contorno em dois pontos:

00
y + λy = 0 y(0) = 0 y(L) = 0 (4.22)

A solução é a mesma mostrada em (4.20)-(4.21) exceto na aplicação da segunda condição de contorno.

Assim, a solução assume a seguinte forma:

9
 y = c1 Cos(µ x) + c2 Sen(µ x)

Aplicando as condições de contorno temos o seguinte:

0 = c1
0 = c2 Sen(µL) =⇒ µL = π, 2π, 3π, 4π, . . .
Das relações anteriores encontramos a seguinte relação:

µ ¶
nπ
µn = (4.23)
L

Portanto, os autovalores assumem a seguinte forma:

n2 π 2
λ = µ2 =⇒ λn =
L2

Finalmente, as autofunções assumem a seguinte forma:

µ ¶
nπ x
yn = c2 Sen(µn x) =⇒ yn = Sen n = 1, 2, 3, . . . (4.24)
L

Os importantes resultados encontrados neste exemplo são usados com frequência nos seguintes tópicos
deste capı́tulo. Os tópicos apresentados são usados para analisar três tipos muito especiais de equações
diferenciais parciais.

4.4. A solução da equação do calor

Três tipos de equações diferenciais parciais são os mais estudados e analisados porque tem aplicação em
várias áreas da engenharia: (i) a equação do calor, (b) a equação de onda e (iii) a equação do potencial.
Todos esses tipos especiais de problemas são resolvidos usando o método de separação de variáveis.
A equação diferencial de segunda ordem parcial, chamada de equação do calor, tem a seguinte estrutura:

∂ 2 u(x, t) ∂u(x, t)
α2 = 0<x<L t>0 (4.25)
∂x2 ∂t

O problema (4.25) representa um conjunto de problemas que aparecem na engenharia sendo o problema
de condução de calor em um corpo sólido o mais representativo e o mais conhecido. Esse tipo de problema
aparece em todos os sistemas em que aparecem a transferência de calor em dispositivos.
Neste caso analisamos o problema especı́fico de condução de calor em uma barra de seção reta constante
como mostrada na figura 4.1. O material é homogêneo e a barra está alinhada com o eixo x. Neste problema
especı́fico α2 é uma constante de difusividade térmica que depende apenas do material do qual é feita a
barra. Assim, a seguinte relação é verdadeira:

k
α2 = (4.26)
ρs

10
 Figura 4.1: Uma barra metálica condutora.

cm2
onde k é a condutividade térmica, ρ é a densidade do material e s é o calor especı́fico. α está em s .

Adicionalmente, não existe transferência de calor pela superficie lateral da barra e os extremos da barra
coincide com x = 0 e x = L. Assim, a equação diferencial parcial de segunda ordem mostrada em (4.25)
assume a seguinte forma:

α2 uxx = ut 0<x<L t>0 (4.27)

onde foi usada notação compacta para indicar a derivada parcial.

Outras suposições particulares para o problema: Para o problema em análise, deve-se ainda especifi-
carse condições iniciais e/ou de contorno. Assim, especificamos o seguinte:

1. Condição inicial:
A distribuição inicial de temperatura na barra satisfaz a relação:

u(x, 0) = f (x) 0≤x≤L (4.28)

onde f (x) é uma função conhecida e especificada em cada problema especı́fico.

2. Condições de contorno:
As temperaturas nos extremos da barra são conhecidas e mantidas fixas. Nesse caso, vamos fixar as
temperaturas nos extremos em zero, isto é, se a temperatura é T1 em x = 0 e T2 em x = L então
T1 = T2 = 0 e temos as seguintes condições de contorno:

u(0, t) = 0 u(L, t) = 0 t>0 (4.29)

Portanto, na equação diferencial parcial que pretendemos resolver relacionado com o problema funda-
mental de condução de calor, pretende-se encontrar a função de distribuição de temperatura u(x, t) que
satisfaz a equação diferencial (4.27), a condição inicial (4.28) e as condições de contorno (4.29).
O problema (4.27)-(4.28)-(4.29) é linear porque u aparece em todas as equações na primeira potência.
Também, o problema é homogêneo em relação às condições de contorno. Assim, uma estratégia para resolver
o problema é encontrar uma solução para a equação diferencial e as condições de contorno. Na sequência
podemos fazer uma superposição para satisfazer a condição inicial.

11
4.4.1. Estratégia de solução

Uma solução da equação diferencial (4.27) que satisfaz as condições de contorno (4.28) pode ser a função
u(x, t) = 0 mas essa proposta de solução não satisfaz a condição inicial exceto para o caso trivial de f (x) = 0.
Para encontrar uma solução não trivial para u(x, t) assumimos a hipótese de que u(x, t) pode ser repre-
sentada como um produto de duas funções, uma que depende apenas de x e outra que depende apenas de
t. Sejam X(x) e T (t) essas funções então temos o seguinte:

u(x, t) = X(x) T (t) (4.30)

Assim temos o seguinte:

∂u(x, t) 0 ∂ 2 u(x, t) 00 ∂u(x, t) 0

= ux = X (x)T (t) =⇒ 2
= uxx = X (x) T (t) = ut = X(x) T (t)
∂x ∂x ∂t

Substituı́ndo as relações anteriores em (4.27) temos o seguinte:

00 0
00 0 X (x) 1 T (t)
α2 X (x) T (t) = X(x) T (t) =⇒ = 2 (4.31)
X(x) α T (t)

A relação (4.31) mostra uma separação de variáveis. A única forma de garantir a igualdade é que ela
seja igual a uma constante que identificaremos por −λ. Assim, temos o seguinte:

00 0
2 00 X (x)
0 1 T (t)
α X (x) T (t) = X(x) T (t) =⇒ = 2 = −λ (4.32)
X(x) α T (t)

Podemos encontrar a solução de cada equação diferencial para qualquer valor de λ resolvendo (4.32) para
encontrar X(x) e T (t). Assim, a solução é o produto dessas funções. Para encontrar a forma matemática de
X(x) e T (t) repassamos as condições de contorno em u(x, t) para X(x).
Para x = 0 temos o seguinte:

u(0, t) = 0 = X(0) T (t) =⇒ X(0) = 0 (4.33)

já que T (t) = 0 levaria para u(x, t) = 0 que é a solução trivial que não interessa.
Para x = L temos o seguinte:

u(L, t) = 0 = X(L) T (t) =⇒ X(L) = 0 (4.34)

Assim, podemos resolver a primeira parcela de (4.32):

00
X (x) 00
= −λ =⇒ X (x) + λ X(x) = 0 X(0) = 0 X(L) = 0 (4.35)
X(x)

A relação (4.35) é uma equação diferencial homogênea com condições de contorno já analisada anterior-
mente (veja exemplo 5). Os autovalores desse problema assumem a seguinte forma:

12
 n2 π 2
λn = n = 1, 2, 3, . . . (4.36)
L2
As autofunções assumem a seguinte forma:

µ ¶
nπx
Xn (x) = Sen n = 1, 2, 3, . . . (4.37)
L

Agora resolvemos a outra equação diferencial de (4.32) para o valor de λ encontrado. Assim, temos o
seguinte:

0
1 T (t) 0 0
= −λ =⇒ T (t) = −λα2 T (t) =⇒ T (t) + λα2 T (t) = 0 =⇒
α2 T (t)

à !
0 α 2 n2 π 2
Tn (t) + Tn (t) = 0 (4.38)
L2

A solução de (4.38) é trivial e assume a seguinte forma:

n2 π 2 α2
Tn (t) = e− L2
t
n = 1, 2, 3, . . . (4.39)

que é obtido usando o seguinte processo:

0 Z 0 Z
T (t) n2 π 2 α2 T (t) n2 π 2 α 2 n2 π 2 α 2 n2 π 2 α2
=− =⇒ dt = − dt =⇒ Ln[T (t)] = − t =⇒ Tn (t) = e− L2 t
T (t) L2 T (t) L 2 L 2

Assim, uma solução do problema assume a seguinte forma:

µ ¶
n2 π 2 α2 nπx
un (x, t) = Xn (x)Tn (t) = e− L2
t
Sen n = 1, 2, 3, . . . (4.40)
L

que satisfaz a equação diferencial parcial (4.27) e as condições de contorno (4.29) para cada n inteiro e
positivo. As funções un (x, t) são chamadas de soluções fundamentais do problema de condução de calor
(4.27)-(4.28)-(4.29). Entretanto, falta satisfazer a condição inicial (4.28), isto é, a relação:

u(x, 0) = f (x) 0≤x≤L (4.41)

Uma solução u(x, t) do problema é uma combinação linear de (4.40) que satisfaz a condição inicial (4.41).
Assim, uma solução geral assume a seguinte forma:

∞
X ∞
X µ ¶
n2 π 2 α2 nπx
u(x, t) = cn un (x, t) = cn e− L2
t
Sen (4.42)
n=1 n=1
L

onde os coeficientes cn ainda não são conhecidos. Lembremos que os termos individuais de (4.42) satisfazem
a equação diferencial parcial (4.27) e as condições de contorno (4.29). Agora, supor que a série infinita (4.42)
converge e satisfaz (4.27) e (4.29). Nesse contexto, para satisfazer a condição inicial (4.28), deve-se cumprir
o seguinte:

13
 ∞
X µ ¶
nπx
u(x, 0) = f (x) = cn Sen (4.43)
n=1
L

Em outras palavras, os coeficientes cn devem ser escolhidos de forma que a série mostrada em (4.43)
represente a função f (x) que representa a função de distribuição inicial de temperatura para 0 ≤ x ≤ L.
Assim, a série mostrada em (4.43) é a série de Fourier em senos de f (x) e seus coeficientes são encontrados
usando a relação:

Z L µ ¶
2 nπx
cn = f (x) Sen dx (4.44)
L 0 L

A solução final é encontrada substituı́ndo os coeficientes cn obtidos em (4.44) na relação geral (4.42).
Exemplo 6: Aplicação da equação do calor.
Encontre a temperatura u(x, t) em qualquer instante em uma barra metálica com 50 cm. de comprimento
e com superficie lateral isolada. A temperatura inicial é de 20o C em toda a barra e as extremidades são
mantidas a 0o C para todo t > 0.
Sabemos que a forma geral da solução assume a seguinte forma:

∞
X µ ¶
n2 π 2 α2 nπx
u(x, t) = cn e− L2
t
Sen
n=1
L

Adicionalmente sabemos que: L = 50 cm. e f (x) = 20 para 0 < x < 50 cm.

Encontramos os elementos de cn da seguinte forma:

Z L µ ¶ Z 50 µ ¶ Z 50 µ ¶
2 nπx 2 nπx 4 nπx
cn = f (x) Sen dx = 20 Sen dx = Sen dx
L 0 L 50 0 50 5 0 50

µ ¶· µ ¶¸50
4 50 nπx 40
cn = −Cos = [1 − Cos(nπ)]
5 nπ 50 0 nπ

Quando n é impar então Cos(nπ) = −1 e portanto (1−Cos(nπ)) = 2 e quando n é par então Cos(nπ) = 1
e portanto (1 − Cos(nπ)) = 0. Assim, cn assume a seguinte forma:

80
cn = para n = 1, 3, 5, . . . cn = 0 para n = 2, 4, 6, . . .
nπ

Finalmente, a temperatura u(x, t) assume a seguinte forma:

∞ µ ¶ µ ¶
80 X 1 − n2 π2 α2 t nπx
u(x, t) = e 2500 Sen
π n=1,3,5,... n 50

∞ µ ¶ µ ¶
80 X 1 (2k+1)2 π 2 α2 (2k + 1)πx
u(x, t) = e− 2500 t Sen
π k=0 2k + 1 50

14
4.5. A solução da equação de onda: Vibrações de uma corda elástica

A equação de onda representa um segundo tipo de equação diferencial parcial importante. Esse tipo
de equação diferencial aparece na análise de fenômenos relacionados com a propagação de ondas em um
meio contı́nuo. Assim, são casos de aplicação especı́ficos a análise de ondas acústicas, ondas sı́smicas, ondas
eletromagnéticas, etc. Neste caso, escolhemos para análise especı́fico a vibração mecânica de uma corda
elástica.
O tipo de problema mais simples de analisar é a vibração mecânica de uma corda elástica de comprimento
L que se encontra ligeiramente esticada entre dois suportes no mesmo nı́vel horizontal de forma que o eixo
x se encontre alinhado com a corda como se mostra na figura 4.2.

Figura 4.2: Uma corda elástica.

Supor que a corda é colocada em movimento no plano vertical por uma ação externa (puxando verti-
calmente a corda de alguma forma). Nesse contexto, a função u(x, t) representa o deslocamento da corda
no ponto x e no instante t. Considerando algumas hipóteses simplificadoras (eliminando a resistência do
ar e para amplitudes de movimento não muito grandes) pode-se provar que u(x, t) cumpre com a seguinte
equação diferencial parcial:

∂ 2 u(x, t) ∂ 2 u(x, t)
a2 = =⇒ a2 uxx = utt 0 < x < L; t>0 (4.45)
∂x2 ∂t2

A equação diferencial parcial (4.45) é chamada de equação de onda. Para o caso particular da vibração
da corda o coeficiente a é dada pela relação:

T
a2 =
ρ

onde T é a tensão da corda e ρ é a massa por unidade de comprimento da corda. Assim, a tem unidade de
velocidade e representa a velocidade de propagação das ondas ao longo da corda.
Outras suposições particulares para o problema: Para o problema em análise, deve-se ainda especifi-
carse condições iniciais e/ou de contorno. Assim, especificamos o seguinte:

1. Condições de contorno:
A corda tem suas extremidades fixas:

15
 u(0, t) = 0 u(L, t) = 0 t≥0 (4.46)

2. Condição iniciais:
Devemos especificar a posição inicial da corda e sua velocidade inicial:

u(x, 0) = f (x) 0≤x≤L (4.47)

∂u(x, 0)
= ut (x, 0) = g(x) 0≤x≤L (4.48)
∂t
onde f (x) e g(x) são funções conhecidas e dadas como parte do problema.

3. Condições adicionais:
Para que as condições impostas sejam consistentes, deve-se cumprir adicionalmente o seguinte:

f (0) = f (L) = 0 g(0) = g(L) = 0 t>0 (4.49)

já que a corda tem seus extremos fixos na posição zero e, estando fixas, a velocidade nos extremos da
corda sempre é nula.

A solução da equação diferencial parcial chamada de equação de onda deve encontrar uma função
u(x, t) (a posição da corda em cada ponto x e para todo t ≥ 0) que satisfaz (4.45), as condições de contorno
(4.46) e as condições iniciais (4.47)-(4.48). Assim como no problema da equação de calor, o problema da
equação de onda é um problema de condições de contorno na variável espacial x e um problema de valor
inicial na variável temporal t. A forma matemática é a mesma para outras aplicações relacionados com
movimentos tais como ondas acústicas e ondas eletromagnéticas.
O problema da equação de onda pode ser generalizado para duas e três dimensões assumindo a seguinte
forma:

a2 (uxx + uyy ) = utt (4.50)

ou

a2 (uxx + uyy + uzz ) = utt (4.51)

e também devemos generalizar as condições de contorno.

A estratégia de solução também é o método de separação de variáveis. Assim, analisamos tipos particu-
lares do problema da equação de onda.

4.5.1. Corda elástica com deslocamento inicial não nulo

Neste caso a corda é deslocada da posição de equilı́brio e solta em t = 0 com velocidade nula. Assim, a
corda pode vibrar livremente. O deslocamento vertical u(x, t) deve satisfazer as seguintes condições:

a2 uxx = utt 0<x<L t>0 (4.52)

16
 u(0, t) = 0 u(L, t) = 0 t≥0 (4.53)

u(x, 0) = f (x) ut (x, 0) = 0 0≤x≤L (4.54)

onde f (x) deve fornecer informação da configuração da corda em t = 0. Usando o método de separação de
variáveis temos o seguinte:

u(x, t) = X(x) T (t)

0 00
ux = X (x) T (t) =⇒ uxx = X (x) T (t)
0 00
ut = X(x) T (t) =⇒ utt = X(x) T (t)

Substituı́ndo as relações anteriores em (4.52) temos o seguinte:

00 00
2 00 X (x) 1 T (t)
00
a X (x) T (t) = X(x) T (t) =⇒ = 2 = −λ (4.55)
X(x) a T (t)

que podem ser separadas nas seguintes equações diferenciais independentes:

00
X (x) + λ X(x) = 0 (4.56)

00
T (t) + a2 λ T (t) = 0 (4.57)

Encontrando uma solução para (4.56):

Inicialmente devemos repassar as condições de contorno de (4.53) para a equação diferencial (4.56) da
seguinte forma:

u(0, t) = 0 =⇒ 0 = X(0) T (t) =⇒ X(0) = 0

(4.58)
u(L, t) = 0 =⇒ 0 = X(L) T (t) =⇒ X(L) = 0

Assim, a solução da equação diferencial (4.56) sujeito às condições de contorno é um tipo de solução que
já conhecemos e essa solução tem os seguintes autovalores:

n2 π 2
λn = n = 1, 2, 3, . . . (4.59)
L2

e as seguintes autofunções (soluções):

µ ¶
nπx
Xn (x) = Sen (4.60)
L

17
Encontrando uma solução para (4.57):
Agora devemos encontrar uma solução para (4.57) sujeito às condições iniciais e a partir da solução
encontrada para (4.56).
Usando a segunda condição inicial, isto é, ut (x, 0) = 0 temos o seguinte:

0 0
ut (x, 0) = X(x) T (0) = 0 =⇒ T (0) = 0 (4.61)

Também, usando o valor de λ encontrado em (4.59) resolvemos o sistema (4.57):

00 n2 π 2 a2
Tn (t) + Tn (t) = 0 (4.62)
L2

A equação caracterı́stica da equação diferencial anterior assume a seguinte forma:

µ ¶
n2 π 2 a2 nπ a
r2 + = 0 =⇒ r = ± i
L2 L

que indica que as raizes são complexas conjugadas e, nesse caso, a solução geral assume a seguinte forma:

µ ¶ µ ¶
nπa nπa
Tn (t) = k1 Cos t + k2 Sen t (4.63)
L L

Impondo a condição inicial (4.61) em (4.63) temos o seguinte:

µ ¶· µ ¶ µ ¶¸
0 nπa nπa nπa
Tn (t) = −k1 Sen t + k2 Cos t
L L L

0 ¡ nπa ¢
Tn (0) = 0 = L k2 =⇒ k2 = 0. Assim, uma solução Tn (t) é proporcional ao seguinte:

µ ¶
nπa
Tn (t) = Cos t
L

Portanto, as funções que satisfazem a equação diferencial (4.52), as condições de contorno (4.53) e a
segunda condição inicial de (4.54) assumem a seguinte forma:

µ ¶ µ ¶
nπx nπa
un (x, t) = Xn (x) Tn (t) = Sen Cos t n = 1, 2, 3, . . . (4.64)
L L

que representam soluções fundamentais do problema. Falta satisfazer a primeira condição inicial de (4.54).
Para satisfazer essa condição inicial vamos considerar uma superposição das soluções fundamentais mostradas
em (4.64). Assim, u(x, t) assume a seguinte forma:

∞
X ∞
X µ ¶ µ ¶
nπx nπa
u(x, t) = cn un (x, t) = cn Sen Cos t (4.65)
n=1 n=1
L L

Para encontrar as constantes cn usamos a seguinte condição inicial:

18
 u(x, 0) = f (x)

na solução mostrada em (4.65). Assim temos o seguinte:

∞
X µ ¶
nπx
f (x) = u(x, 0) = cn Sen (4.66)
n=1
L

Portanto, os coeficientes cn devem ser os coeficientes da série de Fourier em senos de f (x) e com perı́odo
2L. Assim temos o seguinte:

Z L µ ¶
2 nπx
cn = f (x) Sen dx n = 1, 2, 3, . . . (4.67)
L 0 L

Assim, a solução completa do problema (4.52)-(4.53)-(4.54) é a função u(x, t) mostrada em (4.65) com
os coeficientes cn encontrados usando (4.67).
Observação importante:
¡ ¢ ¡ ¢
Para um valor de n especificado a relação Sen nπx
L Cos nπa
L t em (4.64) é periódica no tempo com
2L nπa 2π 2L
perı́odo na (lembremos que L = ω = T =⇒ T = na ). Assim, essa relação representa um movimento
vibratório ¡da corda
¢
com perı́odo 2L nπa
na . Portanto, ωn = L representa as frequências naturais da corda. O
nπx
fator Sen L representa o padrão de deslocamento que ocorre na corda ao vibrar na frequência dada.
Cada padrão de deslocamento é chamado de modo natural de vibração e é periódico na variável espacial
x com perı́odo 2L nπa
n que é chamado de comprimento de onda do modo de frequência natural L . Os três
primeiros modos naturais são mostrados na figura 4.3.

Figura 4.3: Modos naturais na vibração de uma corda.

Exemplo 7: Movimento de uma corda elástica com deslocamento inicial não nulo.
Considere uma corda vibrante de comprimento L = 30 cm. que satisfaz a condição da equação de onda:

4uxx = utt 0 < x < 30 t>0 (4.68)

As extremidades da corda estão fixas e a corda é colocada em movimento sem velocidade inicial e na
seguinte posição inicial:

19
  x

 10 0 ≤ x ≤ 10
u(x, 0) = f (x) = (4.69)

 (30−x)
20 10 < x ≤ 30

Encontrar a equação de deslocamento da corda u(x, t).

Dados conhecidos: a2 = 4 =⇒ a = 2 e também L = 30. Da mesma forma verificamos que f (0) = f (30) =
0 e ut (x, 0) = g(x) = 0. Portanto, a solução geral assume a seguinte forma:

∞
X µ ¶ µ ¶
nπx 2nπ
u(x, t) = cn Sen Cos t
n=1
30 30

Os valores de cn são encontrados usando a relação (4.67):

½Z 10 µ ¶ Z 30 µ ¶ ¾
2 x nπx (30 − x) nπx
cn = Sen dx + Sen dx
30 0 10 30 10 20 30

½Z 10 µ ¶ Z 30 µ ¶ Z 30 µ ¶ ¾
1 x nπx 3 nπx x nπx
cn = Sen dx + Sen dx − Sen dx
15 0 10 30 10 2 30 10 20 30

Precisamos conhecer o resultado da seguinte integral:

Z µ ¶
nπx
x Sen dx
30

que pode ser facilmente integrada por partes fazendo o seguinte:

u = x =⇒ du = dx
¡ nπx ¢ 30
¡ nπx ¢
dv = Sen 30 dx =⇒ v = − nπ Cos 30

Assim temos o seguinte:

Z µ ¶ µ ¶ Z µ ¶
nπx 30 nπx 30 nπx
x Sen dx = − x Cos + Cos dx
30 nπ 30 nπ 30

Z µ ¶ µ ¶ µ ¶2 µ ¶
nπx 30 nπx 30 nπx
x Sen dx = − x Cos + Sen
30 nπ 30 nπ 30

O resultado encontrado pode ser usado para encontrar cn da seguinte forma:

 " µ ¶ µ ¶2 µ ¶#10 · µ ¶¸30
1 1 30 nπx 30 nπx 3 30 nπx
cn = − x Cos + Sen + − Cos +
15  10 nπ 30 nπ 30 0
2 nπ 30 10

" 
µ ¶ µ ¶2 µ ¶#30 
1 30 nπx 30 nπx
− − x Cos + Sen
20 nπ 30 nπ 30 
10

20
 ( " µ ¶ µ ¶2 µ ¶ # · µ ¶¸
1 1 300 nπ 30 nπ 3 30 nπ
cn = − Cos + Sen +0−0 + − Cos(nπ) − Cos( )
15 10 nπ 3 nπ 3 2 nπ 3

" µ ¶2 µ ¶ µ ¶2 µ ¶#)
1 900 30 300 nπ 30 nπ
− − Cos(nπ) + Sen(nπ) + Cos − Sen
20 nπ nπ nπ 3 nπ 3

· µ ¶ µ ¶¸ µ ¶
1 90 nπ 45 nπ 9 nπ
cn = 2
Sen + 2
Sen =⇒ cn = 2
Sen
15 (nπ) 3 (nπ) 3 (nπ) 3

Assim, a solução geral u(x, t) assume a seguinte forma:

∞
X µ ¶ µ ¶ µ ¶
9 nπ nπx nπt
u(x, t) = 2 2
Sen Sen Cos
n=1
n π 3 30 15

µ ¶X
∞ µ ¶ µ ¶ µ ¶
9 1 nπ nπx nπt
u(x, t) = Sen Sen Cos
π2 n=1
n 2 3 30 15

π 2π
A solução encontrada u(x, t) é periódica no tempo com perı́odo T = 30 s. (ω = 15 = T =⇒ T = 30 s.).

4.5.2. Corda elástica com velocidade inicial não nula

Neste caso vamos supor que a corda é colocada em movimento a partir de sua posição de equilı́brio
(f (x) = 0) mas com uma velocidade inicial dada. Assim, a equação de deslocamento vertical u(x, t) deve
satisfazer as seguintes condições:

a2 uxx = utt 0<x<L t>0 (4.70)

com as condições de contorno:

u(0, t) = 0 u(L, t) = 0 t≥0 (4.71)

e as condições iniciais:

u(x, 0) = 0 ut (x, 0) = g(x) 0≤x≤L (4.72)

onde g(x) é a velocidade inicial da corda no ponto x.

Em relação ao problema da corda elástica com deslocamento inicial não nulo, analisado anteriormente,
verificamos que não mudaram as condições de contorno. Entretanto, a condição inicial foi mudada. Assim
temos o seguinte:

u(x, 0) = 0 =⇒ u(x, t) = X(x) T (t) =⇒ u(x, 0) = 0 = X(x) T (0) =⇒ T (0) = 0 (4.73)

21
 Usando a mesma estratégia utilizada para resolver o caso em que a corda tem deslocamento inicial não
nulo verificamos que a solução Xn (x) assume a seguinte forma:
µ ¶
n2 π 2 nπx
λn = 2
=⇒ Xn (x) = Sen n = 1, 2, 3, . . .
L L

Da mesma forma, a solução geral de Tn (t) assume a seguinte forma:

µ ¶ µ ¶
nπa nπa
Tn (t) = k1 Cos t + k2 Sen t
L L

Com a condição inicial T (0) = 0 temos o seguinte:

Tn (0) = 0 = k1 =⇒ k1 = 0

Portanto, a solução de Tn (t) assume a seguinte forma genérica:

µ ¶
nπa
Tn (t) = Sen t
L

A solução fundamental un (x, t) assume a seguinte forma:

un (x, t) = Xn (x) Tn (t)

Para satisfazer a condição inicial (não homogênea) restante vamos supor que u(x, t) pode ser representado
como uma combinação linear das soluções fundamentais. Assim, u(x, t) assume a seguinte forma:

∞
X ∞
X µ ¶ µ ¶
nπ nπa
u(x, t) = kn un (x, t) = kn Sen x Sen t (4.74)
n=1 n=1
L L

Para encontrar os parâmetros kn usamos a segunda condição inicial:

ut (x, 0) = g(x)

∞ µ
X ¶ µ ¶ µ ¶
nπa nπ nπa
ut (x, t) = kn Sen x Cos t
n=1
L L L

Das relações anteriores temos o seguinte:

∞ µ
X ¶ µ ¶
nπa nπ
ut (x, 0) = g(x) = kn Sen x (4.75)
n=1
L L
¡ nπa ¢
Portanto, os termos L kn são os coeficientes da série de Fourier em senos de g(x) com perı́odo 2L.
Assim, temos o seguinte:

µ ¶ Z L µ ¶
nπa 2 nπ
kn = g(x) Sen x dx n = 1, 2, 3, . . . (4.76)
L L 0 L
Portanto, a solução da equação diferencial para a corda elástica com velocidade inicial não nula é a
relação (4.74) com os coeficientes kn obtidos de (4.76).

22
4.5.3. Solução geral para a corda elástica

O problema mais geral da equação diferencial para a corda elástica acontece quando o movimento
é iniciado com a corda deslocada da posição de equilı́brio e com velocidade inicial não nula. Nessa forma
geral devemos resolver o seguinte problema:

a2 uxx = utt 0<x<L t>0 (4.77)

com as condições de contorno:

u(0, t) = 0 u(L, t) = 0 t≥0 (4.78)

e as condições iniciais:

u(x, 0) = f (x) ut (x, 0) = g(x) (4.79)

f (x) define a posição inicial da corda e g(x) a velocidade inicial da corda. Este problema geral também
pode ser resolvido usando o método de separação de variáveis mas pode ser resolvido mais facilmente usando
o princı́pio da superposição.
A solução geral pode ser obtida apenas somando as soluções obtidas em 4.5.1 e 4.5.2 encontradas anteri-
ormente. Para mostrar a validade da proposta consideremos que v(x, t) é a solução do caso do deslocamento
da corda com deslocamento inicial não nulo e seja w(x, t) a solução do caso de deslocamento da corda com
velocidade inicial não nula. Assim, a solução geral para o deslocamento vertical da corda assume a seguinte
forma:

u(x, t) = v(x, t) + w(x, t) (4.80)

Vamos verificar que u(x, t) de (4.80) satisfaz a relação:

a2 uxx = utt =⇒ a2 uxx − utt = 0

∞
X µ ¶ µ ¶
nπ nπa
v(x, t) = cn Sen x Cos t
n=1
L L

∞
X µ ¶ µ ¶
nπ nπa
w(x, t) = kn Sen x Sen t
n=1
L L

a2 uxx − utt = a2 (vxx + wxx ) − (vtt + wtt ) = (a2 vxx − vtt ) + (a2 wxx − wtt ) (4.81)

Derivando v(x, t) temos o seguinte:

∞ µ
X ¶ µ ¶ µ ¶
nπ nπ nπa
vx = cn Cos x Cos t =⇒
n=1
L L L

23
 ∞ µ ¶ µ ¶ µ ¶
X nπ 2 nπ nπa
vxx = − cn Sen x Cos t
n=1
L L L

∞ µ
X ¶ µ ¶ µ ¶
nπa nπ nπa
vt = − cn Sen x Sen t =⇒
n=1
L L L

∞ µ ¶ µ ¶ µ ¶
X nπa 2 nπ nπa
vtt = − cn Sen x Cos t
n=1
L L L

Das relações anteriores podemos verificar facilmente que:

a2 vxx − vtt = 0

Derivando w(x, t) temos o seguinte:

∞ µ
X ¶ µ ¶ µ ¶
nπ nπ nπa
wx = kn Cos x Sen t =⇒
n=1
L L L

∞ µ ¶ µ ¶ µ ¶
X nπ 2 nπ nπa
wxx = − kn Sen x Sen t
n=1
L L L

∞ µ
X ¶ µ ¶ µ ¶
nπa nπ nπa
wt = kn Sen x Cos t =⇒
n=1
L L L

∞ µ ¶ µ ¶ µ ¶
X nπa 2 nπ nπa
wtt = − kn Sen x Sen t
n=1
L L L

Das relações anteriores podemos verificar facilmente que:

a2 wxx − wtt = 0

Assim, provamos que (4.81) assume a seguinte forma:

a2 uxx − utt = (a2 vxx − vtt ) + (a2 wxx − wtt ) = 0

Analisando as condições de contorno temos o seguinte:

u(0, t) = v(0, t) + w(0, t) = 0 + 0 = 0 =⇒ u(0, t) = 0

24
 u(L, t) = v(L, t) + w(L, t) = 0 + 0 = 0 =⇒ u(L, t) = 0

Analisando as condições iniciais temos o seguinte:

u(x, 0) = v(x, 0) + w(x, 0) = f (x) + 0 = f (x) =⇒ u(x, 0) = f (x)

ut (x, 0) = vt (x, 0) + wt (x, 0) = 0 + g(x) = g(x) =⇒ ut (x, 0) = g(x)

Assim, u(x, t) satisfaz todas as condições do problema geral. Portanto, a forma de u(x, t) assume a
seguinte forma:

∞
X µ ¶ µ X∞ ¶ µ ¶ µ ¶
nπ nπa nπ nπa
u(x, t) = cn Sen x Cos t + kn Sen x Sen t (4.82)
n=1
L L n=1
L L

onde os coeficientes cn e kn são encontradas usando as seguintes relações:

Z L µ ¶
2 nπ
cn = f (x) Sen x dx n = 1, 2, 3, . . . (4.83)
L 0 L

µ ¶Z L µ ¶
2 nπ
kn = g(x) Sen x dx n = 1, 2, 3, . . . (4.84)
nπa 0 L

Exemplo 8: Movimento de uma corda elástica com velocidade inicial não nula.
Considere uma corda elástica com comprimento L = 10 cm. e cujas extremidades são mantidas fixas. A
corda é colocada em movimento a partir de sua posição de equilı́brio com a seguinte velocidade inicial:
 x

 5 0≤x≤5
ut (x, 0) = g(x) =

 (10−x)
5 5 < x ≤ 10

Considerando a = 1 encontre a forma matemática de u(x, t).

Para resolver o problema precisamos calcular apenas os valores de kn . Assim, temos o seguinte:

µ ¶ µ ¶Z L µ ¶
nπa 0 2 nπ
kn = kn = g(x) Sen x dx
L L 0 L

½Z 5 µ ¶ Z 10 µ ¶ ¾
0 2 x nπx (10 − x) nπx
kn = Sen dx + Sen dx
10 0 5 10 5 5 10

½ Z 5 µ ¶ Z 10 µ ¶ Z 10 µ ¶ ¾
0 1 1 nπx nπx 1 nπx
kn = x Sen dx + 2 Sen dx − Sen dx
5 5 0 10 5 10 5 5 10

25
 Anteriormente já foi deduzida a seguinte integral:

Z µ ¶ µ ¶ µ ¶ µ ¶2 µ ¶
nπx 10 nπx 10 nπx
x Sen dx = − x Cos + Sen
10 nπ 10 nπ 10

Usando a relação anterior temos o seguinte:

 "
µ ¶ µ ¶ µ ¶2 µ ¶#5 · µ ¶ µ ¶¸10
0 1 1 10 nπx 10 nπx 10 nπx
kn = − x Cos + Sen +2 − Cos +

5 5 nπ 10 nπ 10 nπ 10 5
0

" 
µ ¶ µ ¶ µ ¶2 µ ¶#10 
1 10 nπx 10 nπx
− − x Cos + Sen
5 nπ 10 nπ 10 
5

( " µ ¶ µ ¶2 µ ¶# · µ ¶¸
01 1 50 nπ 10 nπ 10 10 nπ
kn = − Cos + Sen +2 − Cos (nπ) + Cos +
5 5 nπ 2 nπ 2 nπ nπ 2

" µ ¶2 µ ¶ µ ¶2 µ ¶#)
1 100 10 50 nπ 10 nπ
− − Cos (nπ) + Sen (nπ) + Cos − Sen
5 nπ nπ nπ 2 nπ 2

½ µ ¶ µ ¶¾ µ ¶
0 1 20 nπ 20 nπ 0 8 nπ
kn = 2
Sen + 2
Sen =⇒ kn = 2
Sen
5 (nπ) 2 (nπ) 2 (nπ) 2

Agora encontramos o valor de kn :

µ ¶ µ ¶ µ ¶ µ ¶ µ ¶
nπa 0 8 nπ nπ 8 nπ 80 nπ
kn = kn = Sen =⇒ kn = Sen =⇒ kn = Sen
L (nπ)2 2 10 (nπ)2 2 n3 π 3 2

A solução u(x, t) assume a seguinte forma geral:

∞
X µ ¶ µ ¶
nπ nπa
u(x, t) = kn Sen x Sen t
n=1
L L

Assim, substituı́ndo o valor de kn , assim como os valores de a = 1 e L = 10 temos a seguinte solução do

problema:

µ ¶X
∞ µ ¶ µ ¶ µ ¶ µ ¶
80 1 nπ nπx nπt
u(x, t) = Sen Sen Sen
π3 n=1
n3 2 10 10

26
4.6. A equação de Laplace: Equação do potencial

A equação diferencial parcial conhecida como equação de Laplace assume a seguinte forma:

∂ 2 u(x, y) ∂ 2 u(x, y)
+ = 0 =⇒ uxx + uyy = 0 (4.85)
∂x2 ∂y 2

A equação de Laplace pode ainda ser analisada na forma tridimensional:

uxx + uyy + uzz = 0

A equação diferencial de Laplace é um tipo de equação diferencial parcial que aparece em várias aplicações
tais como:

No problema de condução de calor com duas dimensões espaciais e quando a função de temperatura
u(x, y, t) chega a um estado estacionário, isto é, a temperatura não varia mais com o tempo.

Na análise de campos eletrostáticos quando pretendemos analisar o potencial em um meio dielétrico

sem cargas elétricas.

Na análise da função potencial de uma partı́cula livre no espaço sob a ação de forças gravitacionais.

Observações: As seguintes observações são muito importantes:

1. O problema da equação do potencial não depende do tempo e, portanto, não existem condições iniciais
no problema.

2. Devemos especificar as chamadas condições de contorno. Na análise mais geral realmente devemos
especificar uma condição de contorno para cada ponto da fronteira da função u(x, y) da região
onde procuramos uma solução para u(x, y).

3. Nesse contexto, a condição de contorno mais comum ocorre quando é especificado o valor de u
em cada ponto da fronteira. No problema de condução de calor essa exigência equivale a especificar a
temperatura nos extremos da barra metálica. Entretanto, a condição de contorno pode ser especificada
de outras formas tais como a taxa de variação de u (a derivada ou taxa de variação de u na direção
normal à fronteira) ou uma combinação de ambos.

4. Quando a condição de contorno em cada ponto da fronteira é dada especificando os valores de u então
o problema resultante é chamado de problema de Dirichlet. Por outro lado, se são dados os valores da
derivada normal na fronteira então o problema resultante é chamado de problema de Neumann.

4.6.1. O problema de Dirichlet em um retângulo

Neste caso devemos conhecer a condição de contorno em cada ponto da fronteira do retângulo especifi-
cando o valor de u(x, y) em cada ponto dessa fronteira. Pretendemos encontrar u(x, y) que satisfaz a equação
de Laplace:

uxx + uyy = 0 (4.86)

27
no retângulo 0 < x < a, 0 < y < b e as condições de contorno:

u(x, 0) = 0 u(x, b) = 0 0<x<a

(4.87)
u(0, y) = 0 u(a, y) = f (y) 0≤y≤b

onde f (y) é uma função conhecida e dada para o intervalo 0 ≤ y ≤ b. A figura 4.4 mostra o problema em
análise.

y
6
u(x, b) = 0
⇓
b

u(0, y) = 0 uxx + uyy = 0 u(a, y) = f (y)

−→ ←−

-
⇑ a x
u(x, 0) = 0

Figura 4.4: O problema de Dirichlet em um retângulo.

Neste caso também usamos o método de separação de variáveis para resolver a equação diferencial parcial.
Assim, devemos construir um conjunto fundamental de soluções satisfazendo a equação diferencial parcial
e as condições de contorno homogêneas. Depois fazemos a superposição dessas soluções para satisfazer a
condição de contorno restante. Assim temos o seguinte:

u(x, y) = X(x) Y (y) =⇒

0 00
ux = X (x) Y (y) =⇒ uxx = X (x) Y (y)
0 00
uy = X(x) Y (y) =⇒ uyy = X(x) Y (y)

Substituı́ndo a relação anterior na relação (4.86) temos o seguinte:

00 00
00 00 X (x) Y (y)
X (x) Y (y) + X(x) Y (y) = 0 =⇒ =− =λ
X(x) Y (y)

onde λ é a constante de separação. Assim, temos as duas equações diferenciais ordinárias:

28
 00
X (x) − λX(x) = 0 (4.88)

00
Y (y) − λY (y) = 0 (4.89)

As condições de contorno homogêneas são as seguintes:

u(x, 0) = 0 =⇒ X(x) Y (0) = 0 =⇒ Y (0) = 0

(4.90)
u(x, b) = 0 =⇒ X(x) Y (b) = 0 =⇒ Y (b) = 0

u(0, y) = 0 =⇒ X(0) Y (y) = 0 =⇒ X(0) = 0 (4.91)

A continuação devemos encontrar as soluções para Y (y) e X(x):

Solução de Y (y):
De(4.89) e as condições de contorno (4.90) verificamos que a solução de Y (y) assume a seguinte forma:
µ ¶2 µ ¶
nπ nπ
λn = =⇒ Yn (y) = Sen (4.92)
b b
Solução de X(x):
Da solução de Y (y) e de (4.88) temos o seguinte:
A equação caracterı́stica de (4.88) assume a seguinte forma:
µ ¶2 µ ¶
2 2 nπ nπ
r − λ = 0 =⇒ r = =⇒ r = ±
b b

Portanto, a forma matemática de Xn (x) assume a seguinte forma:

1 nπx 1 nπx
Xn (x) = (k1 + k2 )e b + (k1 − k2 )e− b
2 2

· ³ ´¸ · ³ ´¸
1 nπx
− nπx 1 nπx
− nπx
Xn (x) = k1 e b +e b + k2 e b −e b
2 2
µ ¶ µ ¶
nπx nπx
Xn (x) = k1 Cosh + k2 Senh (4.93)
b b
Na relação anterior e usando (4.91) temos o seguinte:

X(0) = 0 =⇒ Xn (0) = k1 =⇒ k1 = 0

Assim uma solução fundamental para X(x) assume a seguinte forma:

µ ¶
nπx
Xn (x) = Senh (4.94)
b

29
 Portanto, uma solução fundamental un (x, y) assume a seguinte forma:

µ ¶ µ ¶
nπx nπy
un (x, y) = Senh Sen n = 1, 2, 3, . . . (4.95)
b b

Usando a mesma estratégia usada para as outras equações diferenciais parciais encontramos uma solução
geral para o problema de Dirichlet em um retângulo. Assim, para satisfazer a condição de contorno não
homogênea em x = a representamos a solução u(x, y) como uma combinação linear das soluções fundamentais
da seguinte forma:

∞
X ∞
X µ ¶ µ ¶
nπx nπy
u(x, y) = cn un (x, y) = cn Senh Sen (4.96)
n=1 n=1
b b

em que os coeficientes cn são obtidos pela condição de contorno:

∞
X µ ¶ µ ¶
nπa nπy
u(a, y) = f (y) = cn Senh Sen (4.97)
n=1
b b
¡ ¢
Portanto, os coeficientes cn Senh nπab devem ser os coeficientes da série de Fourier em senos de f (y) de
perı́odo 2b e são dados pela seguinte relação:

µ ¶ Z b µ ¶
nπa 2 nπy
cn Senh = f (y) Sen dy (4.98)
b b 0 b

Portanto, a solução do problema de Dirichlet em um retângulo é representada por u(x, y) apresentada

em (4.96) e os coeficientes cn são encontradas usando a relação (4.98).
Exemplo 9: Aplicação da equação de Dirichlet em um retângulo.
Resolver a equação de Laplace para os seguintes dados: a = 3, b = 2 e f (y) dada pela relação:


 y 0≤y≤1
f (y) =

 (2 − y) 1≤y≤2

A solução geral tem a seguinte estrutura:

∞
X µ ¶ µ ¶ ∞
X µ ¶ µ ¶
nπx nπy nπx nπy
u(x, y) = cn Senh Sen = cn Senh Sen
n=1
b b n=1
2 2

Encontramos os valores de cn da seguinte forma:

µ ¶ Z b µ ¶
3nπ 0 2 nπy
cn Senh = cn = f (y) Sen dy
2 2 0 2

0
Assim, encontramos inicialmente cn :

Z 1 µ ¶ Z 2 µ ¶
0 nπy nπy
cn = y Sen dy + (2 − y) Sen dy
0 2 1 2

30
 Z 1 µ ¶ Z 2 µ ¶ Z 2 µ ¶
0 nπy nπy nπy
cn = y Sen dy + 2 Sen dy − y Sen dy
0 2 1 2 1 2

Para resolver a integral anterior devemos usar uma relação já analisada anteriormente e que assume a
seguinte forma:

Z µ ¶ µ ¶ µ ¶ µ ¶2 µ ¶
nπy 2 nπy 2 nπy
y Sen dy = − y Cos + Sen
2 nπ 2 nπ 2

0
Usando a relação anterior encontramos a forma matemática de cn :

" µ ¶ µ ¶ µ ¶2 µ ¶#1 · µ ¶ µ ¶¸2

0 2 nπy 2 nπy 2 nπy
cn = − y Cos + Sen +2 − Cos +
nπ 2 nπ 2 0
nπ 2 1

" µ ¶ µ ¶ µ ¶2 µ ¶ #2
2 nπy 2 nπy
− − y Cos + Sen
nπ 2 nπ 2 1

" µ ¶ µ ¶ µ ¶2 µ ¶ # · µ ¶ µ ¶ µ ¶¸
0 2 nπ 2 nπ 2 2 nπ
cn = − Cos + Sen +0−0 +2 − Cos(nπ) + Cos +
nπ 2 nπ 2 nπ nπ 2

" µ ¶ µ ¶2 µ ¶ µ ¶ µ ¶2 µ ¶#
4 2 2 nπ 2 nπ
− − Cos(nπ) + Sen(nπ) + Cos − Sen
nπ nπ nπ 2 nπ 2

µ ¶ µ ¶
0 8 nπ
cn = Sen
n2 π 2 2

0
Com o valor de cn conhecido encontramos o valor de cn :

µ ¶ µ ¶ µ ¶ µ ¶ µ ¶ ¡ nπ ¢
0 3πn 8 nπ 3πn 8 Sen
cn = cn Senh =⇒ Sen = cn Senh =⇒ cn = ³2 ´
2 2
n π 2 2 2 n π2
2 3πn
Senh 2

Finalmente, encontramos a solução geral:

 ¡ ¢ 
µ ¶X
∞ µ ¶ µ ¶ µ ¶
8 1  Sen nπ nπx nπy
u(x, y) = ³ 2 ´  Senh Sen
π2 n=1
n2
Senh 3nπ 2 2
2

31
4.7. Problemas propostos

1. Nos seguintes problemas determine se o método de separação de variáveis pode ser usado para substituir
a equação diferencial parcial dada por um par de equações diferenciais ordinárias. Nesse caso, encontre
as equações:

a) x uxx + ut = 0
b) t uxx + x ut = 0
c) uxx + uxt + ut = 0
d ) [p(x)ux ]x − r(x)utt = 0
e) uxx + (x + y)uyy = 0
f ) uxx + uyy + xu = 0

2. Considere a condução de calor em uma barra com 40 cm. de comprimento cujas extremidades são
mantidas à temperatura 0o C para todo t > 0. Nos seguintes problemas, encontre uma expressão para
a temperatura u(x, t) se a distribuição de temperatura inicial na barra é a função dada. Suponha que
α2 = 1.

a) u(x, 0) = 50, 0 < x < 40

(
x 0 ≤ x < 20
b) u(x, 0) =
(40 − x) 20 ≤ x ≤ 40


 0
0 ≤ x < 10
c) u(x, 0) = 50 10 ≤ x ≤ 30

 0 30 < x ≤ 40
d ) u(x, 0) = x, 0 < x < 40

3. Considere uma corda elástica de comprimento L = 10 cm. cujas extremidades são mantidas fixas. A
corda é colocada em movimento sem velocidade inicial de uma posição inicial u(x, 0) = f (x).
Nos seguintes problemas, com a = 1 encontre o deslocamento u(x, t) para a posição inicial dada:
(
x/5 0≤x<5
a) f (x) =
(10 − x)/5 5 ≤ x ≤ 10


 2x/5 0 ≤ x < 2,5
b) f (x) = 1 2,5 < x < 7,5

 2(10 − x)/5 7,5 ≤ x ≤ 10

c) f (x) = x(10 − x)2 /125

(
1 4≤x<6
d ) f (x) =
0 5caso contrário

32
4. Considere uma corda elástica de comprimento L = 10 cm. cujas extremidades são mantidas fixas. A
corda é colocada em movimento a partir de sua posição de equilı́brio com velocidade inicial
ut (x, 0) = g(x). Nos seguintes problemas, com a = 1 encontre o deslocamento u(x, t) para o g(x) dado:
(
x/5 0≤x<5
a) g(x) =
(10 − x)/5 5 ≤ x ≤ 10


 2x/5 0 ≤ x < 2,5
b) g(x) = 1 2,5 < x < 7,5

 2(10 − x)/5 7,5 ≤ x ≤ 10

c) g(x) = x(10 − x)2 /125

(
1 4≤x<6
d ) g(x) =
0 5caso contrário

5. (a) Encontre a solução u(x, y) da equação de Laplace no retângulo 0 < x < a, 0 < y < b, que satisfaz
as condições de contorno:
u(0, y) = 0, u(a, y) = 0, 0<y<b
u(x, 0) = 0, u(x, b) = g(x), 0≤x≤a
(b) Encontre a solução se:
(
x 0 ≤ x ≤ a/2
g(x) =
(a − x) a/2 ≤ x ≤ a

6. Encontre a solução u(x, y) da equação de Laplace no retângulo 0 < x < a, 0 < y < b, que satisfaz as
condições de contorno:
u(0, y) = 0, u(a, y) = 0, 0<y<b
u(x, 0) = h(x), u(x, b) = 0, 0≤x≤a

7. Encontre a solução u(x, y) da equação de Laplace no retângulo 0 < x < a, 0 < y < b, que satisfaz as
condições de contorno:
u(0, y) = 0, u(a, y) = f (y), 0<y<b
u(x, 0) = h(x), u(x, b) = 0, 0≤x≤a

33