Enunciado

12. 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 y x ,   0 = 0 ,         u x ,   b = g x ,           0 ≤ x ≤ a .

Note que esse não é um problema de Dirichlet nem de Neumann, mas um problema misto no qual u é dada em parte da fronteira e sua derivada normal é dada no resto.

B - Encontre a soluções se

g x = x ,             0 ≤ x ≤ a / 2 a - x ,           a 2 ≤ x ≤ a

Passo 1

L e t r a   A :

Beleza! Vamos fazer, como sempre, a separação de variáveis.

Procuramos uma solução do tipo u x , y = X x . Y ( y ) . Como essa solução deve satisfazer a Equação de Laplace, temos:

u x x + u y y = 0

X ' ' x Y ( y ) + X ( x ) Y ' ' y = 0

Agora nós passamos tudo que tem x para um lado e o que tem y para o outro:

X ' ' x X x = - Y ' ' y Y y = λ

Como essa igualdade vale para todo x e para todo t , essas razões devem ser iguais a uma constante. Separando em duas equações, temos:

X ' ' - λ X = 0

Y ' ' + λ Y = 0

Aplicando a separação de variáveis para as condições homogêneas do problema (a função nós deixamos por último):

u a , y = 0         →         X ( a ) Y y ⏟ ≠ 0 = 0       →       X a = 0

u 0 , y = 0     →         X ( 0 ) Y ( y ) ⏟ ≠ 0 = 0       →       X 0 = 0

Para a outra condição...

u y x ,   0 = 0

Teremos...

u y x ,   0 = 0     →         X ( x ) ⏟ ≠ 0 Y ' ( 0 ) = 0       →       Y ' ( 0 ) = 0

Então, temos o seguinte:

X ' ' - λ X = 0 X 0 = X a = 0 Y ' ' + λ Y = 0 Y ' 0 = 0

Passo 2

Resolver o problema que tem duas condições (de contorno)

Como a função X ( x ) é a que possui duas condições de contorno, vamos resolver primeiro o problema em x , temos:

X ' ' - λ X = 0 X 0 = X a = 0

Temos a seguinte equação característica:

r 2 - λ = 0

r = ± λ

Vamos ver as soluções caso a caso agora.

Caso 1: ∆ > 0 . Ou seja, o que está dentro da raiz é positivo. Nesse caso, λ > 0 . As soluções são do tipo:

X x = c 1 e λ x + c 2 e - λ x

Agora vamos usar as condições de contorno: X 0 = 0 X a = 0

X 0 = c 1 + c 2 = 0

X a = c 1 e a λ + c 2 e - a λ = 0

Isso nos dá

c 2 = - c 1

c 1 e a λ - c 1 e - a λ = 0

c 1 ( e a λ - e - a λ ⏟ ≠ 0 ) = 0

c 1 = 0

c 2 = - c 1 = 0

Substituindo esses valores na nossa solução geral, temos

X x = 0

Mas isso não nos ajuda em nada! Essa é a solução trivial, que não nos interessa. Portanto, esse caso, ∆ > 0 , não nos dá soluções para a EDO.

Caso 2: ∆ = 0 . Nesse caso, λ = 0 . As soluções são do tipo:

X x = c 1 + c 2 x

Agora vamos usar as condições de contorno: X 0 = 0 X a = 0

X 0 = c 1 = 0

X a = c 1 + c 2 a = 0

Isso nos dá

c 2 a = 0

c 2 = 0

Substituindo esses valores na nossa solução geral, temos

X x = 0

Novamente, a solução trivial, que não nos ajuda em nada. Beleza, vamos ver o 3º caso agora!

Caso 3: ∆ < 0 . Ou seja, o que está dentro da raiz é negativo. Nesse caso, λ < 0 . As soluções são do tipo:

X x = c 1 e α x cos ⁡ β x + c 2 e α x sen ⁡ β x

Substituindo r = ± λ = 0 ± i - λ

X x = c 1 cos ⁡ - λ x + c 2 sen ⁡ - λ x

Agora vamos usar as condições de contorno: X 0 = 0 X a = 0

X 0 = c 1 cos ⁡ 0 + c 2 sen ⁡ 0 = 0

c 1 = 0

X a = c 1 cos ⁡ a - λ + c 2 sen ⁡ a - λ = 0

Isso nos dá

c 2 sen ⁡ a λ - = 0

Agora vamos com calma! Se tivermos c 2 = 0 , a solução é a trivial X x = 0 , que não nos interessa. Então, para essa igualdade ai em cima ser verdadeira, precisamos que o seno seja zero.

sen ⁡ a - λ = 0

O seno se anula em π , 2 π , 3 π ... n π , múltiplos inteiros de π . Logo, podemos dizer que:

a - λ = n π             n = 1,2 , 3 …

- λ = n π a

λ n = - n π a 2

As soluções do 3º caso são então do tipo:

X n x = sen ⁡ n π x a                 n = 1,2 , 3 …

Passo 3

Tendo os valores de λ n , resolver a outra EDO

Beleza, agora temos que resolver o problema em y , temos:

Y ' ' + λ Y = 0 Y ' 0 = 0

A equação característica dessa EDO é a seguinte:

r 2 + λ = 0

r = ± - λ

Já vimos que λ < 0 , não é? Quando encontramos seu valor no passo anterior. Então, as soluções dessa EDO são do tipo:

Y y = c 1 e r 1 y + c 2 e r 2 y

Não está errado escrever assim, mas é usual representar a solução dessa EDO (que é a combinação de exponenciais) dessa forma:

Y y = c 1 cosh ⁡ ( r y ) + c 2 s e n h ( r y )

Então, as soluções da EDO em x são:

Y y = c 1 cosh ⁡ n π y a + c 2 s e n h n π y a

Aplicando a condição Y ' 0 = 0 (portando derivando primeiro):

Y ' y = c 1 n π a sinh ⁡ n π y a + c 2 n π a cosh ⁡ n π y a   →

Y ' 0 = c 1 n π a s e n h 0 ⏟ 0 + c 2 n π a cosh ⁡ 0 ⏟ 1   →

Y ' 0 = c 2 n π a = 0   →

c 2 = 0

Chegamos então a:

Y n y = cosh ⁡ n π y a

Então, podemos escrever o seguinte:

u n x , y = s e n n π x a cosh ⁡ n π y a

Passo 4

Escrever u ( x , y ) como série e achar os coeficientes da série

Como as combinações lineares de u n também são soluções do problema, temos que:

u x , y = ∑ n = 1 ∞ b n s e n n π x a cosh ⁡ n π y a  

Agora vamos usar a condição que faltou:   u x , b = g ( x ) .

u x , b = ∑ n = 1 ∞ b n cosh ⁡ n π b a ⏟ c o e f . d a   s é r i e   s e n n π x a = g ( x )

Repare que o coeficiente dos termos dessa série é tudo o que está no colchete! O que depende apenas de n (não de x ). Como essa é uma Série de Fourier de senos, fazendo uma extensão ímpar de g ( x ) , podemos dizer que:

b n cosh ⁡ n π b a = 2 a ∫ 0 a g x s e n   n π x a d x

Repare que o período da extensão é 2 a , não 2 L . Por isso apareceram esses a ’s na integral, tome cuidado. O a veio do denominador do seno.

b n = 2 a   cosh ⁡ n π b a ∫ 0 a g x s e n   n π x a d x

Dessa forma, encontramos os coeficientes b n da série e substituímos em

u x , y = ∑ n = 1 ∞ b n s e n n π x a cosh ⁡ n π b a  

Resposta

u x , y = ∑ n = 1 ∞ b n s e n n π x a cosh ⁡ n π b a  

Onde

b n = 2 a   cosh ⁡ n π b a ∫ 0 a g x s e n   n π x a d x

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

2. Encontre a solução de u x ,y da equação de Laplace no retângulo 0xa , 0yb , que satisfaz as condições de contorno u 0 ,y = 0 , u a ,y = 0 , 0yb , u x , 0 = h x , u x ,b = 0 , 0 ≤ x ≤ a .
Use o problema 43 para mostrar que P n 1 = 1eP n - 1 = - 1 n
Use explicitamente os polinômios de Legendre P0(x), P1(x), P2(x) e P3(x) para calcular ∫ - 1 1 P n 2 x d x ,n = 0 ,1 ,2e3 . Generalize o resultado.
Use a relação de recorrência e P0(x)=1, P1(x)=x para os próximos cinco polinômios de Legendre.
a) Use as equações explícitas y1(x) e y2(x) da equação de Legendre e escolhas apropriadas de c0 e c1 para encontrar os polinômios de Legendre P6 (x) e P7(x). b) Escreva as equações diferenciais explíc