#### Abstract

We consider the initial stage of space-time fractional generalized biological equation in radial symmetry. Dimensionless multiorder fractional nonlinear equation was first given, and approximate solutions were derived in the form of series using the homotopy perturbation method with a new modification. And the influence of fractional derivative is also discussed.

#### 1. Introduction

The problem of biological diffusion is an issue of increasing significance in contemporary ecology [1, 2]. In case of favourable environmental conditions, the alien population may begin to grow and spread over the area and thus the local initial structural perturbation of the native biological community may lead to large-scale dramatic changes in the community structure. Recently, it has turned out that many phenomena in engineering, physics, chemistry, and other sciences [3–5] can be described very successfully by models using mathematical tools fractional calculus [6, 7], such as anomalous transport in disordered systems [8–10], some percolations in porous media, and the diffusion biological population. Mathematical aspects of the biological problem have been considered in many papers [11–13]. El-Sayed et al. [14] studied the fractional-order biological population model in the form using the Adomian decomposition method. Wazwaz and Gorguis [15] gave a detailed study of integer Fisher’s diffusion equation by using Adomian decomposition method. Najeeb et al. [16] studied the time fractional Fisher’s equation and approximate analytical solutions were obtained by using homotopy analysis method. Petrovskii et al. [17] obtained an exact solution of the spatiotemporal dynamics of a predator-prey community by using an approximate change of variables, and the properties of the solution exhibit biologically reasonable dependence on the parameter values. Liu and Xin [18] studied the fractional Lotka-Volterra equations using the homotopy perturbation method.

This paper is devoted to investigating approximate solutions of a generalized fractional nonlinear population diffusion equation in radical symmetry. The structure of the paper is as follows. In Section 2, a brief review of the theory of fractional calculus will be given to fix notation and provide a convenient reference. In Section 3, a mathematical formulation of the generalized multifractional population diffusion model in radical symmetry is given. In Section 4, we extend the homotopy perturbation method and a new reliable modification to the fractional nonlinear population diffusion system and give some properties of this model. Conclusions and prospects will be presented in Section 5.

#### 2. Fractional Calculus

There are several approaches to define the fractional calculus; the Riemann-Liouville and Caputo fractional operators are defined as follows.

*Definition 1. *The Riemann-Liouville fractional integral operator of a function is defined as
where is the well-known gamma function, and some properties of the operator are as follows:

*Definition 2. *The Caputo fractional derivative of a function is defined as
The following are two basic properties of the Caputo fractional derivative:
We have chosen the Caputo fractional derivative because it allows traditional initial and boundary conditions to be included in the formulation of the problem. And some other properties of fractional derivative can be found in [3, 5].

#### 3. Main Equations

It is widely accepted that the spatiotemporal dynamics of a biological community can be qualitatively described by diffusion-reaction equations [2, 19]. A remarkable point is that, in some cases, relatively simple single-species models provide not only qualitative but also quantitative descriptions of the dynamics of a population. In this paper, we first consider a single-species parabolic nonlinear equation arising in the spatial diffusion of biological populations [11, 14] where is the population density at position , and time , coefficient describes the intensity of mixing due to animals self-motion, the term describes multiplication and mortality of a given population, , and parameter is being treated as the carrying capacity for the given population. Let us assume that the initial distribution of the species is symmetrical, with the density of the species depending only on the distance from the origin. Assuming also that the environment is homogeneous, that is, both and not depends explicitly on the position in space, we arrive at the following problem: where and is the typical size of the domain and initial condition promptly approaches zero when . It has turned out that the diffusion of biological population can be described very successfully by fractional calculus. In this paper, we discuss the corresponding fractional equation and the main aim is to solve the nonlinear fractional biological population model in the following form: where is the Caputo derivative. A proper choice of the dimensionless variables, that is, in our case, the choice of scales for the variables , and , is an important point. Coming with the property of the Caputo derivative and using reduced dimensionless variables defined as (8) can be reduced to the respective dimensionless forms (tildes will be omitted hereafter): where .

#### 4. Approximate Solution to the Equation

We now proceed to derive approximate solution to fractional nonlinear population diffusion equation (9).

*Case 1. *In this case, we will examine the following time-space fractional nonlinear population model:
where = constant, corresponding to the *Malthusian* law. According to the homotopy perturbation method, we construct the following simple homotopy:
where is an embedding parameter. In case , (12) is a fractional differential equation, , which is easy to solve, and when , (12) turns out to be the original one (10). The basic assumption is that the solutions can be written as a power series in :
is an initial approximation of (10). The approximate solutions of the original equations can be obtained by setting , that is,
Instituting (13) into (12) and comparing coefficients of terms with identical powers of then applying on both sides of equations yield
and so on, in this manner the rest of components of the solution can be obtained. The solution of (10) in series form is given by

Figure 1 shows the approximate solution for (10) and (11) by using the homotopy perturbation method when choosing . From the figure, it is clear to see the time evolution of nonlinear population diffusion density and we also know that the approximate solution of fractional population model is continuous with the fractional parameter . Figure 2 shows the approximate solution for (10) and (11) when , and the approximate solution of fractional population model is continuous with the fractional parameter . Figures 3 and 4 show the approximate solution for (10) and (11) when the time , from the figures, we also know that the population density changes with the parameters .

*Case 2. *In this case, we will examine the following time-space fractional nonlinear population model:
are constant, corresponding to the *Verhulst* law. According to the homotopy perturbation method, we construct the following simple homotopy:
For this case, it is difficult to solve the multifractional equation. We give a new modification of the homotopy perturbation method; the modified form of the homotopy perturbation method can be established based on the initial condition expressed in the Taylor series. We suggest that be expressed in the Taylor series
Instituting (13) into (19) and comparing coefficients of terms with identical powers of then applying on both sides of equations yield
and so on, in this manner the rest of components of the solution can be obtained. The solution of (17) in series form is given by

*Case 3. *We will consider the following initial value problem of time-space fractional nonlinear diffusion equation:
According to the homotopy perturbation method, we construct the following simple homotopy:
The modified form of the homotopy perturbation method can be established based on the initial condition expressed in the Taylor series; the initial condition is expressed in the Taylor series
Instituting (12) into (23) and equating coefficients of terms with identical powers of
where is the coefficient of in and is the coefficient of in , then applying , the inverse operator of , on both sides of equations, it is obvious that (27) are easy to solve, the components of the homotopy perturbation method can be completely determined, and series solutions are thus entirely determined.

#### 5. Conclusion

Approximate solutions of the multifractional nonlinear diffusion population equations in radial symmetry were derived using the homotopy perturbation method and the new modification of homotopy perturbation method. The solutions are given in the form of series with easily computable terms. The results reveal that the new modified method is very effective for solving nonlinear diffusion equation of multifractional order. This is the first step to study the multifractional nonlinear population diffusion in radical symmetry, and we will make subsequent research, for example, exact solution and self-similar exact solution of these fractional nonlinear system. And we hope that this work is a step in this direction.

#### Acknowledgments

The author expresses thanks to the referees for their fruitful advices and comments. This work was supported by the National Science Foundation of Shandong Province (Grants no. Y2007A06 and ZR2010Al019) and the China Postdoctoral Science Foundation (Grant no. 20100470783.)