Exercícios Resolvidos de Métodos das Diferenças Finitas

Ver Teoria

Enunciado

O PVC abaixo

y ' ' = 4 y - x ,       0 ≤ x ≤ 1 y 0 = 0 y 1 = 2

tem como solução y = e 2 e 4 - 1 - 1 e 2 x - e - 2 x + x . Use o método das diferenças finitas para aproximar a solução e compare com a solução real para:

a h = 0.5

b h = 0.25

Passo 1

a Temos o seguinte:

y ' ' = 4 y - 4 x

E teremos n intervalos de modo que:

n = b - a h = 1 0.5 = 2

Então vai ficar assim:

Beleza, agora vamos substituir as derivadas na nossa EDO:

y ' ' = 4 y - 4 x

y i + 1 - 2 y i + y i - 1 h 2 = 4 y i - 4 x i

Multiplicando os dois lados por h 2 :

y i + 1 - 2 y i + y i - 1 = 4 h 2 y i - 4 h 2 x i

Juntando os termos iguais:

y i + 1 + y i - 2 - 4 h 2 + y i - 1 = - 4 h 2 x i

y i + 1 - 3 y i + y i - 1 = - x i

i = 1 :

y 2 - 3 y 1 + y 0 = - x 1

- 3 y 1 = - y 2 - y 0 - x 1

y 1 = 2.5 3

y 1 = 0.83333

Completando nossa tabela:

Passo 2

Nessa letra o único valor que estimamos foi y 0.5 ≅ 0.83333 . O valor exato seria:

y = e 2 e 4 - 1 - 1 e 2 x - e - 2 x + x

y 0.5 = 0.824027

Logo, houve um erro de:

E A = 0.83333 - 0.824027 = 0.009306

E R = E A 0.824027 × 100 = 1.13 %

Passo 3

b Agora nosso h mudou para 0.25 , de modo que teremos n intervalos:

n = b - a h = 1 0.25 = 4

Então vai ficar assim:

Pegando nossa EDO pronta para o método lá da letra a :

y i + 1 + y i - 2 - 4 h 2 + y i - 1 = - 0.25 x i

y i + 1 - 2.25 y i + y i - 1 = - 0.25 x i

i = 1 :

y 2 - 2.25 y 1 + y 0 = - 0.25 x 1

y 2 - 2.25 y 1 = - 0.0625

i = 2 :

y 3 - 2.25 y 2 + y 1 = - 0.25 x 2

y 3 - 2.25 y 2 + y 1 = - 0.125

i = 3 :

y 4 - 2.25 y 3 + y 2 = - 0.25 x 3

- 2.25 y 3 + y 2 = - 2.1875

Passo 4

O sistema na forma matricial fica:

- 2.25 1 0 1 - 2.25 1 0 1 - 2.25 y 1 y 2 y 3 = - 0.0625 - 0.125 - 2.1875

Fazendo a Eliminação de Gauss:

- 2.25 1 0 1 - 2.25 1 0 1 - 2.25 - 0.0625 - 0.125 - 2.1875 → L 2 ← L 2 + 0.44444 L 1 - 2.25 1 0 0 - 1.80556 1 0 1 - 2.25 - 0.0625 - 0.15278 - 2.1875 → L 3 ← L 3 + 0.55385 L 2 - 2.25 1 0 0 - 1.80556 1 0 0 - 1.69615 - 0.0625 - 0.15278 - 2.27212

Retrosubstituindo:

y 1 = 0.39512

y 2 = 0.82653

y 3 = 1.33957

Completando a tabelinha já com os valores exatos e com o erro ficamos com:

Resposta

a y 0.5 ≅ 0.83333

E A = 0.009306

E R = 1.13 %

b

Exercícios de Livros Relacionados

Uma parede, com 0,12m de espessura e difusividade térmica de 1,5×10 - 6m 2/ s , encontra-se inicialmente a uma temperatura uniforme igual a 85 ° C . Subitamente, uma de suas faces tem a sua temperatur
O critério de estabilidade para o método explícito exige que o coeficiente do termo T m p na equação de diferenças finitas unidimensional seja zero ou positivo. Analise a situação na qual as temperatu
Um fino bastão com diâmetro D está inicialmente em equilíbrio com a sua vizinhança, um grande recipiente, onde há vácuo, a uma temperatura T v i z . De repente, uma corrente elétrica I ( A ) é passada
Um cilindro sólido redondo, feito com um material plástico ( α=6 × 10 - 7 m 2 / s ) , está inicialmente a uma temperatura uniforme de 20 ° C e se encontra isolado ao longo de sua superfície lateral e
Seja o Problema 5.9, exceto pelo fato de que agora o volume combinado do banho de óleo e da esfera é igual a V=1m 3 . O banho de óleo é bem misturado e isolado termicamente. (a) Considerando que as pr