Abstract
This study, considers the fractional order cable model (FCM) in the sense of Riemann–Liouville fractional derivatives (R-LFD). We use a modified implicit finite difference approximation to solve the FCM numerically. The Fourier series approach is used to examine the proposed scheme’s theoretical analysis, including stability and convergence. The scheme is shown to be unconditionally stable, and the approximate solution converges to the exact solution. To demonstrate the application and feasibility of the proposed approach, a numerical example is provided.
1 Introduction
Real-life phenomena have been modeled in a variety of ways, and partial differential equations (PDEs) and ordinary differential equations (ODEs) can be used to model some of these phenomena. For the phenomena that are not sufficiently modeled by PDEs, fractional PDEs have been developed by replacing the non-integers order derivative [1]. Fractional calculus can be applied to every field of science, such as biology, engineering, image processing, wave propagation, rheology, viscoelasticity, etc.
Fractional diffusion equations are a type of fractional differential equation that has sparked a lot of interest due to their various applications. By adding a variable lower limit of integration Rajkovic [2] generalized the notions of fractional the g-integral and g-derivative, and came up with a q-Taylor definition that contains fractional-order q-derivatives of the function. Yakar [3] considered a fractional boundary value problem with a two-part operator. The main problem’s eigenvalues with Eigen functions are the same as the constructed operator’s eigenvalues and corresponding Eigen functions in Hilbert spaces. The non-integers order Cable model is derived from the circuit model based on intracellular and extracellular space [4]. Vitali [5] introduced a Caputo formula as an extension of the FCE, obtained the analytical solution using the Laplace transform, and obtained results in terms of special functions. Yu [6] used the compact difference method and the Fourier method for stability and convergence in his computational treatment of the two-dimensional FCE. Liu et al. [7] discussed the FCE having two fractional temporal derivatives, and proposed implicit schemes with convergence orders of and , respectively. The energy approach was used to investigate the stability and convergence analysis. Lin et al. [8] devised a numerical schema for FCE discretization. They analyzed the proposed schema by providing theoretical and error estimates. The schema was unconditionally stable. Liu et al. [9] used a two-grid approach with the finite element scheme to solve the non-linear FCM, and the stability based on the fully discrete two-grid method was derived. A semi-discrete approach was used for time, and the Galerkin finite element approach was used for space Zhuang [10]. To approximate the time of the FCE involving two fractional temporal derivatives, Nikan [11] proposed a computational scheme for the radial basis function-generated finite difference (RBF-FD). The Grünwald–Letnikov expansion was used to discretize the time domain of the TFCM, and the RBF-FD was used to discretize the spatial derivatives. They demonstrated that the method can easily be implemented on such types of fractional PDEs. The orthogonal spline collocation with a complete theoretical analysis with order was used by Zhang [12]. Quintana-Murillo [13] researched two temporal R-LFD for explicit numerical approaches to solve FCE. The numerical solution was obtained by using the forward difference formula, the Grünwald–Letnikov formula for the first-order derivative and Riemann–Liouville derivatives, respectively, and the three-point centered formula for the spatial derivative. The stability was tested by the von Neumann technique. Baleanu [14] proposed computational schemes for the optimal control problems of fractional order in the R-LFD sense. The approximations were replaced into optimal control equations of fractional order, and an algebraic equation was obtained, which can be solved by a numerical technique. To model the electrodiffusion of ions in nerve cells with anomalous sub-diffusion along and through the nerve cells, Henry [15] introduced fractional Nernst–Planck equations and related FCE. They analyzed fundamental solutions after modeling the sub-diffusion in two different ways, leading to two FCE. The solution approaches the normal non-zero steady state with uniform sub-diffusion along and around the nerve cells, but the approach is delayed by the anomalous sub-diffusion. Realistic electrophysiological studies on actual dendrites may be related to these solutions. Langlands [16] introduced fractional Nernst–Planck equations and derived FCE as macroscopic models for the electro diffusion of ions in nerve cells. They calculated the power lessening along dendrites in response to synaptic inputs of the alpha function. Easy integration and fire variants of the models were also used to calculate action potential firing rates. Tomovski [17] discussed Laplace and Fourier transforms to formulate the Green function of the generalized space-time FCE, and then examined the even moments to demonstrate that it can have a negative sign, indicating that the Green function does not always flow in one direction and that the current can switch directions. Bhrawy [18] used the collocation method in combination with the Shifted-Jacobi operational matrix in the sense of the Caputo fractional derivative. The results of their suggested approach are much more efficient for solving variable-order non-linear PDEs with high accuracy. Liu [19] presented a discrete numerical formula obtained by finite difference and finite element approximation in time and space, respectively, for the FCE. Liang-lian [20] considered the finite volume approach to solving the FCE using an implicit difference scheme. The approach was also convergent and unconditionally stable. Zhang [21] suggested an unconditionally numerical approach for the convection-diffusion of the fractional order problem. A novel shifted version of the Grünwald–Letnikov formula for the fractional order derivatives was used to prove the accuracy, and for theoretical analysis. Hu [22] implemented compact schemes for the FCE, and utilized the energy method to prove that the first scheme is stable and convergent in -norm with the order , while the second one is an inner product. The computed result indicates that both schemes are accurate and effective. Moshtaghi and Saadatmandi [23] researched the cable model of fractional order and solved using the collocation-type approach. They converted the fractional order model into a set of algebraic equations and presented two numerical examples to confirm the accuracy and efficiency. Aslefallah et al. [24] studied the 2D time-fractional order cable model with Dirichlet boundary conditions and implemented the singular boundary method to split the solution of the inhomogeneous governing equation. More studies related to the fractional order differential equation can be seen in [25–35].
The aim of this study is to find out the numerical solution of the fractional-order cable model. The fractional derivative is approximated by the discretized Riemann–Liouville derivative and for the space derivative use the finite difference approximation. For the proposed approach’s complete theoretical analysis as stability and convergence are discussed. The theoretical analysis, confirms the efficiency and effectiveness of the proposed approach.
Suppose the following fractional order cable model [36] as:with initial and boundary conditionswhere are constants and are known functions and the unknown function is to be determined
The is the Riemann–Liouville fractional derivative of fraction order defined by [37]:
The Riemann–Liouville fractional integral can be discretized [38] as:discretizing the equation at the grid point
As by Jumarie property [39] as:
Again, by Jumarie property as:where The same procedure can be followed for .
Lemma 1:
The coefficients
satisfy the following properties [
35]:
(i)
(ii)
(iii) There exists a positive constant , such that
(iv)
2 Methodology
We implement an implicit numerical approximation for the FCE in Eqs 1–3, utilizing the discretization of the Riemann–Liouville integral with backward difference approximation for the partial derivative using central difference approximation. The steps as along x-axis, where and the step where Letting the obtained numerical solution be to , and using Eq. 4 in Eq. 1, we have
Further, applying Eq. 5 in Eq. 6, we can write
Where
By using implicit discretization with respect to time ‘’, we have
Simplifying Eq. 9, we obtained
Withwhere
3 Stability
In this section, we use the Fourier series method to analyze the stability of the implicit numerical scheme. Letting be the approximate solution for Eq. 10, we have
The error is defined as:where satisfies and
The error and initial conditions are
Here,
Here, we need to define grid functions for as the following:
Then, can be written in Fourier series, such as
Here
From the definition of norm and Parseval equality, we have
Supposing thatwhere and by substituting in Eq. 15, we have
Proposition 1: If satisfies Eq. 22, then .
4 Convergence
To investigate the convergence of the proposed implicit scheme. Let be the exact solution represented by Taylor series, then the local truncation error is obtained as
From Eq. 1
Since and are finite, then there exist a positive constant , then we have
The error is defined as
From Eq. 24, as
With error conditions that are
Next, we define the following grid functions for
And
Here, and can be expanded in Fourier series such aswhere
From the definition of norm and the Parseval equality, we have
Based on the above, supposing thatwhere , by using and ( in Eq. 28, we have
Proposition 2: If is the solution of , then there exists a positive constant such that .
Proof: From and Eq. 29 we have
From ( and (, there exists positive constant , such that
Using mathematical induction, for , then from
From Eq. 36
Now suppose
From Eq. 34 and Lemma 1from Eq. 36
Here, because and
The value of is very small, lying between 0 and 1. So, we obtained
5 Numerical tests
In this study, the numerical result of an implicit scheme for one-dimensional FCE are discussed numerically and graphically. The examples are as following.
Example 1Consider the fractional-order cable model [15] with the closed-form solution is given as:with initial and boundary conditions Where and The closed-form solution is
Example 2Consider the 1D fractional Stokes’ first problem for the heated generalized second-grade equation [40].with initial and boundary conditions Where and The closed-form solution is The errors between a numerical and an exact solution are defined as follows:The above problem is solved using the modified implicit scheme. The errors at and for different values of and . The time step is defined by .
6 Results and discussion
The modified implicit finite difference approximation is used to solve the numerical example of fractional order, such as fractional cable model and the fractional order Stokes’ first problem for the heated generalized second grade equation. Numerical results are presented in the form of tables and figures for various values of space and time steps in order to demonstrate the efficiency of the suggested numerical scheme.
In Tables 1–3, the exact and the numerical solution are compared of the given example 1 for fixed values and , and different values of and . The error decreases as the value of increases. Similarly, as the time and space step size and reduce, the errors decrease for a fixed value of . In Table 4, the exact and the numerical solution of example 1 are compared for , , and , and for different values of and . The results show that as we increase the value of N, i.e., reduce the time and space step size and , the errors decrease for different values of . In Table 5, the numerical results are explained for example, 2 of the suggested scheme for the fractional order first problem for a heated generalized second-grade fluid for various values of order , step size , and . Figures 1–3 show the comparison of the numerical and the exact solution of example 1 in Figure 1 at , = 0.5, , = 1/10, and For Figure 2, at , = 0.5, , = 1/20, and . For Figure 3, at , = 0.5, , = 1/40, and . Furthermore, added Figure 4 which shows the graphical representation of example 2 for and , which confirmed our theoretical analysis and demonstrated that the proposed approach is very powerful and efficient.
TABLE 1
| N | |||
|---|---|---|---|
Numerical results for example, 1 of the modified implicit scheme for various values of , and fixed values of .
TABLE 2
| N | |||
|---|---|---|---|
Numerical results example 1 of the modified implicit scheme for various values of , and fixed value of .
TABLE 3
Numerical results example 1 of the modified implicit scheme for various values of , and for fixed value of .
TABLE 4
| N | ||||
|---|---|---|---|---|
Numerical results example 1 of the modified implicit scheme for various values of , , and .
TABLE 5
Numerical results for example 2 of the modified implicit scheme for various values of , , and .
FIGURE 1
FIGURE 2
FIGURE 3
FIGURE 4
7 Conclusion
This paper presented the modified implicit numerical approximation for a fractional one-dimensional linear Cable model. The scheme is convergent and unconditionally stable, as seen by the investigation using the Fourier series method. The time-fractional derivative was calculated using the Riemann–Liouville formula. The outcome of an application to specific examples of fractional order one-dimensional linear Cable model and the fractional order Stokes’ first problem for the heated generalized second-grade equation have been explored graphically and numerically. The scheme is verified through the comparison of the numerical solution with the exact solution, which shows an agreement with the theoretical analysis and the numerical experiment, confirming that the approximate solution converges to the exact solution. This modified approach can also extend to other types of two and three dimensional fractional order differential models.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work under Grant code: 22UQU4310396DSR58.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
ZhuangPLiuF. Finite difference approximation for two-dimensional time fractional diffusion equation. J Algorithms Comput Technol (2007) 1:1–16. 10.1260/174830107780122667
2.
RajkovićPMMarinkovićSDStankovićMS. Fractional integrals, and derivatives in q-calculus. Appl Anal Discret Math (2007) 1:311–23. 10.2298/aadm0701311r
3.
YakarAAkdoganZ. On the fundamental solutions of a discontinuous fractional boundary value problem. Adv Differ Equ (2017) 1:378–15. 10.1186/s13662-017-1433-6
4.
HolmesWR. Cable equation. New York: Springer (2014).
5.
VitaliSCastellaniGMainardiF. Time fractional cable equation, and applications in neurophysiology. Chaos Solitons Fractals (2017) 102:467–72. 10.1016/j.chaos.2017.04.043
6.
YuBJiangX. Numerical identification of the fractional derivatives in the two-dimensional fractional cable equation. J Sci Comput (2016) 68:252–72. 10.1007/s10915-015-0136-y
7.
LiuFYangQTurnerI. Two new implicit numerical methods for the fractional cable equation. J Comput Nonlinear Dyn (2011) 6(1-7):011009. 10.1115/1.4002269
8.
LinYLiXXuC. Finite difference/spectral approximations for the fractional cable equation. Math Comput (2011) 80:1369–96. 10.1090/s0025-5718-2010-02438-x
9.
LiuYDuYWLiHWangJF. A two-grid finite element approximation for a nonlinear time-fractional Cable equation. Nonlinear Dyn (2016) 85:2535–48. 10.1007/s11071-016-2843-9
10.
ZhuangPLiuFTurnerIAnhV. Galerkin finite element method and error analysis for the fractional cable equation. Numer Algorithms (2016) 72:447–66. 10.1007/s11075-015-0055-x
11.
NikanOGolbabaiAMachadoJTNikazadT. Numerical approximation of the time fractional cable model arising in neuronal dynamics. Eng Comput (2020) 1–19.
12.
ZhangHYangXHanX. Discrete-time orthogonal spline collocation method with application to two-dimensional fractional cable equation. Comput Math Appl (2014) 68:1710–22. 10.1016/j.camwa.2014.10.019
13.
Quintana-MurilloJYusteSB. An explicit numerical method for the fractional cable equation. Int J Differ Equ (2011) 2011:1–12. 10.1155/2011/231920
14.
BaleanuDDefterliOAgrawalOP. A central difference numerical scheme for fractional optimal control problems. J Vib Control (2009) 15:583–97. 10.1177/1077546308088565
15.
LanglandsTAMHenryBIWearneSL. Fractional cable equation models for anomalous electrodiffusion in nerve cells: Finite domain solutions. SIAM J Appl Math (2011) 71:1168–203. 10.1137/090775920
16.
LanglandsTAMHenryBIWearneSL. Fractional cable equation models for anomalous electrodiffusion in nerve cells: Infinite domain solutions. J Math Biol (2009) 59:761–808. 10.1007/s00285-009-0251-1
17.
SaxenaRKTomovskiZSandevT. Analytical solution of generalized space-time fractional cable equation. Mathematics (2015) 3:153–70. 10.3390/math3020153
18.
BhrawyAHZakyMA. Numerical simulation for two-dimensional variable-order fractional nonlinear cable equation. Nonlinear Dyn (2015) 80:101–16. 10.1007/s11071-014-1854-7
19.
LiuJLiHLiuY. A new fully discrete finite difference/element approximation for fractional cable equation. J Appl Math Comput (2016) 52:345–61. 10.1007/s12190-015-0944-0
20.
MaLLiuD. An implicit difference approximation for fractional cable equation in high-dimensional case. J Liaoning Tech Univ Nat Sci (2014) 4.
21.
ZhangY. A finite difference method for the fractional partial differential equation. Appl Math Comput (2009) 215:524–9.
22.
HuXZhangL. Implicit compact difference schemes for the fractional cable equation. Appl Math Model (2012) 36:4027–43. 10.1016/j.apm.2011.11.027
23.
MoshtaghiNSaadatmandiA. Numerical solution of time fractional cable equation via the sinc-Bernoulli collocation method. J Appl Comput Mech (2020) 7:1–9.
24.
AslefallahMAbbasbandySShivanianE. Fractional cable problem in the frame of meshless singular boundary method. Eng Anal Bound Elem (2019) 108:124–32. 10.1016/j.enganabound.2019.08.003
25.
AliUSohailMUsmanMAbdullahFAKhan NisarKS. Fourth-order difference approximation for time-fractional modified sub-diffusion equation. Symmetry (2020) 12:691. 10.3390/sym12050691
26.
KumarSChauhanMomaniRPSHadiS. Numerical investigations on COVID‐19 model through singular and non‐singular fractional operators. Numer Methods Partial Differential Equations (2020). 10.1002/num.22707
27.
KumarSKumarASametBDuttaH, A study on fractional host–parasitoid population dynamical model to describe insect species, Numer Methods Partial Differential Equations37 (2), 1673–92. 10.1002/num.226032021). A study on fractional host–parasitoid population dynamical model to describe insect species,
28.
KhanMAUllahSS.K. A robust study on 2019-nCOV outbreaks through non-singular derivative. The Eur Phys J Plus (2021) 136:168.
29.
KumarSKumarROsmanMS. A wavelet based numerical scheme for fractional orderSEIRepidemic of measles by using Genocchi polynomials. Numer Methods Partial Differential Equations (2021) 37(2):1250–68. 10.1002/num.22577
30.
MohammadiHKumarSRezapourSEtemadS, A theoretical study of the Caputo–Fabrizio fractional modeling for hearing loss due to Mumps virus with optimal control, Chaos, Solitons and Fractals (2021) 144, 11066810.1016/j.chaos.2021.110668
31.
KumarSKumarRCattaniCSametB, Chaotic behaviour of fractional predator-prey dynamical system, Chaos Solitons Fractals (2020)135, 10981110.1016/j.chaos.2020.109811
32.
AliUKamalRMohyud-DinST. On nonlinear fractional differential equations. Int J Mod Math Sci (2012) 3(3):116–24.
33.
ZubairTUsmanMAliUMohyud-DinST. Homotopy analysis method for a system of partial differential equations. Int J Mod Eng Sci (2012) 1:67–79.
34.
AliUAbdullahFA. Explicit Saul’yev finite difference approximation for the two-dimensional fractional sub-diffusion equation. In: AIP Conference Proceedings. AIP Publishing LLC (2018). 020111.
35.
AliUSohailMAbdullahFA. An efficient numerical scheme for variable-order fractional sub-diffusion equation. Symmetry (2020) 12:1437. 10.3390/sym12091437
36.
YangQ. Novel analytical and numerical methods for solving fractional dynamical systems. Dr Diss Qld Univ Technol (2010) 1–201.
37.
PovstenkoYOstoja-StarzewskiM. Fractional telegraph equation under moving time-harmonic impact. Int J Heat Mass Transfer (2022) 182:121958. 10.1016/j.ijheatmasstransfer.2021.121958
38.
KhaterMAliUKhanMAMousaAAAttiaRA. A new numerical approach for solving 1D fractional diffusion-wave equation. J Funct Spaces (2021) 2021:1–7. 10.1155/2021/6638597
39.
JumarieG. Modified Riemann-Liouville derivative and fractional Taylor series of non-differentiable functions further results. Comput Math Appl (2006) 51:1367–76. 10.1016/j.camwa.2006.02.001
40.
ChunhongW. Numerical solution for Stocks’ first problem for a heated generalized second grade fluid with fractional derivative. Appl Numer Math (2009) 59:2571–83.
Summary
Keywords
fractional cable equation, implicit approximation, stability, convergence, riemann-liouville fractional derivative
Citation
Ali U, Naeem M, Ganie AH, Fathima D, Salama FM and Abdullah FA (2023) Numerical approach for the fractional order cable model with theoretical analyses. Front. Phys. 11:1160767. doi: 10.3389/fphy.2023.1160767
Received
07 February 2023
Accepted
29 March 2023
Published
25 April 2023
Volume
11 - 2023
Edited by
Alexander Nepomnyashchy, Technion Israel Institute of Technology, Israel
Reviewed by
Yuriy Povstenko, Jan Długosz University, Poland
Sunil Kumar, National Institute of Technology, India
Updates
Copyright
© 2023 Ali, Naeem, Ganie, Fathima, Salama and Abdullah.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Umair Ali, umairkhanmath@gmail.com; Dowlath Fathima, d.fathima@seu.edu.sa
This article was submitted to Mathematical Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.