A hybrid analytical-numerical method for solving advection-dispersion problems on a half-line

•Solution of the advection-dispersion equation on the half-line via the Unified Transform (Fokas Method).•Combines complex analysis with numerics.•Compared to traditional approaches, the unified transform avoids the solution of ordinary differential equations.•Results are successfully compared to ot...

Full description

Saved in:
Bibliographic Details
Published in:International journal of heat and mass transfer Vol. 139; pp. 482 - 491
Main Authors: de Barros, F.P.J., Colbrook, M.J., Fokas, A.S.
Format: Journal Article
Language:English
Published: Oxford Elsevier Ltd 01.08.2019
Elsevier BV
Subjects:
Advection
Advection-dispersion equation
Analytical solution
Complexity
Computation
Convergence
Decay rate
Dispersion
Environmental flows
Fokas method
Integrals
Numerical methods
Ordinary differential equations
Partial differential equations
Transformations (mathematics)
Unified transform
Advection-dispersion equation
Unified transform
Fokas method
Analytical solution
Environmental flows
ISSN:0017-9310, 1879-2189
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:•Solution of the advection-dispersion equation on the half-line via the Unified Transform (Fokas Method).•Combines complex analysis with numerics.•Compared to traditional approaches, the unified transform avoids the solution of ordinary differential equations.•Results are successfully compared to other existing solutions.•Illustrates the advantage of the numerical implementation of the Fokas Method. This paper employs the unified transform, also known as the Fokas method, to solve the advection-dispersion equation on the half-line. This method combines complex analysis with numerics. Compared to classical approaches used to solve linear partial differential equations (PDEs), the unified transform avoids the solution of ordinary differential equations and, more importantly, constructs an integral representation of the solution in the complex plane which is uniformly convergent at the boundaries. As a consequence, such solutions are well suited for numerical computations. Indeed, the numerical evaluation of the solution requires only the computation of a single contour integral involving an integrand which decays exponentially fast for large values of the integration variable. A novel contribution of this paper, with respect to the solution of linear evolution PDEs in general, and the implementation of the unified transform in particular, is the following: using the advection-dispersion equation as a generic example, it is shown that if the transforms of the given data can be computed analytically, then the unified transform yields a fast and accurate method that converges exponentially with the number of evaluations N yet only has complexity O(N). Furthermore, if the transforms are computed numerically using M evaluations, the unified transform gives rise to a method with complexity O(NM). Results are successfully compared to other existing solutions.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0017-9310
1879-2189
DOI:10.1016/j.ijheatmasstransfer.2019.05.018