In this paper, we study the bifurcation and stability of solutions of the extended Fisher-Kolmogorov equation with periodic boundary condition. We prove that the system bifurcates from the trivial solution to an attractor as parameter crosses certain critical value. The topological structure of the attractor is also investigated.
MSC: 35B32, 35K35, 37G35.
Keywords:extended Fisher-Kolmogorov equation; periodic boundary condition; attractor bifurcation; center manifold
In this paper we work with the extended Fisher-Kolmogorov type equation with periodic boundary condition, which reads
The extended Fisher-Kolmogorov (EFK) equation has been proposed as a model for phase transitions and other bistable phenomena [1-3]. It has been extensively studied during past decades. Kalies and van der Vorst  considered the steady-state problem; by analyzing the variational structure, they proved the existence of heteroclinic connections, which are the critical points of a certain functional. Also, by the variational method, Tersian and Chaparova  derived the existence of periodic and homoclinic solutions. Peletier and Troy  were interested in the stationary spatially periodic patterns and showed that the structure of the solutions is enriched by increasing the coefficient of the fourth-order derivative term. The structure of the solution set was also investigated by van den Berg , who enumerated all the possible bounded stationary solutions provided this coefficient is small. Rottschäer and Wayne  showed that for every positive wavespeed there exists a traveling wave. And they also found the critical wavespeed to discriminate the monotonic solution from the oscillatory one. By an iteration procedure, Luo and Zhang  proved that equation (1.1) possesses a global attractor in the Sobolev space for all provided that and p is an odd number. We refer the interested readers to the references in [4-9] for other results on the EFK equation; see also, among others, [10-13].
Returning to problem (1.1), our main interest in the present paper is the bifurcation and stability of solutions. By using a notion of bifurcation called attractor bifurcation developed by Ma and Wang in [14,15], a nonlinear attractor bifurcation theory for this problem is established. Work on the topic of attractor bifurcation also can be seen in [16,17].
The main objectives of this theory include:
(1) existence of attractor bifurcation when the system parameter crosses some critical number,
(2) dynamic stability of bifurcated solutions, and
(3) the topological structure of the bifurcated attractor.
Our main results can be summarized as follows.
Moreover, we apply this theory to a model of the population density for single-species and derive biological results.
This article is organized as follows. The preliminaries are given in Section 2. The mathematical setting is presented in Section 3. The mathematical results are given in Section 4. In Section 5 we apply mathematical results to a model of the population density for single-species and derive biological results. In Section 6 we discuss some existing results and compare them with ours. Finally, Section 7 is devoted to the conclusions.
To prove the main result, we introduce an important theorem.
and the principle of exchange of stabilities holds true:
The following attractor bifurcation theorem can be found in .
(4) Ifis globally asymptotically stable for (2.1) at, then for any bounded open setwith, there is ansuch that, the attractorattractsinH, where Γ is the stable manifold ofwith codimensionm. In particular, if (2.1) has a global attractor for allλnear, then.
Remark 2.1 As and H are infinite dimensional Hilbert spaces, if (2.1) satisfies conditions (2.2)-(2.5) and is a locally (global) asymptotically stable equilibrium point of (2.1) at , then the assertions (1)-(4) of Theorem 2.1 hold; see [14,15].
To get the structure of the bifurcated solutions, we introduce another theorem.
Theorem 2.2 (Theorem 5.10 in )
3 Mathematical setting
4 Mathematical results
As mentioned in the introduction, we study in this manuscript attractor bifurcation of the EFK equation under the periodic boundary condition. Then we have the following bifurcation theorem.
Proof of Theorem 4.1 We shall prove Theorem 4.1 in four steps.
Step 1. In this step, we study the eigenvalue problem of the linearized equation of (3.2) and find the eigenvectors and the critical value of λ.
Consider the eigenvalue problem of the linear equation,
It is not difficult to find that the eigenvalues and the normalized eigenvectors of (4.1) are
under condition that we get the principle of exchange of stabilities
It is easy to get the following inequality:
then condition (2.3) holds true.
Then the reduction equations of (3.2) are as follows:
the first order approximation of (4.3) does not work. Now, we shall find out the second order approximation of (4.3). And the most important of all is to obtain the approximation expression of the center manifold function.
By direct calculation, we have
According to the formula of Theorem 3.8 in  (or Remark 4.1), the center manifold function Φ, in the neighborhood of , can be expressed as
By direct calculation, we have
Putting (4.4) into (4.3), we have the reduction equations
Since the following equality holds true:
Step 4. In the last step, we show that the bifurcated attractor of (3.2) consists of a singularity cycle.
Since the even function space is an invariant subspace of defined by (3.1), we shall consider the bifurcation problem in the even function space and prove that system (1.1) bifurcates from to two steady solutions. For any function v in the even function space can be expressed as follows:
by the Lyapunov-Schmidt reduction method used in Step 3, we can deduce that the reduction equation of (1.1) is as follows:
Since the solutions of (2.1) are translation invariant,
Then near , the center manifold function in Theorem 3.8 in  can be expressed as follows:
In this section, we apply Theorem 4.1 to a model of the population density for single-species as follows:
where μ, α are the diffusion coefficients, v is the population density for single-species, and , , , , . It is easy to see that , and . Inspired by the work of Murray , represents the birth rate, describes the intra specific competition, and stands for the emigration which arises from disease.
we derive the following system:
According to Remark 4.3, if the condition is satisfied, the conclusions of Theorem 4.1 for system (5.3) also hold true. Consequently, from the translation (5.2), we have the following results for (5.1).
(1) If, the steady stateis locally asymptotically stable (Figure 1).
(2) If, system (5.1) bifurcates from the solutionto an attractor. This implies that the stability will switch from the original state (i.e., ) to a new one (i.e., ) (Figure 1).
(3) is homeomorphic toand consists of exactly one cycle of steady solutions of (5.1) (Figure 1).
Figure 1. Bifurcation diagram for the model of the population density for single-species. (1) Bifurcation appears at . (2) Bifurcated attractor is the boundary of the shaded region. (3) The first horizontal solid line from above denotes that the solution is stable, and the horizontal dotted line means this solution is unstable.
Furthermore, Theorem 5.1 and the equality
yield the following biological results.
(1) The population of this single-species is a conservative quantity.
(2) If the birth rate is low, then the population density will keep a uniform spatial distribution (Figure 2(A)).
(3) If the birth rate becomes high enough, then the spatial distribution of the population density will not keep uniform but change periodically with space (Figure 2(B)).
Figure 2. The spatial distribution of the population density. (1) Figure 2(A) shows that the population density keeps a uniform spatial distribution when the birth rate is low. (2) Figure 2(B) shows that the population density changes periodically with space when the birth rate becomes high enough. (3) The area of the shaded regions stands for the population of this single-species. And the area of the shaded region in Figure 2(A) is equal to the area of the shaded region in Figure 2(B).
Taking , , in (1.1), Peletier and Troy  analyzed stationary antisymmetric single-bump periodic solutions. They found that the coefficient of the fourth-order derivative term μ played a role of system parameter. If , the family of periodic solutions is still very similar to that of the Fisher-Kolmogorov equations. However, if , different families of periodic solutions emerged.
Taking , in (1.1), and under hypothesis that , , for , Rottschäer and Wayne  showed that for every positive wavespeed, there exists a traveling wave. And they also found that there exists a critical wavespeed . If , the solution is monotonic; otherwise, the solution is oscillatory.
Unlike the work mentioned above, which focuses on the structure of solutions varying with the system parameter (μ or c), the manuscript presented here investigates the topological structure and the stability of solutions varying with the system parameter, i.e., λ. Firstly, if , the bifurcated attractor consists of the trivial solution; if , the bifurcated attractor consists of only one cycle of steady state solutions and is homeomorphic to . Secondly, if , the trivial solution is locally asymptotically stable. However, if , the stability switches from the trivial solution to the bifurcated attractor.
Since the increment of dimension of spatial domain may lead to much richer bifurcated behavior, further investigation on higher dimension of spatial domain is necessary in the future.
In this article, we first prove the existence of attractor bifurcation when the system parameter crosses critical number , which is the first eigenvalue of the eigenvalue problem of the linearized equation of (1.1). Second, we show that the stability of solutions varies with the system parameter λ. If , the trivial solution is locally asymptotically stable. However, if , the stability switches from to . Third, the topological structure of the attractor is investigated. We prove that the attractor consists of only one cycle of steady state solutions and is homeomorphic to . At last, the expression of bifurcated solution is also obtained.
The authors declare that they have no competing interests.
Both authors read and approved the final manuscript.
The authors are grateful to the anonymous referees whose careful reading of the manuscript and valuable comments were very helpful for revising and improving our work.
Kalies, WD, van der Vorst, RCAM: Multitransition homoclinic and heteroclinic solutions of the extended Fisher-Kolmogorov equation. J. Differ. Equ.. 131, 209–228 (1996). Publisher Full Text
Tersian, S, Chaparova, J: Periodic and homoclinic solutions of extended Fisher-Kolmogorov equations. J. Math. Anal. Appl.. 260, 490–506 (2001). Publisher Full Text
Peletier, LA, Troy, WC: Spatial patterns described by the extended Fisher-Kolmogorov equation: periodic solutions. SIAM J. Math. Anal.. 28, 1317–1353 (1997). Publisher Full Text
van den Berg, JB: The phase-plane picture for a class of fourth-order conservative differential equations. J. Differ. Equ.. 161, 110–153 (2000). Publisher Full Text
Luo, H, Zhang, Q: Regularity of global attractor for the fourth-order reaction-diffusion equation. Commun. Nonlinear Sci. Numer. Simul.. 17, 3824–3831 (2012). Publisher Full Text
Kwapisz, J: Uniqueness of the stationary wave for the extended Fisher-Kolmogorov equation. J. Differ. Equ.. 165, 235–253 (2000). Publisher Full Text
Bartuccelli, MV: On the asymptotic positivity of solutions for the extended Fisher-Kolmogorov equation with nonlinear diffusion. Math. Methods Appl. Sci.. 25, 701–708 (2002). Publisher Full Text
Benguria, RD, Depassier, MC: On the transition from pulled to pushed monotonic fronts of the extended Fisher-Kolmogorov equation. Physica A. 356, 61–65 (2005). Publisher Full Text
Zhang, YD, Song, LY, Axia, W: Dynamical bifurcation for the Kuramoto-Sivashinsky equation. Nonlinear Anal.. 74(4), 1155–1163 (2011). Publisher Full Text