Research Article | Open Access
Global Existence to an Attraction-Repulsion Chemotaxis Model with Fast Diffusion and Nonlinear Source
This paper deals with the global existence of solutions to a strongly coupled parabolic-parabolic system of chemotaxis arising from the theory of reinforced random walks. More specifically, we investigate the attraction-repulsion chemotaxis model with fast diffusive term and nonlinear source subject to the Neumann boundary conditions. Such fast diffusion guarantees the global existence of solutions for any given initial value in a bounded domain. Our main results are based on the method of energy estimates, where the key estimates are obtained by a technique originating from Moser’s iterations. Moreover, we notice that the cell density goes to the maximum value when the diffusion coefficient of the cell density tends to infinity.
Chemotaxis is known as the active orientation of moving organisms along the chemical gradient. It is observed in many natural systems. For example, myxobacteria produce so-called slime trails on which their cohorts can move more readily. The mathematical models of chemotaxis were introduced by Patlak in  and Keller and Segel in . During the past four decades, chemotaxis models have been studied extensively (see, e.g., [3–10] and the rich references therein). For instance, Othmer and Stevens in  modeled myxobacteria as individual random walkers and proposed the microscopic model based on the velocity jump process. By taking the parabolic limit of microscopic model, we can obtain the macroscopic chemotaxis model which is the well-known Keller-Segel system:where is a bounded connected domain with a smooth boundary , , , and are positive constants, and . The function denotes the cell density and represents the chemotactic concentration, for example, the oxygen. The constant is called the chemosensitive coefficient, and the sign of corresponds to chemoattraction if , and chemorepulsion if . The function and the constant are the diffusion coefficient of the cell motility and the chemical, respectively. The function represents the kinetic function describing production and degradation of chemicals, and is commonly referred to the chemotactic potential function.
In the absence of logistic source (i.e., ), there have been extensive studies to system (1). The main feature of solution to the Keller-Segel model is the possibility of blowup in finite time in [3, 8, 11, 12]. For instance, the first result on finite-time blowup for a radially symmetric solutions was shown in [13, 14] when is a ball in under the condition that and is sufficiently small. For a general domain, Nagai in  further showed the finite-time blowup of nonradial solutions provided that is sufficiently small and if or if with . Winkler in  studied finite-time blowup of radially symmetric solutions to the full parabolic system with logistic sources in a ball with parameters , , and . Moreover, Winkler in  gave the set of blowup to the system (1) by enforcing initial data with respect to the topology of for any , where is a ball in with . In summary, the solution of system (1) never blows up when , whereas there is finite-time or infinite-time blowup when . Moreover, recent results in  confirmed that the attraction-repulsion is a plausible mechanism to regularize the classical Keller-Segel system (2) with whose solutions may blow up in higher dimensions. Authors in [5, 9, 10, 18, 19] proved the global existence of solutions.
However, in many biology progresses, the cells usually interact with not only the attractive combination but also repulsive signalling. Therefore, it is necessary to study the attraction-repulsion chemotaxis model:System (2) with was proposed in  to describe the aggregation of microglia observed in Alzheimer’s disease. Here represents the cell density and denotes the concentration of the chemoattractant and is the concentration of the chemorepulsion. The constants satisfy and , where and measure the strength of chemotactic signal of attraction and repulsion, respectively. The constants , , , and are positive, and they denote the production and degradation rates of the two chemicals, respectively. The first cross-diffusive term in the first equation of (2) means that the orientation movement of the cell is directed to the chemorepulsion, whereas the second cross-diffusive term implies that cells move down the chemoattraction. The second and third equations in (2) elaborate that the two chemicals of chemoattraction and chemorepulsion are released by cells and go through decay. For the case of , the theorem of competing effects has been established in  with . Moreover, the global existence, asymptotic behavior, and steady states of classical solution were studied in  for one-dimensional case with .
In this paper, we consider the following attraction-repulsion chemotaxis system including three parabolic equations:with the initial-boundary value conditions:where is a bounded domain for some and denotes the unit outer normal vector field on . The function of the cell density is the fraction of volume occupied by cells, whereas the fast diffusion coefficient of the cell is described asIt is easy to see that the function is a monotonically increasing function of which guarantees that the solution will not blow up. The functions and denote the concentration of the chemoattraction and chemorepulsion, respectively. The constants and are the diffusion coefficients of the chemoattraction and chemorepulsion, respectively. and denote the chemosensitive coefficients, and describes the kinetic function.
In particular, the global existence of solutions was obtained by increasing the diffusion coefficient : for example, , in , and in . Liu and Tao in  established global boundedness of classical solutions to the parabolic attraction-repulsion chemotaxis system (2) when . For the case of including logistic source, global existence and boundedness of classical solution were studied in [23, 24]. For the case , the existence of families of traveling impulses and fronts was analyzed in .
In this paper, we will prove the existence of global classical solutions to the generalized system (3) with the initial-boundary conditions and the no-flux boundary condition (4). In addition, we need the following assumptions.(A1)Let . Suppose that functions are of class, and for and for .(A2)Assume that functions are of class and that , , and for all , where .(A3)Let be of class, where .(A4), , and for and for .The main result is the following.
Throughout this paper, we introduce some notations. is a bounded open interval in and denotes a general constant which may have different values in different place. is the usual Lebesgue space with the norm for and When , we denote for convenience. is the th-order Sobolev space with the norm For notational convenience, we write and as and , respectively. Moreover, we denote
The rest of this paper is organized as follows. In Section 2, we establish the local existence and give some preliminary lemmas. Some necessary priori estimates will be established in Section 3. We will complete the proof of Theorem 1 in Section 4.
3. Local Existence
Theorem 2 (local existence). Let assumptions (A1)–(A4) hold. Then(1)There is a positive constant depending on initial value such that system (3) has a unique maximal solution in the space with .(2)There is a global classical solution of system (3) if is bounded away from 1 for each finite time , which means .
Proof. Let . System (3) can be written aswhereSystem (9) satisfies the initial valueFor some and given initial conditions , it is clear that the eigenvalues of matrix are positive at . Thus, system (9) is normally parabolic and there exists local solution by Theorem 7.3 in ; that is, there is a such that the unique solution .
Next we rewrite the first equation of (3) as Recalling that assumption (A2) ensures that for all and for all , this means that for all and , where and . Hence, the source term is zero when at some . By applying the maximum principle to , therefore, we obtain that whenever and . Similarly, we infer (or ) from the second equation of (3) (or the third of (3)) whenever . The proof of (1) in Theorem 2 is finished.
The proof of (2) in Theorem 2 is completed by applying the theorem of quasilinear parabolic equations in  since system (9) is an upper triangular system.
4. A Priori Estimates
Lemma 3. Let be a bounded open domain satisfying the uniform cone property in with :Then there exists a positive constant , which depends on , , and , such that for all : where and .
Lemma 4. Let the conditions in Lemma 3 hold: and. Then there exists a constant independent of such that
Proof. Recalling for and employing Lemma 3 with and , we obtain ; thus,where we have used the inequality for .
Lemma 5. Let a nonnegative numerical sequence satisfying with and for . Suppose and . Then one has
Proof. By the definition of the sequence , we obtain that for
Lemma 6. Let conditions (A1)–(A4) hold. Then there exists a unique global solution to system (3) such that , , and are in Furthermore, there exists two positive constants and satisfying
Proof. After appropriate scaling, the system is limited in a region satisfying . Next, we introduce the auxiliary scalar equation:with the initial-boundary condition:where , , and are given functions. To complete the proof of Lemma 6, we will establish the following lemma.
Lemma 7. Suppose for all . LetLet be a classical solution to system (20)-(21) satisfying for . Then, for any , there exists a constant such that for all , where the constant depends only on , , , and .
Proof. We take satisfying for all and notice that by scaling and in the governing equation. To complete the proof, we have to divide the proof into three steps.
Step 1. We will firstly prove that for any . Multiplying both sides of the first equation of (3) by and integrating the resultant equation, we yieldwhere the last inequality follows from (A3); then (23) is rewritten asNoticing , we obtain the fact thatSubstituting (24) into (25), we havewhere we have applied Young’s inequality with .
Combing with the assumption of Lemma 6, and , we obtain the following inequality:where .
For convenience, we define . Therefore, for some constant , we have from which we yield that for some positive constant which is independent of . Thus, for any , we conclude that which means that for any and . We notice that depends only on and and .
Step 2. Choose a satisfying and . Then we obtain from (25) thatRecalling -Young’s inequality , , we havewhere and Combining (31) with (32), we yieldwhere . We notice that Hence, we choose a sufficiently small such thatFurthermore, we have the estimatewhose proof is similar to Lemma 4. Thus, we omit it. Squaring (36), we obtain that Thus, That is,Substituting the above inequality into (39), we havefor some constants .
Step 3. We notice that , and define We also see that is bounded and depends only on , , , and . It is easy to observe that is nondecreasing in with and is a nondecreasing function with respect to . Combining the above analysis, we yield That is,Integrating (43) in , we haveAssuming by the definition of , we getwhereFrom (44)–(46), we obtain thatBy a similar argument as in , we immediately yield that . Thus, the proof of Lemma 7 is completed.
5. Proof of Theorem 1
In this section, we will complete the proof of Theorem 1 by the local existence and some priori estimates as given below.
Proof. Suppose there exists a solution for the maximal time . Then according to Lemma 6 we know . Thus, if we treat as a source term in the second and third equations of (9), we can obtain thatBy employing the estimate for parabolic equation in  when goes to infinity, we obtain the norm bound of , , and , . So we have being infinite from the second conclusion in Theorem 1 which contradicts our original assumption. That is, the maximal time .
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors are very grateful to Professor Yong Li for providing valuable advice. Moreover, they very thankful to the anonymous referees for the careful reading and various corrections which greatly improved the exposition of the paper. This work partially supported by the NSF of China (under Grants nos. 11171350 and 11204019), the Scientific and Technological Research Project of Jilin Province’s Education Department (nos. 2013287 and 2014312), the Youth Project of Jilin Province’s Science and Technology Department (nos. 20130522099JH and 201201140), and the Twelfth Five-Year Plan Project of Jilin Province’s Educational Science (no. ZD2014078).
- C. S. Patlak, “Random walk with persistence and external bias,” Bulletin of Mathematical Biophysics, vol. 15, pp. 311–338, 1953.
- E. F. Keller and L. A. Segel, “Traveling bands of chemotactic bacteria: a theoretical analysis,” Journal of Theoretical Biology, vol. 30, no. 2, pp. 235–248, 1971.
- H. G. Othmer and A. Stevens, “Aggregation, blowup, and collapse: the ABCs of taxis in rein-forced random walks,” SIAM Journal on Applied Mathematics, vol. 57, no. 4, pp. 1044–1081, 1997.
- T. Hillen and K. J. Painter, “A user's guide to PDE models for chemotaxis,” Journal of Mathematical Biology, vol. 58, no. 1-2, pp. 183–217, 2009.
- M. Alber, R. Gejji, and B. Kazmierczak, “Existence of global solutions of a macroscopic model of cellular motion in a chemotactic field,” Applied Mathematics Letters, vol. 22, no. 11, pp. 1645–1648, 2009.
- B. D. Sleeman, M. J. Ward, and J. C. Wei, “The existence and stability of spike patterns in a chemotaxis model,” SIAM Journal on Applied Mathematics, vol. 65, no. 3, pp. 790–817, 2005.
- T. Hillen and H. A. Levine, “Blow-up and pattern formation in hyperbolic models for chemotaxis in 1-D,” ZAMP—Zeitschrift für Angewandte Mathematik und Physik, vol. 54, pp. 839–868, 2003.
- M. Winkler, “Finite-time blow-up in the higher-dimensional parabolic-parabolic Keller-Segel system,” Journal de Mathématiques Pures et Appliquées, vol. 100, no. 5, pp. 748–767, 2013.
- C. Mu, L. Wang, P. Zheng, and Q. Zhang, “Global existence and boundedness of classical solutions to a parabolic-parabolic chemotaxis system,” Nonlinear Analysis: Real World Applications, vol. 14, no. 3, pp. 1634–1642, 2013.
- Y. Tao and M. Winkler, “Eventual smoothness and stabilization of large-data solutions in a three-dimensional chemotaxis system with consumption of chemoattractant,” Journal of Differential Equations, vol. 252, no. 3, pp. 2520–2543, 2012.
- Y. Sugiyama and Y. Yahagi, “Extinction, decay and blow-up for Keller-Segel systems of fast diffusion type,” Journal of Differential Equations, vol. 250, no. 7, pp. 3047–3087, 2011.
- T. Cieślak and C. Stinner, “Finite-time blowup and global-in-time unbounded solutions to a parabolic-parabolic quasilinear Keller-Segel system in higher dimensions,” Journal of Differential Equations, vol. 252, no. 10, pp. 5832–5851, 2012.
- W. Jáger and S. Luckhaus, “On explosions of solutions to a system of partial differential equations modelling chemotaxis,” Transactions of the American Mathematical Society, vol. 329, no. 2, pp. 819–824, 1992.
- T. Nagai, “Blow-up of radially symmetric solutions to a chemotaxis system,” Advances in Mathematical Sciences and Applications, vol. 5, no. 2, pp. 581–601, 1995.
- T. Nagai, “Blowup of nonradial solutions to parabolic-elliptic systems modeling chemotaxis in two-dimensional domains,” Journal of Inequalities and Applications, vol. 6, no. 1, pp. 37–55, 2001.
- M. Winkler, “Blow-up in a higher-dimensional chemotaxis system despite logistic growth restriction,” Journal of Mathematical Analysis and Applications, vol. 384, no. 2, pp. 261–272, 2011.
- Y. S. Tao and Z.-A. Wang, “Competing effects of attraction vs. repulsion in chemotaxis,” Mathematical Models and Methods in Applied Sciences, vol. 23, no. 1, pp. 1–36, 2013.
- Y.-S. Choi and Z.-A. Wang, “Prevention of blow-up by fast diffusion in chemotaxis,” Journal of Mathematical Analysis and Applications, vol. 362, no. 2, pp. 553–564, 2010.
- Y. S. Tao and M. Winkler, “Boundedness in a quasilinear parabolic-parabolic Keller-Segel system with subcritical sensitivity,” Journal of Differential Equations, vol. 252, no. 1, pp. 692–715, 2012.
- M. Luca, A. Chavez-Ross, L. Edelstein-Keshet, and A. Mogilner, “Chemotactic signaling, microglia, and Alzheimer's disease senile plaques: is there a connection?” Bulletin of Mathematical Biology, vol. 65, no. 4, pp. 693–730, 2003.
- J. Liu and Z.-A. Wang, “Classical solutions and steady states of an attraction-repulsion chemotaxis in one dimension,” Journal of Biological Dynamics, vol. 6, pp. 31–41, 2012.
- D.-M. Liu and Y.-S. Tao, “Global boundedness in a fully parabolic attraction-repulsion chemotaxis model,” Mathematical Methods in the Applied Sciences, 2015.
- L. C. Wang, C. L. Mu, and P. Zheng, “On a quasilinear parabolic-elliptic chemotaxis system with logistic source,” Journal of Differential Equations, vol. 256, no. 5, pp. 1847–1872, 2014.
- X. R. Cao and S. N. Zheng, “Boundedness of solutions to a quasilinear parabolic-elliptic Keller-Segel system with logistic source,” Mathematical Methods in the Applied Sciences, vol. 37, no. 15, pp. 2326–2330, 2014.
- F. S. Berezovskaya, A. S. Novozhilov, and G. P. Karev, “Families of traveling impulses and fronts in some models with cross-diffusion,” Nonlinear Analysis: Real World Applications, vol. 9, no. 5, pp. 1866–1881, 2008.
- H. Amann, “Dynamic theory of quasilinear parabolic equations. II: reaction-diffusion systems,” Differential and Integral Equations, vol. 3, no. 1, pp. 13–75, 1990.
- O. A. Ladyzhenskaya, V. A. Solonnik, and N. N. Ural'ceva, Linear and Quasi-Linear Equations of Parabolic Type, Translations of Mathematical Monographs, American Mathematical Society, 1968.
- L. Nirenberg, “An extended interpolation inequality,” Annali Della Scuola Normale Superiore di Pisa-Classe di Scienze, vol. 20, pp. 733–737, 1966.
Copyright © 2015 Yingjie Zhu and Fuzhong Cong. 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.