skip to main content


Title: A fast‐high order compact difference method for the fractional cable equation

The Cable equation is one of the most fundamental equations for modeling neuronal dynamics. In this article, we consider a high order compact finite difference numerical solution for the fractional Cable equation, which is a generalization of the classical Cable equation by taking into account the anomalous diffusion in the movement of the ions in neuronal system. The resulting finite difference scheme is unconditionally stable and converges with the convergence order ofin maximum norm, 1‐norm and 2‐norm. Furthermore, we present a fast solution technique to accelerate Toeplitz matrix‐vector multiplications arising from finite difference discretization. This fast solution technique is based on a fast Fourier transform and depends on the special structure of coefficient matrices, and it helps to reduce the computational work fromrequired by traditional methods towithout using any lossy compression, whereandτis the size of time step,andhis the size of space step. Moreover, we give a compact finite difference scheme and consider its stability analysis for two‐dimensional fractional Cable equation. The applicability and accuracy of the scheme are demonstrated by numerical experiments to support our theoretical analysis.

 
more » « less
PAR ID:
10060870
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Numerical Methods for Partial Differential Equations
Volume:
34
Issue:
6
ISSN:
0749-159X
Page Range / eLocation ID:
p. 2237-2266
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In this paper, the author derives an‐superconvergence for the piecewise linear Ritz‐Galerkin finite element approximations for the second‐order elliptic equationequipped with Dirichlet boundary conditions. This superconvergence error estimate is established between the finite element solution and the usual Lagrange nodal point interpolation of the exact solution, and thus the superconvergence at the nodal points of each element. The result is based on a condition for the finite element partition characterized by the coefficient tensorand the usual shape functions on each element, called‐equilateral assumption in this paper. Several examples are presented for the coefficient tensorand finite element triangulations which satisfy the conditions necessary for superconvergence. Some numerical experiments are conducted to confirm this new theory of superconvergence.

     
    more » « less
  2. Abstract

    We prove that a WLD subspace of the spaceconsisting of all bounded, countably supported functions on a set Γ embeds isomorphically intoif and only if it does not contain isometric copies of. Moreover, a subspace ofis constructed that has an unconditional basis, does not embed into, and whose every weakly compact subset is separable (in particular, it cannot contain any isomorphic copies of).

     
    more » « less
  3. Abstract

    In this paper, we are interested in the following question: given an arbitrary Steiner triple systemonvertices and any 3‐uniform hypertreeonvertices, is it necessary thatcontainsas a subgraph provided? We show the answer is positive for a class of hypertrees and conjecture that the answer is always positive.

     
    more » « less
  4. Abstract

    Letbe the random directed graph onnvertices where each of thepossible arcs is present independently with probabilityp. A celebrated result of Frieze shows that ifthentypically has a directed Hamilton cycle, and this is best possible. In this paper, we obtain a strengthening of this result, showing that under the same condition, the number of directed Hamilton cycles inis typically. We also prove a hitting‐time version of this statement, showing that in the random directed graph process, as soon as every vertex has in‐/out‐degrees at least 1, there are typicallydirected Hamilton cycles.

     
    more » « less
  5. Abstract

    It is proved that for every countable structureand a computable successor ordinal α there is a countable structurewhich is‐least among all countable structuressuch thatis Σ‐definable in the αth jump. We also show that this result does not hold for the limit ordinal. Moreover, we prove that there is no countable structurewith the degree spectrumfor.

     
    more » « less