一类趋化消耗模型解的整体有界性
Global Boundedness of Solution to a Chemotaxis-Consumption Model
DOI: 10.12677/PM.2023.1311325, PDF, HTML, XML, 下载: 85  浏览: 129 
作者: 陈 越, 牛 聪:辽宁师范大学数学学院,辽宁 大连
关键词: 趋化消耗次线性敏感Logistic源整体有界性Chemotaxis-Consumption Sublinear Sensitivity Logistic Source Global Boundedness
摘要: 趋化描述了生物有机体受化学信号刺激所产生的偏向性运动,它在生物学、化学、医学等学科领域有着广泛的应用。本文研究齐次Neumann边界条件下的具次线性敏感及带logistic源的趋化消耗模型:,其中μ,χ > 0,r∈ℝ,α∈(0,1)和k > 1。在一维情形下,模型存在整体有界的古典解。
Abstract: Chemotaxis describes the biased movement of biological organisms stimulated by chemical signals. It is widely used in biology, chemistry, medicine and other fields. In this paper, a sublinear sensitive chemotactic-consumption model with logistic source under homogeneous Neumann boundary conditions is studied:, where μ,χ > 0, r∈ℝ, α∈(0,1) and k > 1. In the one-dimensional case, the model has a globally bounded classical solution.
文章引用:陈越, 牛聪. 一类趋化消耗模型解的整体有界性[J]. 理论数学, 2023, 13(11): 3139-3145. https://doi.org/10.12677/PM.2023.1311325

1. 前言

上世纪70年代,Keller和Segel提出了生物趋化模型,它简单而又直观地刻画了生物有机体受化学信号刺激所产生的趋向性运动,在生物学、医学、环境科学等领域有着重要的应用。对于如下带logistic源的趋化消耗模型:

{ u t = Δ u χ ( u v α v ) + r u μ u 2 , x Ω , t > 0 , v t = Δ v u v , x Ω , t > 0 , u ν = v ν = 0 , x Ω , t > 0 , u ( x , 0 ) = u 0 ( x ) , v ( x , 0 ) = v 0 ( x ) , x Ω , (1.1)

其中 Ω n ( n 1 ) 为有界光滑区域, χ > 0 代表趋化强度, α [ 0 , 1 ] ,并且logistic源中 r μ > 0 u , v 分别表示细胞密度和氧气的浓度。另外,初始数据 u 0 C 0 ( Ω ¯ ) v 0 W 1 , ( Ω ) 非负。

现回顾模型(1.1)的相关结论。若 α = 0 ,文献 [1] 证明了当 r > 0 , 0 < v 0 L ( Ω ) 1 6 ( n + 1 ) χ 时,模型存在全局有界的古典解;当 μ 充分大时文献 [2] 建立了模型古典解的整体有界性,并且 r > 0 时收敛到 ( r μ , 0 )

α = 1 ,当 n = 2 , r 0 , μ > μ * ( χ ) > 0 r > 0 且对 μ 再进行限制,文献 [3] 得到了模型(1.1)古典解在时间上的有界性和长时间行为;当 r 0 ,文献 [4] 证明了 n = 1 n 2 , χ < 2 / n μ > ( n 2 ) / 2 n 古典解的整体有界性。另外,将(1.1)中第一个方程变成拟线性形式: u t = ( D ( u ) u ) ( u S ( x , u , v ) v ) ,其中

D ( u ) k D u m 1 , | S ( x , u , v ) | S 0 ( v ) v α ( S 0 是非减的),并且 k D > 0 , m > 1 , α 0 ,文献 [5] 证明了当 m > 3 n 2 2 n , α [ 0 , 1 ) 时,系统存在整体有界的古典解。

本文受以上结论启发,研究一维情形下带logistic源及具次线性敏感的趋化消耗模型:

{ u t = u x x χ ( u v α v x ) x + r u μ u k , x Ω , t > 0 , v t = v x x u v , x Ω , t > 0 , u x = v x = 0 , x Ω , t > 0 , u ( x , 0 ) = u 0 ( x ) , v ( x , 0 ) = v 0 ( x ) , x Ω , (1.2)

其中 μ χ > 0 r α ( 0 , 1 ) k > 1 Ω 为有界区间。

根据以上假设,可以得到系统(1.2)古典解的整体有界性。

定理1.1 设 χ > 0 r α , β ( 0 , 1 ) k > 1 。有

u L ( Ω ) + v W 1 , ( Ω ) + v x v β L ( Ω ) C , t > 0 , (1.3)

其中 C > 0

注 文献 [4] 证明了一维情形下,当 r 0 , α = 1 , k = 2 时,模型古典解的整体有界性。与文献 [4] 不同,若 r ,定理1.1证明了 α ( 0 , 1 ) k > 1 时模型古典解的整体有界性。

2. 预备知识

根据Banach不动点理论,可以得到如下古典解局部存在性,参见 [4] 。

引理2.1 假设 u 0 C 0 ( Ω ¯ ) v 0 W 1 , ( Ω ) 是非负的。若 μ χ > 0 k > 1 r α ( 0 , 1 ) ,则存

T max ( 0 , ] 及唯一非负函数 ( u , v ) ( C 0 ( Ω ¯ × [ 0 , T max ) ) C 2 , 1 ( Ω ¯ × ( 0 , T max ) ) ) 2 在古典解意义下满足系统(1.2)。另外,当 T max < 时, lim t T max u L ( Ω ) =

引理2.2 设 μ χ > 0 k > 1 r α ( 0 , 1 ) ,则

Ω u ( , t ) d x max { Ω u 0 d x , | Ω | ( | r | μ ) 1 k 1 } , t ( 0 , T max ) , (2.1)

(2.2)

证明 根据(1.2)的第一个方程,知

Ω u t d x = r Ω u d x μ Ω u k d x , t ( 0 , T max )

由Hölder不等式有

Ω u t d x r Ω u d x μ | Ω | k 1 ( Ω u d x ) k , t ( 0 , T max ) .

应用伯努利不等式即得(2.1)。由抛物方程的极值原理易得(2.2)。

3. 主要结果

引理3.1 设 μ χ > 0 k > 1 r α ( 0 , 1 ) 。若 p > 1 ,则

d d t Ω u p d x + 2 ( p 1 ) p Ω | ( u p 2 ) x | 2 d x p ( p 1 ) χ 2 2 Ω u p v 2 α v x 2 d x + p r Ω u p d x p μ Ω u k + p 1 d x , t ( 0 , T max ) . (3.1)

证明 根据(1.2)的第一个方程,利用分部积分及Young不等式,知

d d t Ω u p d x = p ( p 1 ) Ω u p 2 u x 2 d x + p ( p 1 ) χ Ω u p 1 v α u x v x d x + p r Ω u p d x p μ Ω u k + p 1 d x 2 ( p 1 ) p Ω | ( u p 2 ) x | 2 d x + p ( p 1 ) χ 2 2 Ω u p v 2 α v x 2 d x + p r Ω u p d x p μ Ω u k + p 1 d x , t ( 0 , T max ) .

将上式进行移项即得(3.1)。

根据文献 [5] 中引理3.3,3.4和3.5,有

引理3.2 设 μ χ > 0 k > 1 r α ( 0 , 1 ) 。若 q 2 ,则存在 c i ( q ) > 0 ( i = 1 , 4 ) ,使得

d d t Ω v q + 1 v x q d x + c 1 Ω v q + 1 v x q 2 | v x x | 2 d x + c 2 Ω v q 1 v x q + 2 d x c 3 Ω v d x + c 4 Ω u v q + 2 v x q 2 | v x x | d x , t ( 0 , T max ) . (3.2)

引理3.3 设 μ χ > 0 k > 1 r α ( 0 , 1 ) 。定义

F = Ω u p d x + Ω v q + 1 v x q d x .

p > 1 q 2 ,则存在 Γ = Γ ( p , q , χ ) > 0 ,使得

F + F + 2 ( p 1 ) p Ω | ( u p 2 ) x | 2 d x Γ Ω u p ( q + 2 ) q v 2 ( q + 1 ) 2 ( q + 2 ) α q d x + Γ Ω u q + 2 2 v d x + ( c 3 + ( 3 c 2 ) q 2 ) Ω v d x + c 5 , t ( 0 , T max ) . (3.3)

证明 将(3.1)与(3.2)相加,有

F + F + 2 ( p 1 ) p Ω | ( u p 2 ) x | 2 d x + c 1 Ω v q + 1 v x q 2 | v x x | 2 d x + c 2 Ω v q 1 v x q + 2 d x p ( p 1 ) χ 2 2 Ω u p v 2 α v x 2 d x + ( p r + 1 ) Ω u p d x p μ Ω u k + p 1 d x + c 3 Ω v d x + c 4 Ω u v q + 2 v x q 2 | v x x | d x + Ω v q + 1 v x q d x , t ( 0 , T max ) (3.4)

由Young不等式,存在 c 5 > 0 ,有

p ( p 1 ) χ 2 2 Ω u p v 2 α v x 2 d x ( 3 c 2 ) 2 q ( p ( p 1 ) χ 2 2 ) q + 2 q Ω u p ( q + 2 ) q v 2 ( q + 1 ) 2 ( q + 2 ) α q d x + c 2 3 Ω v q 1 v x q + 2 d x , (3.5)

c 4 Ω u v q + 2 v x q 2 | v x x | d x c 1 Ω v q + 1 v x q 2 | v x x | 2 d x + c 4 2 4 c 1 Ω u 2 v q + 3 v x q 2 d x c 1 Ω v q + 1 v x q 2 | v x x | 2 d x + c 2 3 Ω v q 1 v x q + 2 d x + ( 3 c 2 ) q 2 4 ( c 4 2 4 c 1 ) q + 2 4 Ω u q + 2 2 v d x , (3.6)

Ω v q + 1 v x q d x c 2 3 Ω v q 1 v x q + 2 d x + ( 3 c 2 ) q 2 Ω v d x , (3.7)

( p r + 1 ) Ω u p d x p μ Ω u k + p 1 d x + c 5 . (3.8)

Γ = max { ( 3 c 2 ) 2 q ( p ( p 1 ) χ 2 2 ) q + 2 q , ( 3 c 2 ) q 2 4 ( c 4 2 4 c 1 ) q + 2 4 } ,将(3.5),(3.6),(3.7),(3.8)代入到(3.4)中可得(3.3)。

引理3.4 设 k > 1 μ , χ > 0 r α ( 0 , 1 ) 。对于 p , q > 1 ,存在 C 0 > 0 ,使得

Ω u p d x + Ω v q + 1 v x q d x C 0 , t ( 0 , T max ) . (3.9)

证明 取 p > max { 1 , 4 α 3 2 2 α } max { 2 α 1 1 α , p } < q < 2 + 2 p 。根据引理2.2与Gagliardo-Nirenberg不等

式及Young不等式,有

Γ Ω u p ( q + 2 ) q v 2 ( q + 1 ) 2 ( q + 2 ) α q d x Γ v 0 L ( Ω ) 2 ( q + 1 ) 2 ( q + 2 ) α q Ω u p ( q + 2 ) q d x = Γ v 0 L ( Ω ) 2 ( q + 1 ) 2 ( q + 2 ) α q u p 2 L 2 ( q + 2 ) q ( Ω ) 2 ( q + 2 ) q c 6 ( u p 2 ) x L 2 ( Ω ) 2 ( q + 2 ) q a u p 2 L 2 p ( Ω ) 2 ( q + 2 ) q ( 1 a ) + c 6 u p 2 L 2 p ( Ω ) 2 ( q + 2 ) q c 7 ( u p 2 ) x L 2 ( Ω ) 2 p q + 4 p 2 q q + p q + c 7 p 1 p ( u p 2 ) x L 2 ( Ω ) 2 + c 8 , (3.10)

Γ Ω u q + 2 2 v d x Γ v 0 L ( Ω ) Ω u q + 2 2 d x = Γ v 0 L ( Ω ) u p 2 L q + 2 p ( Ω ) q + 2 p c 9 ( u p 2 ) x L 2 ( Ω ) q + 2 p b u p 2 L 2 p ( Ω ) q + 2 p ( 1 b ) + c 9 u p 2 L 2 p ( Ω ) q + 2 p c 10 ( u p 2 ) x L 2 ( Ω ) p q p + p 2 + c 10 p 1 p ( u p 2 ) x L 2 ( Ω ) 2 + c 11 , (3.11)

其中 a = p q q + 2 1 + p ( 0 , 1 ) , b = p 2 p q + 2 1 + p ( 0 , 1 ) , c j > 0 ( j = 6 , , 11 )

将(3.10),(3.11)代入(3.3),并应用引理2.2,存在 c 12 > 0 ,有

F + F c 12 , t ( 0 , T max ) .

应用常数变易法得到

F max { c 12 , Ω u 0 p d x + Ω v 0 q + 1 v 0 x q d x } , t ( 0 , T max ) .

定理1.1 证明 由热半群表示,知

u ( , t ) = e t 2 x 2 u 0 χ o t e ( t s ) 2 x 2 ( u v α v x ) x d s + r o t e ( t s ) 2 x 2 u d s μ o t e ( t s ) 2 x 2 u k d s , (3.12)

t ( 0 , T max ) 。利用文献 [6] 中Neumann热半群的 L P L q 估计,对 p > 1

u L ( Ω ) e t 2 x 2 u 0 L ( Ω ) + χ o t e ( t s ) 2 x 2 ( u v α v x ) x L ( Ω ) d s + | r | o t e ( t s ) 2 x 2 u L ( Ω ) d s u 0 L ( Ω ) + c 13 o t ( 1 + ( t s ) 3 4 ) e λ ( t s ) u v α v x L 2 ( Ω ) d s + c 13 o t ( 1 + ( t s ) 1 2 p ) e λ ( t s ) u L p ( Ω ) d s , t ( 0 , T max ) , (3.13)

其中 c 13 > 0 。由引理3.4,有

Ω u 4 + 1 2 ( 1 α ) d x + Ω v ( 8 + 1 ( 1 α ) ) + 1 | v x | 8 + 1 ( 1 α ) d x c 14 ,

其中 c 14 > 0 。进一步,利用Young不等式及(2.2),存在 c 15 > 0 ,使得

Ω ( u v α v x ) 2 d x Ω u 4 d x + Ω | v α v x | 4 d x Ω u 4 + 1 2 ( 1 α ) d x + Ω | v α v x | 8 + 1 ( 1 α ) d x + 2 | Ω | Ω u 4 + 1 2 ( 1 α ) d x + 2 | Ω | + v 0 L ( Ω ) 8 ( 1 α ) Ω v ( 8 + 1 ( 1 α ) ) + 1 | v x | 8 + 1 ( 1 α ) d x c 15 , t ( 0 , T max ) . (3.14)

类似地,根据半群理论可以证得 v W 1 , ( Ω ) C 。进而由引理2.1知 T max =

w = v 1 β β ( 0 , 1 )

w t = ( 1 β ) v β v t = w x x + β 1 β w x 2 w ( 1 β ) w u , t > 0.

通过热半群表示,知

w ( , t ) = e t 2 x 2 w 0 + β 1 β o t e ( t s ) 2 x 2 w x 2 w d s ( 1 β ) o t e ( t s ) 2 x 2 w u d s , t > 0. (3.15)

应用Neumann热半群估计,对 q > max { 1 1 β , 2 + β 2 } ,有

( 1 β ) v x v β L ( Ω ) = w x L ( Ω ) ( e t 2 x 2 w 0 ) x L ( Ω ) + β 1 β o t ( e ( t s ) 2 x 2 w x 2 w ) x L ( Ω ) d s + ( 1 β ) o t ( e ( t s ) 2 x 2 w u ) x L ( Ω ) d s = ( e t 2 x 2 v 0 1 β ) x L ( Ω ) d s + β ( 1 β ) o t ( e ( t s ) 2 x 2 v x 2 v 1 + β ) x L ( Ω ) d s + ( 1 β ) o t ( e ( t s ) 2 x 2 u v 1 β ) x L ( Ω ) d s ( v 0 1 β ) x L ( Ω ) + c 16 o t ( 1 + ( t s ) 1 2 1 + β 2 ( 2 q 1 ) ) e λ ( t s ) v x 2 v 1 + β L 2 q 1 1 + β ( Ω ) d s + c 16 o t ( 1 + ( t s ) 1 2 1 2 q ) e λ ( t s ) u v 1 β L q ( Ω ) d s , t > 0 , (3.16)

其中 c 16 > 0 。对(3.16)应用引理3.3及(2.2),可得 v x v β L ( Ω ) C , t > 0

参考文献

[1] Wang, L., Khan, U.D. and Khan, U.D. (2013) Boundedness in a Chemotaxis System with Consumption of Chemoattractant and Logistic Source. Electronic Journal of Differential Equations, 2013, 172-203.
[2] Lankeit, J. and Wang, Y. (2017) Global Existence, Boundedness and Stabilization in a High-Dimensional Chemotaxis System with Consumption. Discrete & Continuous Dynamical Systems, 37, 6099-6121.
https://doi.org/10.3934/dcds.2017262
[3] Wang, W. (2019) The Logistic Chemotaxis System with Singular Sensitivity and Signal Absorption in Dimension Two. Nonlinear Analysis: Real World Applications, 50, 532-561.
https://doi.org/10.1016/j.nonrwa.2019.06.001
[4] Lankeit, E. and Lankeit, J. (2019) Classical Solutions to a Lo-gistic Chemotaxis Model with Singular Sensitivity and Signal Absorption. Nonlinear Analysis: Real World Applications, 46, 421-445.
https://doi.org/10.1016/j.nonrwa.2018.09.012
[5] Winkler, M. (2022) Approaching Logarithmic Singularities in Quasilinear Chemotaxis-Consumption Systems with Signal-Dependent Sensitivities. Discrete and Continuous Dynam-ical Systems-B, 27, 6565-6587.
https://doi.org/10.3934/dcdsb.2022009
[6] Winkler, M. (2010) Aggregation vs. Global Diffusive Behavior in the Higher-Dimensional Keller-Segel Model. Journal of Differential Equations, 248, 2889-2905.
https://doi.org/10.1016/j.jde.2010.02.008