Complexity
Volume 2020 (2020), Article ID 2371019, 15 pages
https://doi.org/10.1155/2020/2371019

Research Article

On the Fractional View Analysis of Keller–Segel Equations with Sensitivity Functions

Haobin Liu,1,2 Hassan Khan,³ Rasool Shah,³ A. A. Alderremy,⁴ Shaban Aly,⁵ and Dumitru Baleanu6,7

¹Data Recovery Key Laboratory of Sichuan Province, Neijiang Normal University, Neijiang 641000, Sichuan, China
²School of Mathematics and Information Sciences, Neijiang Normal University, Neijiang 641000, Sichuan, China
³Department of Mathematics, Abdul Wali Khan University Mardan (AWKUM), Mardan, Pakistan
⁴Department of Mathematics, Faculty of Science, King Khalid University, Abha 61413, Saudi Arabia
⁵Department of Mathematics, Faculty of Science, AL-Azhar University, Assiu, 71516, Egypt
⁶Department of Mathematics, Faculty of Arts and Sciences, Cankaya University, Ankara 06530, Turkey
⁷Institute of Space Sciences, Magurele-Bucharest, Romania

Correspondence should be addressed to Hassan Khan; hassanmath@awkum.edu.pk

Received 2 April 2020; Revised 22 June 2020; Accepted 1 July 2020; Published 28 July 2020

Academic Editor: Dimitri Volchenkov

Copyright © 2020 Haobin Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In this paper, the fractional view analysis of the Keller–Segal equations with sensitivity functions is presented. The Caputo operator has been used to pursue the present research work. The natural transform is combined with the homotopy perturbation method, and a new scheme for implementation is derived. The modified established method is named as the homotopy perturbation transform technique. The derived results are compared with the solution of the Laplace Adomian decomposition technique by using the systems of fractional Keller–Segal equations. The solution graphs and the table have shown that the obtained results coincide with the solution of the Laplace Adomian decomposition method. Fractional-order solutions are determined to confirm the reliability of the current method. It is observed that the solutions at various fractional orders are convergent to an integer-order solution of the problems. The suggested procedure is very attractive and straight forward and therefore can be modified to solve high nonlinear fractional partial differential equations and their systems.

1. Introduction

Fractional differential equations (FDEs) are the generalizations of the standard integer-order differential equations. FDEs have gained much attention in the recent decades, as they are broadly used in several areas to analyze different processes, such as the processing of the signal, control theory, identification of the system, fluid flow, biomathematics, and other fields [1–3]. FDEs have been implemented in the modeling of the description of unification of diffusion, fractional random walk, and systems of both diffusive and subdiffusive. It is also an investigation that in many physical phenomena, experiments have proved that the solution of fractional-order DEs is in good agreement with experimental data of any physical phenomena to integer-order DEs. For example, the noninteger differential equations (DEs) can model some physical phenomena more effectively such as delineates memory, heredity effects, properties of different materials, and internal friction process [4–9].

The importance of FDEs is found in the literature because it can model most of the physical phenomena in science and engineering more accurately as compared to integer-order models [10–12], and therefore the researchers have shown much interest to study fractional calculus and FDEs during the last decades. Because of the significance of FDEs, the mathematicians are working to develop many useful and powerful numerical and analytical techniques to determine the actual or approximate solutions of the targeted problems [13–15]. Many other powerful and efficient techniques have been suggested to obtain the exact or analytical solutions for FDEs. For instance, these techniques include the Laplace Adomian decomposition method (LADM) [16–18], Chebyshev wavelet methods (CWM) [19], collocation-shooting method [20], power series methods (PSM) [21], fractional Bernstein polynomials along with shooting method [22], fractional-Legendre spectral Galerkin method [23], variational iterative transform method (VITM) [24], homotopy perturbation transform method (HPTM) [25–27], homotopy analysis transform method (HATM) [28, 29], reduced differential transform method (RDTM) [30, 31], finite element technique (FET) [32], finite difference technique (FDT) [33], and q-homotopy analysis transform method (q-HATM) [34]. Based on these techniques, a wide range of FDEs have been analyzed.

The most important cell motion which can be used as a response to the gradient of a chemical compound is known as chemotaxis. It has an important contribution in the population cell number and growth of biology. E. Keller and L. Segal suggested the first mathematical equation of chemotaxis in 1970. They displayed parabolic schemes to describe the method of aggregating a cell slime mold depending on a molecule attraction [35]. In this article, we considered the coupled time-fractional Keller–Segel equation [36–39]:

With initial conditions,

where the unidentified function defines the concentrations of amoebae, the chemical concentration is denoted by , and the chemotactic is defined by the partial derivative of , is the sensitivity function. Some other types of sensitivity functions can be observed in the coupled time-fractional K-S chemotactic model in [40] and also in the K-S model with a logic sensitivity function and small diffusivity [41]. Also, singular sensitivity can be seen in the system of the two-dimensional K-S system in [36, 42, 43]. The K-S model has been extensively analyzed recently. For example, Kamel et.al. have used LADM for the solution of the S-K equation [37]. Bournaveas [38] gave the one-dimensional Keller–Segel model with a fractional cell diffusion. Zayernouri [39] formed a fractional class of implicit Adams–Moulton and explicit Adams–Bashforth techniques.

The HPM is mixed with the natural transform method to generate a highly effective method to handle the solution of several nonlinear problems and is known as the homotopy perturbation transform method (HPTM). HPTM generates a convergent series form solution that converges to the exact solution of the problems and provides closed-form solutions. The proposed method can combine two important techniques to achieve an effective solution for nonlinear equations. Gorbhani first proposed the use of He’s polynomial in the nonlinear terms [44, 45]. It should be remembered that the HPTM is implemented without any flexibility or restrictive assumptions or transforms and round off free error. Many authors have solved fractional-order diffusion equations [46], partial differential equations [47, 48], and wave-like equations [49] by using linear and nonlinear problems with the help of HPTM.

In the current research paper, HPTM is implemented to solve fractional-order Keller–Segel equations. The solutions achieved through the suggested technique are straightforward and simple. The quality of the current method is appropriate to provide the analytical results to the given examples. The HPTM solutions are shown to be in close contact with the solutions of other existing techniques.

2. Preliminary Concepts

Definition 1. Let and , and then the Riemann–Liouville integral of order is given by [1–3]

Definition 2. Let and , and then the Caputo definition of the fractional derivative of order if with is described as [1–3]

Remark 1. Some basic properties are given below [1–3]:

Definition 3. The natural transform of the function for is defined as [50]

Theorem 1. The natural transform of the fractional derivative of the function in the Riemann–Liouville sense is given as [50]

Theorem 2. The natural transform of the fractional derivative of the function in the Caputo sense is given as [50]

3. Implementation of HPTM [51]

To explain the fundamental concept of this technique, we consider a particular fractional-order nonlinear partial differential nonhomogeneous equation:

With initial conditions,

where the fractional derivative in equation (9) is defined in the Caputo sense. The operators and describe the linear and nonlinear operators, respectively, and is the source term.

Using natural transformation in equation (9), we get [51]

With the help of the fractional derivative natural property, we have

Taking the inverse natural transformation of equation (12), we get

Using the HPTM procedure, the solution is expressed as

The nonlinear term can be decomposed as

Few He’s polynomials [35, 40] are described by

Utilizing equations (14) and (15) in equation (13), we have

Using He’s polynomials in HPTM and comparing the coefficient of power p, we get

The remaining of the components can be totally obtained, and the sequence result is thus fully determined.

Finally, we calculate the analytical result of :

Generally, the solutions of the above series converge very quickly.

4. Numerical Examples

Example 1. The Keller–Segal equations with a fractional derivative are given as [37–39]

Subjecting to the initial solutions, we get

Taking natural transformation of equation (20), we get

Simplifying the above equation, we get

Using inverse natural transform, we have

Then, the chemotactic term .
Now, implementing HPM, we get

By comparing the coefficients of powers p, we get

The series form of equation (20) is given by

The HPTM solution is given as

The series solutions are obtained by using the numerical values , , , , , and :

Example 2. The Keller–Segal equations with a fractional derivative are given as [37–39]

With the initial conditions:

Applying natural transformation to equation (30), we get

Simplifying the above equation, we get

Using inverse natural transform, we have

Now, implementing HPM, we get

where and are He’s polynomials that show the nonlinear terms. So, representing He’s polynomials for , we find that

For , we find that

The series solutions are obtained by using the numerical values , , , , , and .
By comparing the coefficients of powers p, we get

The series form of equation (30) is given by

The HPTM solution is given as

5. Results and Discussion

In Figures 1 and 2, the comparative study of HPTM and LADM solutions has been made for variables and at for Example 1. The graphical representation has shown the close relation between HPTM and LADM solutions. In Figures 3 and 4, the graphs of and verses are plotted for fixed for both fractional and integer orders of Example 2. The graphical representation has confirmed the convergence of fractional-order solutions towards integer-order solutions. In Table 1, the LADM and HPTM solutions of Example 1 are compared at and . It is observed that LADM and HPTM solutions are identical and justify the reliability of the proposed techniques. In Figures 5 and 6, the HPTM solutions, and in Figure 7, the HPTM solutions at fractional orders , and 1, are discussed for variables and . It is investigated that both the methods have a higher degree of accuracy and provide the closed-form solution to Example 2.

(a)

(b)

Figure 1: HPTM solutions of and at .

(a)

(b)

Figure 2: LADM solutions of and at .

(a)

(b)

Figure 3: HPTM solutions of and at .

(a)

(b)

Figure 4: HPTM solutions of and at different fractional orders and .

Table 1: Comparison of HDM [36], LADM [37], and HPTM at and .

HDM
HDM
LADM
LADM
HPTM
HPTM

0.2
2.70
8.0
2.70
8.0
2.70
8.0

0.4
3.0
2.0
3.0
2.0
3.0
2.0

0.6
5.90
4.0
5.90
4.0
5.90
4.0

0.8
8.50
7.0
8.50
7.0
8.50
7.0

1.00
1.50
1.20
1.50
1.20
1.50
1.20

(a)

(b)

Figure 5: HPTM solutions of and at .

(a)

(b)

Figure 6: HPTM solutions of and at .

(a)

(b)

Figure 7: HPTM solutions of and at different fractional orders and .

6. Conclusion

In the current article, an effective technique which is known as the Laplace homotopy transform method is implemented to solve the systems of Keller–Segal equations for both fractional and integer orders of the derivatives within the Caputo operator. Two numerical examples of fractional Keller–Segal equations are presented to verify the reliability of the suggested method. The graphical and tabular representation has confirmed that the derived results are in close agreement with the solution of the Laplace Adomian decomposition method. Moreover, the current technique needs very small calculations and has a higher degree of accuracy for the targeted problems. In conclusion, the present technique is found to be an accurate and effective analytical technique to solve high nonlinear fractional systems of partial differential equations.

Data Availability

The numerical data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University, Saudi Arabia, for funding this work through the research groups program under the grant number (R.G.P2./99/41).

References

K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, Wiley, Hoboken, NJ, USA, 1993.

R. Hilfer, Applications of Fractional Calculus in Physics, World Science Publishing, River Edge, NJ, USA, 2000.

I. Podlubny, “Fractional differential equations,” in An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, Elsevier, Amsterdam, Netherlands, 1998.

L. Debnath, “Recent applications of fractional calculus to science and engineering,” International Journal of Mathematics and Mathematical Sciences, vol. 2003, no. 54, pp. 3413–3442, 2003.

M. Ortigueira, “An introduction to the fractional continuous-time linear systems: the 21st century systems,” IEEE Circuits and Systems Magazine, vol. 8, no. 3, pp. 19–26, 2008.

T. Kaczorek, “Positive linear systems consisting of $n$ subsystems with different fractional orders,” IEEE Transactions on Circuits and Systems I: Regular Papers, vol. 58, no. 6, pp. 1203–1210, 2011.

R. Shah, H. Khan, S. Mustafa, P. Kumam, and M. Arif, “Analytical solutions of fractional-order diffusion equations by natural transform decomposition method,” Entropy, vol. 21, no. 6, p. 557, 2019.

R. Jan and Y. Xiao, “Effect of partial immunity on transmission dynamics of dengue disease with optimal control,” Mathematical Methods in the Applied Sciences, vol. 42, no. 6, pp. 1967–1983, 2019.

R. Jan and Y. Xiao, “Effect of pulse vaccination on dynamics of dengue with periodic transmission functions,” Advances in Difference Equations, vol. 1, p. 368, 2019.

A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and applications of fractional differential equations, vol. 204, Elsevier, Amsterdam, Netherlands, 2006.

V. Lakshmikantham and A. S. Vatsala, “Basic theory of fractional differential equations,” Nonlinear Analysis: Theory, Methods & Applications, vol. 69, no. 8, pp. 2677–2682, 2008.

K. B. Oldham, “Fractional differential equations in electrochemistry,” Advances in Engineering Software, vol. 41, no. 1, pp. 9–12, 2010.

B. Lu, “Bäcklund transformation of fractional Riccati equation and its applications to nonlinear fractional partial differential equations,” Physics Letters A, vol. 376, no. 28-29, pp. 2045–2048, 2012.

Z. Odibat and S. Momani, “The variational iteration method: an efficient scheme for handling fractional partial differential equations in fluid mechanics,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2199–2208, 2009.

H. Khan, R. Shah, P. Kumam, and M. Arif, “Analytical solutions of fractional-order heat and wave equations by the natural transform decomposition method,” Entropy, vol. 21, no. 6, p. 597, 2019.

S. Mahmood, R. Shah, H. khan, and M. Arif, “Laplace adomian decomposition method for multi dimensional time fractional model of Navier-Stokes equation,” Symmetry, vol. 11, no. 2, p. 149, 2019.

R. Shah, H. Khan, D. Baleanu, P. Kumam, and M. Arif, “A novel method for the analytical solution of fractional ZakharovKuznetsov equations,” Advances in Difference Equations, vol. 2019, no. 1, pp. 1–14, 2019.

H. M. Srivastava, R. Shah, H. Khan, and M. Arif, “Some analytical and numerical investigation of a family of fractional‐order Helmholtz equations in two space dimensions,” Mathematical Methods in the Applied Sciences, vol. 43, no. 1, pp. 199–212, 2020.

H. Khan, R. Shah, M. Arif, and S. Bushnaq, “The Chebyshev wavelet method (CWM) for the numerical solution of fractional HIV infection of CD4 + T cells model,” International Journal of Applied and Computational Mathematics, vol. 6, no. 2, pp. 1–17, 2020.

Q. M. Al-Mdallal, M. I. Syam, and M. N. Anwar, “A collocation-shooting method for solving fractional boundary value problems,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 12, pp. 3814–3822, 2010.

M. A. Hajji, Q. M. Al-Mdallal, and F. M. Allan, “An efficient algorithm for solving higher-order fractional Sturm-Liouville eigenvalue problems,” Journal of Computational Physics, vol. 272, pp. 550–558, 2014.

Q. M. Al-Mdallal and M. A. Hajji, “A convergent algorithm for solving higher-order nonlinear fractional boundary value problems,” Fractional Calculus and Applied Analysis, vol. 18, no. 6, p. 1423, 2015.

Q. M. Al-Mdallal, “On fractional-Legendre spectral Galerkin method for fractional Sturm-Liouville problems,” Chaos, Solitons & Fractals, vol. 116, pp. 261–267, 2018.

R. Shah, H. Khan, D. Baleanu, P. Kumam, and M. Arif, “A semi-analytical method to solve family of Kuramoto-Sivashinsky equations,” Journal of Taibah University for Science, vol. 14, no. 1, pp. 402–411, 2020.

M. Madani, M. Fathizadeh, Y. Khan, and A. Yildirim, “On the coupling of the homotopy perturbation method and Laplace transformation,” Mathematical and Computer Modelling, vol. 53, no. 9-10, pp. 1937–1945, 2011.

Y. Khan and Q. Wu, “Homotopy perturbation transform method for nonlinear equations using He’s polynomials,” Computers & Mathematics with Applications, vol. 61, no. 8, pp. 1963–1967, 2011.

S. J. Johnston, H. Jafari, S. P. Moshokoa, V. M. Ariyan, and D. Baleanu, “Laplace homotopy perturbation method for Burgers equation with space- and time-fractional order,” Open Physics, vol. 14, no. 1, pp. 247–252, 2016.

V. F. Morales-Delgado, J. F. Gómez-Aguilar, H. Yépez-Martínez, D. Baleanu, R. F. Escobar-Jimenez, and V. H. Olivares-Peregrino, “Laplace homotopy analysis method for solving linear partial differential equations using a fractional derivative with and without kernel singular,” Advances in Difference Equations, vol. 2016, no. 1, pp. 1–17, 2016.

Y. Li, B. T. Nohara, and S. Liao, “Series solutions of coupled Van der Pol equation by means of homotopy analysis method,” Journal of Mathematical Physics, vol. 51, no. 6, Article ID 063517, 2010.

Y. Keskin and G. Oturanc, “Reduced differential transform method for partial differential equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, no. 6, pp. 741–750, 2009.

P. K. Gupta, “Approximate analytical solutions of fractional Benney-Lin equation by reduced differential transform method and the homotopy perturbation method,” Computers & Mathematics with Applications, vol. 61, no. 9, pp. 2829–2842, 2011.

K. H. Huebner, D. L. Dewhirst, D. E. Smith, and T. G. Byrom, The Finite Element Method for Engineers, John Wiley & Sons, Hoboken, NJ, USA, 2001.

G. D. Smith, G. D. Smith, and G. D. S. Smith, Numerical Solution of Partial Differential Equations: Finite Difference Methods, Oxford University Press, Oxford, UK, 1985.

R. Shah, H. Khan, and D. Baleanu, “Fractional whitham-broer-kaup equations within modified analytical approaches,” Axioms, vol. 8, no. 4, p. 125, 2019.

E. F. Keller and L. A. Segel, “Initiation of slime mold aggregation viewed as an instability,” Journal of Theoretical Biology, vol. 26, no. 3, pp. 399–415, 1970.

J. Lankeit and M. Winkler, A Generalized Solution Concept for the Keller-Segel System with Logarithmic Sensitivity: Global Solvability for Large Nonradial Data, 2017, https://arxiv.org/abs/1701.07391.

F. Haq, K. Shah, Q. M. Al-Mdallal, and F. Jarad, “Application of a hybrid method for systems of fractional order partial differential equations arising in the model of the one-dimensional Keller-Segel equation,” The European Physical Journal Plus, vol. 134, no. 9, p. 461, 2019.

N. Bournaveas and V. Calvez, “The one-dimensional Keller-Segel model with fractional diffusion of cells,” Nonlinearity, vol. 23, no. 4, pp. 923–935, 2010.

M. Zayernouri and A. Matzavinos, “Fractional Adams-Bashforth/Moulton methods: an application to the fractional Keller-Segel chemotaxis system,” Journal of Computational Physics, vol. 317, pp. 1–14, 2016.

V. F. Morales-Delgado, J. F. Gómez-Aguilar, S. Kumar, and M. A. Taneco-Hernández, “Analytical solutions of the Keller-Segel chemotaxis model involving fractional operators without singular kernel,” The European Physical Journal Plus, vol. 133, no. 5, p. 200, 2018.

M. Winkler, “The two-dimensional Keller-Segel system with singular sensitivity and signal absorption: global large-data solutions and their relaxation properties,” Mathematical Models and Methods in Applied Sciences, vol. 26, no. 5, pp. 987–1024, 2016.

Y. Dolak and C. Schmeiser, “The keller--segel model with logistic sensitivity function and small diffusivity,” SIAM Journal on Applied Mathematics, vol. 66, no. 1, pp. 286–308, 2005.

S. Kumar, A. Kumar, and I. K. Argyros, “A new analysis for the Keller-Segel model of fractional order,” Numerical Algorithms, vol. 75, no. 1, pp. 213–228, 2017.

A. Ghorbani and J. Saberi-Nadjafi, “He’s homotopy perturbation method for calculating adomian polynomials,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 229–232, 2007.

A. Ghorbani, “Beyond Adomian polynomials: he polynomials,” Chaos, Solitons & Fractals, vol. 39, no. 3, pp. 1486–1492, 2009.

H. Khan, R. Shah, D. Baleanu, P. Kumam, and M. Arif, “Analytical solution of fractional-order hyperbolic telegraph equation, using natural transform decomposition method,” Electronics, vol. 8, no. 9, p. 1015, 2019.

S. Maitama, “A hybrid natural transform homotopy perturbation method for solving fractional partial differential equations,” International Journal of Differential Equations, vol. 2016, Article ID 9207869, 7 pages, 2016.

S. Maitama, “Local fractional natural homotopy perturbation method for solving partial differential equations with local fractional derivative,” Progress in Fractional Differentiation and Applications, vol. 4, no. 3, pp. 219–228, 2018.

A. Khalouta and A. Kadem, “Comparison of new iterative method and natural homotopy perturbation method for solving nonlinear time-fractional wave-like equations with variable coefficients,” Nonlinear Dynamics and Systems Theory, vol. 19, pp. 160–169, 2019.

Z. H. Khan and W. A. Khan, “N-transform-properties and applications,” NUST Journal of Engineering Sciences, vol. 1, no. 1, pp. 127–133, 2008.

K. Shah, H. Khalil, and R. A. Khan, “Analytical solutions of fractional order diffusion equations by natural transform method,” Iranian Journal of Science and Technology, Transactions A: Science, vol. 42, no. 3, pp. 1479–1490, 2018.