This paper focuses on the numerical modeling of drug diffusion in biological tissues using fractional time-dependent parabolic equations with non-local boundary conditions. The model includes a Caputo fractional derivative to capture the non-local effects and memory inherent in biological processes, such as drug absorption and transport. The theoretical framework of the problem is based on the work of Alhazzani, et al.,which demonstrates the solution’s goodness, existence, and uniqueness. Building on this foundation, we present a robust numerical method designed to deal with the complexity of fractional derivatives and non-local interactions at the boundaries of biological tissues. Numerical simulations reveal how fractal order and non-local boundary conditions affect the drug concentration distribution over time, providing valuable insights into drug delivery dynamics in biological systems. The results underscore the potential of fractal models to accurately represent diffusion processes in heterogeneous and complex biological environments.
Ghezal et al. (Mon,) studied this question.