Next Article in Journal
Heavy-Ball-Based Hard Thresholding Pursuit for Sparse Phase Retrieval Problems
Next Article in Special Issue
Lump-Type Solutions, Lump Solutions, and Mixed Rogue Waves for Coupled Nonlinear Generalized Zakharov Equations
Previous Article in Journal
The Study of Bicomplex-Valued Controlled Metric Spaces with Applications to Fractional Differential Equations
Previous Article in Special Issue
Analytical Solutions to the Chavy-Waddy–Kolokolnikov Model of Bacterial Aggregates in Phototaxis by Three Integration Schemes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Boundedness in a Logarithmic Keller–Segel System

1
School of Mathematics, Harbin Institute of Technology, Harbin 150001, China
2
School of Statistics, Chengdu University of Information Technology, Chengdu 610103, China
3
School of Economics and Management, Zhejiang Sci-Tech University, Hangzhou 310018, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2743; https://doi.org/10.3390/math11122743
Submission received: 24 March 2023 / Revised: 7 June 2023 / Accepted: 12 June 2023 / Published: 16 June 2023
(This article belongs to the Special Issue Applications of Partial Differential Equations)

Abstract

:
In this paper, we propose a user-friendly integral inequality to study the coupled parabolic chemotaxis system with singular sensitivity under the Neumann boundary condition. Under a low diffusion rate, the classical solution of this system is uniformly bounded. Our proof replies on the construction of the energy functional containing Ω | v | 4 v 2 with v > 0 . It is noteworthy that the inequality used in the paper may be applied to study other chemotaxis systems.

1. Introduction

Our work considers the coupled parabolic chemotaxis system with singular sensitivity
u t = · ( u χ u v v ) , x Ω , t > 0 , v t = k Δ v v + u , x Ω , t > 0 , ν u = ν v = 0 , x Ω , t > 0 , u ( x , 0 ) = u 0 , v ( x , 0 ) = v 0 , x Ω ,
for parameters χ , k > 0 with the Neumann boundary condition, where Ω R 2 is a bounded domain with a smooth boundary. u and v are the cell density and concentration of chemical stimulus with respect to time t and x, respectively. k represents the diffusion rate of the chemical signal. The initial functions u 0 C 0 ( Ω ¯ ) and v 0 W 1 , ( Ω ) satisfy u 0 0 and v 0 > 0 .
In 1970, Keller and Segel [1] originally introduced the system
u t = · ( u u χ ( v ) v ) , x Ω , t > 0 , τ v t = k Δ v α v + β u , x Ω , t > 0 ,
to describe chemotaxis, the oriented movement of cells in response to the concentration of chemical signal produced by themselves and self-diffusion, where τ , k , α , β > 0 are parameters. The chemical signal experiences random diffusion and decay. Particular cases and derivatives of chemotaixs models have been developed extensively, such as the parabolic–elliptic case [2,3,4,5], the fully parabolic case [6,7,8,9,10] and other extensive versions [11,12,13]. Some studies have focused on the problem of whether the solution to the respective model undergoes a chemotactic collapse in the sense that the cell density becomes unbounded in finite or infinite time [3,6,7,12]. Given the initial conditions u 0 0 , v 0 > 0 and the Neumann boundary conditions, others have concentrated on the aggregation effect of the chemotactic sensitivity χ ( v ) .
If χ ( v ) = χ with τ = k = α = β = 1 , Osaki and Yagi [14] showed the global boundedness of solutions to (2) for n = 1 and Nagai et al. [15] proved the results if Ω u 0 < 4 π for n = 2 . For n 3 , if u 0 L n 2 ( Ω ) is small enough, there exist global weak solutions [16]. Another form of sensitivity function is
χ ( v ) = χ 0 ( c + α v ) k
for c , χ 0 > 0 , k > 1 and α > 0 , which is non-singular. In this case, the global existence is established for k = 2 , c = 1 by [17] and for k = 1 , α = 1 by [12]. Furthermore, if χ ( v ) = χ 0 v k for k > 1 , χ 0 > 0 , there exist global classical solutions to (2) [18].
The logarithmic sensitivity function χ ( v ) = χ v with χ > 0 is commonly considered because it is in compliance with the Weber–Fechner law [19]. Taking this form with τ = k = α = β = 1 , the chemotaxis model becomes the classical version:
u t = · ( u χ u v v ) , x Ω , t > 0 , v t = Δ v v + u , x Ω , t > 0 , ν u = ν v = 0 , x Ω , t > 0 , u ( x , 0 ) = u 0 , v ( x , 0 ) = v 0 , x Ω .
Global bounded solutions to (3) are provided by Osaki and Yagi [14] in a one-dimensional case. As for n = 2 , Lankeit [7] introduced an energy functional and proved that the solutions are uniform bounded in a convex domain with the range of χ extending to slightly more than one. Moreover, Winkler [20] proved that there exist global classical solutions if 0 < χ < 2 n , and Fujie [6] showed the solutions are uniformly time bounded. In [21], global bounded solutions are constructed under the the condition of χ 4 n with Ω R n being the convex domain. Furthermore, (3) employs global weak solutions when χ < n + 2 3 n 4 [20]. In the radially symmetric setting, weak solutions are constructed by [22] under the condition χ < n n 2 . These results imply that there is a balance between χ and dimension n for the establishment of global solutions to classic models (3). The work to extend both χ and n is laborious without giving any condition of (3). Lankeit and Winkler [23] extended the range of χ to
χ < if n = 2 8 if n = 3 n n 2 if n 4 ,
under the definition of the generalized solution, which is constructed on the basis of the global weak solution.
There are also other results established on the changing of parameters, referring to [9,24]. Indeed, the parameters in (2) have an impact on the aggregation of cell density. Xiangdong [25] constructed global solutions to (1) with n 8 under some conditions, where the relationship between k and χ is established. However, if n = 2 , the diffusion rate of the concentration of chemicals k does not work, since χ is still less than one, as in [25].
In [26], the estimates containing Ω | v | 2 are established to study the system where the chemotactic sensitivity is a constant and the source of the signal is modeled by v. In the work of Winkler [27], the only evident global quasi-dissipative structure involving Ω | v | 2 v 2 , ( v > 0 ) is established to address the difficulty brought about by the nonlinear source of signal. However, the system with logarithmic sensitivity presents more challenges, and the structure of Ω f ( v ) | v | n (n is even) is essential to the estimates. Hence, motivated by Lankeit [7] and Nagai [15], we establish an energy-type functional containing
Ω | v | 4 v 2 .
The fractional term of v in the energy-type functional may alleviate the difficulty of preventing the aggregation caused by nonlinear kinetics in some derivate systems such as [27,28], where the source of the signal is modeled by u v .
In this paper, the global existence and uniform boundedness of the classical solutions of (1) are established as follows:
Theorem 1.
Let Ω be a bounded domain with a smooth boundary Ω on R 2 , initial data v 0 > 0 and u 0 , 0 in Ω with u 0 C 0 ( Ω ¯ ) and v 0 W 1 , ( Ω ) . For all χ > 0 , there exists a constant C k that depends on u 0 , v 0 , Ω and χ, such that whenever
k C k ,
then (1) admits a unique classical solution ( u , v ) C 0 ( Ω ¯ × [ 0 , ) ) C 2 , 1 ( Ω ¯ × ( 0 , ) ) . Moreover, there exist constants δ , C > 0 such that δ v < C and 0 u < C for all t ( 0 , ) .
Intuitively, this shows that the large diffusion rate of chemical signals can prevent the aggregation of cell density resulting from a large χ .
In the paper, we first demonstrate the local existence of and recall some inequalities in the preliminaries. Then, we prove our key integral inequality in the Section 3 and give some useful a priori estimates in the Section 4. Finally, we prove the uniform boundedness of the solutions.

2. Preliminaries

2.1. Local Existence

The local existence of classical solutions to chemotaxis systems has been well-established using the methods of standard parabolic regularity theory and an appropriate fixed-point framework, which is shown in the following. Details of proof can be seen in Theorem 2.1 of [7] or [20].
Proposition 1.
Let Ω R n be a bounded domain with a smooth boundary, and u 0 C 0 ( Ω ¯ ) and v 0 W 1 , q ( Ω ) , q > n 1 are non-negative; then, for any k , χ > 0 , there exists T max ( 0 , ] and a pair of unique non-negative solutions satisfying
u C 0 ( Ω ¯ × [ 0 , T max ) ) C 2 , 1 ( Ω ¯ × ( 0 , T max ) ) , v C 0 ( Ω ¯ × [ 0 , T max ) ) C 2 , 1 ( Ω ¯ × ( 0 , T max ) ) L l o c ( [ 0 , T max ) ; W 1 , q ( Ω ) ) ,
such that ( u , v ) solves (1) classically in Ω × [ 0 , T max ) and, moreover, if T max < , then lim t T max u ( · , t ) L ( Ω ) + v ( · , t ) W 1 , q ( Ω ) = .

2.2. The Positive Lower Boundedness of v

In order to prove the lower boundedness of v in (1), we first prove the boundedness of u L 1 and v L 1 . Integrating the first and the second PDE in (1), we have the mass identities
Ω u = Ω u 0 = : m , t > 0
and
Ω v = Ω u 0 + Ω v 0 Ω u 0 · e t , t > 0 .
Based on these facts, one can deduce the non-negative lower boundedness of v from the abstract representation formula of the v equation. Copying Lemma 2.2 of [7], we write it as follows:
Lemma 1.
Let ( u , v ) satisfy Proposition 1; then, there exists T max > 0 and a positive constant δ depending on v 0 such that
v ( x , t ) δ > 0 , ( x , t ) Ω ¯ × [ 0 , T max ) .
Proof. 
Firstly, by the comparison principle and the fact of v 0 > 0 on Ω ¯ , we have for a small t
v ( x , t ) min x Ω ¯ v 0 · e t > 0 .
Let us fix τ = τ ( u 0 , v 0 ) . Then, it follows that
v ( x , t ) min x Ω ¯ v 0 · e τ : = δ 1 > 0 , t [ 0 , τ ) .
Now, from the well-known Neumann heat semigroup estimate for e t Δ (see Lemma 1.3 in [29] and Lemma 2.2 in [20]), we denote by d the diameter of the Ω and have for Ω R 2 that
( e t Δ ω ) 1 4 π t e d 2 4 t · Ω ω > 0 , ω C 0 ( Ω ) .
Then, the abstract representation formula of v shows
v ( · , t ) = e t ( Δ 1 ) v 0 + Ω e ( t s ) ( Δ 1 ) u ( · , t ) d s 0 t 1 4 π ( t s ) e ( ( t s ) + d 2 4 ( t s ) ) ( Ω u ( · , t ) ) d s m 0 t 1 4 π r e ( r + d 2 4 r ) d r : = δ 2 > 0 , t [ τ , ) ,
where r : = t s . Choosing δ = min { δ 1 , δ 2 } , we deduce (4). □

2.3. Recall of Useful Theorems

The well-known general Young’s inequality [30] is recalled.
Lemma 2.
Let f , g 0 be the continuous function with p , q > 0 satisfying 1 p + 1 q = 1 , then
f g ϵ f p + 1 q ( ϵ p ) q p g q
holds for all ϵ > 0 . Moreover, for continuous h > 0 and any ϵ 1 , ϵ 2 > 0 , taking p = 2 , q = 3 , r = 6 such that 1 p + 1 q + 1 r = 1 , we have
f g h ϵ 1 f 2 + ϵ 2 4 ϵ 1 g 3 + 6 36 ϵ 1 ϵ 2 h 6 .
Proof. 
In (7) is given the result of the straightforward calculation of the well-known inequality (6). □
Lemma 3.
Let Ω R n , n 1 be a smooth bounded domain. Any function f C 2 ( Ω ) satisfies
i . | f | 2 = 2 f · D 2 f ,
i i . ( Δ f ) 2 n | D 2 f | 2 ,
i i i . f · Δ f = 1 2 Δ | f | 2 | D 2 f | 2 .
All the identities and inequalities in the above lemma can be obtained from straightforward calculation. One can see [7,31] and Lemma 3.1 in [8] for their application. We could not find a precise reference in the literature that covers all that is necessary for our purpose; therefore, we conclude with a short lemma here.

3. A User-Friendly Integral Inequality

The proof of Theorem 1 is based on the extension and application of an integral inequality, which is generated within one dimension by Q. Wang [28]. The following theorem has a multidimensional form. It is worth noting that the integral inequality connects the fraction of the gradient and the second derivative. A similar inequality can be found in [7]. Furthermore, the explicit coefficient in the integral inequality is easy to use for readers.
Theorem 2.
Let Ω R n be a smooth bounded domain with w > 0 satisfying w C 2 ( Ω ¯ ) and w ν = 0 on Ω . Then,
Ω | w | 2 p + 2 w q + 2 n + 4 p ϵ 2 q + 1 p ϵ Ω | D 2 w | 2 | w | 2 p 2 w q
for all p 1 , q > 1 2 and ϵ > p 2 q + 1 > 0 .
Proof. 
Let J : = Ω | Δ log w | 2 | w | 2 p 2 w q 2 > 0 for p 1 . Directly calculating | Δ log w | 2 leads to
J = Ω | Δ w | 2 | w | 2 p 2 w q 2 Ω | w | 2 p Δ w w q + 1 J 0 + Ω | w | 2 p + 2 w q + 2 .
Since w ν = 0 on Ω , integration by parts gives
J 0 = 2 Ω | w | 2 p · w w q + 1 2 ( q + 1 ) Ω | w | 2 p + 2 w q + 2 = 2 p Ω | w | 2 p 2 | w | 2 · w w q + 1 2 ( q + 1 ) Ω | w | 2 p + 2 w q + 2 .
By (8) of Lemma 3 and (6), we have for ϵ > 0 that
J 0 = 4 p Ω | w | 2 p · D 2 w w q + 1 2 ( q + 1 ) Ω | w | 2 p + 2 w q + 2 4 p ϵ Ω | w | 2 p 2 | D 2 w | 2 w q 2 ( q + 1 ) p ϵ Ω | w | 2 p + 2 w q + 2 .
By (9), substituting (13) into (12) gives
J ( n + 4 p ϵ ) Ω | w | 2 p 2 | D 2 w | 2 w q ( 2 q + 1 ) p ϵ Ω | w | 2 p + 2 w q + 2 .
Due to q > 1 2 , ϵ > p 2 q + 1 > 0 ; thus, ( 2 q + 1 ) p ϵ > 0 , and we conclude with (11). □
Remark 1.
Letting Ω R 2 and taking q = p = 2 , ϵ > 2 5 , then n + 4 p ϵ 2 q + 1 p ϵ = 2 + 8 ϵ 5 2 ϵ . Note that 2 + 8 ϵ 5 2 ϵ achieves its global minimum over ( 2 5 , ) at ϵ = 4 + 26 10 ( 0.9099) with the value 2 21 4 26 ( 3.3117). Therefore,
Ω | w | 6 w 4 2 21 4 26 Ω | D 2 w | 2 | w | 2 w 2 .

4. Some Useful A Priori Estimates

Let us first give an inequality to estimate the boundary integration.
Lemma 4.
Let Ω R 2 be a bounded smooth domain. If v C 2 ( Ω ¯ ) satisfies v ν = 0 , the for any ϵ ^ > 0 , there exists C ( ϵ ^ ) depending on Ω such that
Ω | v | 2 v 2 ( | v | 2 ) ν ϵ ^ Ω | v | 2 | D 2 v | 2 v 2 + C ( ϵ ^ )
for all t ( 0 , T max ) and n 1 .
Proof. 
Firstly, we show that
Ω | v | 2 v 2 ( | v | 2 ) ν = 16 Ω | v | 2 ( | v | 2 ) ν .
From the Neumann boundary condition, we calculate the right-hand side, obtaining
16 Ω | v | 2 ( | v | 2 ) ν = Ω | v | 2 v | v | 2 v ν = Ω | v | 2 v | v | 2 ν v v ν | v | 2 v 2 = Ω | v | 2 v 2 ( | v | 2 ) ν
for all t ( 0 , T max ) .
Now, according to (3.17) in [11], we have for any ε > 0 and constant C ε > 0 depending on Ω that
Ω | v | 2 ( | v | 2 ) ν ε Ω | v | 2 2 + C ε for all t ( 0 , T max ) .
By straightforward calculation, we have
Ω | v | 2 2 = 1 16 Ω | v | 2 v 2 = 1 16 Ω ( | v | 2 ) v v | v | 2 v 2 2 = 1 16 Ω | v | 2 2 v 2 2 | v | 2 · v | v | 2 v 3 + | v | 6 v 4 1 16 Ω ( 2 ε 1 + 1 ) | v | 2 2 v 2 + ( 1 2 ε 1 + 1 ) | v | 6 v 4
for ε 1 > 0 . Then, we have from (14) that
Ω | v | 2 2 1 16 2 ε 1 + 1 + ϵ ˜ 8 ε 1 + ϵ ˜ 4 Ω | v | 2 2 v 2
for all t ( 0 , T max ) , where ϵ ˜ = 2 21 4 26 for simplicity. Combining (18) and (17) with (16), we can obtain that
Ω | v | 2 v 2 ( | v | 2 ) ν ε 2 ε 1 + 1 + ϵ ˜ 8 ε 1 + ϵ ˜ 4 Ω | v | 2 2 v 2 + 16 C ε .
Denoting ϵ ^ = ε ( 2 ε 1 + 1 + ϵ ˜ 8 ε 1 + ϵ ˜ 4 ) and C ( ϵ ^ ) = 16 C ε , we prove (15) for any ϵ ^ > 0 . □
In preparation for the construction and estimation of energy-type functionals, some important a p r i o r i estimates are provided and collected into two lemmas in the following.
Lemma 5.
Let k > 0 and ( u , v ) be the solutions of (1) satisfying Proposition 1. Then, we have for any ϵ ^ > 0 that
d d t Ω | v | 4 v 2 + Ω | v | 4 v 2 ( 4 k 3 2 k ϵ ^ ) Ω | v | 2 | D 2 v | 2 v 2 2 Ω | v | 4 u v 3 + 4 Ω | v | 2 v · u v 2 + C ( ϵ ^ ) .
Proof. 
Through straightforward calculation, we can show
d d t Ω | v | 4 v 2 = 4 Ω | v | 2 v · v t v 2 2 Ω | v | 4 v t v 3 = 4 k Ω | v | 2 v · Δ v v 2 I 1 2 Ω | v | 4 v 2 + 4 Ω | v | 2 v u v 2 2 k Ω | v | 4 Δ v v 3 I 2 2 Ω | v | 4 u v 3 .
In light of (10), we have from (15) that
I 1 = 2 k Ω | v | 2 Δ | v | 2 v 2 4 k Ω | v | 2 | D 2 v | 2 v 2 = 2 k Ω | v | 2 v 2 ( | v | 2 ) ν 2 k Ω | v | 2 v 2 · | v | 2 4 k Ω | v | 2 | D 2 v | 2 v 2 = 2 k Ω | v | 2 v 2 ( | v | 2 ) ν 2 k Ω ( | v | 2 ) 2 v 2 + 4 k Ω | v | 2 v · | v | 2 v 3 I 3 4 k Ω | v | 2 | D 2 v | 2 v 2 ( 12 k 2 k ϵ ^ ) Ω | v | 2 | D 2 v | 2 v 2 + 4 k Ω v | v | 2 · | v | 2 v 3 I 3 + C ( ϵ ^ ) .
Similarly, we calculate that
I 2 = 2 k Ω | v | 4 v 3 v = 4 k Ω | v | 2 v · | v | 2 v 3 6 k Ω | v | 6 v 4 .
Given by the sum of I 2 and I 3 and taking ϵ = 1 3 , (6) implies that
I 2 + I 3 = 8 k Ω | v | 2 v · | v | 2 v 3 6 k Ω | v | 6 v 4 8 k ϵ Ω ( | v | 2 ) 2 v 2 + ( 2 k ϵ 6 k ) Ω | v | 6 v 4 = 32 k 3 Ω | v | 2 | D 2 v | 2 v 2 .
Substituting (22) and (21) into (20), we finish the proof by taking the first identity of Lemma 3. □
Lemma 6.
Supposing that ( u , v ) solves (1) and all conditions of Proposition 1 hold, then there exist small ϵ 1 , ϵ 2 > 0 and δ > 0 such that
1 2 d d t Ω u 2 ( 1 χ ϵ 1 δ * ) Ω | u | 2 + χ ϵ 2 4 ϵ 1 Ω u 3 + χ 6 36 ϵ 1 ϵ 2 Ω | v | 6 v 4 .
Proof. 
In light of the u equation of (1) and integration by parts, we can show that
1 2 d d t Ω u 2 = Ω u · ( u χ v v u ) = Ω | u | 2 + Ω χ u u · v v .
The employment of (7) implies
Ω χ u u · v v χ ϵ 1 Ω | u | 2 v 2 3 + χ ϵ 2 4 ϵ 1 Ω u 3 + χ 6 36 ϵ 1 ϵ 2 Ω | v | 6 v 4
for small ϵ 1 , ϵ 2 > 0 . Note that v has the lower bound for any t > 0 . Let δ * : = δ 2 3 be the upper bound of v 2 3 and substitute (25) into (24) to obtain (23). □

5. Uniform Boundedness

In this section, we shall finish the proof of Theorem 1. Firstly, we construct the energy functional and prove that each item of the functional is uniform bounded.
Theorem 3.
For α > 0 , let F α ( u , v ) take the following form:
F α ( u , v ) = α Ω u 2 + Ω | v | 4 v 2 .
Then, for Ω R 2 and any χ > 0 , there exists a constant C k ( u 0 , v 0 , Ω , χ ) > 0 such that if k > C k ( u 0 , v 0 , Ω , χ ) , then for some C > 0
d d t F α ( u , v ) + F α ( u , v ) < C for   all   t ( 0 , T max ) .
Proof. 
Combining (19) and (23), we achieve
d d t ( α Ω u 2 + Ω | v | 4 v 2 ) + ( α Ω u 2 + Ω | v | 4 v 2 ) ( 2 α 2 α χ ϵ 1 δ * ) Ω | u | 2 + α χ ϵ 2 2 ϵ 1 Ω u 3 I 1 + α χ 6 18 ϵ 1 ϵ 2 Ω | v | 6 v 4 I 2 ( 4 k 3 2 k ϵ ^ ) Ω | v | 2 | D 2 v | 2 v 2 2 Ω | v | 4 u v 3 + 4 Ω | v | 2 v · u v 2 I 3 + α Ω u 2 + C ( ϵ ^ ) .
The Gagliardo–Nirenberg inequality and the boundedness of u L 1 ( Ω ) imply that there exists C > 0 depending on u 0 L 1 ( Ω ) , Ω such that
Ω u 2 η Ω | u | 2 + C
for some small η > 0 and
I 1 = α χ ϵ 2 2 ϵ 1 u L 3 ( Ω ) 3 α χ ϵ 2 2 ϵ 1 ( C 1 u L 2 ( Ω ) 2 + C 2 ) ,
where C 1 , C 2 > 0 , depending on u 0 L 1 ( Ω ) and Ω . For I 3 , we employ Lemma 2 to obtain
I 3 1 ϵ 3 Ω | v | 6 v 4 + 4 ϵ 3 Ω | u | 2
for any ϵ 3 > 0 . Combining the first item of (29) with I 2 and employing (14), we have
( α χ 6 18 ϵ 1 ϵ 2 + 1 ϵ 3 ) Ω | v | 6 v 4 ( α χ 6 18 ϵ 1 ϵ 2 + 1 ϵ 3 ) ϵ ˜ Ω | v | 2 | D 2 v | 2 v 2 ,
where we denote ϵ ˜ = 2 21 4 26 for simplicity. Thus, substituting (28)–(30) into (27) gives
d d t ( α Ω u 2 + Ω | v | 4 v 2 ) + ( α Ω u 2 + Ω | v | 4 v 2 ) ( 2 α ( 1 χ ϵ 1 δ * χ ϵ 2 4 ϵ 1 C 1 η ) 4 ϵ 3 ) κ 1 Ω | u | 2 ( 4 k 3 2 k ϵ ^ α χ 6 18 ϵ 1 ϵ 2 ϵ ˜ 1 ϵ 3 ϵ ˜ ) κ 2 Ω | v | 2 | D 2 v | 2 v 2 + C .
Let ϵ 1 , ϵ 2 and ϵ 3 be small, such that κ 1 = 0 . Then, taking a small ϵ ^ such that 2 3 > ϵ ^ , we denote
C ( ϵ 1 , ϵ 2 , ϵ 3 ) : = α χ 6 6 ϵ 1 ϵ 2 ( 4 6 ϵ ^ ) ϵ ˜ + 3 ϵ 3 ( 4 6 ϵ ^ ) ϵ ˜ > 0 ,
and let C k depending on u 0 , v 0 , Ω , χ be the lower bound of C ( ϵ 1 , ϵ 2 , ϵ 3 ) provided κ 1 = 0 . Therefore, for any k C k > 0 , we have κ 2 0 and can then deduce (26). □
Theorem 4.
Let (u,v) be the solutions of (1) satisfying all conditions in Proposition 1. Then,
Ω u 2 ( · , t ) C and Ω | v ( · , t ) | 2 C ,
with t ( 0 , T max ) .
Proof. 
According to (26), there is C > 0 such that
Ω u 2 ( · , t ) C and Ω | v | 4 v 2 ( · , t ) C ,
for all t ( 0 , T max ) . From Young’s inequality and the Gagliardo–Nirenberg inequality, there exist ϵ 4 , ϵ G N > 0 and C > 0 such that
Ω | v ( · , t ) | 2 C ϵ 4 Ω | v | 4 v 2 ( · , t ) + ϵ 4 Ω v 2 ( · , t ) C ϵ 4 Ω | v | 4 v 2 ( · , t ) + ϵ 4 ϵ G N Ω | v ( · , t ) | 2 + C
for all t ( 0 , T max ) . Taking ϵ 4 < 1 2 ϵ G N , then we have ϵ 4 ϵ G N < 1 2 and prove (32). □
Proof of Theorem 1. 
Using the well-known Moser’s technique [32], the L boundedness of u follows from Theorem 4. Indeed, one can follow the estimates of Nagai [15] or directly employ Lemma 2.3 in [7] to prove the theorem. □

6. Conclusions

Our paper proves the uniform boundedness of solutions of the chemotaxis system with singular sensitivity under a small diffusion rate of the chemical signal. We prove a user-friendly inequality that has certain parameters, and construct a new energy functional that is applicable to the double Keller–Segel model with nonlinear sources.

Author Contributions

Conceptualization, J.L., B.T. and D.W.; Methodology, D.W.; Validation, J.L., B.T. and J.T.; Formal analysis, J.L., J.T. and Y.W.; Writing—original draft, B.T. and D.W.; Writing—review & editing, J.L., J.T. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Keller, E.F.; Segel, L.A. Initiation of slime mold aggregation viewed as an instability. J. Theor. Biol. 1970, 26, 399–415. [Google Scholar] [CrossRef] [PubMed]
  2. Biler, P. Global solutions to some parabolic-elliptic systems of chemotaxis. Adv. Math. Sci. Appl. 1999, 9, 347–359. [Google Scholar]
  3. Fujie, K.; Winkler, M.; Yokota, T. Blow-up prevention by logistic sources in a parabolic–elliptic Keller-Segel system with singular sensitivity. Nonlinear Anal. Theory Methods Appl. 2014, 109, 56–71. [Google Scholar] [CrossRef]
  4. Fujie, K.; Winkler, M.; Yokota, T. Boundedness of solutions to parabolic–elliptic Keller-Segel systems with signal-dependent sensitivity. Math. Methods Appl. Sci. 2015, 38, 1212–1224. [Google Scholar] [CrossRef]
  5. Zhigun, A. Generalised supersolutions with mass control for the Keller–Segel system with logarithmic sensitivity. J. Math. Anal. Appl. 2018, 467, 1270–1286. [Google Scholar] [CrossRef] [Green Version]
  6. Fujie, K. Boundedness in a fully parabolic chemotaxis system with singular sensitivity. J. Math. Anal. Appl. 2015, 424, 675–684. [Google Scholar] [CrossRef]
  7. Lankeit, J. A new approach toward boundedness in a two-dimensional parabolic chemotaxis system with singular sensitivity. Math. Methods Appl. Sci. 2015, 39, 394–404. [Google Scholar] [CrossRef] [Green Version]
  8. Marras, M.; Viglialoro, G. Boundedness in a fully parabolic chemotaxis-consumption system with nonlinear diffusion and sensitivity, and logistic source. Math. Nachrichten 2018, 291, 2318–2333. [Google Scholar] [CrossRef]
  9. Zhang, Q. Global bounded solutions to a Keller-Segel system with singular sensitivity. Appl. Math. Lett. 2020, 107, 106397. [Google Scholar] [CrossRef]
  10. Liujie, G.; Fei, G.; Hui, Z. Existence, uniqueness and L∞ -bound for weak solutions of a time fractional Keller-Segel system. Chaos Solitons Fractals 2022, 160, 112–185. [Google Scholar]
  11. Liu, D.; Tao, Y. Global boundedness in a fully parabolic attraction–repulsion chemotaxis model. Math. Methods Appl. Sci. 2015, 38, 2537–2546. [Google Scholar] [CrossRef]
  12. Wang, Q. Global solutions of a Keller-Segel system with saturated logarithmic sensitivity function. Commun. Pure Appl. Anal. 2015, 14, 383–396. [Google Scholar] [CrossRef]
  13. Xu, J.; Liu, Z.; Shi, S. Large time behavior of solutions for the attraction–repulsion Keller-Segel system with large initial data. Appl. Math. Lett. 2019, 87, 13–19. [Google Scholar] [CrossRef]
  14. Osaki, K.; Yagi, A. Finite dimensional attractor for one-dimensional Keller-Segel equations. Funkc. Ekvacioj 2001, 44, 441–469. [Google Scholar]
  15. Nagai, T.; Senba, T.; Yoshida, K. Application of the Trudinger-Moser inequality to a parabolic system of chemotaxis. Funkc. Ekvacioj 1997, 40, 411–434. [Google Scholar]
  16. Corrias, L.; Perthame, B.; Zaag, H. Global solutions of some chemotaxis and angiogenesis systems in high space dimensions. Milan J. Math. 2004, 72, 1–28. [Google Scholar] [CrossRef]
  17. Winkler, M. Absence of collapse in a parabolic chemotaxis system with signal-dependent sensitivity. Math. Nachrichten 2010, 283, 1664–1673. [Google Scholar] [CrossRef]
  18. Fujie, K.; Yokota, T. Boundedness in a fully parabolic chemotaxis system with strongly singular sensitivity. J. Math. Anal. Appl. 2014, 38, 140–143. [Google Scholar] [CrossRef]
  19. Keller, E.F.; Segel, L.A. Traveling bands of chemotactic bacteria: A theoretical analysis. J. Theor. Biol. 1971, 30, 235–248. [Google Scholar] [CrossRef]
  20. Winkler, M. Global solutions in a fully parabolic chemotaxis system with singular sensitivity. Math. Methods Appl. Sci. 2011, 34, 176–190. [Google Scholar] [CrossRef]
  21. Ahn, J.; Kang, K.; Lee, J. Global well-posedness of logarithmic Keller-Segel type systems. J. Differ. Equ. 2021, 287, 185–211. [Google Scholar] [CrossRef]
  22. Stinner, C.; Winkler, M. Global weak solutions in a chemotaxis system with large singular sensitivity. Nonlinear Anal. Real World Appl. 2011, 12, 3727–3740. [Google Scholar] [CrossRef]
  23. Lankeit, J.; Winkler, M. A generalized solution concept for the Keller-Segel system with logarithmic sensitivity: Global solvability for large nonradial data. Nonlinear Differ. Equ. Appl. 2017, 24, 49–73. [Google Scholar] [CrossRef] [Green Version]
  24. Fujie, K.; Senba, T. Global existence and boundedness of radial solutions to a two dimensional fully parabolic chemotaxis system with general sensitivity. Nonlinearity 2016, 29, 2417–2450. [Google Scholar] [CrossRef]
  25. Zhao, X.; Zheng, S. Global boundedness of solutions in a parabolic-parabolic chemotaxis system with singular sensitivity. J. Math. Anal. Appl. 2016, 443, 445–452. [Google Scholar] [CrossRef] [Green Version]
  26. Horstmann, D.; Winkler, M. Boundedness vs. blow-up in a chemotaxis system. J. Differ. Equ. 2005, 215, 52–107. [Google Scholar] [CrossRef] [Green Version]
  27. Winkler, M. The two-dimensional Keller–Segel system with singular sensitivity and signal absorption: Global large-data solutions and their relaxation properties. Math. Model. Methods Appl. Sci. 2016, 26, 987–1024. [Google Scholar] [CrossRef] [Green Version]
  28. Wang, Q.; Wang, D.; Feng, Y. Global well-posedness and uniform boundedness of urban crime models: One-dimensional case. J. Differ. Equ. 2020, 269, 6216–6235. [Google Scholar] [CrossRef]
  29. Winkler, M. Aggregation vs. global diffusive behavior in the higher-dimensional Keller-Segel model. J. Differ. Equ. 2010, 248, 2889–2905. [Google Scholar] [CrossRef] [Green Version]
  30. Young, W.H. On classes of summable functions and their Fourier series. Proc. R. Soc. Lond. Ser. A Contain. Pap. Math. Phys. Character 1912, 87, 225–229. [Google Scholar]
  31. Winkler, M. Global Large-Data Solutions in a Chemotaxis-(Navier–)Stokes System Modeling Cellular Swimming in Fluid Drops. Commun. Partial. Differ. Equ. 2012, 37, 319–351. [Google Scholar] [CrossRef]
  32. Alikakos, N.D. Lp bounds of solutions of reaction-diffusion equations. Commun. Partial. Differ. Equ. 1979, 4, 827–868. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, J.; Tian, B.; Wang, D.; Tang, J.; Wu, Y. Global Boundedness in a Logarithmic Keller–Segel System. Mathematics 2023, 11, 2743. https://doi.org/10.3390/math11122743

AMA Style

Liu J, Tian B, Wang D, Tang J, Wu Y. Global Boundedness in a Logarithmic Keller–Segel System. Mathematics. 2023; 11(12):2743. https://doi.org/10.3390/math11122743

Chicago/Turabian Style

Liu, Jinyang, Boping Tian, Deqi Wang, Jiaxin Tang, and Yujin Wu. 2023. "Global Boundedness in a Logarithmic Keller–Segel System" Mathematics 11, no. 12: 2743. https://doi.org/10.3390/math11122743

APA Style

Liu, J., Tian, B., Wang, D., Tang, J., & Wu, Y. (2023). Global Boundedness in a Logarithmic Keller–Segel System. Mathematics, 11(12), 2743. https://doi.org/10.3390/math11122743

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop