Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 10 June 2020
Sec. Statistical and Computational Physics
This article is part of the Research Topic New Trends in Fractional Differential Equations with Real-World Applications in Physics View all 16 articles

Some Effective Numerical Techniques for Chaotic Systems Involving Fractal-Fractional Derivatives With Different Laws

  • 1Department of Engineering Science, Kermanshah University of Technology, Kermanshah, Iran
  • 2Department of Mathematics, Faculty of Engineering and Natural Sciences, Bahçeşehir University, Istanbul, Turkey
  • 3Department of Mathematics, College of Arts and Sciences, Prince Sattam Bin Abdulaziz University, Wadi Aldawaser, Saudi Arabia

Chaotic systems are dynamical systems that are highly sensitive to initial conditions. Such systems are used to model many real-world phenomena in science and engineering. The main purpose of this paper is to present several efficient numerical treatments for chaotic systems involving fractal-fractional operators. Several numerical examples test the performance of the proposed methods. Simulations with different values of the fractional and fractal parameters are also conducted. It is demonstrated that the fractal-fractional derivative enables one to capture all the useful information from the history of the phenomena under consideration. The numerical schemes can also be implemented for other chaotic systems with fractal-fractional operators.

1. Introduction

In recent decades numerical methods have been recognized as powerful mathematical techniques for solving nonlinear equations that model real-world problems [15]. Numerical techniques have been used to solve different classes of differential equations, including those of arbitrary order. Due to the outstanding contribution of nonlinear models to human understanding of many phenomena and the prediction of the future behavior of systems, many researchers are devoting their attention to developing new and reliable numerical techniques that could be applied to more complex cases. Chaotic behaviors are among the natural phenomena that have attracted the attention of many researchers, who aim to replicate and predict those behaviors. To achieve this, differential operators are often used as mathematical tools to construct the underlying models. Recently a new class of differential operators was introduced, which are convolutions of fractal derivatives and fractional kernels with different forms such as power law, exponential decay and the Mittag-Leffler function [6]. These differential operators are able to represent complexities that cannot be described with classical fractional differentiation and integration. One of the strengths of these operators is their two orders, where one is considered a fractional order and the other is the fractal dimension [79]. With these efficient operators, a new class of nonlinear differential and integral equations can be constructed, and existing chaotic models can be extended. Nevertheless, to verify the effectiveness of such fractal-fractional operators in modeling chaotic attractors, one needs to solve the models numerically as their exact solutions cannot be easily obtained by existing analytical methods. So far, a few numerical methods have been used to discretize such models, and some numerical results for chaotic attractors have been obtained. A different numerical scheme was suggested very recently [1027] and was found to be efficient for solving nonlinear differential equations. Some articles have examined the analysis of errors and the determination of error bands in the context of the possible deficit differential equations [2830]. Since the operators defined in this paper are novel in the field, the numerical methods associated with them are also very limited. One of the main motivations for this article was to introduce methods for solving fractal-fractional problems that have not been considered before. By using the proposed methods, approximate solutions to these problems can be determined more easily and with higher accuracy. In addition, the methods can be applied to real-world problems. In this work, we present applications of such numerical schemes in solving chaotic models that involve the new class of differential operators. We consider some well-known chaotic models with fractal-fractional differential operators, so that we can compare our results with those in the literature. The article is organized as follows. In section 2, we give a brief overview of some basic definitions of fractional differential calculus. Two efficient and effective numerical methods for determining approximate solutions to fractal-fractional problems are presented in section 3. The first is for the Caputo derivative and the second is for the Atangana-Baleanu-Caputo derivative. The kernels used in these two definitions are singular and non-singular, respectively. Several numerical simulations for chaotic systems are described in section 4. The results obtained are accurate, interesting, and meaningful. Finally, we present our overall conclusions.

2. Preliminary Definitions

In this section, we give a brief review of some existing definitions of fractal-fractional operators. These operators result from the combination of two important concepts: fractional differentiation and fractal derivatives. Most of the definitions given here are taken from Atangana [6] and Atangana and Qureshi [31].

Definition 1. If g(t) is a differentiable function on a finite open interval, then we define the fractal-fractional derivative of g(t) in the Caputo sense as

D0         FFCtρ,τg(t)=1Γ(kρ1)0tdg(ω)dtτ(tω)kρ1dω,                                      k1<ρk, 0<k1<τk,    (1)

where

dg(ω)dtτ=limtyg(t)g(ω)tβωβ.    (2)

Definition 2. If g(t) is a differentiable function on a finite open interval, then we define the fractal-fractional derivative of g(t) in the Caputo-Fabrizio sense as

D0           FFCFtρ,τg(t)=L(ρ)1ρ0tdg(ω)dtτexp[ρ1ρ(tω)]dω,                                       n1<ρ,τn,    (3)

provided K(0) = K(1) = 1.

Definition 3. If g(t) is a differentiable function on a finite open interval, then we define the fractal-fractional derivative of g(t) in the Atangana-Baleanu sense as

D0            FFABtρ,τg(t)=AB(ρ)1ρ0tdg(ω)dtτEρ[ρ1ρ(tω)]dω,                                         n1<ρ,τn,    (4)

where Eρ(.) is the Mittag-Leffler function, defined by

Eρ(ω)=k=0ωkΓ(ρk+1),   ρ>0.    (5)

This function is an essential function in the modeling of physical processes using fractional calculus concepts. We know that Equation (5) reduces to the exponential function ex if one takes ρ = 1. Also, AB(·) is a function used for normalization, and it satisfies the property AB(0) = AB(1) = 1. One of the most popular definitions of AB(·) is

AB(ρ)=1-ρ+ρΓ(ρ).

Definition 4. If g(t) is a differentiable function on a finite open interval, then we define the fractal-fractional integral of g(t) in the Caputo sense [31] as

I0        FFCtρ,τg(t)=τΓ(ρ)0tωτ1g(ω)dω(tω)1ρ.    (6)

Definition 5. If g(t) is a differentiable function on a finite open interval, then we define the fractal-fractional integral of g(t) in the Caputo-Fabrizio sense [31] as

I0           FFCFtρ,τg(t)=τρM(ρ)0tg(ω)dωω1ρ+τ(1ρ)tτ1g(t)M(ρ).    (7)

Definition 6. If g(t) is a differentiable function on a finite open interval, then we define the fractal-fractional integral of g(t) in the Atangana-Baleanu sense [31] as

I0            FFABtρ,τg(t)=τρAB(ρ)0tωτ1g(ω)dω(tω)1ρ+τ(1ρ)tτ1g(t)AB(ρ).    (8)

3. The Proposed Numerical Methods

In what follows, the main aim is to construct two equations involving fractal-fractional derivatives,

D0  FFtρ,τξ(t)=N(t,ξ(t)),   t[t0,T],    (9)

with the initial condition ξ(0) = ξ0.

3.1. The Caputo Fractal-Fractional Derivative

From the results in Atangana and Qureshi [31], Equation (9) reduces to the following Caputo fractional representation:

D0Ctρξ(t)=τtτ1N(t,ξ(t)).    (10)

Now, taking into account the fundamental theorem of calculus, one gets

ξ(t)-ξ(t0)=τΓ(ρ)0tωτ-1N(ω,ξ(ω))dω(t-ω)1-ρ.    (11)

Then, inserting t = tn = t0 + nΔt into (17) leads to the following expression:

ξ(tn)=ξ(t0)+τΓ(ρ)i=0n1titi+1ωτ1N(ω,ξ(ω))dω(tnω)1ρ,              1nN,    (12)

where Δt=T-t0N is the time step for the discretization points.

Next, taking the linear Lagrange interpolation into account for the function of f(ω) = ωτ−1ℕ(ω, ξ(ω)), we obtain

ωτ-1N(ω,ξ(ω))ti+1τ-1N(ti+1,ξi+1)+ω-ti+1Δt(ti+1τ-1N(ti+1,ξi+1)-tiτ-1N(ti,ξi)),ω[ti,ti+1],    (13)

where ξi = ξ(ti). At this point, the result obtained in formula (13) can be used in relation (12). This substitution results in the following relationship:

ξn=ξ0+τΔtρ(ξnt0τ-1N(t0,ξ0)+i=0nðn-itiτ-1N(ti,ξi)),    (14)

where

ξn=(n1)ρ+1nρ(nρ1)Γ(ρ+2),ðm={1Γ(ρ+2),j=0,(m1)ρ+12mρ+1+(m+1)ρ+1Γ(ρ+2),m=1,2,,n1.    (15)

3.2. The Atangana-Baleanu Fractal-Fractional Derivative

In this case we can convert Equation (9) to the following Atangana-Baleanu fractional form [31]:

D0ABtρξ(t)=τtτ1N(t,ξ(t)).    (16)

Upon applying the integral operator to both sides of Equation (16), the following Volterra integral equation is constructed:

ξ(t)ξ(t0)=ττρAB(ρ)tτ1N(t,ξ(t))                          +τρAB(ρ)Γ(ρ)0t(tω)ρ1ωτ1N(ω,ξ(ω))dω.    (17)

Taking t = tn = t0 + nΔt in (17), we have

ξ(tn)=ξ(t0)+ττρAB(ρ)tnτ1N(tn,ξ(tn))             +τρAB(ρ)Γ(ρ)i=0n1titi+1(tnω)ρ1ωτ1N(ω,ξ(ω))dω.    (18)

Now, substituting (13) into (18), we get the following implicit Atangana-Baleanu-Caputo scheme:

ξn=ξ0+ττρAB(ρ)tnτ1N(tn,ξn)      +τρΔtρAB(ρ)(ξnt0τ1N(t0,ξ0)+i=0nðnitiτ1N(ti,ξi)),    (19)

where ξn and dj are the coefficients defined in (15). A closer look at relations (14) and (19) shows that these expressions are implicit equations for determining yn. The Newton iteration method is one of the most popular and efficient techniques for solving such problems. In this article, as in Garrappa [32], Ghanbari and Kumar [33], and Ghanbari et al. [34], we will use the Newton method to solve these equations. By solving these equations, approximate solutions to the original problem will be determined.

4. Numerical Simulations for Some Chaotic Systems

In this section, to show the validity of the proposed numerical technique, three chaotic systems involving Atangana-Baleanu-Caputo fractional derivatives are considered.

Example 1. Consider the following modified cyclically symmetric Thomas attractor given by [35]

D0FFtρ,τx(t)=sin(exp(y(t)))Bx(t),D0FFtρ,τy(t)=cos(sin(z(t)))By(t),D0FFtρ,τz(t)=exp(cos(x(t)))Bz(t).    (20)

Now we take into account the iterative methods given by (14) and (19) to solve (20). Taking (14) into account, we get the iterative structure

xn=x0+τΔtρ(ξnt0τ1[sin(exp(y0))Bx0]       +i=0nðnitiτ1[sin(exp(yi))Bxi]),yn=y0+τΔtρ(ξnt0τ1[cos(sin(sin(z0))By0]       +i=0nðnitiτ1[cos(sin(zi)Byi]),zn=z0+τΔtρ(ξnt0τ1[exp(cos(x0))Bz0]       +i=0nðnitiτ1[exp(cos(xi))Bzi]).    (21)

Moreover, from (19) we get

xn=x0+ττρAB(ρ)tnτ1[sin(exp(yn))Bxn]         +τρΔtρAB(ρ)(ξnt0τ1[sin(exp(y0))Bx0]      +i=0nðnitiτ1[sin(exp(yi))Bxi]),yn=y0+ττρAB(ρ)tnτ1[cos(sin(zn))Byn]         +τρΔtρAB(ρ)(ξnt0τ1[cos(sin(z0))By0]       +i=0nðnitiτ1[cos(sin(zi))Byi]),zn=z0+ττρAB(ρ)tnτ1[exp(cos(xn))Bzn]         +τρΔtρAB(ρ)(ξnt0τ1[exp(cos(x0))Bz0]       +i=0nðnitiτ1[exp(cos(xi))Bzi]).    (22)

In Figures 14 we plot the results of numerical simulations using the two iterative schemes (21) and (22) with B = 0.2. For the numerical implementations we took (x0, y0, z0) = (0, 3, 11) as the initial condition. The results of the iterative scheme (21) are plotted in Figures 58. In each figure we have taken a fixed value of τ and different values of ρ. In these experiments we set Δt = 10−3 and T = 800.

FIGURE 1
www.frontiersin.org

Figure 1. Evolution of the chaotic system (20) for τ = 0.93 and different values of ρ, using the iterative scheme (21).

FIGURE 2
www.frontiersin.org

Figure 2. Evolution of the chaotic system (20) for τ = 0.97 and different values of ρ, using the iterative scheme (21).

FIGURE 3
www.frontiersin.org

Figure 3. Evolution of the chaotic system (20) for τ = 0.93 and different values of ρ, using the iterative scheme (22).

FIGURE 4
www.frontiersin.org

Figure 4. Evolution of the chaotic system (20) for τ = 0.97 and different values of ρ, using the iterative scheme (22).

FIGURE 5
www.frontiersin.org

Figure 5. Evolution of the chaotic system (23) for τ = 0.96 and different values of ρ, using the iterative scheme (24).

FIGURE 6
www.frontiersin.org

Figure 6. Evolution of the chaotic system (23) for τ = 0.98 and different values of ρ, using the iterative scheme (24).

FIGURE 7
www.frontiersin.org

Figure 7. Evolution of the chaotic system (23) for τ = 0.96 and different values of ρ, using the iterative scheme (25).

FIGURE 8
www.frontiersin.org

Figure 8. Evolution of the chaotic system (23) for τ = 0.98 and different values of ρ, using the iterative scheme (25).

Example 2. In this example we consider the chaotic system

D0FFtρ,τx(t)=x(t)z(t)Bx(t)Dy(t),D0FFtρ,τy(t)=Dx(t)+y(t)z(t)By(t),D0FFtρ,τz(t)=C+Az(t)z(t)33x(t)2+Ez(t)x(t)3,    (23)

To solve this chaotic system, we use the following iterative scheme obtained from (14):

xn=x0+τΔtρ(ξnt0τ1[x0z0Bx0Dy0]       +i=0nðnitiτ1[xiziBxiDyi]),yn=y0+τΔtρ(ξnt0τ1[Dx0+y0z0By0]       +i=0nðnitiτ1[Dxi+yiziByi]),zn=z0+τΔtρ(ξnt0τ1[C+Az0z033x02+Ez0x03]                 +i=0nðnitiτ1[C+Azizi33xi2+Ezixi3]).    (24)

Using (19) one gets

xn=x0+ττρAB(ρ)tnτ1[xnznBxnDyn]         +τρΔtρAB(ρ)(ξnt0τ1[x0z0Bx0Dy0]      +i=0nðnitiτ1[xiziBxiDyi]),yn=y0+ττρAB(ρ)tnτ1[Dxn+ynznByn]         +τρΔtρAB(ρ)(ξnt0τ1[Dx0+y0z0By0]       +i=0nðnitiτ1[Dxi+yiziByi]),zn=z0+ττρAB(p)tnτ1[C+Aznzn33xn2+Eznxn3]         +τρΔtρAB(ρ)(ξnt0τ1[C+Az0z033x02+Ez0x03]                      +i=0nðnitiτ1[C+Azizi33xi2+Ezixi3]).    (25)

The approximate solutions obtained from (24) are plotted in Figures 5, 6, and those obtained from (25) are plotted in Figures 7, 8. The parameter values used in the model are A = 0.95, B = 0.7, C = 0.6, D = 3.5, and E = 0.1. We used the initial guess (x0, y0, z0) = (0.1, 0, 0) for different fractional orders of ρ and τ.

Example 3. Consider the chaotic system given by [31]

D0FFtρ,τx(t)=A(y(t)x(t)),D0FFtρ,τy(t)=(CA)x(t)[S0z(t)S1sin(z(t))]x(t)+Cy(t),D0FFtρ,τz(t)=x(t)y(t)Bz(t),    (26)

For this system, from (14) we obtain the following iterative scheme:

xn=x0+τΔtρ(ξnt0τ1[A(y0x0)]+i=0nðnitiτ1[A(yixi)]),yn=y0+τΔtρ(ξnt0τ1[(CA)x0[S0z0S1sin(z0)]x0+Cy0]                 +i=0nðnitiτ1[(CA)x0[S0ziS1sin(zi)]xi+Cyi]),zn=z0+τΔtρ(ξnt0τ1[x0y0Bz0]+i=0nðnitiτ1[xiyiBzi]).    (27)

Moreover, from (19), the approximate solution is obtained as

xn=x0+ττρAB(ρ)tnτ1[A(ynxn)]      +τρΔtρAB(ρ)(ξnt0τ1[A(y0x0)]+i=0nðnitiτ1[A(yixi)]),yn=y0+ττρAB(ρ)tnτ1[(CA)x0[S0znS1sin(zn)]xn+Cyn]     +τρΔtρAB(ρ)(ξnt0τ1[(CA)x0[S0z0S1sin(z0)]x0+Cy0]           +i=0nðnitiτ1[(CA)x0[S0ziS1sin(zi)]xi+Cyi]),zn=z0+ττρAB(ρ)tnτ1[xnynBzn]      +τρΔtρAB(ρ)(ξnt0τ1[x0y0Bz0]+i=0nðnitiτ1[xiyiBzi]).    (28)

Figures 912 show the portraits corresponding to the chaotic system in (26) obtained using (27) and (28). The parameters for the model are A = 50, B = 2, C = 30, S1 = 1, and S2 = 20. We used the initial guess (x0, y0, z0) = (2, 1, 1) for different fractional orders of ρ.

FIGURE 9
www.frontiersin.org

Figure 9. Evolution of the chaotic system (26) for τ = 0.7 and different values of ρ, using the iterative scheme (27).

FIGURE 10
www.frontiersin.org

Figure 10. Evolution of the chaotic system (26) for τ = 0.9 and different values of ρ, using the iterative scheme (27).

FIGURE 11
www.frontiersin.org

Figure 11. Evolution of the chaotic system (26) for τ = 0.7 and different values of ρ, using the iterative scheme (28).

FIGURE 12
www.frontiersin.org

Figure 12. Evolution of the chaotic system (26) for τ = 0.9 and different values of ρ, using the iterative scheme (28).

5. Conclusion

Although many numerical methods are available, the development of new efficient numerical schemes has always been one of the most important concerns in applied mathematics and engineering. A foremost reason for the widespread interest in new numerical methods is that they may reveal new facts about real-world phenomena. This paper has presented some efficient approximate methods for solving chaotic systems that use new definitions for the derivative, called fractal-fractional derivatives. The concept of memory is one of the most important features of these types of derivatives. With this valuable feature, the evolution of the phenomena modeled by such systems can be more accurately predicted. The proposed new techniques have been tested by using them to solve several important practical problems. Application of the methods to these problems revealed very interesting behaviors of the systems that have meaningful interpretations. The numerical methods presented in this article have the potential to be used for solving similar models. Since any new numerical method should be validated in terms of convergence, stability and consistency of solutions, these are important research directions left to future work.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

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.

References

1. Kilbas AAA, Srivastava HM, Trujillo JJ. Theory and Applications of Fractional Differential Equations. Vol. 204. Amsterdam; Boston, MA: Elsevier Science Limited (2006).

Google Scholar

2. Caputo M, Fabrizio M. A new definition of fractional derivative without singular kernel. Progr Fract Differ Appl. (2015) 1:1–13. doi: 10.18576/pfda/020101

CrossRef Full Text | Google Scholar

3. Atangana A, Baleanu D. New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model. arXiv preprint arXiv:160203408. (2016). doi: 10.2298/TSCI160111018A

CrossRef Full Text | Google Scholar

4. Gao W, Ghanbari B, Baskonus HM. New numerical simulations for some real world problems with Atangana-Baleanu fractional derivative. Chaos Solit Fract. (2019) 128:34–43. doi: 10.1016/j.chaos.2019.07.037

CrossRef Full Text | Google Scholar

5. Ghanbari B, Gómez-Aguilar J. Modeling the dynamics of nutrient-phytoplankton-zooplankton system with variable-order fractional derivatives. Chaos Solit Fract. (2018) 116:114–20. doi: 10.1016/j.chaos.2018.09.026

CrossRef Full Text | Google Scholar

6. Atangana A. Fractal-fractional differentiation and integration: connecting fractal calculus and fractional calculus to predict complex system. Chaos Solit Fract. (2017) 102:396–406. doi: 10.1016/j.chaos.2017.04.027

CrossRef Full Text | Google Scholar

7. Losada J, Nieto JJ. Properties of a new fractional derivative without singular kernel. Progr Fract Differ Appl. (2015) 1:87–92.

Google Scholar

8. Chen W, Sun H, Zhang X, Korošak D. Anomalous diffusion modeling by fractal and fractional derivatives. Comput Math Appl. (2010) 59:1754–8. doi: 10.1016/j.camwa.2009.08.020

CrossRef Full Text | Google Scholar

9. Kanno R. Representation of random walk in fractal space-time. Phys A. (1998) 248:165–75. doi: 10.1016/S0378-4371(97)00422-6

CrossRef Full Text | Google Scholar

10. Singh J, Kumar D, Baleanu D, Rathore S. An efficient numerical algorithm for the fractional Drinfeld-Sokolov-Wilson equation. Appl Math Comput. (2018) 335:12–24. doi: 10.1016/j.amc.2018.04.025

CrossRef Full Text | Google Scholar

11. Qureshi S, Yusuf A, Shaikh AA, Inc M, Baleanu D. Fractional modeling of blood ethanol concentration system with real data application. Chaos. (2019) 29:013143. doi: 10.1063/1.5082907

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Gill V, Singh J, Singh Y. Analytical solution of generalized space-time fractional advection-dispersion equation via coupling of Sumudu and Fourier transforms. Front Phys. (2019) 6:151. doi: 10.3389/fphy.2018.00151

CrossRef Full Text | Google Scholar

13. Qureshi S, Atangana A, Shaikh AA. Strange chaotic attractors under fractal-fractional operators using newly proposed numerical methods. Eur Phys J Plus. (2019) 134:523. doi: 10.1140/epjp/i2019-13003-7

CrossRef Full Text | Google Scholar

14. Kumar S, Kumar R, Singh J, Nisar K, Kumar D. An efficient numerical scheme for fractional model of HIV-1 infection of CD4+ T-cells with the effect of antiviral drug therapy. Alexandria Eng J. (2020). doi: 10.1016/j.aej.2019.12.046. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Berhe HW, Qureshi S, Shaikh AA. Deterministic modeling of dysentery diarrhea epidemic under fractional Caputo differential operator via real statistical analysis. Chaos Solit Fract. (2020) 131:109536. doi: 10.1016/j.chaos.2019.109536

CrossRef Full Text | Google Scholar

16. Singh J. A new analysis for fractional rumor spreading dynamical model in a social network with Mittag-Leffler law. Chaos. (2019) 29:013137. doi: 10.1063/1.5080691

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Qureshi S, Rangaig NA, Baleanu D. New numerical aspects of Caputo-Fabrizio fractional derivative operator. Mathematics. (2019) 7:374. doi: 10.3390/math7040374

CrossRef Full Text | Google Scholar

18. Srivastava H, Kumar D, Singh J. An efficient analytical technique for fractional model of vibration equation. Appl Math Model. (2017) 45:192–204. doi: 10.1016/j.apm.2016.12.008

CrossRef Full Text | Google Scholar

19. Qureshi S. Effects of vaccination on measles dynamics under fractional conformable derivative with Liouville-Caputo operator. Eur Phys J Plus. (2020) 135:63. doi: 10.1140/epjp/s13360-020-00133-0

CrossRef Full Text | Google Scholar

20. Solís-Pérez J, Gómez-Aguilar J, Atangana A. Novel numerical method for solving variable-order fractional differential equations with power, exponential and Mittag-Leffler laws. Chaos Solit Fract. (2018) 114:175–85. doi: 10.1016/j.chaos.2018.06.032

CrossRef Full Text | Google Scholar

21. Owolabi KM. Computational study of noninteger order system of predation. Chaos. (2019) 29:013120. doi: 10.1063/1.5079616

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Yavuz M, Bonyah E. New approaches to the fractional dynamics of schistosomiasis disease model. Phys A. (2019) 525:373–93. doi: 10.1016/j.physa.2019.03.069

CrossRef Full Text | Google Scholar

23. Jajarmi A, Ghanbari B, Baleanu D. A new and efficient numerical method for the fractional modeling and optimal control of diabetes and tuberculosis co-existence. Chaos. (2019) 29:093111. doi: 10.1063/1.5112177

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Ghanbari B, Gómez-Aguilar J. Analysis of two avian influenza epidemic models involving fractal-fractional derivatives with power and Mittag-Leffler memories. Chaos. (2019) 29:123113. doi: 10.1063/1.5117285

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Ganji R, Jafari H, Baleanu D. A new approach for solving multi variable orders differential equations with Mittag-Leffler kernel. Chaos Solit Fract. (2020) 130:109405. doi: 10.1016/j.chaos.2019.109405

CrossRef Full Text | Google Scholar

26. Kadkhoda N, Jafari H. An analytical approach to obtain exact solutions of some space-time conformable fractional differential equations. Adv Differ Equat. (2019) 2019:428. doi: 10.1186/s13662-019-2349-0

CrossRef Full Text | Google Scholar

27. Jafari H, Babaei A, Banihashemi S. A novel approach for solving an inverse reaction-diffusion-convection problem. J Optim Theory Appl. (2019) 183:688–704. doi: 10.1007/s10957-019-01576-x

CrossRef Full Text | Google Scholar

28. Liu B, Hill DJ, Chen G. Synchronization errors and uniform synchronization with an error bound for chaotic systems. Int J Bifur Chaos. (2008) 18:3341–54. doi: 10.1142/S021812740802241X

CrossRef Full Text | Google Scholar

29. Sorrentino F, Barlev G, Cohen AB, Ott E. The stability of adaptive synchronization of chaotic systems. Chaos. (2010) 20:013103. doi: 10.1063/1.3279646

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kazemi BF, Jafari H. Error estimate of the MQ-RBF collocation method for fractional differential equations with Caputo-Fabrizio derivative. Math Sci. (2017) 11:297–305. doi: 10.1007/s40096-017-0232-2

CrossRef Full Text | Google Scholar

31. Atangana A, Qureshi S. Modeling attractors of chaotic dynamical systems with fractal-fractional operators. Chaos Solit Fract. (2019) 123:320–337. doi: 10.1016/j.chaos.2019.04.020

CrossRef Full Text | Google Scholar

32. Garrappa R. Numerical solution of fractional differential equations: a survey and a software tutorial. Mathematics. (2018) 6:16. doi: 10.3390/math6020016

CrossRef Full Text | Google Scholar

33. Ghanbari B, Kumar D. Numerical solution of predator-prey model with Beddington-DeAngelis functional response and fractional derivatives with Mittag-Leffler kernel. Chaos. (2019) 29:063103. doi: 10.1063/1.5094546

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Ghanbari B, Kumar S, Kumar R. A study of behaviour for immune and tumor cells in immunogenetic tumour model with non-singular fractional derivative. Chaos Solit Fract. (2020) 133:109619. doi: 10.1016/j.chaos.2020.109619

CrossRef Full Text | Google Scholar

35. Chepyzhov VV, Vishik MI. Attractors for Equations of Mathematical Physics. Vol. 49. American Mathematical Society (2002). doi: 10.1090/coll/049

CrossRef Full Text | Google Scholar

Keywords: chaotic attractors, computational efficiency, fractal-fractional operators, thomas attractor, Newton's method, product integration rule

Citation: Ghanbari B and Nisar KS (2020) Some Effective Numerical Techniques for Chaotic Systems Involving Fractal-Fractional Derivatives With Different Laws. Front. Phys. 8:192. doi: 10.3389/fphy.2020.00192

Received: 19 March 2020; Accepted: 30 April 2020;
Published: 10 June 2020.

Edited by:

Jordan Yankov Hristov, University of Chemical Technology and Metallurgy, Bulgaria

Reviewed by:

Sania Qureshi, Mehran University of Engineering and Technology, Pakistan
Hossein Jaferi, University of South Africa, South Africa

Copyright © 2020 Ghanbari and Nisar. 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: Kottakkaran Sooppy Nisar, n.sooppy@psau.edu.sa; ksnisar1@gmail.com

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.