A time dependent fractional order diffusion equation with constant diffusivity matrix
This work addresses numerical challenges in modeling composite materials with anisotropic diffusion, but the contribution is incremental as it applies known spectral and backward Euler methods to a specific equation.
The paper develops a numerical approximation scheme for a time-dependent fractional diffusion equation with a constant diffusivity matrix, capturing boundary behavior accurately. Numerical experiments demonstrate the effects of nonhomogeneous domains and support the theoretical error analysis.
Of primary interest in this paper is the numerical approximation of a time dependent fractional, in space, diffusion equation where the domain is assumed to be nonhomogeneous, having different axial diffusion coefficients. This work is motivated from the consideration of composite material which can exhibit different material properties along, and perpendicular to, internal planar structures. Careful attention is paid to accurately capture the boundary behavior of the solution. A spectral approximation scheme is used for the spatial discretization and a backward Euler approximation used for the temporal discretization. Following an error analysis for the approximation scheme, numerical experiments are given to demonstrate the effects of the nonhomogeneous domain and to support the theoretical analysis.