Discrete Dynamics in Nature and Society

Volume 2013 (2013), Article ID 412343, 7 pages

http://dx.doi.org/10.1155/2013/412343

Research Article

An Algebraic Method on the Eigenvalues and Stability of Delayed Reaction-Diffusion Systems

Jian Ma1,2 and Baodong Zheng²

¹Department of Mathematics, Northeast Forestry University, Harbin 150040, China

²Department of Mathematics, Harbin Institute of Technology, Harbin 150001, China

Received 2 April 2013; Accepted 8 July 2013

Academic Editor: Xiang Ping Yan

Copyright © 2013 Jian Ma and Baodong Zheng. 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.

The eigenvalues and stability of the delayed reaction-diffusion systems are considered using the algebraic methods. Firstly, new algebraic criteria to determine the pure imaginary eigenvalues are derived by applying the matrix pencil and the linear operator methods. Secondly, a practical checkable criteria for the asymptotic stability are introduced.

1. Introduction

The reaction-diffusion system is a semilinear partial differential equation, which has been used for the study of morphogenesis, population dynamics, and autocatalytic oxidation reactions. The early reaction-diffusion models were constructed to describe the process of chemical reaction. For example, the Brusselator reaction-diffusion system and its improved systems had been researched by many scholars [1, 2]. Besides, during the past years, many results about the stability, steady state bifurcation, and the Hopf bifurcation on the reaction-diffusion systems had been derived [3–8]. From those results, we see that most of them had been mainly focused on two methods: the analytical methods and the numerical methods. The analytical methods, such as the center manifold theorem, the normal form theory, and the Laplace transformation, were mainly used to research the analytic solutions and dynamical property. The numerical methods, such as the Runge-Kutta methods and linear multistep methods, paid close attention to checking the changes and convergence of solutions. For example, Wu determined the direction and stability of periodic solutions occurring through the Hopf bifurcation by the center manifold theory and the normal form theory, which are the classic analytical methods in functional differential systems [1]. Wei et al. in papers [3–5] demonstrated the bifurcation and stability of different reaction-diffusion systems.

In this paper, we will introduce a new algebraic method. In fact, in the past recent years, many algebraic methods were applied to the functional differential systems, such as matrix theory, polynomial theory, the Lie group, and algebraic system. In particular for the study on the more complex systems, such as the high dimensional functional differential systems, the algebraic methods have their advancement. So in the following, we will discuss the delayed reaction-diffusion systems by the algebraic methods. We will study a general delayed reaction-diffusion system with spatial domain in the following form:

where the spatial domain is restricted to one-dimensional interval , that is, . Also    stands for the concentrations of reactants or densities of species. denotes the Laplace operator and .    denotes the diffusion coefficients of , respectively. This paper will consider (1) subject to the origin condition and the homogeneous Neumann boundary condition on spatial domain :

For the system (1) and the determined condition (2), we suppose that , for any , and is , . That is to say, the origin point is an equilibrium of system (1) and (2). So the matrix form of system (1) can be written as

where

The linearization of system (3) at origin point is

Furthermore, it is well known that the eigenvalue problem

has eigenvalues , , with corresponding eigenfunctions . Then eigenvalues of are , . Let

Then is an eigenfunction for (5) with eigenvalue . By straightforward analysis, we obtain

where . So we can find all of the eigenvalues of system (5) from (8). Certainly, the stability, and Hopf bifurcation of system (1) can be described by studying the linear system (5).

The remaining parts of the paper are structured in the following way. In Section 2, we demonstrated the critical condition on the delay of the system (1) and got the algebraic criteria for determining the pure imaginary eigenvalues. In Section 3, we researched the stability and the Hopf bifurcation of the delayed reaction-diffusion equation (39) with Neumann boundary condition and derived the corresponding algebraic criteria. At last, we described a specific reaction-diffusion equation and simulated the results by MATLAB.

2. Algebraic Criteria for Determining the Pure Imaginary Eigenvalues

Firstly, we research an ordinary differential equation

where , . Let denote the vector space , let denote the identity operator on , and let operator on be given by

With , the system (9) can be written as

Suppose that is a matrix solution of the system (9). We have

For any complex , let be the operator . Then

For any complex , let , satisfying

For any complex , let , satisfying

By simple computations, we can get

By the operator language, that is,

In the following, we will convert matrix ordinary differential equation (10) to vector form. Let be the elementary transform, . That is,

Let , and

Using the property of the Kronecker product, we have

where . So (11) can be written as

Similarly, by denoting , , and as follows:

we have

Lemma 1. For all complex , , and so .

Proof. Let , where , , , and . Noting the Kronecker product identities , we have , . By regularity of , we know that is nonsingular for enough complex . So by the property of the polynomial and easy computations, we can get that for any complex . Elsewhere,

so

Thus

Theorem 2. Any pure imaginary eigenvalue of the system (5) or the system (1) is a zero point of and thus also one of the eigenvalues of the matrix pencil .

Proof . For (8), let , . Then we have

Let be a pure imaginary eigenvalue of the system (5), let be associated eigenvector, and let . We have . By conjugating and transforming, we can get

Via the elementary transform , we get

That is , . We know that , and so,

From Theorem 2, we know that all of the pure imaginary eigenvalues of the system (1) are zero points of the algebraic equation

So the pure imaginary eigenvalues of the system (1) or (5) can be computed via the algebraic equation (31). In fact, (31) is usually called a polynomial eigenvalue problem. The classical and most widely used approach to research the polynomial eigenvalue problems is linearization, where the polynomial is converted into a larger matrix pencil with the same eigenvalues. There are many forms for the linearization, but the companion form is most typically commission [9]. So we have the following results.

Theorem 3. The pure imaginary eigenvalues of the system (5) are the general eigenvalues of matrix pencil , where

, , and .

Proof . Equation (31) can be rewritten as

where . For , then

Let , , and ; then we have

Equation (35) is a quadratic eigenvalue problem, and its linearization is

where . So let

Then has nonzero solutions if and only if

So is the general eigenvalue of matrix pencil .

Suppose that the system (1) is stable at the initial time. Then the stability of the system will change as the positive real part root of characteristic equation emerged for the different parameter values of system (1). So the first pure imaginary eigenvalue is the critical condition and plays an important role in researching the stability and the Hopf bifurcation of the delayed reaction-diffusion equation. From the previous results, we can get all of the pure imaginary eigenvalues for a delayed reaction-diffusion equation. Here let the time delay be the parameter. Then we can get the stable or unstable condition of the delayed reaction-diffusion equation for different parameter values.

3. The Stability and the Hopf Bifurcation of a Delayed Reaction-Diffusion Equation

In the recent years, scholars discussed many different reaction-diffusion systems [10–12]. For different systems, they expounded different results to justify the stability. Here we will study a general high-dimensional reaction-diffusion equation with single delay. The considered reaction-diffusion system is

where and is a bounded connected domain in with smooth boundary . The functions are continuously differentiable, and is the out normal derivative of . We assume that such that

That is to say, is the positive equilibrium of system (39). The linearization of system (39) is

where and . For simplicity we suppose that the eigenvalues of the Laplace operator by the Neumann boundary condition are simple. From Section 1, we know that the eigenvalues of in with the Neumann boundary condition are , . Let denote an eigenvalue of system (39). Then there exists a nonzero vector such that

The stability of the system (39) is determined by

We know that the system (39) is asymptotically stable when all of the eigenvalues of (42) are in the left part of the imaginary axis. Apparently, for nondiffusive case, that is, , the system (39) becomes a delay differential equation, and the stability has been studied by many scholars. Here we mainly discuss the diffusive system. Let

Firstly, when , the system (39) becomes a partial differential equation with no delay. So we can get the stability in the following.

Lemma 4. For , let be any element of the set . That is, .

(i)The system (39) is eigenmode -independently stable if , for all ;(ii)else, the system (39) is eigen-mode -dependently stable if , for a fixed .

Next we research the delay-independent and the delay-dependent stabilities of the system (39). Let denote the set which contains all the general eigenvalues of matrix pencil . From (38) in Section 2, we have the following results.

Theorem 5. For any and a fixed eigen-mode , if the coefficient matrices of the system (39) satisfy the following conditions:

(i), (ii),

then the system (39) is asymptotically stable for all with fixed eigen-mode ; for example, the stability is delay-independent and eigenmode -dependent. On the other hand, if, for all eigen-mode , the stability is delay independent, then stability of system (39) is delay independent and eigenmode -independent.

Lemma 6. Given a fixed eigen-mode , let , . For all the pure imaginary roots of (31) (or general pure imaginary eigenvalues of matrix pencil ), one can get the critical value of delay parameter for the fixed eigen-mode .

Conclusion. Suppose that the system (39) is eigen-mode -independently stable at , and for multiple eigen mode , the critical values of delay parameter are , . Let . One can get the minimal critical value of delay parameter . When , the system (39) is asymptotically stable, and when , the stability of system (39) changes. That is to say, the positive equilibrium of system (39) is delay dependently stable, and the system (39) first generates bifurcation at .

From the above theorems, we can discuss the stability and the Hopf bifurcation of system (39). So the dynamical character of a generic class of nonlinear reaction-diffusion system with the Neumann boundary condition can be studied by the algebraic methods. It is well known that reaction and diffusion of physical chemistry and chemical or biochemical species can produce all kinds of spatial patterns. Next we will give a general example.

Example 7. Consider a delayed reaction-diffusion equation with the Neumann boundary condition:

where , , , , and especially denote the decay rate. This system is a biochemically mathematical model, such as the improved Gierer gene model. Apparently, is the unique constant positive equilibrium of system (45). Here we will only consider the numerical case of , , , , , , and . Hence the positive equilibrium . Then (43) has no root in the right of the imaginary axis for , has two roots in the imaginary axis and for , and has two roots in the imaginary axis and for . The parameter has four sequences , , , and , . So the smallest Hopf bifurcation value is . Then the unique positive equilibrium of system (45) is asymptotically stable when and unstable when . The simulation results are shown in Figures 1 and 2.

(a)

(b)

Figure 1: For system (45), when , the equilibrium is locally asymptotically stable.

(a)

(b)

Figure 2: For system (45), when , the equilibrium is unstable.

4. Conclusion

In this paper, we consider a general high-dimensional delayed reaction-diffusion system. By the algebraic methods, such as the matrix pencil and the linear operator, we discussed the eigenvalues and the stability of the delayed reaction-diffusion system (1) and (39). In fact, we only find the pure imaginary eigenvalues, which are a few parts of the infinite eigenvalues. So we can easily compute pure imaginary eigenvalues from the algebraic equation (31) by MATLAB. Certainly, applying the algebraic methods to analyze the dynamical properties of the reaction-diffusion system with delays is still a new and immature field, so we believe that the algebraic methods used to research the stability of the delayed reaction-diffusion systems would be of more interest in the future.

Acknowledgment

The research is funded by the Fundamental Research Funds for the Central Universities (DL12BB24).

References

J. Wu, Theory and Applications of Partial Functional-Differential Equations, Springer, 1996.

M. Ghergu, “Non-constant steady-state solutions for Brusselator type systems,” Nonlinearity, vol. 21, no. 10, pp. 2331–2345, 2008.

J. Fang, J. Wei, and X.-Q. Zhao, “Spatial dynamics of a nonlocal and time-delayed reaction-diffusion system,” Journal of Differential Equations, vol. 245, no. 10, pp. 2749–2770, 2008.

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

Y. You, “Global dynamics of the Brusselator equations,” Dynamics of Partial Differential Equations, vol. 4, no. 2, pp. 167–196, 2007.

C. C. Travis and G. F. Webb, “Existence, stability, and compactness in the -norm for partial functional differential equations,” Transactions of the American Mathematical Society, vol. 240, pp. 129–143, 1978.

M. C. Memory, “Stable and unstable manifolds for partial functional differential equations,” Transactions of the American Mathematical Society, vol. 200, pp. 395–418, 1974.

D. S. Mackey, N. Mackey, C. Mehl, and V. Mehrmann, “Vector spaces of linearizations for matrix polynomials,” SIAM Journal on Matrix Analysis and Applications, vol. 28, no. 4, pp. 971–1004, 2006.

J. S. Dos Santos and M. A. Bená, “The delay effect on reaction-diffusion equations,” Applicable Analysis, vol. 83, no. 8, pp. 807–824, 2004.

X.-P. Yan and W.-T. Li, “Stability of bifurcating periodic solutions in a delayed reaction-diffusion population model,” Nonlinearity, vol. 23, no. 6, pp. 1413–1431, 2010.

X. D. Lin, J. W.-H. So, and J. H. Wu, “Centre manifolds for partial differential equations with delays,” Proceedings of the Royal Society of Edinburgh A, vol. 122, no. 3-4, pp. 237–254, 1992.