Algoritmo de Stehfest

• Aplicação – Inversão numérica de funções no domínio

de Laplace para o domínio do tempo.

• Algoritmo – F(t) é a função no domínio do tempo

– é a função resultante da aplicação da transformada de Laplace à função F(t)

– F(t0) é o valor aproximado da função F(t) no ponto t=t0 obtido com o algoritmo de Stehfest.

– Onde Vi é dado por:

)(sF

• Condições – O número n tem que ser par e os melhores

resultados são obtidos com 6 < n < 16. Tipicamente o valor n=8 apresenta bons resultados para problemas associados a equação da difusividade hidráulica.

– A somatória dos componente do vetor Vi é sempre igual a zero para qualquer valor de n par. Exemplo para n=10

• LimitaçãoAlgoritmo – O algoritmo de Stehfest não proporciona

bons resultados para funções com grande variações bruscas no domínio do tempo.