Ian W. Turner

Megan E. Farquhar, Timothy J. Moroney, Qianqian Yang et al.

Heart failure is one of the most common causes of death in the western world. Many heart problems are linked to disturbances in cardiac electrical activity, such as wave re-entry caused by ischaemia. In terms of mathematical modelling, the monodomain equation is widely used to model electrical activity in the heart. Recently, Bueno-Orovio et al. [J. R. Soc. Interface 11: 20140352, 2014] pioneered the use of a fractional Laplacian operator in the monodomain equation to account for the complex heterogeneous structures in heart tissue. In this work we consider how to extend this approach to apply to hearts with regions of damaged tissue. This requires the use of a fractional Laplacian operator whose fractional order varies spatially. We develop efficient numerical methods capable of solving this challenging problem on domains ranging from simple one-dimensional intervals with uniform meshes, through to full three-dimensional geometries on unstructured meshes. Results are presented for several test problems in one dimension, demonstrating the effects of different fractional orders in regions of healthy and damaged tissue. Then we showcase some new results for a three-dimensional fractional monodomain equation with a Beeler-Reuter ionic current model on a rabbit heart mesh. These simulation results are found to exhibit wave re-entry behaviour, brought about only by varying the value of the fractional order in a region representing damaged tissue.

Kevin Burrage, Pamela M. Burrage, Ian W. Turner et al.

In this paper we study the class of mixed-index time fractional differential equations in which different components of the problem have different time fractional derivatives on the left hand side. We prove a theorem on the solution of the linear system of equations, which collapses to the well-known Mittag-Leffler solution in the case the indices are the same, and also generalises the solution of the so-called linear sequential class of time fractional problems. We also investigate the asymptotic stability properties of this class of problems using Laplace transforms and show how Laplace transforms can be used to write solutions as linear combinations of generalised Mittag-Leffler functions in some cases. Finally we illustrate our results with some numerical simulations.