#### 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 getTaking natural transformation of equation (20), we getSimplifying the above equation, we getUsing inverse natural transform, we haveThen, the chemotactic term .

Now, implementing HPM, we getBy comparing the coefficients of powers *p*, we getThe series form of equation (20) is given byThe HPTM solution is given asThe 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 getSimplifying the above equation, we getUsing inverse natural transform, we haveNow, implementing HPM, we getwhere and are He’s polynomials that show the nonlinear terms. So, representing He’s polynomials for , we find thatFor , we find thatThe series solutions are obtained by using the numerical values , , , , , and .

By comparing the coefficients of powers *p*, we getThe series form of equation (30) is given byThe 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)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

#### 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).