Prova Metodos Numericos Aplicados 1

Fazer download em pdf ou txt
Fazer download em pdf ou txt
Você está na página 1de 3

FURG Universidade Federal do Rio Grande

IMEF Instituto de Matemtica, Estatstica e Fsica


Prova de Mtodos Numricos Aplicados
Prof. Dr. Wiliam Correa Marques
DATA DE ENTREGA: 31 de Julho de 2014

Nome:

Total

Matrcula:

AVISOS
Questes sem a apresentao de desenvolvimento sero consideradas incorretas.
A data de entrega da prova NO ser alterada. Se o material no for entregue at a data marcada, o
discente ter sua prova cancelada e sua nota ser considerada ZERO (0).
As provas podem ser realizadas individualmente ou em duplas. Os problemas devem ser resolvidos e
os programas utilizados devem ser entregues junto com um breve relato dos resultados.
A resoluo da prova pode ser escrita (letra legvel) ou digitada. A prova deve ser assinada e datada
pelo discente.
Questo 1. Problema de ejeo de fluidos: No projeto de uma estao de tratamento de esgotos
a Equao Diferencial Ordinria (EDO) apresentada abaixo utilizada para estimar o nvel (H(t)) de um
fluido em uma cmara de ejeo como funo do tempo (t):
2 H(t)
+ (19.42)
(77.7)
t2

H(t)
t

2
+ H(t) 60 = 0

Com condies iniciais:

H(0) =

H(0)
=0
t

Utilizar os mtodos de Euler e Runge-Kutta 2 para estimar a variao do nvel de fluido como funo do
tempo (H(t)) no intervalo de tempo t = [0, 5] com h = 0.05.

a) Apresentar a descrio e montagem do sistema de equaes de primeira ordem que sero utilizadas nos
algoritmos.
b) Montar um algoritmo utilizando sintaxe em: Matlab, Python, Fortran, etc. importante ressaltar
que o algoritmo tem que ser construido sem a utilizao de funes especficas destes programas para a
soluo de EDOs.
c) Apresentar uma descrio suscinta na forma de um relato para resultado obtido.

Questo 2. Espcies em competio: Considere que pi (t) indique, respectivamente, as populaes de


trs espcies competindo Si , i = 1, 2, 3. Suponha que tenham as mesmas taxas de crescimento, e a populao
mxima que o habitat pode admitir a mesma para cada espcie (Considerando que ela seja uma unidade).
Suponha tambm que a vantagem competitiva que S1 tem sobre S2 a mesma de S2 sobre S3 e de S3 sobre
S1 . Esta situao modelada pelo sistema de EDOs dado abaixo:

p1 (t)

= p1 (t) (1 p1 (t) ap2 (t) bp3 (t))

p2 (t)
= p2 (t) (1 bp1 (t) p2 (t) ap3 (t))

p3 (t) = p (t) (1 ap (t) bp (t) p (t))


3
1
2
3
t
Sendo: a e b constantes positivas. Para demonstrar a dinmica da populao deste sistema quando a = b =
0.5 use o mtodo de Runge-Kutta 2 para o sistema com h = 0.01. Neste sentido, aproxime as populaes
pi (t) pelo intervalo de tempo t = [0, 10] utilizando cada uma das seguintes condies iniciais:
(a) p1 (0) = 1.0 , p2 (0) = 0.1, p3 (0) = 0.1
(b) p1 (0) = 0.1 , p2 (0) = 1.0, p3 (0) = 0.1
(c) p1 (0) = 0.1 , p2 (0) = 0.1, p3 (0) = 1.0
Com base nos resultados dos itens (a) - (c), decida o que voc acredita que acontecer com essas populaes
quando t +. Monte um algoritmo utilizando sintaxe em: Matlab, Python, Fortran, etc, e apresente
uma descrio suscinta na forma de um relato para os resultados obtidos. O algoritmo deve ser construido
sem a utilizao de funes especficas destes programas para a soluo de EDOs.

Questo 3. Soluo da equao de transferncia de calor com efeito advectivo: Considere a


Equao Diferencial Parcial (EDP) que pode ser utilizada para representar a transferncia de calor (u(x, t))
em uma barra unidimensional de comprimento L com carcatersticas fsicas que podem ser desprezadas de
forma a simplificar a soluo do problema. Esta situao pode ser modelada pela EDP dada abaixo, durante
um tempo total de siulao T , onde a velocidade de adveco e o coeficiente de difusividade so representados
por c e , respectivamente:
u(x, t)
2 u(x, t)
u(x, t)
= c
+
t
x
x2
Com condio inicial:
u(x, 0) = 500K, 0 x L
Com condies de contorno:
u(0, t) = 270K, 0 t T
u(L, t) = 270K, 0 t T

Utilizar o mtodo de diferenas finitas (explcito ou implcito) para resolver o problema considerando os
seguintes parmetros:
4x = 0.1m
4t = 0.0025s
=1

m2
s

T = 0.25s
L = 1m

a) Montar um algoritmo utilizando sintaxe em: Matlab, Python, Fortran, etc. importante ressaltar que o
algoritmo tem que ser construido sem a utilizao de funes especficas destes programas para a soluo de
EDPs.
b) Realizar duas simules considerando c = 2.5

m
m
e c = 2.5 .
s
s

c) Apresentar uma descrio suscinta na forma de um relato para resultado obtido.

Você também pode gostar