Abstract

A new semi-analytical solution to the advection-dispersion-reaction equation for modelling solute transport in layered porous media is derived using the Laplace transform. Our solution approach involves introducing unknown functions representing the dispersive flux at the interfaces between adjacent layers, allowing the multilayer problem to be solved separately on each layer in the Laplace domain before being numerically inverted back to the time domain. The derived solution is applicable to the most general form of linear advection-dispersion-reaction equation, a finite medium comprising an arbitrary number of layers, continuity of concentration and dispersive flux at the interfaces between adjacent layers and transient boundary conditions of arbitrary type at the inlet and outlet. The derived semi-analytical solution extends and addresses deficiencies of existing analytical solutions in a layered medium, which consider analogous processes such as diffusion or reaction-diffusion only and/or require the solution of complicated nonlinear transcendental equations to evaluate the solution expressions. Code implementing our semi-analytical solution is supplied and applied to a selection of test cases, with the reported results in excellent agreement with a standard numerical solution and other analytical results available in the literature.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.