BioMed Research International

Volume 2014 (2014), Article ID 815690, 9 pages

http://dx.doi.org/10.1155/2014/815690

Research Article

Modelling the Aggregation Process of Cellular Slime Mold by the Chemical Attraction

Abdon Atangana and P. D. Vermeulen

Institute for Groundwater Studies, Faculty of Natural and Agricultural Sciences, University of the Free State, Bloemfontein 9300, South Africa

Received 30 May 2014; Accepted 3 June 2014; Published 8 July 2014

Academic Editor: Suares Clovis Oukouomi Noutchie

Copyright © 2014 Abdon Atangana and P. D. Vermeulen. 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.

We put into exercise a comparatively innovative analytical modus operandi, the homotopy decomposition method (HDM), for solving a system of nonlinear partial differential equations arising in an attractor one-dimensional Keller-Segel dynamics system. Numerical solutions are given and some properties show evidence of biologically practical reliance on the parameter values. The reliability of HDM and the reduction in computations give HDM a wider applicability.

1. Introduction

In 1970, Keller and Segel have offered parabolic systems to illustrate the aggregation process of cellular slime mold by the chemical attraction [1]. The system of a simplified form in the one-dimensional case is written as

subject to the boundaries conditions

and initial conditions

where   is a bounded open interval and , , , and are positive constants. The unknown functions and denote the concentration of amoebae and the concentration of chemical substance, respectively, in . The chemotactic term indicates that the cells are sensitive to the chemicals and are attracted by them. called the sensitivity function is a smooth function of which describes cell’s perception and response to the chemical stimulus . Several normalized forms have been suggested , , , and , and so forth (see [2, 3]). Recently, the Keller-Segel (KS) equations attracted interests of many mathematicians. Since the modelling of chemotaxis has developed into a large and diverse discipline, one model which is widely used is the Keller-Segel model of chemotaxis; it is important to recall that chemotaxis describes the movement of single or multicellular organisms when they move up or down a chemical gradient [4]. This movement allows the organism to explore its extracellular environment. Organisms move randomly, away from repellents and towards attractants. Questions have arisen on how organisms can detect small changes in their extracellular environment [5]. Usually the organism will undergo a random walk, consisting of smooth swimming and brief direction changes (tumbles). By increasing the attractant, the tumbling is suppressed, which leads to a biased random walk [4]. The organism will then accumulate in areas of high attractant concentration. This type of movement is referred to as runs [5]. A combination of tumbles and runs allows the organism to explore and respond to changes in its extracellular environment [4].

The local solutions were studied by the second author [6]. It was also suggested in [6] that, in the one-dimensional case, (KS) possesses a global solution and that, in the two-dimensional case, when ( being a positive constant) is a linear function, (KS) possesses a global solution for any sufficiently small initial function . Horstmann and Wang [7] showed more strongly that the global solution exists if the norm is smaller than a specific number, which is given from the coefficients of the equations. Recently, in the same case, the asymptotic behaviour of the global solutions was studied in [8]. On the other hand, Herrero and Velázquez [9] showed that when is linear and the domain is a circular disc, there exist radial local solutions which blow up in a finite time. The blowup of nonradial local solutions was shown recently by [10, 11]. For the study of stationary solutions, we refer to [12–14]. In the field of dynamical systems theory some work has been carried out although the suggested models and algorithms are still in an introductory platform of establishing. It is perhaps important to notice that several analytical methods have been proposed to deal with nonlinear equations, but there exist a lot of nonlinear ordinary differential equations and nonlinear partial differential equations for which exact analytical solution cannot be found. There is no exact solution of (4) in the literature. To solve these problems, some eminent scholars have proposed some powerful iteration methods to deal with this class of nonlinear equation.

As V.M. Alexandrov wrote in the introduction of a well-liked science book Asymtotology: Ideas, Methods, and Applications [15, 16], asymptotic methods belong to the, perhaps, most romantic area of modern mathematics [15–19]. Though computer science is growing very fast and numerical simulation is applied everywhere, nonnumerical issues will still play a large role [16, 20–22]. There exist some alternative analytical asymptotic approaches such as the nonperturbative method, modified Lindstedt-Poincare method [21], variational iteration method [22], Adomian decomposition method [23], homotopy perturbation method [17, 24], and bookkeeping artificial parameter perturbation method [18].

The purpose of this paper is to derive analytical solutions of attractor one-dimensional Keller-Segel equations (1) via the relatively new analytical method the modified homotopy perturbation method. The HDM was recently used in [19, 25–27]. This method displays some advantages over existing methods.

The paper is prearranged as follows: in Section 2, we present the basic idea of the HDM for solving high orders differential equations. We present the application of the HDM for attractor one-dimensional Keller-Segel equations and numerical results in Section 3. In Section 4 we present the discussions. The conclusions are then given in Section 5.

2. Basic Properties of Homotopy Decomposition Method

With the purpose of making the fundamental possessions of the homotopy decomposition method [28] clear, we think about a universal nonlinear nonhomogeneous partial differential equation with initial conditions of the following form:

Subject to the initial condition

is the order of the derivative.

Where is a known function, is the general nonlinear differential operator and represents a linear differential operator, and is the order of the derivative. The method’s first step is to apply the inverse operator of on both sides of (4) to obtain

The multi-integral in (4) can be transformed to

so that (4) can be reformulated as

Using the homotopy scheme the solution of the above integral equation is given in series form as

and the nonlinear term can be decomposed as

where is an embedding parameter. is the He’s polynomials that can be generated by

The modified homotopy perturbation method is obtained by the coupling of decomposition method with Abel integral and is given by

with

Comparing the terms of same powers of produces solutions of various orders. The initial guess of the approximation is [26, 27]. This is actually the Taylor series of the exact solution of order . Note that this initial guess insures the uniqueness of the series decompositions [26, 27].

3. Application

In this section we apply this method for solving coupled attractor one-dimensional Keller-Segel equations.

Example 1. Consider the following Keller Segel equation with the sensitivity function .

Then the chemotactic term

Subject to the initial conditions

In the view of the HDM, we obtain the following equation:

Now comparing the terms of same power of , we obtained the following integral equations:

The following solutions are obtained straightforwardly:

Using the iterative formula, the remaining terms can be obtained. But here, only few terms of the series solutions are considered and the asymptotic solution is given as

The following figures show the biological behaviour of the coupled solutions for the following set of theoretical parameters: , , , , , and , first for a fixed time and secondly for a fixed distance .

Figures 1, 2, and 3 show the behaviour of the solution of the system of (15) describing the concentrations of the chemical substance and the amoebae in the human body. While on one hand, Figure 1 shows the behaviour as function of space, Figures 2 and 3 show the behaviour of these solutions as function of time. From the above figures, one can see that the concentration of amoebae reduces in space as the concentration of the chemical substance reduces. This simply implies that if the concentration of the chemical substance introduced in the human system to combat the spread of the disease is not sufficient enough, the amoebae will spread all over and the patient will certainly die. However, if this concentration is sufficient enough, the amoebae will decrease in space. It is observed from the graphical representation that the approximate solutions obtained here display the behaviour of the real world problem.

(a)

(b)

Figure 1: Biological behaviour of concentrations of the chemical substance and amoebae as function of space.

(a)

(b)

Figure 2: Biological behaviour of concentrations as function of time.

Figure 3: The behaviour of the coupled solutions.

Example 2. Consider the following Keller-Segel equation with the sensitivity function .

With the chemotactic term = + ,

subject to the initial conditions

In the view of the homotopy decomposition method, we arrive at the following set of integral equations that are very easy to handle:

Here we will consider two cases. Case one: we suppose that and . The following series solutions are obtained:

The remaining terms can be obtained by using the iteration formula. But here, only few terms of the series solutions are considered and the asymptotic solution is given as

Therefore when tends to infinity the concentration of amoebae converges to

The following figures show the biological behaviour of the coupled solution for the following set of theoretical parameters: , , , , , and , first for a fixed time and secondly for a fixed distance .

Second case, we suppose that

Following the homotopy decomposition steps, we arrived at the following series solutions:

Using the iteration formulas, the remaining terms can be obtained. But here, only few terms of the series solutions are considered and the asymptotic solution is given as

The following figures show the biological behaviour of the coupled solutions for the following set of theoretical parameters: , , , , , and , first for a fixed time and secondly for a fixed distance .

The above figures show the behaviour of the solution of the system of (20) with initial conditions in (21) and (26). These solutions are describing the concentrations of the chemical substance and the amoebae in the human body for a given set of theoretical parameters chosen according to the literatures. While on one hand, Figure 4 shows the behaviour as function of space, Figures 5, 6, and 7 show the behaviour of these solutions as function of time. From the above figures, one can deduce that the concentration of amoebae reduces in space as the concentration of the chemical substance reduces. It is observed from the graphical representation that the approximate solutions obtained display the behaviour of the real world problem.

Figure 4: Coupled solutions.

Figure 5: Coupled solutions.

Figure 6: Coupled solutions.

Figure 7: Coupled solutions.

Example 3. Consider the following Keller-Segel equation with the sensitivity function .

With the chemotactic term = + ,

subject to the initial conditions

Following the homotopy decomposition steps, we arrive at the following integral equations that are very easy to solve:

We will suppose that and . The following series solutions are obtained:

Using the iteration formulas, the remaining terms can be obtained. But here, only few terms of the series solutions are considered and the asymptotic solution is given as

4. Discussion

The homotopy decomposition method is chosen to solve this kind of nonlinear problem. Because of the following advantages that, the HDM has over the exiting methods. The method does not require the linearization or assumptions of weak nonlinearity [29, 30]. The solutions are not generated in the form of general solution as in Adomian decomposition method [29, 31, 32]. No Lagrange multiplier and correction functional are required as in the case of the variational iteration method [22, 30, 31, 34]. It is more realistic compared to the method of simplifying the physical problems. If the exact solution of the partial differential equation exists, the approximated solution via the method converges to the exact solution [26]. A construction of a homotopy is not needed as in the case of the homotopy perturbation method, because in this case one needs first to continuously deform a difficult problem into another one, which is easy to solve [17, 18, 24, 33]. HDM provides us with a convenient way to control the convergence of approximation series without adapting , as in the case of [24] which is a fundamental qualitative difference in analysis between HDM and other methods [29–32, 34–38].

5. Conclusion

An interesting biological problem describing theaggregation process of cellular slime mold by the chemical attraction was investigated in this paper. We made use of the efficient method called homotopy decomposition method to derive the solution of the mathematical equation underpinning this problem. Analysis and results of nonlinear system of attractor one-dimensional Keller-Segel equation indicate that the model matches the regular biological diffusion behaviour observed in the field.

Conflict of Interests

The authors declare that they have no conflict of interests.

Authors’ Contribution

Abdon Atangana wrote the first draft, P. D. Vermeulen revised the paper, and all the authors corrected the last version.

References

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.

R. Lapidus and M. Levandowsky, “Modeling chemosensory responses of swimming eukaryotes,” in Biological Growth and Spread (Proc. Conf., Heidelberg, 1979), vol. 38 of Lecture Notes in Biomathematics, pp. 388–396, Springer, Berlin, Germany, 1979.

R. Schaaf, “Stationary solutions of chemotaxis systems,” Transactions of the American Mathematical Society, vol. 292, no. 2, pp. 531–556, 1985.

M. J. Tindall, P. K. Maini, S. L. Porter, and J. P. Armitage, “Overview of mathematical approaches used to model bacterial chemotaxis II: bacterial populations,” Applied Numerical Mathematics, vol. 70, no. 6, pp. 1570–1607, 2009.

E. Hildebrand and U. B. Kaupp, “Sperm chemotaxis: a primer,” Annals of the New York Academy of Sciences, vol. 1061, pp. 221–225, 2005.

A. Yagi, “Norm homotopy of solutions to a parabolic system of chemotaxis,” Mathematica Japonica, vol. 45, no. 2, pp. 241–265, 1997.

D. Horstmann and G. Wang, “Blow-up in a chemotaxis model without symmetry assumptions,” European Journal of Applied Mathematics, vol. 12, no. 2, pp. 159–177, 2001.

H. Gajewski and K. Zacharias, “Global omotopy of a reaction-diffusion system omotopy chemotaxis,” Mathematische Nachrichten, vol. 195, pp. 77–114, 1998.

M. A. Herrero and J. J. L. Velázquez, “A blow-up mechanism for a chemotaxis model,” Annali della Scuola Normale Superiore di Pisa. Classe di Scienze. Serie IV, vol. 24, no. 4, pp. 633–683, 1997.

T. Nagai, T. Senba, and K. Yoshida, “Application of the Trudinger-Moser inequality to a parabolic system of chemotaxis,” Funkcialaj Ekvacioj, vol. 40, no. 3, pp. 411–433, 1997.

T. Nagai, T. Senba, and T. Suzuki, “Chemotactic collapse in a parabolic system of mathematical biology,” Hiroshima Mathematical Journal, vol. 30, no. 3, pp. 463–497, 2000.

W. Ni and I. Takagi, “Locating the peaks of least-energy solutions to a semilinear Neumann problem,” Duke Mathematical Journal, vol. 70, no. 2, pp. 247–281, 1993.

T. Senba and T. Suzuki, “Some structures of the solution set for a stationary system of chemotaxis,” Advances in Mathematical Sciences and Applications, vol. 10, no. 1, pp. 191–224, 2000.

H. Wiebers, “-shaped bifurcation curves of nonlinear elliptic boundary value problems,” Mathematische Annalen, vol. 270, no. 4, pp. 555–570, 1985.

I. V. Andrianov and L. I. Manevitch, Asymptotology : Ideas, Methods, and Applications, Kluwer Academic, 2003.

I. Andrianov and J. Awrejcewicz, “Construction of periodic solutions to partial differential equations with non-linear boundary conditions,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 1, no. 4, pp. 327–332, 2000.

J. He, “A coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000.

J. He, “Bookkeeping parameter in perturbation methods,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 2, no. 3, pp. 257–264, 2001.

A. Atangana and J. F. Botha, “Analytical solution of the groundwater flow equation obtained via homotopy decomposition method,” Journal of Earth Science & Climatic Change, vol. 3, no. 2, p. 115, 2012.

C. M. Bender, K. A. Milton, S. S. Pinsky, and L. M. Simmons Jr., “A new perturbative approach to nonlinear problems,” Journal of Mathematical Physics, vol. 30, no. 7, pp. 1447–1455, 1989.

B. Delamotte, “Nonperturbative method for solving differential equations and finding limit cycles,” Physical Review Letters, vol. 70, no. 22, pp. 3361–3364, 1993.

J. He, “Variational iteration method: a kind of non-linear analytical technique: some examples,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 699–708, 1999.

G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988.

J. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.

Y. Tan and S. Abbasbandy, “Homotopy analysis method for quadratic Riccati differential equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 3, pp. 539–546, 2008.

A. Atangana and A. Secer, “The time-fractional coupled-Korteweg-de-Vries equations,” Abstract and Applied Analysis, vol. 2013, Article ID 947986, 8 pages, 2013.

A. Atangana, O. A. Ahmed, and N. Bıldık, “A generalized version of a low velocity impact between a rigid sphere and a transversely isotropic strain-hardening plate supported by a rigid substrate using the concept of noninteger derivatives,” Abstract and Applied Analysis, vol. 2013, Article ID 671321, 9 pages, 2013.

J. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006.

J. He, “Asymptotic methods for solitary solutions and compactons,” Abstract and Applied Analysis, vol. 2012, Article ID 916793, 130 pages, 2012.

G. Wu, “Challenge in the variational iteration method—a new approach to identification of the Lagrange multipliers,” Journal of King Saud University-Science, vol. 25, no. 2, pp. 175–178, 2013.

G. Wu and D. Baleanu, “Variational iteration method for the Burgers' flow with fractional derivatives—new Lagrange multipliers,” Applied Mathematical Modelling, vol. 37, no. 9, pp. 6183–6190, 2013.

J. S. Duan, R. Rach, D. Buleanu, and A. M. Wazwaz, “A review of the Adomian decomposition method and its applications to fractional differential equations,” Communications in Fractional Calculus, vol. 3, no. 2, pp. 73–99, 2012.

M. Matinfar and M. Ghanbari, “The application of the modified variational iteration method on the generalized Fisher's equation,” Journal of Applied Mathematics and Computing, vol. 31, no. 1-2, pp. 165–175, 2009.

M. Javidi and M. A. Raji, “Combinaison of Laplace transform and homotopy perturbation method to solve the parabolic partial differential equations,” Communications in Fractional Calculus, vol. 3, pp. 10–19, 2012.

L. Xu, “Variational principles for coupled nonlinear Schrödinger equations,” Physics Letters A, vol. 359, no. 6, pp. 627–629, 2006.

L. Xu, “Variational approach to solitons of nonlinear dispersive equations,” Chaos, Solitons & Fractals, vol. 37, no. 1, pp. 137–143, 2008.

L. Xu, “Application of He's parameter-expansion method to an oscillation of a mass attached to a stretched elastic wire,” Physics Letters A, vol. 368, no. 3-4, pp. 259–262, 2007.

L. Xu, “Determination of limit cycle by He's parameter-expanding method for strongly nonlinear oscillators,” Journal of Sound and Vibration, vol. 302, no. 1-2, pp. 178–184, 2007.