Numerical techniques for simulating a fractional mathematical model of epidermal wound healing

A number of mathematical models investigating certain aspects of the complicated process of wound healing are reported in the literature in recent years. However, effective numerical methods and supporting error analysis for the fractional equations which describe the process of wound healing are still limited.In this paper, we consider the numerical simulation of a fractional mathematical model of epidermal wound healing (FMM-EWH), which is based on the coupled advection-diffusion equations for cell and chemical concentration in a polar coordinate system. The space fractional derivatives are defined in the Left and Right Riemann-Liouville sense. Fractional orders in the advection and diffusion terms belong to the intervals (0,1) or (1,2], respectively. Some numerical techniques will be used. Firstly, the coupled advection-diffusion equations are decoupled to a single space fractional advection-diffusion equation in a polar coordinate system. Secondly, we propose a new implicit difference method for simulating this equation by using the equivalent of Riemann-Liouville and Grünwald-Letnikov fractional derivative definitions. Thirdly, its stability and convergence are discussed, respectively. Finally, some numerical results are given to demonstrate the theoretical analysis.

[1]  Fawang Liu,et al.  The fundamental solution and numerical solution of the Riesz fractional advection–dispersion equation , 2008 .

[2]  Fawang Liu,et al.  Numerical methods and analysis for a class of fractional advection-dispersion models , 2012, Comput. Math. Appl..

[3]  Cem Çelik,et al.  Crank-Nicolson method for the fractional diffusion equation with the Riesz fractional derivative , 2012, J. Comput. Phys..

[4]  J. P. Roop.  Computational aspects of FEM approximation of fractional advection dispersion equations on bounded domains in R 2 , 2006 .

[5]  Avner Friedman,et al.  Wound angiogenesis as a function of tissue oxygen tension: A mathematical model , 2008, Proceedings of the National Academy of Sciences.

[6]  Mihály Kovács,et al.  Numerical solutions for fractional reaction-diffusion equations , 2008, Comput. Math. Appl..

[7]  Fawang Liu,et al.  Numerical Methods for the Variable-Order Fractional Advection-Diffusion Equation with a Nonlinear Source Term , 2009, SIAM J. Numer. Anal..

[8]  Yoshiyuki Yamashita,et al.  Hard x-ray photoelectron spectroscopy of buried Heusler compounds , 2009 .

[9]  I. Turner,et al.  Numerical methods for fractional partial differential equations with Riesz space fractional derivatives , 2010 .

[10]  M. Caputo,et al.  The memory formalism in the diffusion of drugs through skin membrane , 2009 .

[11]  Fawang Liu,et al.  Numerical approximations and solution techniques for the space-time Riesz–Caputo fractional advection-diffusion equation , 2011, Numerical Algorithms.

[12]  I. Podlubny.  Fractional differential equations , 1998 .

[13]  Sophia Maggelakis,et al.  A mathematical model of tissue replacement during epidermal wound healing , 2003 .

[14]  M. Meerschaert,et al.  Finite difference approximations for fractional advection-dispersion flow equations , 2004 .

[15]  Fawang Liu,et al.  THE FUNDAMENTAL AND NUMERICAL SOLUTIONS OF THE RIESZ SPACE-FRACTIONAL REACTION–DISPERSION EQUATION , 2008, The ANZIAM Journal.

[16]  Fawang Liu,et al.  Numerical solution of the space fractional Fokker-Planck equation , 2004 .