628
Views
5
CrossRef citations to date
0
Altmetric
Research papers

Numerical solution of the advection-diffusion equation using Laplace transform finite analytical method

Pages 177-188 | Received 20 Jan 2012, Accepted 22 Mar 2012, Published online: 01 May 2012
 

Abstract

In this study, a numerical method has been investigated and developed to solve the one-dimensional advection-diffusion equation to predict the quality of water in rivers. In this method, time variable is eliminated first by the Laplace transformation, and then a finite analytical method is applied in space. Both are based on a local element of the discretized domain in a finite-volume method. Since the Laplace transformation has been used for temporal approximation, an efficient and accurate inverse Laplace transform method of De Hoog 1982 [An improved method for numerical inversion of Laplace transform. SIAM, Journal of Scientific and Statistical Computing, 3 (3), 357–366] is employed to obtain the solution in real time. The proposed method is compared against analytical solutions and two finite-difference methods. The present computations and comparisons show that the proposed method is superior to the finite-difference methods. The results of the proposed method also agree with analytical solutions without numerical oscillation or diffusion. The present method is applied to steady and unsteady flows and it also provides flexibility for uniform and non-uniform grid spacing and for a wide range of Péclet numbers. It takes less computational effort than finite-difference methods.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.