Home A numerical study of anomalous electro-diffusion cells in cable sense with a non-singular kernel
Article Open Access

A numerical study of anomalous electro-diffusion cells in cable sense with a non-singular kernel

  • Azhar Iqbal EMAIL logo and Tayyaba Akram EMAIL logo
Published/Copyright: September 23, 2022
Become an author with De Gruyter Brill

Abstract

The time-fractional cable model is solved using an extended cubic B-spline (ECBS) collocation strategy. The B-spline function was used for space partitioning, while the Caputo-Fabrizio (CF) was used for temporal discretization. The finite difference technique was used to discretize the CF operator. For the first time in cable modeling, the CF operator has been used. In terms of time, the convergence of order τ . An ECBS collocation approach is investigated by numerical example at different values, and comparisons with published work are made. The numerical results show that the scheme performed well, and the graphical representations show that the results are very close to exact values. The Von Neumann technique is applied to investigate the stability of the proposed scheme.

MSC 2010: 35R11; 26A33; 65N12; 65N06

1 Introduction

Anomalous diffusion in cells has been observed in several research experiments, and this anomalous diffusion is primarily exhibited by fractional derivatives. Brownian motion can be used to explain the process of signal transmission to the insides of cells. As a consequence of the irregularity, the diffusion in the cells can become aberrant. The voltage transmission in axons is affected in this situation [1]. The fractional diffusion model has caught the interest of many academicians in recent years because of its numerous applications in various areas. The special case of time-fractional cable model [2] is one of these models that we will look at in this article.

(1) U ( y , t ) t + D 0 , t γ C U ( y , t ) 2 U ( y , t ) y 2 + U 2 ( y , t ) = F ( y , t ) , y Ω , t > 0 ,

with initial and boundary value conditions, respectively,

(2) U ( y , 0 ) = U 0 ( y ) , y Ω , U ( y , t ) y Ω = 0 , t 0 ,

where γ > 0 is a constant, Ω = [ a , b ] , U 0 F , is the smooth function. Also, D 0 , t γ C is the γ th order Caputo-Fabrizio (CF) operator and can be defined as [3]:

(3) γ G ( t ) t γ = N ( γ ) 1 γ 0 t G ( τ ) exp γ ( t τ ) 1 γ d τ , 0 < γ < 1 ,

where N ( γ ) is a normalizing function, such that N ( 0 ) = N ( 1 ) = 1 .

If G ( t ) is a constant function, then the CF derivative of F ( t ) is zero, which is the same as the Caputo derivative. As a result, no singularity exists in the kernel.

The CF derivative has added a new dimension to the study of fractional differential equations (FDEs). The Caputo, Riemann-Liouville, and other singular kernel operators have kernels for power law that have limitations in representing physical situations [4]. It is produced using the exponential function and regular derivative convolution and it retains the same intrinsic inspirational qualities of heterogeneity and configuration for many scales as the Caputo and Riemann-Liouville fractional derivatives [5,6, 7,8]. This derivative has already been mentioned in a number of articles, such as fractional Maxwell fluid [9], fluid flows [10], diffusive transport system [6], non-linear Fisher’s diffusion equation [11], and electric circuits [5]. To explain the electrotonic characteristics of spiny neural dendrites, Henry et al. [12] presented cable equations with fractional order temporal operators. The anomalous electrodiffusion can be alternatively modeled by a fractal-fractional model in [13,14]. He’s fractional derivative for the evolution equation has been adopted by Wang and Yao [15].

The time-fractional cable model is equivalent to the classic cable equation other than that the order of derivative is fractional with regard to time. Many fractional models have no analytical solutions due to the non-locality of fractional derivatives. Hence, numerical solutions to fractional differential model are important in terms of both science and practical application. As a result, several investigators have recently taken fractional differential models as a challenge and discovered numerical solutions to such problems. Two implicit compact difference techniques for solving fractional cable model were developed by Hu and Zhang [16]. Yu and Jiang [17] studied time-fractional cable model using fourth-order compact difference method. The local discontinuous Galerkin technique has been developed by Zheng and Zhao [18]. Yang et al. [19] proposed a temporal-spacial spectral-tau technique for solving fractional cable and its inverse model. Liu et al. [20,21] solved the fractional cable equations via finite element approximations and Grünwald difference approximations. Zhu et al. [22] presented the Galerkin finite element and convolution quadrature methods for the solution of fractional cable equation. Various researchers have worked on the solution of the FDEs by using numerical technique based on spline function [23,24, 25,26]. Li [27] has studied the Wavelet collocation approach based on B-spline functions for solving the FDEs. B splines helped to resolve the unknowns in the equation. It was shown that the B-spline finite element is useful for the numerical solutions of FDE, particularly when continuity of the solutions was essential. Saad et al. [28] used a Atangana-Baleanu fractional derivative with spectral collocation methods to approximate the solution of the fractional Fisher’s type equations. The results indicate that the proposed methodology is a simple and effective tool for studying nonlinear equations with local and non-local singular kernels. Akram et al. [29,30] investigated time fractional cable equation and sub-diffusion equation using B-spline methods, Riemann-Liouville, and Caputo derivatives. The stability and convergence of the proposed methods have also been examined. Moshtaghi and Saadatmandi [31] developed the numerical technique for time fractional cable problem using Sinc-Bernoulli and Riemann-Liouville methods. A Galerkin method with spline function for solving time fractional diffusion equation has been proposed by Pezza and Pitolli [32]. Tabriz et al. [33] developed the operational matrix method based on Caputo derivative to solve fractional optimal control problem.

Despite the fact that B-spline has been used to solve fractional partial differential equations in several articles, we are not aware of any research that has used B-spline to solve time-fractional cable model equations. The B-spline function has very suitable properties for designing such convex hull and continuity properties [34]. Thus, the main objective of this research is to use the extended cubic B-spline (ECBS) to develop the collocation approach for the time-fractional cable model. The strategy’s main feature is that it converts such problems into algebraic system of equations that can be used in programming. The following is the organization of this research article: Section 2 presents the extended cubic B-spline (ECBS) functions and time approximation. The ECBS collocation method is presented in Section 3. The stability analysis is proved in Section 4. Sections 5 and 6 discuss the numerical experiment as well as the conclusion.

2 Preliminaries

2.1 ECBS functions

Assume that { y j } is a division of equal length depending on the given interval with j Z . The assumed interval at the endpoints is split into N equivalent subintervals using y j = y 0 + j σ , where σ is the grid length. The basis function at the nodes y j over the given interval is described as follows:

(4) E i ( y , δ ) = 1 24 σ 4 4 σ ( 1 δ ) ( y y j 2 ) 3 + 3 δ ( y y j 2 ) 4 , y [ y j 2 , y j 1 ) , ( 4 δ ) σ 4 + 12 σ 3 ( y y j 1 ) + 6 σ 2 ( 2 + δ ) ( y y j 1 ) 2 12 σ ( y y j 1 ) 3 3 δ ( y y j 1 ) 4 , y [ y j 1 , y j ) , ( 4 δ ) σ 4 + 12 σ 3 ( y j + 1 y ) + 6 σ 2 ( 2 + δ ) ( y j + 1 y ) 2 12 σ ( y j + 1 y ) 3 3 δ ( y j + 1 y ) 4 , y [ y j , y j + 1 ) , 4 σ ( 1 δ ) ( y j + 2 y ) 3 + 3 δ ( y j + 2 y ) 4 , y [ y j + 1 , y j + 2 ) , 0 , otherwise ,

where j = 1 ( 1 ) N + 1 , δ R represent the free parameter in the [ 8 , 1 ] . y R represents a variable. The extended and cubic B-spline functions have the same features for δ [ 8 , 1 ] . The ECBS changes to a cubic B-spline when δ = 0 . For a function U ( y , t ) there is a unique U ˆ ( y , t ) , which assures the given requirements, as

(5) U ˆ ( y , t ) = j = k 1 k + 1 d j ( t ) E j ( y , δ ) .

The ECBS functions (4) and equation (5) yield the equations as follows:

(6) U ˆ ( y , t ) = j = k 1 k + 1 d j ( t ) E j ( y , δ ) = 4 δ 24 d k 1 + 8 + δ 12 d k + 4 δ 24 d k + 1 ,

(7) U ˆ ( y , t ) = j = k 1 k + 1 d j ( t ) E j ( y , δ ) = 1 2 σ d k 1 + 1 2 σ d k + 1 ,

(8) U ˆ ( y , t ) = j = k 1 k + 1 d j ( t ) E j ( y , δ ) = 2 + δ 2 σ 2 d k 1 + 4 + 2 δ 2 σ 2 d k + 2 + δ 2 σ 2 d k + 1 .

2.2 Difference method for the CF

In this section, we look at using CF to discretize fractional time derivatives. Assume that t l = t 0 + l τ , l = 0 , 1 , , L , where τ = T L is the grid size in time direction. The difference method is applied to discretize CF in the bounded domain. The CF can be expressed as follows using equation (3):

U γ ( y , t l + 1 ) t γ = 1 1 γ 0 t l + 1 U ( y , η ) η exp γ 1 γ ( t l + 1 η ) d η .

The discretized formulation [23] is as follows:

(9) U γ ( y , t l + 1 ) t γ = 1 τ γ m = 0 l ε m [ U ( y , t l m + 1 ) U ( y , t l m ) ] 1 exp γ 1 γ τ + G τ l + 1 .

The truncation error is [23] as follows:

(10) G τ l + 1 R τ 2 ,

where ε m = exp γ 1 γ τ m ; it is simple to demonstrate the characteristics of ε m [23]:

  • ε 0 = 1 ,

  • ε 0 > ε 1 > ε 2 > > ε m , ε m 0 as m ,

  • ε m > 0 for m = 0 , 1 , , l ,

  • m = 0 l ( ε m ε m + 1 ) + ε l + 1 = ( 1 ε 1 ) + m = 1 l 1 ( ε m ε m + 1 ) + ε l = 1 .

3 Derivation of the ECBS

To develop the numerical strategy for solving the cable model, we use the finite difference method, the CF, and the EBSM in this section. By plugging equations (9) and (4) in equation (1), we obtain

(11) U ˆ ( y , t l + 1 ) U ˆ ( y , t l ) τ + 1 τ γ m = 0 l ε m [ U ˆ ( y , t l m + 1 ) U ˆ ( y , t l m ) ] 1 exp γ 1 γ τ U ˆ ( y , t l + 1 ) + U 2 ˆ ( y , t l + 1 ) = F ( y , t l + 1 ) .

The nonlinear term is linearized by Taylor’s series expansion as follows:

( U 2 ) l + 1 = ( U 2 ) l + τ [ ( U 2 ) t ] l + O ( τ ) 2 = ( U 2 ) l + τ [ 2 U U t ] l + O ( τ ) 2 = ( U 2 ) l + τ 2 U l U l + 1 U l τ + O ( τ ) + O ( τ ) 2 = ( U 2 ) l + 2 U l U l + 1 2 U l U l + O ( τ ) 2 = 2 U l U l + 1 U l U l + O ( τ ) 2 ,

which implies that

(12) ( U 2 ) l + 1 = 2 U l U l + 1 ( U 2 ) l .

Using equation (12) in (11), we obtain

U ˆ ( y , t l + 1 ) U ˆ ( y , t l ) τ + 1 τ γ m = 0 l ε m [ U ˆ ( y , t l m + 1 ) U ˆ ( y , t l m ) ] 1 exp γ 1 γ τ U ˆ ( y , t l + 1 ) + 2 U ˆ ( y , t l + 1 ) U ˆ ( y , t l ) U 2 ˆ ( y , t l ) = F ( y , t l + 1 ) .

The aforementioned equation can be written in another form as:

s U ˆ l + 1 + s 1 τ U ˆ l + 1 s τ U ˆ y y l + 1 + 2 τ s U ˆ l U ˆ l + 1 = s U ˆ l + s 1 τ U ˆ l + s τ ( U 2 ˆ ) l s 1 τ m = 1 l ε m [ U ˆ l m + 1 U ˆ l m ] + s τ F ˆ l + 1 ,

where s 1 = 1 exp γ 1 γ τ and s = γ τ .

(13) s j = k 1 k + 1 d j l + 1 E j + s 1 τ j = k 1 k + 1 d j l + 1 E j s τ j = k 1 k + 1 d j l + 1 E j + 2 τ s j = k 1 k + 1 d j l E j j = k 1 k + 1 d j l + 1 E j = s j = k 1 k + 1 d j l E j + s 1 τ j = k 1 k + 1 d j l + 1 E j + s τ j = k 1 k + 1 d j l E j 2 + s τ F ˆ l + 1 s 1 τ m = 1 l ε m j = k 1 k + 1 d j l m + 1 E j j = k 1 k + 1 d j l m E j .

The aforementioned system of equations have the order ( N + 1 ) × ( N + 3 ) . In order to obtain a unique solution, two equations are required which can be derived from the boundary conditions. The following two boundary conditions are needed to solve the aforementioned system uniquely.

4 δ 24 d 1 + 8 + δ 12 d 0 + 4 δ 24 d 1 = U ( y 0 , t ) ,

4 δ 24 d N 1 + 8 + δ 12 d N + 4 δ 24 d N + 1 = U ( y N , t ) .

Hence, we have the ( N + 3 ) × ( N + 3 ) order system. Now, we will employ the initial vector from the initial conditions to start the iteration on system (13).

3.1 Algorithm

  • Describe the model.

  • Discretize the modeled problem using the CF operator in time direction.

  • Linearize the nonlinear terms.

  • Discretize the problem using ECBS formulation in space direction.

  • Apply the two boundary conditions.

  • Resulting system of order ( N + 3 ) × ( N + 3 ) is obtained.

  • Solve the aforementioned system using Mathematica.

4 Stability and convergence analysis

The stability concept is linked to computation approach errors that do not increase as the procedure continues. The Von Neumann stability technique also known as Fourier stability is a method utilized to analyze the stability of difference techniques. Here, we will use the Von Neumann technique to investigate the stability of the proposed method. Assume that η l is the growth factor in the Fourier mode form, and that η ˆ l is the obtained value. As a result, we defined the error concept at the l th time step as

(14) ζ l = η l η ˆ l .

To linearize the non-linear term in (11), we use U 2 = α U , while α is a positive constant, we obtain

U ˆ ( y , t l + 1 ) U ˆ ( y , t l ) + τ s m = 0 l ε m [ U ˆ ( y , t l m + 1 ) U ˆ ( y , t l m ) ] s 1 τ U ˆ ( y , t l + 1 ) + τ α U ˆ ( y , t l + 1 ) = τ F ˆ ( y , t l + 1 ) .

The aforementioned equation becomes

(15) s U ˆ l + 1 + s 1 τ U ˆ l + 1 s τ U ˆ y y l + 1 + τ α U ˆ l + 1 = s U ˆ l + s 1 τ U ˆ l s 1 τ m = 1 l ε m [ U ˆ l m + 1 U ˆ l m ] + τ F ˆ l + 1 .

By inserting equation (14) in equation (15), we obtain the following error equation:

(16) s ζ l + 1 + s 1 τ ζ l + 1 s τ ζ y y l + 1 + τ α ζ l + 1 = s ζ l + s 1 τ ζ l s 1 τ m = 1 l ε m [ ζ l m + 1 ζ l m ] .

Consider the error solution for the B-spline can be written in one Fourier mode as

(17) ζ j l = Φ l e i j σ ρ ,

where Φ , σ , ρ and i = 1 are the Fourier coefficient, element size, and mode number, respectively. By combining equation (17) and B-spline functions in (16), we obtain the following result:

s Φ l + 1 ( E 1 e i ( j 1 ) σ ρ + E 2 e i j σ ρ + E 1 e i ( j + 1 ) σ ρ ) + s 1 τ Φ l + 1 ( E 1 e i ( j 1 ) σ ρ + E 2 e i j σ ρ + E 1 e i ( j + 1 ) σ ρ ) s τ ( E 4 e i ( j 1 ) σ ρ + E 5 e i j σ ρ + E 4 e i ( j + 1 ) σ ρ ) Φ l + 1 + τ α ( E 1 e i ( j 1 ) σ ρ + E 2 e i j σ ρ + E 1 e i ( j + 1 ) σ ρ ) Φ l + 1 = s Φ l ( E 1 e i ( j 1 ) σ ρ + E 2 e i j σ ρ + E 1 e i ( j + 1 ) σ ρ ) + s 1 τ Φ ˆ l ( E 1 e i ( j 1 ) σ ρ + E 2 e i j σ ρ + E 1 e i ( j + 1 ) σ ρ ) s 1 τ m = 1 l ε m ( E 1 e i ( j 1 ) σ ρ + E 2 e i j σ ρ + E 1 e i ( j + 1 ) σ ρ ) [ Φ l m + 1 Φ l m ] ,

where E 1 = 4 δ 24 , E 2 = 8 + δ 12 , E 3 = 1 2 σ , E 4 = 2 + δ 2 σ 2 , and E 5 = 2 + δ σ 2 . By calculating and rearranging, we obtain

( s + s 1 τ + τ α ) Φ l + 1 ( E 2 + 2 E 1 cos ( ε σ ) ) s τ ( E 5 + 2 E 4 cos ( ε σ ) ) Φ l + 1 = s Φ l ( E 2 + 2 E 1 cos ( ε σ ) ) + s 1 τ ε l ( E 2 + 2 E 1 cos ( ε σ ) ) Φ 0 + s 1 τ ( E 2 + 2 E 1 cos ( ε σ ) ) m = 0 l 1 [ ε m ε m + 1 ] Φ l m .

Taking the part that is identical on both sides and dividing it by E 2 + 2 E 1 cos ( ε σ ) , we obtain

(18) [ ( s + s 1 τ + τ α ) + s τ β ] Φ l + 1 = s Φ l + s 1 τ ε l Φ 0 + s 1 τ m = 0 l 1 [ ε m ε m + 1 ] Φ l m ,

where β = 6 ( 2 + δ ) sin 2 ρ σ 2 σ 2 ( 6 + ( δ 4 ) sin 2 ρ σ 2 ) > 0 , δ 2 .

Proposition 4.1

Consider Φ l , l = 0 , 1 , , L be the solution of model problem (1), we write

(19) Φ l Φ 0 , l = 0 , 1 , , L .

Proof

Here, we use mathematical induction to validate this result. The following equation can be obtained when we choose l = 0 in equation (18):

[ ( s + s 1 τ + τ α ) + s τ β ] Φ 1 = s Φ 0 + s 1 τ ε 0 Φ 0 = ( s + s 1 τ ε 0 ) Φ 0 Φ 1 = s + s 1 τ ε 0 ( s + s 1 τ + τ α ) + s τ β Φ 0 Φ 1 Φ 0 . ε 0 = 1 .

Suppose that Φ l Φ 0 is true for l = 0 , 1 , , L 1 . For l + 1 , we have

[ ( s + s 1 τ + τ α ) + s τ β ] Φ l + 1 = s Φ l + s 1 τ ε l Φ 0 + s 1 τ m = 0 l 1 [ ε m ε m + 1 ] Φ l m , [ ( s + s 1 τ + τ α ) + s τ β ] Φ l + 1 s Φ l + s 1 τ ε l Φ 0 + s 1 τ m = 0 l 1 [ ε m ε m + 1 ] Φ l m , s Φ 0 + s 1 τ ε l + m = 0 l 1 [ ε m ε m + 1 ] Φ 0 , Φ l + 1 Φ 0 .

Thus, Φ l + 1 = ζ j l + 1 Φ 0 = ζ j 0 , so that ζ j l + 1 2 Φ 0 2 . As a result, the presented approach is unconditionally stable.□

Theorem 1

Let U ( y j , t l ) be the exact solution and U ˆ l be the time discrete result of the given equation, then

(20) ϒ l + 1 B τ ,

where ϒ l + 1 = U ( y j , t l + 1 ) U ˆ l + 1 .

Proof

The non-linear term in (11) can be linearized as U 2 = α U , where α is a positive constant, we obtain the difference equation of time discrete solution and exact solution as follows:

(21) s ϒ l + 1 + s 1 τ ϒ l + 1 s τ ϒ y y l + 1 + τ α ϒ l + 1 = s ϒ l + s 1 τ ϒ 0 + s 1 τ m = 1 l 1 [ ε m ε m + 1 ] ϒ l m + τ R l + 1 .

For l = 0 , we obtain

s ϒ 1 + s 1 τ ϒ 1 s τ ϒ y y 1 + τ α ϒ 1 = s ϒ 0 + s 1 τ ϒ 0 + τ R 1 , s ϒ 1 + s 1 τ ϒ 1 s τ ϒ y y 1 + τ α ϒ 1 = τ R 1 , ϒ = 0 .

Taking inner product of the aforementioned equation with ϒ , we obtain

s ϒ 1 , ϒ 1 + s 1 τ ϒ 1 , ϒ 1 s τ ϒ y y 1 , ϒ 1 + τ α ϒ 1 , ϒ 1 = τ R 1 , ϒ 1 .

For u y y , u = u y , u y , u , u = u 2 , and u , u 1 u u 1 , the following result is obtained:

s ϒ 1 , ϒ 1 + s 1 τ ϒ 1 , ϒ 1 + τ α ϒ 1 , ϒ 1 = τ R 1 , ϒ 1 s τ ϒ y 1 , ϒ y 1 , s ϒ 1 2 + s 1 τ ϒ 1 2 + τ α ϒ 1 2 τ R 1 ϒ 1 , ϒ 1 τ s + s 1 τ + τ α R 1 τ s + s 1 τ + τ α R 1 τ 2 ,

which implies that

(22) ϒ 1 B 1 τ ,

where B 1 is a constant free of τ . Assume that (22) is true when l = 0 , 1 , , L 1 . Taking the inner product of (21) with ϒ l + 1 , we obtain

s ϒ l + 1 , ϒ l + 1 + s 1 τ ϒ l + 1 , ϒ l + 1 s τ ϒ y y l + 1 , ϒ l + 1 + τ α ϒ l + 1 , ϒ l + 1 = s ϒ l , ϒ l + 1 + s 1 τ ϒ 0 , ϒ l + 1 + s 1 τ m = 1 l 1 [ ε m ε m + 1 ] ϒ l m , ϒ l + 1 + τ R l + 1 , ϒ l + 1 .

Using ϒ 0 = 0 , u y y , u = u y , u y , u , u = u 2 , and u , u 1 u u 1 , we obtain

s ϒ l + 1 , ϒ l + 1 + s 1 τ ϒ l + 1 , ϒ l + 1 + τ α ϒ l + 1 , ϒ l + 1 = s τ ϒ y l + 1 , ϒ y l + 1 + s ϒ l , ϒ l + 1 + s 1 τ m = 1 l 1 [ ε m ε m + 1 ] ϒ l m , ϒ l + 1 + τ R l + 1 , ϒ l + 1 ,

which implies that

s ϒ l + 1 , ϒ l + 1 + s 1 τ ϒ l + 1 , ϒ l + 1 + τ α ϒ l + 1 , ϒ l + 1 s ϒ l , ϒ l + 1 + s 1 τ m = 1 l 1 [ ε m ε m + 1 ] ϒ l m , ϒ l + 1 + τ R l + 1 , ϒ l + 1 ,

s ϒ l + 1 2 + s 1 τ ϒ l + 1 2 + τ α ϒ l + 1 2 s ϒ l ϒ l + 1 + s 1 τ m = 1 l 1 [ ε m ε m + 1 ] ϒ l m ϒ l + 1 + τ R l + 1 ϒ l + 1 .

By dividing throughout with ϒ l + 1 , we obtain

s ϒ l + 1 + s 1 τ ϒ l + 1 + τ α ϒ l + 1 s ϒ l + s 1 τ m = 0 l 1 [ ε m ε m + 1 ] ϒ l m + τ R l + 1 , s + s 1 τ + α τ s 1 τ ϒ l + 1 s s 1 τ ϒ l + m = 0 l 1 [ ε m ε m + 1 ] ϒ l m + R l + 1 s 1 .

By Gronwall’s inequality, we obtain

s + s 1 τ + α τ s 1 τ ϒ l + 1 s s 1 τ ϒ l + R l + 1 s 1 exp m = 0 l 1 [ ε m ε m + 1 ] , ϒ l + 1 s ϒ l + τ R l + 1 s + s 1 τ + α τ exp ( 1 ε l ) , ϒ l + 1 s R l τ 2 + τ R l + 1 τ 2 s + s 1 τ + α τ exp ( 1 ε l ) , ϒ l + 1 B τ .

5 Numerical results

In this section, we incorporate the simulation findings of given problem by using the provided scheme. The numerical findings are used to test the claimed method validity and efficiency. The main goal of these examples is to examine the degree of convergence of the estimated results for γ .

Example 1

U ( y , t ) t + C D 0 , t γ U ( y , t ) 2 U ( y , t ) y 2 + U 2 ( y , t ) = F ( y , t ) , y Ω = [ 0 , 2 π ] , t [ 0 , 1 ] ,

with

U ( y , 0 ) = U 0 = 0 ,

and

F ( y , t ) = 2 t + t γ ( 1 γ ) γ 1 e γ 1 γ sin y + t 2 sin y ( 1 + t 2 sin y ) .

The exact solution of the problem is:

U ( y , t ) = t 2 sin y .

By taking γ = 0.5 , σ = 0.05 , in Table 1, we show the L 2 , L and convergence order at different time-steps. The comparison of L 2 error norm has also presented at T = 1 . In Table 2, we demonstrate the L 2 , L and convergence order and comparison of L 2 error for fixed σ = L 20 , γ = 0.01 when T = 1 . Figure 1 depicts the comparison of approximated and exact results with γ = 0.5 , L 50 when T = 1 . Figures 2 and 3 show that the absolute error graph of Example 1 for σ = 2 π 50 , γ = 0.2 , 0.4 , 0.5 , and 0.6 at T = 1 . Figure 4 displays the contour plot for γ = 0.5 , σ = 2 π 25 , L = 25 at T = 1 (Table 3).

Table 1

L , L 2 erorrs and convergence rate of (1) for γ = 0.5 , σ = 2 π 20 , at T = 1

τ γ = 0.5
L 2 [2] Proposed method L 2 Proposed method L Order
1/5 1.4927 × 1 0 1 3.94746 × 1 0 3 3.64701 × 1 0 2
1/10 4.0466 × 1 0 2 2.16374 × 1 0 3 1.79683 × 1 0 2 1.02126
1/20 1.0758 × 1 0 2 1.30405 × 1 0 3 8.44512 × 1 0 3 1.08926
1/40 2.7868 × 1 0 3 4.66346 × 1 0 4 4.16857 × 1 0 3 1.01856
Table 2

L , L 2 errors of (1) for γ = 0.01 , σ = 2 π 20 , at T = 1

τ γ = 0.01
L 2 [2] Proposed method L 2 Proposed method L Order
1/5 1.5542 × 1 0 1 8.99805 × 1 0 3 5.85797 × 1 0 2
1/10 4.3922 × 1 0 2 3.37174 × 1 0 3 2.84031 × 1 0 2 1.04435
1/20 1.1966 × 1 0 2 1.72413 × 1 0 3 1.36493 × 1 0 3 1.05722
1/40 3.1359 × 1 0 3 9.17447 × 1 0 4 6.59511 × 1 0 3 1.04936
Figure 1 
               Exact and approximated solution of Example 1 for value of 
                     
                        
                        
                           γ
                           =
                           0.5
                        
                        \gamma =0.5
                     
                  , at 
                     
                        
                        
                           T
                           =
                           1
                        
                        T=1
                     
                  .
Figure 1

Exact and approximated solution of Example 1 for value of γ = 0.5 , at T = 1 .

Figure 2 
               Space-time absolute error plot for value of 
                     
                        
                        
                           γ
                           =
                           0.2
                           ,
                           0.4
                        
                        \gamma =0.2,0.4
                     
                  , respectively, at 
                     
                        
                        
                           T
                           =
                           1
                        
                        T=1
                     
                  .
Figure 2

Space-time absolute error plot for value of γ = 0.2 , 0.4 , respectively, at T = 1 .

Figure 3 
               Space-time absolute error plot for value of 
                     
                        
                        
                           γ
                           =
                           0.5
                           ,
                           0.6
                        
                        \gamma =0.5,0.6
                     
                  , respectively, at 
                     
                        
                        
                           T
                           =
                           1
                        
                        T=1
                     
                  .
Figure 3

Space-time absolute error plot for value of γ = 0.5 , 0.6 , respectively, at T = 1 .

Figure 4 
               Contour plot for value of 
                     
                        
                        
                           γ
                           =
                           0.5
                        
                        \gamma =0.5
                     
                   at 
                     
                        
                        
                           T
                           =
                           1
                        
                        T=1
                     
                  .
Figure 4

Contour plot for value of γ = 0.5 at T = 1 .

Table 3

Absolute errors for different γ at various knots

y γ = 0.2 γ = 0.4 γ = 0.5 γ = 0.6
0.1 0.000505826 0.00093832 0.00137872 0.00183893
0.2 0.000839704 0.00160334 0.00239615 0.00340647
0.3 0.000852242 0.00161478 0.00240241 0.00340672
0.4 0.000587011 0.00103535 0.00145486 0.00185383
0.5 0.000342879 0.00053744 0.00058313 0.00040270
0.6 0.000419648 0.00080647 0.00090600 0.00103866
0.7 0.000677912 0.00148961 0.00193893 0.00299914
0.8 0.000658044 0.00147090 0.00192816 0.00299863
0.9 0.000319919 0.00068692 0.00081109 0.00101909

6 Conclusion

This study describes an ECBS collocation strategy for solving the time-fractional cable model. For space partitioning, the ECBS was used, whereas for temporal discretization, the CF was applied. The CF operator was discretized using the finite difference technique. The CF operator has been utilized for the first time in the cable model. The technique has order τ in time direction. Therefore, the ECBS with CF has reliable computational outcomes; the theoretical and numerical results indicate its compatibility with each other. The graphical representations also show that the results are very near to exact values.

Acknowledgments

The authors would like to thank Prince Mohammad Bin Fahd Center for Futuristic Studies and World Futures Studies Federation, Saudi Arabia, for their financial assistance.

  1. Author contributions: All authors contributed equally. All authors read and approved the final manuscript.

  2. Conflict of interest: The authors state no conflict of interest.

References

[1] E. J. López-Sáchez, J. M. Romero, and Y. Martínez, Fractional cable equation for general geometry: A model of axons with swellings and anomalous diffusion, Phys. Rev. E. 96 (2017), 032411, https://doi.org/10.1103/PhysRevE.96.032411. Search in Google Scholar PubMed

[2] Z. Wang, On Caputo type cable equation: Analysis and computation, Comput. Model. Eng. Sci. 123 (2020), 353–376, https://doi.org/10.32604/cmes.2020.08776. Search in Google Scholar

[3] J. Losada and J. J. Nieto, Properties of a new fractional derivative without singular Kernel, Prog. Fract. Differ. Appl. 1 (2015), no. 2, 87–92, https://doi.org/10.12785/pfda/010202. Search in Google Scholar

[4] M. Caputo and M. Fabrizio, A new definition of fractional derivative without singular Kernel, Prog. Fract. Differ. Appl. 1 (2015), no. 3, 73–85, https://doi.org/10.12785/pfda/010201. Search in Google Scholar

[5] A. Atangana and B. S. T. Alkahtani, Extension of the resistance inductance, capacitance electrical circuit of fractional derivative without singular kernel, Adv. Mech. Eng. 7 (2015), no. 6, 1–6, https://doi.org/10.1177/1687814015591937. Search in Google Scholar

[6] J. F. Gómez-Aguilar, M. G. Lopez-Lopez, V. M. Alvarado-Martınez, J. Reyes-Reyes, and M. Adam-Medina, Modeling diffusive transport with a fractional derivative without singular kernel, Phys. A 447 (2016), 467–481, https://doi.org/10.1016/j.physa.2015.12.066. Search in Google Scholar

[7] S. Djennadi, N. Shawagfeh, M. S. Osman, J. F. Gómez-Aguilar, and O. A. Abu Arqub, The Tikhonov regularization method for the inverse source problem of time fractional heat equation in the view of ABC-fractional technique, Phys. Scr. 96 (2021), 094006, https://doi.org/10.1088/1402-4896/ac0867. Search in Google Scholar

[8] S. Momani, B. Maayah, and O. AbuArqub, The reproducing kernel algorithm for numerical solution of Van der Poldamping model in view of the Atangana-Baleanu fractional approach, Fractals 28 (2020), no. 8, 2040010, https://doi.org/10.1142/S0218348X20400101. Search in Google Scholar

[9] X. J. Yang, Z. Z. Zhang, and H. M. Srivastava, Some new applications for heat and fluid flows via fractional derivatives without singular kernel, Therm. Sci. 20 (2016), no. 3, 833–839, https://doi.org/10.2298/TSCI16S3833Y. Search in Google Scholar

[10] O. AbuArqub, Numerical simulation of time-fractional partial differential equations arising in fluid flows via reproducing Kernel method, Int. J. Numer. Meth. Heat. Fluid 30 (2020), no. 11, 4711–4733, https://doi.org/10.1108/HFF-10-2017-0394. Search in Google Scholar

[11] A. Atangana, On the new fractional derivative and application to nonlinear Fisher’s reaction-diffusion equation, Appl. Math. Comput. 273 (2016), 948–956, https://doi.org/10.1016/j.amc.2015.10.021. Search in Google Scholar

[12] B. I. Henry, T. A. M. Langlands, and S. L. Wearne, Fractional cable models for spiny neuronal dendrites, Phys. Rev. Lett. 100 (2008), no. 12, 128103, https://doi.org/10.1103/PhysRevLett.100.128103. Search in Google Scholar PubMed

[13] X. Y. Liu, Y. P. Liu, and Z. W. Wu, Optimization of a fractal electrode-level charge transport model, Therm. Sci. 25 (2021), no. 3B, 2213–2220, https://doi.org/10.2298/TSCI200301108L. Search in Google Scholar

[14] D. D. Dai, T. T. Ban, and Y. L. Wang, The piecewise reproducing kernel method for the time variable fractional order advection-reaction-diffusion equations, Therm. Sci. 25 (2021), no. 2B, 1261–1268, https://doi.org/10.2298/TSCI200302021D. Search in Google Scholar

[15] K. L. Wang and S. W. Yao, He’s fractional derivative for the evolution equation, Therm. Sci. 24 (2020), no. 4, 2507–2513, https://doi.org/10.2298/TSCI2004507W. Search in Google Scholar

[16] X. L. Hu and L. M. Zhang, Implicit compact difference schemes for the fractional cable equation, Appl. Math. Model. 36 (2012), no. 9, 4027–4043, https://doi.org/10.1016/j.apm.2011.11.027. Search in Google Scholar

[17] B. Yu and X. Y. Jiang, Numerical identification of the fractional derivatives in the two-dimensional fractional cable equation, J. Sci. Comput. 68 (2016), 252–272, https://doi.org/10.1007/s10915-015-0136-y. Search in Google Scholar

[18] Y. Y. Zheng and Z. G. Zhao, The discontinuous Galerkin finite element method for fractional cable equation, Appl. Numer. Math. 115 (2017), 32–41, https://doi.org/10.1016/j.apnum.2016.12.006. Search in Google Scholar

[19] X. Yang, X. Y. Jiang, and H. Zhang, A time-space spectral tau method for the time fractional cable equation and its inverse problem, Appl. Numer. Math. 130 (2018), 95–111, https://doi.org/10.1016/j.apnum.2018.03.016. Search in Google Scholar

[20] Y. Liu, Y. W. Du, H. Li, and J. F. Wang, A two-grid finite element approximation for a nonlinear time-fractional cable equation, Nonlinear Dyn. 85 (2016), 2535–2548, https://doi.org/10.1007/s11071-016-2843-9. Search in Google Scholar

[21] Y. Liu, Y. W. Du, H. Li, and J. F. Wang, Some second-order θ schemes combined with finite element method for nonlinear fractional cable equation, Numer. Algorithms 80 (2019), 533–555, https://doi.org/10.1007/s11075-018-0496-0. Search in Google Scholar

[22] P. Zhu, S. Xie, and X. Wang, Non-smooth data error estimates for FEM approximations of the time fractional cable equation, Appl. Numer. Math. 121 (2017), 170–184, https://doi.org/10.1016/j.apnum.2017.07.005. Search in Google Scholar

[23] T. Akram, M. Abbas, A. Ali, A. Iqbal, and D. Baleanu, A numerical approach of a time fractional reaction-diffusion model with a non-singular kernel, Symmetry 12 (2020), no. 10, 1653, https://doi.org/10.3390/sym12101653. Search in Google Scholar

[24] T. Akram, M. Abbas, M. B. Riaz, A. I. Ismail,and N. M. Ali, Development and analysis of new approximation of extended cubic B-spline to the non-linear time fractional Klein-Gordon equation, Fractals 28 (2020), no. 8, 2040039, https://doi.org/10.1142/S0218348X20400393 Search in Google Scholar

[25] K. Zeynab and S. Habibollah, B-spline wavelet operational method for numerical solution of time-space fractional partial differential equations, Int. J. Wavelets Multiresolut Inf. Process 15 (2017), no. 4, 1750034, https://doi.org/10.1142/S0219691317500345. Search in Google Scholar

[26] B. P. Moghaddam and J. A.T. Machado, A stable three-level explicit spline finite difference scheme for a class of nonlinear time variable order fractional partial differential equations, Comput. Math. Appl. 36 (2017), no. 6, 1262–1269, https://doi.org/10.1016/j.camwa.2016.07.010. Search in Google Scholar

[27] X. Li, Numerical solution of fractional differential equations using cubic B-spline wavelet collocation method, Commun. Nonlinear Sci. Numer. Simul. 17 (2012), no. 10, 3934–3946, https://doi.org/10.1016/j.cnsns.2012.02.009. Search in Google Scholar

[28] K. M. Saad, M. M. Khader, J. F. Gómez-Aguilar, and D. Baleanu, Numerical solutions of the fractional Fisher’s type equations with Atangana-Baleanu fractional derivative by using spectral collocation methods, Chaos, 29 (2019), 023116, https://doi.org/10.1063/1.5086771. Search in Google Scholar PubMed

[29] T. Akram, M. Abbas, and A. I. Ismail, An extended cubic B-spline collocation scheme for time fractional sub-diffusion equation, AIP Confer. Proc. 2184 (2019), 060017, https://doi.org/10.1063/1.5136449. Search in Google Scholar

[30] T. Akram, M. Abbas, and A. I. Ismail, Numerical solution of fractional cable equation via extended cubic B-spline, AIP Confer. Proc. 2138 (2019), 030004, https://doi.org/10.1063/1.5121041. Search in Google Scholar

[31] N. Moshtaghi and A. Saadatmandi, Numerical solution of time fractional cable equation via the Sinc-Bernoulli collocation method, J. Appl. Comput. Mech. 7 (2021), no. 4, 1916–1924, https://doi.org/10.22055/JACM.2020.31923.1940. Search in Google Scholar

[32] L. Pezza and F. Pittoli, A fractional spline collocation-Galerkin method for the time-fractional diffusion equation, Commun Appl. Ind. Math. 9 (2018), 104–120, https://doi.org/10.1515/caim-2018-0007. Search in Google Scholar

[33] E. Y. Tabriz, M. Lakestani, and M. Razzaghi, Study of B-spline collocation method for solving fractional optimal control problems, Trans. Ins. Meas. Cont. 43 (2021), 2425–2437, https://doi.org/10.1177/0142331220987537 Search in Google Scholar

[34] C. De Boor, A Practical Guide to Splines, Springer-Verlag, New York, vol. 27, 1978, p. 325. 10.1007/978-1-4612-6333-3Search in Google Scholar

Received: 2022-04-08
Revised: 2022-06-08
Accepted: 2022-07-07
Published Online: 2022-09-23

© 2022 Azhar Iqbal and Tayyaba Akram, published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. On some summation formulas
  3. A study of a meromorphic perturbation of the sine family
  4. Asymptotic behavior of even-order noncanonical neutral differential equations
  5. Unconditionally positive NSFD and classical finite difference schemes for biofilm formation on medical implant using Allen-Cahn equation
  6. Starlike and convexity properties of q-Bessel-Struve functions
  7. Mathematical modeling and optimal control of the impact of rumors on the banking crisis
  8. On linear chaos in function spaces
  9. Convergence of generalized sampling series in weighted spaces
  10. Persistence landscapes of affine fractals
  11. Inertial iterative method with self-adaptive step size for finite family of split monotone variational inclusion and fixed point problems in Banach spaces
  12. Various notions of module amenability on weighted semigroup algebras
  13. Regularity and normality in hereditary bi m-spaces
  14. On a first-order differential system with initial and nonlocal boundary conditions
  15. On solving pseudomonotone equilibrium problems via two new extragradient-type methods under convex constraints
  16. Local linear approach: Conditional density estimate for functional and censored data
  17. Some properties of graded generalized 2-absorbing submodules
  18. Eigenvalue inclusion sets for linear response eigenvalue problems
  19. Some integral inequalities for generalized left and right log convex interval-valued functions based upon the pseudo-order relation
  20. More properties of generalized open sets in generalized topological spaces
  21. An extragradient inertial algorithm for solving split fixed-point problems of demicontractive mappings, with equilibrium and variational inequality problems
  22. An accurate and efficient local one-dimensional method for the 3D acoustic wave equation
  23. On a weighted elliptic equation of N-Kirchhoff type with double exponential growth
  24. On split feasibility problem for finite families of equilibrium and fixed point problems in Banach spaces
  25. Entire and meromorphic solutions for systems of the differential difference equations
  26. Multiplication operators on the Banach algebra of bounded Φ-variation functions on compact subsets of ℂ
  27. Mannheim curves and their partner curves in Minkowski 3-space E13
  28. Characterizations of the group invertibility of a matrix revisited
  29. Iterates of q-Bernstein operators on triangular domain with all curved sides
  30. Data analysis-based time series forecast for managing household electricity consumption
  31. A robust study of the transmission dynamics of zoonotic infection through non-integer derivative
  32. A Dai-Liao-type projection method for monotone nonlinear equations and signal processing
  33. Review Article
  34. Remarks on some variants of minimal point theorem and Ekeland variational principle with applications
  35. Special Issue on Recent Methods in Approximation Theory - Part I
  36. Coupled fixed point theorems under new coupled implicit relation in Hilbert spaces
  37. Approximation of integrable functions by general linear matrix operators of their Fourier series
  38. Sharp sufficient condition for the convergence of greedy expansions with errors in coefficient computation
  39. Approximation of conic sections by weighted Lupaş post-quantum Bézier curves
  40. On the generalized growth and approximation of entire solutions of certain elliptic partial differential equation
  41. Existence results for ABC-fractional BVP via new fixed point results of F-Lipschitzian mappings
  42. Linear barycentric rational collocation method for solving biharmonic equation
  43. A note on the convergence of Phillips operators by the sequence of functions via q-calculus
  44. Taylor’s series expansions for real powers of two functions containing squares of inverse cosine function, closed-form formula for specific partial Bell polynomials, and series representations for real powers of Pi
  45. Special Issue on Recent Advances in Fractional Calculus and Nonlinear Fractional Evaluation Equations - Part I
  46. Positive solutions for fractional differential equation at resonance under integral boundary conditions
  47. Source term model for elasticity system with nonlinear dissipative term in a thin domain
  48. A numerical study of anomalous electro-diffusion cells in cable sense with a non-singular kernel
  49. On Opial-type inequality for a generalized fractional integral operator
  50. Special Issue on Advances in Integral Transforms and Analysis of Differential Equations with Applications
  51. Mathematical analysis of a MERS-Cov coronavirus model
  52. Rapid exponential stabilization of nonlinear continuous systems via event-triggered impulsive control
  53. Novel soliton solutions for the fractional three-wave resonant interaction equations
  54. The multistep Laplace optimized decomposition method for solving fractional-order coronavirus disease model (COVID-19) via the Caputo fractional approach
  55. Special Issue on Problems, Methods and Applications of Nonlinear Analysis
  56. Some recent results on singular p-Laplacian equations
  57. Infinitely many solutions for quasilinear Schrödinger equations with sign-changing nonlinearity without the aid of 4-superlinear at infinity
  58. Special Issue on Recent Advances for Computational and Mathematical Methods in Scientific Problems
  59. Existence of solutions for a nonlinear problem at resonance
  60. Asymptotic stability of solutions for a diffusive epidemic model
  61. Special Issue on Computational and Numerical Methods for Special Functions - Part I
  62. Fully degenerate Bernoulli numbers and polynomials
  63. Wigner-Ville distribution and ambiguity function associated with the quaternion offset linear canonical transform
  64. Some identities related to degenerate Stirling numbers of the second kind
  65. Two identities and closed-form formulas for the Bernoulli numbers in terms of central factorial numbers of the second kind
  66. λ-q-Sheffer sequence and its applications
  67. Special Issue on Fixed Point Theory and Applications to Various Differential/Integral Equations - Part I
  68. General decay for a nonlinear pseudo-parabolic equation with viscoelastic term
  69. Generalized common fixed point theorem for generalized hybrid mappings in Hilbert spaces
  70. Computation of solution of integral equations via fixed point results
  71. Characterizations of quasi-metric and G-metric completeness involving w-distances and fixed points
  72. Notes on continuity result for conformable diffusion equation on the sphere: The linear case
Downloaded on 25.10.2025 from https://www.degruyterbrill.com/document/doi/10.1515/dema-2022-0155/html
Scroll to top button