Aula Teórica
description
Transcript of Aula Teórica
![Page 1: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/1.jpg)
Aula Teórica
Séries de Taylor e resolução numérica da equação de advecção - difusão
![Page 2: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/2.jpg)
Equações que vamos resolver• Conservação da massa:
• Num modelo Hidrodinâmico também a equação de Transporte de Quantidade de Movimento:
)( kkjjj
kjk PFxc
xxcu
tc
)( GravidadePressãoxu
xxuu
tu
j
i
jj
ij
i
![Page 3: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/3.jpg)
Como se resolvem as equações
• Métodos Numéricos:• Diferenças finitas/Volumes finitos• Elementos Finitos/Elementos de fronteira.
• Como se constrói o método das diferenças finitas?
• Série de Taylor:t
in
nnt
i
t
i
t
i
ti
tti t
cnt
tct
tct
tctcc
!....
!3!2 3
33
2
22
![Page 4: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/4.jpg)
O que representa a série de Taylor?t
in
nnt
i
t
i
t
i
ti
tti t
cnt
tct
tct
tctcc
!....
!3!2 3
33
2
22
t1 t1+Δt
Δt
Δc
Outras derivadas Δc 1ª Derivada: Δc/ Δt
t
c
![Page 5: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/5.jpg)
Como usar para calcular as derivadas?t
in
nnt
i
t
i
t
i
ti
tti t
cnt
tct
tct
tctcc
!....
!3!2 3
33
2
22
)(
)( 2
ttcc
tc
ttctcc
ti
tti
t
i
t
i
ti
tti
Método Explícito: A derivada é calculada à esquerda “em t” e tem precisão de 1ª ordem, ou seja, as derivadas que foram ignoradas estão multiplicadas por )( t
Isto significa que o erro do cálculo aumenta quando o passo de tempo aumenta.
![Page 6: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/6.jpg)
Mas poderia ter feito calculado a derivada à direita
tt
in
nntt
i
tt
i
tt
i
tti
ti t
cnt
tct
tct
tctcc
!
....!3!2 3
33
2
22
)(
)( 2
ttcc
tc
ttctcc
ti
tti
tt
i
tt
i
tti
ti
Método Implícito: A derivada é calculada à direita “em t+dt” e tem precisão de 1ª ordem, ou seja, todas as derivadas que foram ignoradas estão multiplicadas por )( tIsto significa que o erro do cálculo aumenta quando o passo de tempo aumenta. Os métodos implícitos e explícitos têm a mesma precisão.
![Page 7: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/7.jpg)
Para calcular a derivada no centro do intervalo teria que calcular os valores nos extremos a partir daquele
2/2/
3
332/
2
222/2/
!2/....
!32/
!22/2/
tt
in
nntt
i
tt
i
tt
i
tti
tti t
cnt
tct
tct
tctcc
Subtraindo uma da outra:
2/2/
3
332/
2
222/2/
!2/....
!32/
!22/2/
tt
in
nntt
i
tt
i
tt
i
tti
ti t
cnt
tct
tct
tctcc
22/
32/
2/
2/
ttcc
tc
ttctcc
ti
tti
tt
i
tt
i
ti
tti
Neste método a derivada é calculada no centro do intervalo de tempo e tem precisão de 2ª ordem. Dá a solução exacta até uma evolução parabólica. As derivadas ignoradas estão multiplicadas por 22/t
![Page 8: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/8.jpg)
O que representa a série de Taylor?t
in
nnt
i
t
i
t
i
ti
tti t
cnt
tct
tct
tctcc
!....
!3!2 3
33
2
22
t
c
t1 t1+Δt
Δt
Δc
Outras derivadas
1ª Derivada: Δc/ Δt
Método ImplícitoMétodo Explícito
Método Diferenças Centrais
![Page 9: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/9.jpg)
Derivadas espaciaist
in
nnt
i
t
i
t
i
ti
txi x
cnx
xcx
xcx
xcxcc
!....
!3!2 3
33
2
22
)(
)( 2
xxcc
xc
xtcxcc
ti
txi
t
i
t
i
ti
txi
Derivada à direita, Método downwind, se velocidade positiva
Neste método a derivada espacial num ponto é calculada a partir da informação no ponto e da informação à direita.
Veremos mais adiante que este cálculo cria problemas se esta derivada for usada para calcular o termo advectivo quando a velocidade é positiva.
![Page 10: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/10.jpg)
Derivadas espaciais
Derivada à esquerda, Método upwind se velocidade positiva.
t
in
nnt
i
t
i
t
i
ti
txi x
cnx
xcx
xcx
xcxcc
!....
!3!2 3
33
2
22
)(
)( 2
xxcc
xc
xtcxcc
txi
ti
t
i
t
i
ti
txi
Neste método a derivada espacial num ponto é calculada a partir da informação no ponto e da informação à esquerda.
Veremos mais adiante que este cálculo cria problemas se esta derivada for usada para calcular o termo advectivo quando a velocidade é positiva.
![Page 11: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/11.jpg)
Subtraindo uma equação da outra
)(2
)(2
2
3
xxcc
xc
xtcxcc
txi
txi
t
i
t
i
txi
txi
Diferenças Centrais
t
in
nnt
i
t
i
t
i
ti
txi x
cnx
xcx
xcx
xcxcc
!....
!3!2 3
33
2
22
t
in
nnt
i
t
i
t
i
ti
txi x
cnx
xcx
xcx
xcxcc
!....
!3!2 3
33
2
22
![Page 12: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/12.jpg)
2ª Derivada**
3
33*
2
22***
!....
!3!2i
n
nn
iiiixi x
cnx
xcx
xcx
xcxcc
**
3
33*
2
22***
!....
!3!2 in
nn
iiiixi x
cnx
xcx
xcx
xcxcc
)(2
)(2
22
****
2
2
4*
2
22***
xx
cccxc
xtcxccc
xiixi
i
iixixi
Adicionando:
![Page 13: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/13.jpg)
Equações Algébricas• Obtêm-se substituindo as derivadas pelas
aproximações:
• Explícito, diferenças centrais. Precisão de 2ª ordem no espaço e 1ª no tempo.
• Semi-implícito (Crank-Nicholson) diferenças centrais espaço. Precisão de 2ª ordem no tempo e no espaço.
222 2
2x
xcccx
xccut
tcc t
xxtx
txx
txx
txx
txx
tt
22
2/2/2/2
2/2/2 2
2x
xcccx
xccut
tcc tt
xxtt
xttxx
ttxx
ttxx
txx
tt
O que se paga pela precisão de 2ª ordem no tempo?
![Page 14: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/14.jpg)
Como se obtém o valor em (t+Δt/2) ?Fazendo a média…..
• Adicionando as equações!
2/2/
3
332/
2
222/2/
!2/....
!32/
!22/2/
tt
in
nntt
i
tt
i
tt
i
tti
tti t
cnt
tct
tct
tctcc
2/2/
3
332/
2
222/2/
!2/....
!32/
!22/2/
tt
in
nntt
i
tt
i
tt
i
tti
ti t
cnt
tct
tct
tctcc
22/
2/
2
222/
2/2
.....2/2
tccc
tctccc
tti
titt
i
tt
i
tti
ti
tti
• Substituindo estes termos nas equações obtém-se a equação a resolver
![Page 15: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/15.jpg)
Explícito Upwind
• Precisão de 1ª ordem no tempo e no espaço para advecção. Segunda ordem para difusão.
• Esta equação pode ser organizada na forma:
222 2 x
xcccx
xccut
tcc t
xxtx
txx
txx
tx
tx
ttx
)(1 11 PFcfcecdc tii
tii
tii
tti
ti
ti
ti
tti c
xtc
xt
xtuc
xt
xtuc 12212
21
![Page 16: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/16.jpg)
Forma geral da Equação )(11111 1111 PFcfkcekcdkckfckeckd t
iitii
tii
ttii
ttii
ttii
Explicito, upwind:
xtuCr
2ºxtDifN
Números de Courant e de Difusão
ti
ti
ti
tti c
xtc
xt
xtuc
xt
xtuc 12212
21
K=1=> implícito. K=0 => Explicito, k=0.5=> Crank-Nicholson:
![Page 17: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/17.jpg)
Sobre a precisão do cálculo• No cálculo implícito e no cálculo explícito as derivadas são calculadas
nos extremos do intervalo de tempo. Estes métodos ignoram todas as derivadas a partir da primeira: têm precisão de primeira ordem ou “até à primeira ordem”.
• Os termos da série de Taylor ignorados estão multiplicados por• Quando a derivada é calculada no centro do intervalo de tempo as
derivadas só são ignoradas a partir da segunda. São métodos com precisão de 2ª ordem, ou “até à 2ª ordem”. Se a função for uma recta ou uma parábola o cálculo da derivada é exacto.
• Os termos da série de Taylor ignorados estão multiplicados por• Mas >1 então quanto maior é a ordem de precisão do cálculo,
maior é o coeficiente dos termos ignorados. Porque é que a precisão do cálculo aumenta?
)( t
22/t)( t
![Page 18: Aula Teórica](https://reader036.fdocumentos.com/reader036/viewer/2022082521/56814ac9550346895db7df3a/html5/thumbnails/18.jpg)
Porque aumenta a precisão com o expoente de ? )( t
2/1
!2/
tt
in
nn
tc
nt
Porque os termos ignorados são da forma:
O cálculo da derivada faz aparecer em denominador o intervalo de tempo elevado n e o coeficiente está elevado a (n-1) e por isso o produto é proporcional a ou seja à primeira derivada multiplicada pelo inverso do factorial de n e por isso quanto maior é o valor do expoente do intervalo de tempo, menor é o valor dos temos desprezados.
Esta conclusão é consistente como facto de as derivadas perderem importância à medida que a ordem aumenta.
)/()( tc