Journal of Applied Mathematics
Volume 2015 (2015), Article ID 657307, 12 pages
http://dx.doi.org/10.1155/2015/657307

Research Article

Normal Forms of Hopf Bifurcation for a Reaction-Diffusion System Subject to Neumann Boundary Condition

Cun-Hua Zhang and Xiang-Ping Yan

Department of Mathematics, Lanzhou Jiaotong University, Lanzhou 730070, China

Received 14 May 2015; Accepted 12 July 2015

Academic Editor: Alain Miranville

Copyright © 2015 Cun-Hua Zhang and Xiang-Ping Yan. 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.

A reaction-diffusion system coupled by two equations subject to homogeneous Neumann boundary condition on one-dimensional spatial domain with is considered. According to the normal form method and the center manifold theorem for reaction-diffusion equations, the explicit formulas determining the properties of Hopf bifurcation of spatially homogeneous and nonhomogeneous periodic solutions of system near the constant steady state are obtained.

1. Introduction

As an important dynamic bifurcation phenomenon in dynamical systems, Hopf bifurcation of periodic solutions has attracted great interest of many authors in the last several decades [1–8]. In general, the study of Hopf bifurcation includes the existence and the properties such as the direction of bifurcation and the stability of bifurcating periodic solutions. In application, however, it is more difficult to determine the properties of Hopf bifurcation than to find the existence of a Hopf bifurcation. An approach applied to determine the properties of Hopf bifurcation is to derive the projected equation of original equations on the associated center manifold, that is, the so-called normal form. Then one may explore the local dynamical behaviors of a higher dimensional or even infinitely dimensional dynamical system near a certain nonhyperbolic steady state according to the normal form obtained. The normal form of Hopf bifurcation in ordinary differential equations (ODEs) with or without delays has been established well [1, 3, 5] since in this case the equilibrium is always constant and there are also no effects of spatial diffusion.

Under some certain conditions, the reaction-diffusion equations under the homogeneous Neumann boundary condition may have the constant steady state and thus one can study the Hopf bifurcation of system at this constant steady state. Compared with the ODEs, it is more difficult to derive the normal form of Hopf bifurcation for reaction-diffusion equations at the constant steady state. Although Hassard et al. [3] established the method computing the normal form of Hopf bifurcation in reaction-diffusion equations with the homogeneous Neumann boundary condition and also considered the Hopf bifurcation of spatially homogeneous periodic solutions in Brusselator system, using the same method, Jin et al. [9] and Ruan [10] as well as Yi et al. [11, 12] considered the Hopf bifurcation of spatially homogeneous periodic solutions for Gierer-Meinhardt system and CIMA reaction, respectively. There are few results regarding Hopf bifurcation of spatially nonhomogeneous periodic solutions for spatially homogeneous reaction-diffusion equations [7].

Based on the reason mentioned above, in this paper we consider the normal form of Hopf bifurcation of reaction-diffusion equations at the constant steady state following the idea in [3]. In order to have a clearer structure, we are concerned with the following general reaction-diffusion system coupled by two equations defined on one-dimensional spatial domain with and subject to Neumann boundary conditions; that is,

in which are the diffusion coefficients, is the parameter, and are    functions with    for any . Although Yi et al. [7] described the algorithm determining the properties of Hopf bifurcation of spatially homogeneous and nonhomogeneous periodic solutions for (1) at and also considered the Hopf bifurcation of a diffusive predator-prey system with Holling type-II functional response and subject to the homogeneous Neumann boundary condition, they did not give the normal form of Hopf bifurcation of spatially homogeneous and nonhomogeneous periodic solutions of the general reaction-diffusion system (1) at .

This paper is organized as follows. In the next section, following the abstract method according to [3], we describe the algorithm determining the properties of Hopf bifurcation of spatially homogeneous and nonhomogeneous periodic solutions for system (1) at the constant steady state . In Section 3, the explicit formulas determining the properties of Hopf bifurcation of spatially homogeneous periodic solutions for system (1) at are obtained. The explicit formulas determining the properties of Hopf bifurcation of spatially nonhomogeneous periodic solutions for (1) at are also derived in Section 4.

2. Algorithm Determining the Properties of Hopf Bifurcation

In this section, we will describe the explicit algorithm determining the direction of Hopf bifurcation and the stability of the bifurcating periodic solutions of system (1) at .

Define the real-valued Sobolev space by

In terms of , the complex-valued Sobolev space is given by

and the inner product on is defined by

Let , , , and and define the linear operator with the domain by

Assume that, for some , the following condition holds:

(H)There exists a neighborhood of such that, for , has a pair of simple and continuously differentiable eigenvalues with , , and . In addition, all other eigenvalues of have nonzero real parts for .

Then from [3, 7] we know that system (1) undergoes a Hopf bifurcation at when crosses through .

Define the second-order matrix sequence by

Then the characteristic equation of is

where

The eigenvalues of can be determined by the eigenvalues of    and we have the following conclusion.

Lemma 1. If is an eigenvalue of the operator , then there exists some such that is the eigenvalue of and vice versa.

Proof. It is well known that the eigenvalue problem

has eigenvalues    with eigenfunctions . Assume that is an eigenvalue of the operator and the corresponding eigenfunction is ; that is,

Notice that can be represented as

where   . Then (10) can be written into

From the orthogonality of the function sequence , one can get from (12) that, for each ,

Since is the eigenfunction of corresponding to the eigenvalue , and so there must be some such that . Therefore, is the eigenvalue of the matrix .
If is the eigenvalue of some matrix , then there exists a nonzero vector such that (13) holds. Let

Then and

This demonstrates that is an eigenvalue of and thus the proof is complete.

Lemma 1 shows that, under assumption (H), there is a unique such that are purely imaginary eigenvalues of ; that is, and . Furthermore, it is easy to see that for any . Therefore,    has eigenvalues with zero real parts if and only if . Assume that is the eigenvalue of for sufficiently approaching . Then by the smoothness of    we know that is also the eigenvalue of ; namely, satisfies the following equation:

Under the assumption (H), differentiating the above equation with respect to at yields

Based on the above discussion, condition (H) has the following equivalent form:

Then we know that and cannot be equal to zero simultaneously when the hypothesis (H) is satisfied. Therefore, the eigenvector of corresponding to the eigenvalue can be chosen as

and thus the eigenfunction of corresponding to the eigenvalue has the form

Let the linear operator with the domain be defined by

Then is the adjoint operator of the operator such that with . Similar to the choice of the eigenfunction of the operator corresponding to the eigenvalue , we can choose

such that

Define and by and , respectively. Then can be decomposed as the direct sum of and ; that is, . Thus, for any , there exists and such that

Define by

Then system (1) can be rewritten into the following abstract form:

When , system (26) is reduced to

where . In terms of (23) and decomposition (24), system (27) can be transformed into the following system in coordinates:

where

For , , and , define the symmetric multilinear forms and , respectively, by

Then, for , we have

For the simplicity of notations, we will use and to denote and , respectively.

Let

Then from (29) and (32), one can get

From the center manifold theorem in [3], we can rewrite in the form

The second equation of (28), (33), and (35) yields

Substituting (35) into the first equation of (28) gives the equation of reaction-diffusion system (1) restricted on the center manifold at as

where , , , and

The dynamics of (28) can be determined by the dynamics of (37).

In addition, it can be observed from [3] that when approaches sufficiently , the Poincaré normal form of (26) has the form

where is a complex variable, , and are complex-valued coefficients with

The direction of Hopf bifurcation and the stability of the bifurcating periodic solutions of (1) at can be determined by the sign of and we have the following conclusion.

Theorem 2. Assume that condition (H) (or equivalently (18)) holds. Then system (1) undergoes a supercritical (or subcritical) Hopf bifurcation at when if

In addition, if all other eigenvalues of have negative real parts, then the bifurcating periodic solutions are stable (resp., unstable) when ().

3. Spatially Homogeneous Hopf Bifurcation

From the description in the previous section we know that Hopf bifurcation of (1) at is spatially homogeneous if condition (18) holds when . In the present section, we compute in (40) in order to determine the direction of spatially homogeneous Hopf bifurcation and the stability of bifurcating periodic solutions of (1) at following the algorithm described in this pervious section.

Lemma 3. If condition (18) is satisfied when , then .

Proof. From (20) and (22) one can see

where

Let all the partial derivatives of    be evaluated at , and let , , and    be defined, respectively, by

Then from (30) and (31) we can get

Therefore,

From (34) and (46), one can obtain

Notice from (42) that

The conclusion follows by substituting (48) into (47).

Lemma 3 and (36) imply when in (18) and thus we have

It follows from (40) that

We represent , , and by , and , respectively, for the simplicity of notations and under assumption (18) with , substituting in (42) into (44) yields that, for ,

From (42), (50), and (51) one can derive

where

Thus we have the following result.

Theorem 4. Assume that condition (18) is satisfied when and is defined by (53). Then the spatially homogeneous Hopf bifurcation of system (1) at is supercritical (resp., subcritical) if

Moreover, if each eigenvalue of has negative real parts for all , then the above spatially homogeneous bifurcating periodic solutions are stable (resp., unstable) when

4. Spatially Nonhomogeneous Hopf Bifurcation

Notice that the spatially nonhomogeneous periodic solutions of (1) at from Hopf bifurcation are unstable. Accordingly, in this section we will calculate in (40) in order to determine the direction of Hopf bifurcation of spatially nonhomogeneous periodic solutions of system (1) at . To this end, we always assume that in (18) throughout this section and still represent , , and by , and , respectively. Thus defined in (22) has the form

where

Since when ,

one can obtain

Thus, in order to calculate , it remains to compute

Let all the second- and third-order partial derivatives of    with respect to and be evaluated at and let

Then from (30) and (31), one can observe

In view of (34), (59), and (62), we have

Equalities (63) show that the calculation of and will be restricted on the subspaces spanned by eigenmodes and .

Let

Then, under condition (18), one can derive

From (36) and (61), we have

Similarly, we can get

where

From (30) we have

with

Notice that, for ,

It follows from (69) that

Substituting into (61) gives

Then from (70) and (73), we have

Since and , one can get

In addition, it follows from and (62) that

Therefore,

Now, by (40), we have

Thus we have the following result.

Theorem 5. Assume that condition (18) holds for . Then the spatially nonhomogeneous Hopf bifurcation of system (1) at is supercritical (resp., subcritical) if

here is given by (78).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Cun-Hua Zhang was supported by the Natural Science Foundation of Gansu Province (1212RJZA065). Xiang-Ping Yan was supported by the National Natural Science Foundation of China (11261028), Gansu Province Natural Science Foundation (145RJZA216), and China Scholarship Council.

References

A. Andronov, E. A. Leontovich, I. I. Gordon, and A. G. Maier, Theory of Bifurcations of Dynamical Systems on a Plane, Israel Program for Scientific Translations, Jerusalem, Israel, 1971.

J. Hale, Theory of Functional Differential Equations, Springer, New York, NY, USA, 1977.

B. D. Hassard, N. D. Kazarinoff, and Y. H. Wan, Theory and Applications of Hopf Bifurcation, Cambridge University Press, Cambridge, UK, 1981.

H. Kielhófer, Bifurcation Theory: An Introduction with Applications to PDEs, vol. 156 of Applied Mathematical Sciences, Springer, New York, NY, USA, 2004.

Y. A. Kuznetsov, Elements of Applied Bifurcation Theory, Applied Mathematical Sciences, Springer, New York, NY, USA, 1998.

J.-H. Wu, Theory and Applications of Partial Functional Differential Equations, Springer, New York, NY, USA, 1996.

F. Yi, J. Wei, and J. Shi, “Bifurcation and spatiotemporal patterns in a homogeneous diffusive predator-prey system,” Journal of Differential Equations, vol. 246, no. 5, pp. 1944–1977, 2009.

L.-J. Yang and X.-W. Zeng, “Stability of singular Hopf bifurcations,” Journal of Differential Equations, vol. 206, no. 1, pp. 30–54, 2004.

J. Y. Jin, J. P. Shi, J. J. Wei, and F. Q. Yi, “Bifurcations of patterned solutions in the diffusive Lengyel-Epstein system of CIMA chemical reactions,” Rocky Mountain Journal of Mathematics, vol. 43, no. 5, pp. 1637–1674, 2013.

S. G. Ruan, “Diffusion-driven instability in the Gierer-Meinhardt model of morphogenesis,” Natural Resource Modeling, vol. 11, no. 2, pp. 131–142, 1998.

F. Yi, J. Wei, and J. Shi, “Diffusion-driven instability and bifurcation in the Lengyel-Epstein system,” Nonlinear Analysis: Real World Applications, vol. 9, no. 3, pp. 1038–1051, 2008.

F. Q. Yi, J. J. Wei, and J. P. Shi, “Global asymptotical behavior of the Lengyel-Epstein reaction-diffusion system,” Applied Mathematics Letters, vol. 22, no. 1, pp. 52–55, 2009.