1204.4754 (Kristopher L. Kuhlman)
Kristopher L. Kuhlman
A boundary element method (BEM) simulation is used to compare the efficiency of numerical inverse Laplace transform strategies, focusing on requirements of Laplace-space numerical approaches in general. The two-dimensional BEM solution is used to solve the Laplace-transformed diffusion equation, producing a time-domain solution after a numerical Laplace transform inversion. Motivated by the needs of numerical methods posed in Laplace-transformed space, we compare five inverse Laplace transform algorithms and discuss implementation techniques to minimize the number of Laplace-space function evaluations. We investigate the ability to calculate a sequence of time domain values using the fewest Laplace-space model evaluations possible. We find Fourier series based inversion algorithms work for common time behaviors, are the most robust with respect to free parameters, and allow for image function evaluation re-use across at least a log cycle of time.
View original:
http://arxiv.org/abs/1204.4754
No comments:
Post a Comment