Semi-analytical solution of multilayer diffusion problems with time-varying boundary conditions and general interface conditions

Abstract

We develop a new semi-analytical method for solving multilayer diffusion problems with time-varying external boundary conditions and general internal boundary conditions at the interfaces between adjacent layers. The convergence rate of the semi-analytical method, relative to the number of eigenvalues, is investigated and the effect of varying the interface conditions on the solution behaviour is explored. Numerical experiments demonstrate that solutions can be computed using the new semi-analytical method that are more accurate and more efficient than the unified transform method of Sheils [Appl. Math. Model., 46:450-464, 2017]. Furthermore, unlike classical analytical solutions and the unified transform method, only the new semi-analytical method is able to correctly treat problems with both time-varying external boundary conditions and a large number of layers. The paper is concluded by replicating solutions to several important industrial, environmental and biological applications previously reported in the literature, demonstrating the wide applicability of the work.

24 pages, 8 figures, accepted version of paper published in Applied Mathematics and Computation

Cited by in corpus (10)

• New semi-analytical solutions for advection-dispersion equations in multilayer porous media
• Characteristic time scales for diffusion processes through layers and across interfaces
• Calculating how long it takes for a diffusion process to effectively reach steady state without computing the transient solution
• Algorithms for Brownian dynamics across discontinuities
• A fast algorithm for semi-analytically solving the homogenization boundary value problem for block locally-isotropic heterogeneous media
• A Langevin dynamics approach for multi-layer mass transfer problems
• Multilayer heat equations and their solutions via oscillating integral transforms
• Mean exit time in irregularly-shaped annular and composite disc domains
• Multilayer heat equations: application to finance
• Generalized integral transform method for solving multilayer diffusion problems