Discrete Dynamics in Nature and Society
Volume 2019 (2019), Article ID 8140380, 6 pages
https://doi.org/10.1155/2019/8140380

Research Article

Existence of a Nontrivial Steady-State Solution to a Parabolic-Parabolic Chemotaxis System with Singular Sensitivity

Yingjie Zhu1,2,3

¹Northeast Institute of Geography and Agroecology, Chinese Academy of Science, Changchun, 130012, China
²Institute of Science, Changchun University, Changchun, 130022, China
³Department of Mathematics and Statistics, University of Ottawa, Ottawa, ON, K1N 6N5, Canada

Correspondence should be addressed to Yingjie Zhu; yingjiezhu2012@163.com

Received 23 August 2018; Accepted 25 October 2018; Published 1 January 2019

Academic Editor: Rodica Luca

Copyright © 2019 Yingjie Zhu. 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.

This paper establishes the existence of a nontrivial steady-state solution to a parabolic-parabolic coupled system with singular (or logarithmic) sensitivity and nonlinear source arising from chemotaxis. The proofs mainly rely on the maximum principle, the implicit function theorem, and the Hopf bifurcation theorem.

1. Introduction

Chemotaxis is the biased movement of cells toward the concentration gradient of a chemical. It plays a critical role in a wide range of biological phenomena. For example, cells migrate toward resources of food and stay away from harmful substances. The first mathematical model of chemotaxis was introduced by Patlak in [1] and Keller and Segel in [2]. There are numerous works dedicated to the analysis of chemotaxis models. For example, Othmer and Stevens in [3] modeled myxobacteria as individual random walkers and proposed a microscopic model based on a velocity jump process. By taking the parabolic limit of the microscopic model, they obtain the macroscopic chemotaxis model, which is the well-known Keller-Segel system

where is a bounded connected domain with a smooth boundary . The function denotes the cell density and represents the chemical concentration, for example, oxygen. The constant is called the chemotactic coefficient, and the sign of corresponds to chemoattraction if and chemorepulsion if . Parameters and are the diffusion coefficients of the cells and the chemical, respectively. The function represents the kinetic function describing production and degradation of cells, and is commonly referred to as the chemotactic potential function. Function describes the production and degradation of the chemical.

The existence of global solutions, blow-up, and traveling wave solutions to the chemotactic system (1) were extensively studied during the past four decades (see, e.g., [4–12] and references therein). The authors studied the roles of growth, death and random in promoting population persistence through band popation in [4], and the work was related to the significance of cell motility and chemotaxis in microbial ecology.

X.F. Wang addressed the trivial and nontrivial steady states with small and of the quasi-linear system (1) in [13]. In the absence of population dynamics (i.e., ), there have been extensive studies. The main feature of solutions to the Keller-Segel model is the possibility of blow-up in finite time in [9, 14, 15]. Moreover, recent results in [6, 10, 11, 16, 17] proved the global existence of solutions under some conditions. Moreover, the global existence, asymptotic behavior, and steady states of classical solutions were studied in [18] for the one-dimensional case.

The purpose of this paper is to study the existence of the nontrivial steady-state solution to a parabolic-parabolic coupled system arising from chemotaxis with singular sensitivity. We consider the following system

with the initial-boundary value conditions

where and “”. The function denotes the density of the cells and denotes the concentration of the chemical. The parameters , and are all positive constants, with being the diffusion coefficient of cells and the chemotactic sensitivity coefficient as above. The function is assumed to be smooth.

The steady-state problem corresponding to system (2)-(3) is

where . The linear chemotactic potential was considered in [19], and X.F. Wang’s model in [13] is similar to ours.

In this paper, we establish the existence of nontrivial steady-state solutions to the parabolic-parabolic coupled chemotactic system (2)-(3) with singular sensitivity. Here we assume that the function satisfies

Obviously, and are the trivial solutions to system (4).

Our main results in this paper are presented in the following theorem.

Theorem 1. Consider system (4) with condition (5). Then the following alternative holds:

(1)If , then system (4) has a pair of nonnegative steady-state solutions and (2)If , then there exists at least a pair of nontrivial steady-state solutions .

The rest of this paper is organized as follows. In Section 2, we give some preliminary lemmas. In Section 3, we complete the proof of Theorem 1.

2. Lemmas

To prove Theorem 1, we need to establish boundary estimates of solutions to system (4). Since the state variables represent densities, we only consider nonnegative solutions, that is, and on . Firstly, we state one result concerning the estimate of a solution .

Lemma 2. Assume the solution is a nonnegative steady state. Then we obtain the following:

(i), and .(ii).

Proof.

(i)By the second equation of (4), we have and, integrating once, we obtain . Hence, and .(ii)Integrating the first equation of (4) yields

Lemma 3. Let solve system (2) with the initial-boundary value condition of (3). Then we have the lower-bound estimate

and is defined in (11) below.

Proof. From the second equation in (4), we obtain

Integrating the above equation once, we have

Integrating by part and applying the condition , we rearrange the above equation

From the boundary value condition of in (3), we have . Define the function

From the definition of the function , we easily get

Applying the implicit function theorem, we have the estimate (7).

Lemma 4. There are two positive constants and such that a nonnegative steady state of system (4) satisfies

Proof. By the steady-state system (4), we know , and for every , we have

Applying the condition , we obtain the following equation

Similarly, we also get

On the other hand, integrating once the first equation of (4), we get

which is equal to

Then, applying the boundary condition and rewriting (18), we have the following estimate

Next, we review the condition , which implies that . Combined with the condition of in Lemma 3, inequality (19) can be rewritten as

Integrating (20) once, we have

Here and . Hence, we obtain

which implies that

Define the function . Then we have . Combining the above inequalities, we obtain the following bound

By Gronwall’s inequality, we get

Substituting the definition of into inequality (25), we obtain

Combining this result with in Lemma 2, we obtain estimate (13).

Next, we suppose that for any . Given , we consider the equation

Then, we prove the following lemma.

Lemma 5. For any small , we pick a sufficiently small such that . Next, we choose

and set

Finally, we denote by a nonnegative solution to system (4) when . Then we have .

Proof. If not, there exists some such that . According to Lemma 2 and the definition of , we have

From (13), we can easily get

However, .
In fact, if not we have

which contradicts the fact that .
Next, we consider

which is a contradiction. In other words, we have .

Lemma 6.

(i)For each , system (27) has a unique solution .(ii) is continuous.(iii)If , then .

Proof. For any given , it is obvious that and are a pair of sup-sub-solutions. According to the standard comparison theorem, we easily obtain the estimate of the solution to system (27). If and are solutions to system (27), then we have

The difference between the first and second equation of (34) is

with the notation .
By the maximum principle, we obtain . Moreover, is uniquely defined and .
Next we give the proof of (ii) and (iii) in Lemma 6, respectively.

(i)Assume that there exists a sequence , which converges to . Then we only point out that there exists a subsequence of such that on . By regularity theory, we know that is the solution to system (34) and , which contradicts the uniqueness of .(ii)We can directly prove (iii) from the maximum principle.

3. Proof of Theorem 1

We dedicate this section to the proof of Theorem 1.

Assume there exists such that satisfies the following system

It is easy to see that is the nonnegative solution to system (36) by the definition of and .

Define to be the unique solution to the system

where is a constant satisfying . For any given , we define to be the unique solution to the following system

Then we rewrite (36) and obtain the following operator equation

We can reformulate (39) as the following operator equation

Next, we introduce some notation:

It is easy to see that the operator is compact and linear.

Define the operator as the map

Similar to the proof of Lemma 5(i), we can show that the operator is continuous and bounded from to . We also have

Hence, is compact operator from to with the norm . Moreover, the operator is also compact from to and .

From the above lemmas, we can show the proof of Theorem 1.

Proof of Theorem 1. Assume that is the nonnegative solution to system (2). Then we have

However, from the proof of Lemma 3, we have

which is a contradiction. Thus, we obtain .
Assume that there is a constant such that .
Next, we choose a variable as the parameter of bifurcation. Furthermore, is a single-value characteristic root. Thus we know that system produces a bifurcation in a small neighbourhood of the point according to the classical bifurcation theorem. Then there is a bifurcating solution , represented as

where and is sufficiently small. Furthermore, and . Therefore, when is sufficient small, the bifurcation solution is in the set of

where is the characteristic eigenvector of characteristic eigenvalue and is a positive constant. There is a nontrivial connected set satisfying the following alternative conclusion:

(i) connects the two points and , where is another characteristic eigenvalue of .(ii) is unbounded and .

We firstly prove . Otherwise, there exists a point such that and

Choose a sufficient large constant . Then by (36), we can obtain

However, which means that is equal to zero at some point. We can directly have by the maximum principle. That is, is a bifurcation point of (40).
Denote . Because is a compact operator, there exists a convergent subsequence of . For the convenience, we still write it as which satisfies . Next, we divide by both sides of (40) and let go to infinity. Then we obtain

where . By the definition of we know that . However, has a unique characteristic eigenvalue with positive characteristic function. Thus, we have which contradicts the assumption. So conclusion (i) is impossible. In other words, is unbounded and .
According to (i) in Theorem 1, we know that if , then , which is equivalent to the condition . However, from Lemma 5, we have that is bounded for any . By Lemma 4, we have that is bounded. Furthermore, because the set is unbounded and connected, we can obtain that , which means that when , system (4) has nontrivial steady-state solution .
Finally, we show that system (2)-(3) has a positive solution when . In fact, we know that because . Next we prove that . If is zero at some point, then the initial value problem has a solution , which contradicts the boundary condition .

Data Availability

No data were used to support this study.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

This work is partially supported by the National Natural Science Foundation of China (No. 41701054), the China Postdoctoral Science Foundation (No. 2018M631890), the Scientific and Technological Research Project of Jilin Province’s Education Department (No. 2016285), Teaching Reform of Higher Education Research Project of Jilin Province (No. SJYB17-01). Finally, the author is grateful for the support by the China Scholarship Council (No. 201707535005).

References

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.

D. Lauffenburger, C. R. Kennedy, and R. Aris, “Traveling bands of chemotactic bacteria in the context of population growth,” Bulletin of Mathematical Biology, vol. 46, no. 1, pp. 19–40, 1984.

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,” Discrete Dynamics in Nature and Society, vol. 54, no. 5, 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.

S. Ai and Z. Wang, “Traveling bands for the Keller-Segel model with population growth,” Mathematical Biosciences and Engineering, vol. 12, no. 4, pp. 717–737, 2015.

X. Wang, “Qualitative behavior of solutions of chemotactic diffusion systems: effects of motility and chemotaxis and dynamics,” SIAM Journal on Mathematical Analysis, vol. 31, no. 3, pp. 535–560, 2000.

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.

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.

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.

K. Deep and M. Thakur, “A real coded multi parent genetic algorithms for function optimization,” Applied Mathematics and Computation, vol. 1, no. 2, pp. 67–83, 2008.