Trabalho Final de EDO Finalizado PDF

UNIVERSIDADE FEDERAL DO CEARÁ
CENTRO DE TECNOLOGIA – CT
EQUAÇÕES DIFERENCIAIS APLICADAS AS ENGENHARIAS DE ENERGIA E

MEIO AMBIENTE.
PROFESSOR: FRANCISCO MURILO TAVARES DE LUNA
TRABALHO #1
ATIVIDADE 6
EQUIPE: ANDRÊSA DA SILVA OLIVEIRA - 427802
ESTER FERREIRA RODRIGUES - 415219
FRANCISCA PATRÍCIA ROSA MOREIRA - 405696
LINO SOUZA FILHO – 403968
LÍVIA MARIA DA SILVA MELLO - 427850
FORTALEZA – 2020
1
ÍNDICE
1. INTRODUÇÃO....................................................................................................3
2. MÉTODO ANALÍTICO E COMPUTACIONAL
2.1 MÉTODO ANALÍTICO ................................................................................ 4
2.2 EQUAÇÕES PARA RUNGE-KUTTA 2° ORDEM. .................................... 8
2.3 EQUAÇÕES PARA RUNGE-KUTTA 4° ORDEM. .................................... 8
3. CONSIDERAÇÕES FINAIS ................................................................................ 11
4. REFERÊNCIAS BIBLIOGRÁFICAS .................................................................. 12
2
1. INTRODUÇÃO
Alguns problemas de engenharias podem ser modelados através da matemática,

portanto sendo uma representação ideal e simplificada daquele problema. Entretanto,
muitas vezes não é possível justificar algumas hipóteses que mudam a natureza desse
problema, visto que as soluções precisam preservar os aspectos do problema para que se
tenha uma solução exata. Os métodos numéricos para edos, equações diferenciais
ordinárias, são de grande importância pois alguns problemas não possuem solução
analítica. Mas, nesses casos pode-se obter solução numérica para o caso.
O método de Runge-Kutta é um dos métodos mais populares e o de quarta ordem é
um dos mais utilizados para achar soluções aproximadas de valor inicial. No método de
Runge-Kutta, faz-se uma comparação do polinômio de Taylor para eliminar o cálculo das
derivadas, dessa forma é feita várias avaliações d a função a cada passo.
O objetivo desse trabalho é a resolução de problemas de valor inicial, ou seja, a
equação diferencial ordinária possui uma condição inicial conhecida, portanto, através do
método de Runge-Kutta seremos capazes de determinar a concentração da mistura a cada
momento. Ao final, será realizada uma comparação entre os métodos computacionais
testados e o método analítico, baseadas nos erros percentuais encontrados e utilizando
gráficos e tabelas para deixar nosso estudo mais lúdico e se fazer melhor compreendido.
3
2. MÉTODO ANALÍTICO E COMPUTACIONAL
Problema #1:
Uma bebida contendo 6% de álcool por litro é bombeada em um tonel que
inicialmente contém 400 litros de bebida com 3% de álcool. A taxa de bombeamento é
de 3 litros por minuto, enquanto o líquido misturado é drenado a uma taxa de 4 litros por
minuto. Encontre a solução analítica que representa a concentração de álcool na saída do
tanque em qualquer instante. Utilize os Métodos de Runge-Kutta de segunda e quarta
ordem com passo de 0.1, para apresentar a comparação dos perfis de concentração de
álcool em função do tempo, a partir da solução analítica e das soluções numéricas. Calcule
os erros dos métodos numéricos utilizados.
RESULTADOS:
Solução analítica
Considerando A (t) a quantidade presente em qualquer momento, então 𝑑𝐴 é a taxa na
𝑑𝑡
qual a quantidade de álcool que está mudando:
Resolvendo E.D.O:
𝑑𝐴
= 𝐸𝑛𝑡𝑟𝑎𝑑𝑎 − 𝑠𝑎í𝑑𝑎 𝝁(𝒕) = 𝑬−𝟒 (𝒍𝒏(𝟒𝟎𝟎−𝒕))
𝑑𝑡
6𝑙/𝐴𝑙𝑐𝑜𝑜𝑙 3𝑙/𝑆𝑜𝑙𝑢çã𝑜 𝝁(𝒕) = (𝟒𝟎𝟎 − 𝒕)−𝟒
𝑅1 = ( )( )= 0,13 l/min
100 𝑙/𝑠𝑜𝑙𝑢çã𝑜 1 𝑚𝑖𝑛
𝒅[𝟒𝟎𝟎 − 𝒕]−𝟒 𝑨 = 𝟎, 𝟏𝟑 (𝟒𝟎𝟎 − 𝒕)−𝟒
𝐴𝑙/𝐴𝑙𝑐𝑜 𝑜 𝑙 4𝑙/𝑆𝑜𝑙𝑢çã𝑜 4𝐴
𝑅2 = ( )( )= l/min (𝟒𝟎𝟎 − 𝒕)−𝟒 𝑨 = 𝟎, 𝟎𝟔 (𝟒𝟎𝟎 − 𝒕)−𝟑 + 𝑪
400−𝑡𝑙/𝑠𝑜 𝑙𝑢çã𝑜 1 𝑚𝑖𝑛 400−𝑡
𝑑𝐴 = 0,18 −
4𝐴
=> 𝑑𝐴 4𝐴 𝑨 = 𝟎, 𝟎𝟔 (𝟒𝟎𝟎 − 𝒕) + 𝑪 (𝟒𝟎𝟎 − 𝒕)𝟒
+ = 0,18
𝑑𝑡 400−𝑇 𝑑𝑡 400−𝑡
Temos que A(0) = 12, então:
12 = 0,06 (𝟒𝟎𝟎 − 𝟎)+ 𝑪 (𝟒𝟎𝟎 − 𝟎)𝟒
C= -4,6875x𝟏𝟎−𝟏𝟎
Logo, o número de litros de álcool em qualquer momento é:
𝑨(𝒕) = 𝟎, 𝟎𝟔 (𝟒𝟎𝟎 − 𝒕) − 𝟒, 𝟔𝟖𝟕𝟓𝐱𝟏𝟎 −𝟏𝟎 (𝟒𝟎𝟎 − 𝒕)𝟒
Com essa fórmula substituímos o valor do passo (0,1) de 0 a 1, obtivemos 10 valores que
representam a concentração de álcool em função do tempo nesse problema, o cálculo
analítico está abaixo demonstrado.
4
5
6
Utilizando a ferramenta Excel como base para o método numérico de Runge-Kutta
de 2° e 4° ordem e como base de cálculos as fórmulas das figuras 1 e 3 obtivemos os
resultados tabelados na figura 2 e 4.
Nas figuras abaixo (2 e 4), é possível observar os valores de concentração obtidos
quando simulados de acordo com o método acima citado e os resultados corroboram para
que consideremos a equação de 4° ordem como alguns décimos mais eficaz, apesar de as
duas terem mostrado valores similares aos obtidos de forma analítica. Os valores para o
tempo demonstrados nas tabelas, um total de 10 partições foram suficientes para mostrar
que existe uma constância nos valores de concentração encontrados com o auxílio das
fórmulas para ordem 2 e 4. Para comprovar isso de forma súbita calculamos o erro
percentual dos dois métodos numéricos com relação ao método analítico, esses valores
também estão dispostos nas duas tabelas, por fim, para melhor observar o comportamento
da concentração em função do tempo plotamos 2 gráficos para as equações de 2° e 4°
ordem.
7
2.2 Equações para Runge-Kutta 2° ordem:
Figura 01
Fonte: https://www.youtube.com/watch?v=mHAjIW_gVPw
Figura 03 – Fonte autoral

2.3 Equações para Runge-Kutta 4° ordem:
Figura 03
Fonte: https://pt.wikipedia.org/wiki/M%C3%A9todo_de_Runge-Kutta
8
Figura 04 – Fonte autoral
9
• Gráficos da concentração em função do tempo.
2° ORDEM
12,07
12,06
12,05
Concentração
12,04
12,03
12,02
12,01
12
11,99
0 0,2 0,4 0,6 0,8 1 1,2
Tempo
Para a 2° ordem podemos observar um crescimento linear e quase não se percebe

uma variação nos valores de concentração.
4° ordem
12,060
12,050
12,040
concentração
12,030
12,020
12,010
12,000
11,990
0 0,2 0,4 0,6 0,8 1 1,2
tempo
Já para 4° ordem percebemos uma variação bem significativa nos valores de

concentração, no entanto a função permanece crescendo, apesar de não tão linear como a
anterior, essa equação por tanto é mais realista ao demonstrar como realmente se
comporta a concentração em relação ao tempo desse problema.
10
3. CONSIDERAÇÕES FINAIS
No trabalho apresentado, os métodos de Runge-Kutta de segunda e quarta ordem

foram utilizados comparando-os com a solução analítica calculada. Existem situações em
que é preferível um método numérico ao método analítico ainda que este exista, por
exemplo se a solução para um problema envolve muitos cálculos. A maior parte dos
problemas concretos são, em geral, complexos e envolvem fenômenos não lineares, pelo
que é comum de se encontrar numa situação em que os conhecimentos de matemática não
são suficientes para a descoberta de uma solução exata para um problema real.
De acordo com o problema apresentado neste trabalho, nota-se que o método de
Runge-Kutta de quarta ordem foi mais eficaz em relação ao método de Runge-Kutta de
segunda ordem, o que era esperado, pois o método de Runge-Kutta é um método que
apresenta erro de quarta ordem. E comparando os dois métodos computacionais utilizados
neste trabalho com o método analítico, os valores obtidos pelas duas equações mostram-
se bastantes similares ao encontrado pelo método tradicional, no entanto ainda assim os
valores obtidos pelo método computacional Runge-Kutta de quarta ordem traduzem-se
mais precisos para a resolução do problema. Uma evidencia para reforçar essa afirmação,
são os valores dos erros de cada método quando comparados com o método analítico,
visto que, o erro para os dois métodos de Runge-Kutta de segunda e quarta ordem com o
tempo de 0,2 segundos foi de 0,1997 % e 0,1996 % respectivamente, isso reafirma nossa
premissa de que os dois métodos são bons, mas que o Runge-Kutta de quarta ordem se
sobressai ao de segunda.
11
4. REFERÊNCIAS BIBLIOGRÁFICAS
• STERZA, R. L.; BRANDI, A. C. Comparação entre métodos numéricos: Runge-

Kutta de quarta ordem e previsor-corretor. C.Q.D.– Revista Eletrônica Paulista de
Matemática, Bauru, v. 7, p. 12-22, dez. 2016. Edição ERMAC. DOI:
10.21167/cqdvol7ermac201623169664rlsacb1222 - Disponível em:
http://www.fc.unesp.br/#!/departamentos/matematica/revista-cqd/ 2
• https://www.youtube.com/watch?v=mHAjIW_gVPw
Data e hora de acesso: 20/10/2020 às 21:55.
• https://pt.wikipedia.org/wiki/M%C3%A9todo_de_Runge-Kutta
Data e hora de acesso: 20/10/2020 às 22:30.
12

Trabalho Final de EDO Finalizado PDF

Enviado por

Direitos autorais:

Formatos disponíveis

Trabalho Final de EDO Finalizado PDF

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Trabalho Final de EDO Finalizado PDF

Enviado por

Direitos autorais:

Formatos disponíveis

UNIVERSIDADE FEDERAL DO CEARÁ

EQUAÇÕES DIFERENCIAIS APLICADAS AS ENGENHARIAS DE ENERGIA E

PROFESSOR: FRANCISCO MURILO TAVARES DE LUNA

EQUIPE: ANDRÊSA DA SILVA OLIVEIRA - 427802

ESTER FERREIRA RODRIGUES - 415219

FRANCISCA PATRÍCIA ROSA MOREIRA - 405696

LINO SOUZA FILHO – 403968

LÍVIA MARIA DA SILVA MELLO - 427850

3. CONSIDERAÇÕES FINAIS ................................................................................ 11

4. REFERÊNCIAS BIBLIOGRÁFICAS .................................................................. 12

Alguns problemas de engenharias podem ser modelados através da matemática,

Logo, o número de litros de álcool em qualquer momento é:

𝑨(𝒕) = 𝟎, 𝟎𝟔 (𝟒𝟎𝟎 − 𝒕) − 𝟒, 𝟔𝟖𝟕𝟓𝐱𝟏𝟎 −𝟏𝟎 (𝟒𝟎𝟎 − 𝒕)𝟒

Figura 03 – Fonte autoral

Para a 2° ordem podemos observar um crescimento linear e quase não se percebe

Já para 4° ordem percebemos uma variação bem significativa nos valores de

No trabalho apresentado, os métodos de Runge-Kutta de segunda e quarta ordem

• STERZA, R. L.; BRANDI, A. C. Comparação entre métodos numéricos: Runge-

Você também pode gostar