A numerical inversion algorithm for Laplace transforms that is capable of handling rapid changes in the computed function is applied to the Laplace transform solution to the problem of convergent radial dispersion in a homogeneous aquifer. Prior attempts by the author to invert this solution were unsuccessful for highly advective systems where the Peclet number was relatively large. The algorithm used in this note allows for rapid and accurate inversion of the solution for all Peclet numbers of practical interest, and beyond. Dimensionless breakthrough curves are illustrated for tracer input in the form of a step function, a Dirac impulse, or a rectangular input.
Additional publication details
|Publication Subtype||Journal Article|
|Title||Convergent radial dispersion: A note on evaluation of the Laplace transform solution|
|Series title||Water Resources Research|
|Publisher||American Geophysical Union|
|Contributing office(s)||Toxic Substances Hydrology Program|
|Google Analytic Metrics||Metrics page|