Time and Space Fractional Diffusion in Finite Systems

Spatiotemporal nonlocal diffusion in a bounded system is addressed by considering fractional diffusion in a linear, composite system. By considering limiting conditions, solutions for combinations of Neumann and Dirichlet boundary conditions (either zero or nonzero) at the ends of a finite system are derived in terms of Mittag–Leffler functions by the Laplace transformation. Computational viability is demonstrated by inverting the solutions numerically and comparing resulting calculations with asymptotic solutions. Time and space fractional derivatives, defined by variables \(\alpha \) and \(\beta \), respectively, are employed in the Caputo sense; a single-sided, asymmetric space derivative is used. Inspection of the asymptotic solutions leads to insights on the structure of the solutions that may not be available otherwise; the resulting deductions are verified through the numerical inversions. For pure superdiffusion, characteristics of some of the solutions presented here are very similar to those of classical diffusion but combined effects for the corresponding situation result in power-law behaviors. Incidentally, to our knowledge, the pressure distribution for space fractional diffusion at long enough times in a finite system is derived based on first principles for the first time.