一类二维Keller-Segel系统整体解的存在性
The Global Existence of Solutions to a Class of Two-Dimensional Keller-Segel Systems
摘要: 该文研究一类Keller-Segel方程组的Cauchy问题: { u t =div( uχuv ), x 2 , t( 0, ), κ v t =Δvλv+u, x 2 , t( 0, ), u( x,0 )= u 0 ( x )0,v( x,0 )= v 0 ( x )0, x 2 , 其中 κ0 χ>1 λ>0 ,且 u 0 L 1 ( 2 ) L ( 2 ),  v 0 L 1 ( 2 ) H 1 ( 2 ) ,通过构造更一般的能量泛函和能量不等式,证明得到若细胞或细菌的初始质量满足 2 u 0 dx<min{ 8π χ , 1 2( κ ( χ 1 ) 2 4 ε * + ( χ 1 ) 2 )max{ C GN 2 ,9 } }, 其中 ε * ( 0, 1 2 ( χ 1 ) 2 ) C GN 为Gagliardo-Nirenberg常数,则该模型的解整体存在。
Abstract: This paper investigates a class of Keller-Segel systems described by the Cauchy problem: { u t =div( uχuv ), x 2 , t( 0, ), κ v t =Δvλv+u, x 2 , t( 0, ), u( x,0 )= u 0 ( x )0,v( x,0 )= v 0 ( x )0, x 2 , where κ0 , χ>1 , and λ>0 . By constructing a more general energy functional together with a corresponding energy inequality, we prove global-in-time existence provided the initial cell (bacterial) mass satisfies 2 u 0 dx<min{ 8π χ , 1 2( κ ( χ 1 ) 2 4 ε * + ( χ 1 ) 2 )max{ C GN 2 ,9 } }, for some ε * ( 0, 1 2 ( χ 1 ) 2 ) . Here C GN denotes the Gagliardo-Nirenberg constant. Under this smallness condition on the initial mass, the model admits a global solution. Then the global existence of classical solutions to this system can be established.
文章引用:李鉴峰, 韩永杰. 一类二维Keller-Segel系统整体解的存在性[J]. 应用数学进展, 2025, 14(12): 385-398. https://doi.org/10.12677/aam.2025.1412515

1. 引言

趋化性是生物体或细胞对化学信号刺激产生的定向移动性质,是一种广泛存在于自然界的适应性机制。从微观的细胞运动到宏观的生态系统调节,趋化性在生物体的生存、繁殖和环境适应中起到了关键作用。例如,细胞或细菌通过感知化学物质的浓度梯度,能够向目标方向移动以获取营养、修复损伤或进行免疫防御。同时,趋化性异常与癌症转移等重大疾病的发生密切相关,其研究对于疾病的防治具有重要意义[1]-[12]

为了揭示趋化行为背后的动力学机制,研究者们于1970年提出了描述趋化现象的Keller-Segel系统 [1] [2]

{ u t =div( uχuv ), xΩ, t( 0, ), κ v t =Δvλv+u, xΩ, t( 0, ), (1.1)

其中 u( x,t ) 表示细胞或细菌密度, v( x,t ) 表示化学物质浓度,信号扩散系数 κ0 ,趋化系数 χ>0 ,信号衰减系数 λ>0 。当 κ>0 时,该系统为抛物–抛物型方程组。同时考虑到当化学物质的扩散速度要比细胞扩散速度快得多时,该Keller-Segel系统可简化为 κ=0 时的抛物——椭圆型Keller-Segel系统。

系统(1.1)的初边值问题有很多丰富的数学研究结果,特别是当空间维数是二维时,相应的初边值问题存在临界质量现象,即细胞的初始质量 Ω u ( x,0 )dx 的大小是区分相应初边值问题的解是否整体存在或发生爆破特性的指标。具体而言,当 Ω 2 是具有光滑边界的有界区域, κ=χ=λ=1 时,若初始质量 Ω u ( x,0 )dx<8π ,系统(1.1)的齐次Neumann初边值问题存在径向对称的整体经典解[13];当初始质量 Ω u ( x,0 )dx>8π 时,相应初边值问题的解会在有限时刻爆破[14]

另一方面,当 κ=0 χ=λ=1 时,系统(1.1)的齐次Neumann初边值问题的相应结果与 κ=1 时的情况类似。即,如果 Ω u ( x,0 )dx<8π ,则该初边值问题的径向对称解整体存在[3] [9] [15];如果 Ω u ( x,0 )dx>8π ,则可以构造初值使得该初边值问题的解在有限时刻发生爆破[16]-[18]。同时系统(1.1)在二维全空间上的Cauchy问题也有如下相应的结果。Nagai与Ogawa证明了当 κ=0 χ=1 λ>0 时,若 2 u ( x,0 )dx<8π ,则该初值问题的解整体存在;当 κ>0 χ=1 λ>0 时,若 2 u ( x,0 )dx<4π ,该初值问题的解整体存在[19]。随后,Mizoguchi于2013年证明了当 κ=1 χ=1 λ>0 时,若 2 u ( x,0 )dx<8π 则相应的初值问题存在整体解[20]

因此一个自然的问题是:对任意的 κ0 χ>1 λ>0 ,系统(1.1)在二维全空间上的Cauchy问题是否存在整体解。鉴于此,本文在 κ0 χ>1 λ>0 的情况下,通过建立新的积分不等式,证明了系统(1.1)在二维全空间上的Cauchy问题解的整体存在性。本文的研究不仅进一步统一和推广了现有理论结果,还为二维全空间下更一般的Keller-Segel模型提供了新的研究工具和理论支撑。

2. 主要结果

本文研究如下一类Keller-Segel系统解的整体存在性。

{ u t =div( uχuv ), x 2 ,t( 0, ), κ v t =Δvλv+u, x 2 ,t( 0, ), u( x,0 )= u 0 ( x )0,v( x,0 )= v 0 ( x )0, x 2 (2.1)

其中 χ κ λ 为常数,满足 κ0 χ>1 λ>0

本文的主要结论如下:

定理2.1 κ0 χ>1 λ>0 ( u 0 , v 0 )( L 1 ( 2 ) L ( 2 ) )×( L 1 ( 2 ) H 1 ( 2 ) ) ,若初始质量条件满足

M:= 2 u 0 dx<min{ 8π χ , 1 2( κ ( χ 1 ) 2 4 ε * + ( χ 1 ) 2 )max{ C GN 2 ,9 } },

其中 ε ( 0, 1 2 ( χ 1 ) 2 ), C GN 为Gagliardo-Nirenberg常数,则系统(2.1)的解整体存在。

如引言所述,Nagai、Ogawa及Mizoguchi分别在不同条件下探讨了二维全空间上的Keller-Segel模型解的整体存在性,但Nagai、Ogawa的证明对 χ>1 的情况不适用;Mizoguchi的证明只考虑了 κ=χ=1 的情况。为了统一并推广现有理论结果,本文将引入一个新的积分不等式,用于建立 u L 2 ( 2 ) 2 的一个有效的上界估计,实现了参数 κ0 χ>1 的任意性和对应 κ0 阈值的一致性。

3. 预备知识

在研究Keller-Segel系统解的整体存在性之前,我们首先需要建立系统(2.1)解的局部存在性。下面的引理表明在适当的初始条件下,系统(2.1)存在唯一的非负局部光滑解。由于Yagi在其研究中证明了该结论(见Nagai等人的引述[13]),我们可以直接得到下面引理。

引理3.1 κ0 χ>0 λ>0 ( u 0 , v 0 )( L 1 ( 2 ) L ( 2 ) )×( L 1 ( 2 ) H 1 ( 2 ) ) ,则存在 T( 0, ] ,使得初值问题(2.1)在 2 ×( 0,T ) 上存在唯一的非负局部光滑解 ( u,v ) ,满足

( u,v )C( [ 0,T ): L 1 ( 2 ) H 1 ( 2 ) ) .

下面引理可参见文献[13]的引理3.2。

引理3.2 κ0 χ>0 λ>0 ( u 0 , v 0 )( L 1 ( 2 ) L ( 2 ) )×( L 1 ( 2 ) H 1 ( 2 ) ) ,则存在 T( 0, ] ,使得对任意的 t( 0,T ) ,系统(2.1)的局部光滑解 ( u,v ) 满足

(i) u( t ) L 1 ( 2 ) = u 0 L 1 ( 2 )

(ii) { v( t ) L 1 ( 2 ) = 1 λ u 0 L 1 ( 2 ) ,                                         κ=0, v( t ) L 1 ( 2 ) = e λ κ t v 0 L 1 ( 2 ) + 1 λ ( 1 e λ κ t ) u 0 L 1 ( 2 ) ,  κ>0.

为了进一步探讨系统(2.1)解的整体存在性,我们需要构造细致的能量泛函与更一般化的广义能量等式。

引理3.3 κ0 χ>0 λ>0 存在 T( 0, ] ,使得对任意的 t( 0,T ) ,系统(2.1)的局部光滑解 ( u,v ) 满足

d dt W( t )+ 2 | ( log( 1+u( t ) )+ χ2 χ 2 v( t ) ) | 2 dx +κ t v( t ) L 2 ( 2 ) 2 + 2 u ( t ) | ( log( 1+u( t ) ) χ v( t ) ) | 2 dx = ( χ 1 ) 2 2 u( t )v( t )dx + ( 2 χ χ ) 2 4 2 | v | 2 dx ,

其中

W( t )= 2 ( 1+u( t ) )( log( 1+u( t ) ) )dx + 1 2 v( t ) L 2 ( 2 ) 2 + λ 2 v( t ) L 2 ( 2 ) 2 2 u ( t )v( t )dx. (3.1)

证明 对于系统(2.1)的局部光滑解 ( u,v ) ,注意到对任意的 t( 0,T )

d dt ( 1+u( t ) )log( 1+u( t ) )= t u( t )log( 1+u( t ) )+ t u( t )

从而对系统(2.1)的第一个方程乘以 log( 1+u( t ) ) 并进行分部积分可以得到

d dt 2 ( 1+u( t ) )log( 1+u( t ) )dx d dt 2 u ( t )dx+ 2 1 1+u( t ) | u( t ) | 2 dx χ 2 u ( t )log( 1+u( t ) )v( t )dx=0. (3.2)

对系统(2.1)的第二个方程乘以 t v( t ) ,得到

d dt ( 1 2 v( t ) L 2 ( 2 ) 2 + λ 2 v( t ) L 2 ( 2 ) 2 )+κ t v( t ) L 2 ( 2 ) 2 2 u ( t ) t v( t )dx=0, (3.3)

同时对系统(2.1)的第一个方程乘以 v 并进行分部积分

2 t u( t )v( t )dx + 2 u( t )v( t )dxχ 2 u ( t ) | v( t ) | 2 dx=0, (3.4)

结合(3.3),(3.4)式,对任意的 t( 0,T ) ,有

d dt ( 1 2 v( t ) L 2 ( 2 ) 2 + λ 2 v( t ) L 2 ( 2 ) 2 2 u ( t )v( t )dx )+κ t v( t ) L 2 ( 2 ) 2 = 2 u( t )v( t )dxχ 2 u ( t ) | v( t ) | 2 dx. (3.5)

结合(3.2),(3.5)式和 d dt 2 u ( t )dx=0 ,可知对任意的 t( 0,T )

d dt { 2 ( 1+u( t ) ) log( 1+u( t ) )dx+ 1 2 v( t ) L 2 ( 2 ) 2 + λ 2 v( t ) L 2 ( 2 ) 2 2 u ( t )v( t )dx } +κ t v( t ) L 2 ( 2 ) 2 + 2 1 1+u( t ) | u( t ) | 2 dxχ 2 u ( t )log( 1+u( t ) )v( t )dx = 2 u( t )v( t )dxχ 2 u ( t ) | v( t ) | 2 dx. (3.6)

直接计算得到如下三个等式,

2 u ( t ) | ( log( 1+u( t ) ) χ v( t ) ) | 2 dx = 2 u ( log( 1+u ) ) 2 dx 2 χ 2 uvlog( 1+u )dx+χ 2 u | v | 2 dx ,

2 | ( log( 1+u( t ) ) 2 χ χ 2 v( t ) ) | 2 dx = 2 | log( 1+u( t ) ) | 2 dx ( 2 χ χ ) 2 log( 1+u( t ) )v( t )dx + ( 2 χ χ ) 2 4 2 | v( t ) | 2 dx ,

2 1 1+u( t ) | u( t ) | 2 dx= 2 u( t ) | log( 1+u( t ) ) | 2 dx + 2 | log( 1+u( t ) ) | 2 dx ,

从而对任意的 t( 0,T ) ,有

d dt { 2 ( 1+u( t )log( 1+u( t ) ) )dx + 1 2 v( t ) L 2 ( 2 ) 2 + λ 2 v( t ) L 2 ( 2 ) 2 2 u ( t )v( t )dx } +κ t v( t ) L 2 ( 2 ) 2 + 2 u ( t ) | ( log( 1+u( t ) ) χ v( t ) ) | 2 dx + 2 | ( log( 1+u( t ) ) 2 χ χ 2 v( t ) ) | 2 dx +( 2 χ χ1 ) 2 u( t )v( t )dx = ( 2 χ χ ) 2 4 2 | v | 2 dx.

在研究系统(2.1)的过程中,利用Trudinger-Moser不等式来控制(2.1)中非线性项的增长是关键步骤。为了证得本文的结果,我们需要建立下述Trudinger-Moser不等式的变体。

引理3.4 G 为二维平面 2 的一个有界区域, | G | 表示 G 的Lebesgue测度,则存在仅依赖于 G 的常数 C G | G |>0 ,使得对任意的非负函数 f h H 0 1 ( G ) ,有

G f hdx G f log( f )dx+M{ log( C G )+ 1 16π h L 2 (G) 2 }MlogM,

其中 M= G f dx.

证明[21]的引理5.4,下列积分不等式成立

G f hdx G f log( f )dx+Mlog( G e h dx )MlogM.

其中 M= G f dx 。再回顾由[22][23]得出的不等式并参考[13] [20] [24]的不等式变体,可知存在仅依赖于 G 的常数 C G | G |>0 ,使得对任意的 h H 0 1 ( G )

log G e |h| dxlog( C G )+ 1 16π h L 2 (G) 2 .

结合这两个不等式可以得到

G f hdx G f log( f )dx+M{ log( C G )+ 1 16π h L 2 (G) 2 }MlogM.

下面引理出自文献[19]的引理2.1。

引理3.5对任意 ε>0 ,存在常数 C( ε )>0 ,使得对 f L 3 ( 2 ) ,满足

f L 3 ( 2 ) ε ( 1+f )log( 1+f ) L 1 ( 2 ) 1/3 f L 2 ( 2 ) 2/3 +C( ε ) f L 1 ( 2 ) 1/3 ,

其中 C( ε )=O( ( e ε 3 1 ) 2/3 ) ε0

在后续建立引理4.2过程中,需要利用下面引理中的积分不等式来建立 u L 2 ( 2 ) 2 的一个上界估计,从而将估计式(4.9)右端 u L 2 ( 2 ) 2 有效地转化为能量泛函在初始时刻的值 W( 0 ) 和对数增长项的组合形式。

引理3.6 对任意非负函数 u L 2 ( 2 ) ,定义 M:= u L 1 ( 2 ) <

I( u ):= 2 1 1+u | u | 2 dx

则存在常数 C=max{ C GN 2 ,9 } ,有

u L 2 ( 2 ) 2 CMI( u )+CM

证明 定义 E 0 ={ x:u( x )<1 },  E 1 ={ x:u( x )1 } ,易知在 E 0 u 2 u ,故有

E 0 u 2 dx E 0 udx M .

ϕ= 1+u 10 ,显然有

ϕ 2 = ( 1+u 1 ) 2 =u+22 1+u u,  ϕ L 2 ( 2 ) 2 M.

又由 ϕ= u 2 1+u ϕ L 2 ( 2 ) 2 = 1 4 2 1 1+u | u | 2 dx= 1 4 I( u )

再由Gagliard-Nirenberg不等式可知,存在常数 C 1 >0 ,有

ϕ L 4 ( 2 ) 4 C 1 ϕ L 2 ( 2 ) 2 ϕ L 2 ( 2 ) 2 ,

从而有

ϕ L 4 ( 2 ) 4 C GN 4 MI( u ).

对任意的 t0, φ= 1+t 1 ,则 t= φ 2 +2φ ,从而

t 2 =( φ 2 +2φ )= φ 4 +4 φ 3 +4 φ 2     φ 4 + φ 4 +4 φ 2 +4 φ 2    2 φ 4 +8t.

t=u( x ) 并在 E 1 上积分,有

E 1 u 2 dx 2 E 1 ϕ 4 dx +8 E 1 udx            2 ϕ L 4 ( 2 ) 4 +8M.

综上所述,可知存在常数 C=max{ C GN 2 ,9 } ,有

u L 2 ( 2 ) 2 = E 0 u 2 dx + E 1 u 2 dx             M+ C 1 2 MI( u )+8M             CMI( u )+CM.

4. 解的整体存在性

本节证明系统(2.1)解的整体存在性。为此,我们首先利用引理3.4建立能量泛函 W( t ) 的一个下界估计。

引理4.1 κ0 χ>0 λ>0 ( u 0 , v 0 )( L 1 ( 2 ) L ( 2 ) )×( L 1 ( 2 ) H 1 ( 2 ) ) ,假设 ( u,v ) 是模型(2.1)在 2 ×( 0,T ) 上的解,其中 0<T<+ ,且满足 M< 8π χ ,则存在依赖于 u 0 L 1 ( 2 ) v 0 L 1 ( 2 ) 的常数 Q>0 和依赖于 u 0 L 1 ( 2 ) 的常数 δ>0 ,使得对任意的 t( 0,T ) ,有

W( t )δ 2 ( u( t )+1 ) log( u( t )+1 )dx+ λ 2 v( t ) L 2 ( 2 ) 2 Q. (4.1)

证明 根据引理3.2可知,对任意的 t( 0,T )

v( t ) L 1 ( 2 ) v 0 L 1 ( 2 ) + 1 λ u 0 L 1 ( 2 ) , (4.2)

选取足够大的常数 s>0 ,定义

Ω( t ):={ x 2 :v( x,t )>s },  v ˜ ( x,t ):=max{ v( x,t )s,0 }.

结合(4.2)式知,对任意的 t( 0,T ) ,可以估计

| Ω( t ) |< 1 s ( v 0 L 1 ( 2 ) + 1 λ u 0 L 1 ( 2 ) ). (4.3)

结合 u( t ) L 1 ( 2 ) = u 0 L 1 ( 2 ) ,对任意的 t( 0,T ) ,有

2 u( t )v( t )dx = Ω(t) u( t )v( t )dx + 2 \Ω(t) u( t )v( t )dx Ω(t) u( t )( v ˜ ( t )+s )dx +s 2 \Ω(t) u( t )dx Ω(t) u( t ) v ˜ ( t )dx +s u 0 L 1 ( 2 ) .  (4.4)

对任意的 t( 0,T ) ,若 Ω( t )= ,则引理4.1得证;若 Ω( t ) ,则将 Ω( t ) 划分为至多可列个互不相交的开区域的并,即 Ω( t )= i=1 m(t) Ω i ( t ) ,其中 m( t ){ + }

首先由(4.4)式可知,对任意的 t( 0,T )

2 u( t )v( t )dx s u 0 L 1 ( 2 ) Ω(t) ( u( t )+1 ) v ˜ ( t )dx . 

根据引理3.4得到,取足够小的 δ>0 ,对任意的 t( 0,T )

Ω i (t) ( u( t )+1 ) v ˜ ( t )dx = Ω i (t) ( 1δ )( u( t )+1 ) v ˜ ( t ) ( 1δ ) dx Ω i (t) ( 1δ )( u( t )+1 ) log( ( 1δ )( u( t )+1 ) )dx + Ω i (t) ( 1δ )( u( t )+1 )dx ( log C Ω i (t) + 1 16π v ˜ ( t ) 1δ 2 ) Ω i (t) ( 1δ )( u( t )+1 )dx log( Ω i (t) ( 1δ )( u( t )+1 )dx ) ( 1δ ) Ω i (t) ( u( t )+1 ) log( u( t )+1 )dx + 1 16π( 1δ ) ( M i ( t )+| Ω i ( t ) | ) Ω i (t) | v( t ) | 2 dx +( 1δ )log( C Ω i (t) )( M i ( t )+| Ω i ( t ) | ),

其中 M i ( t )= Ω i (t) u ( t )dx i=1,2,,m( t )

根据积分的可列可加性得到,对任意的 t( 0,T )

Ω(t) ( u( t )+1 ) v ˜ ( t )dx = i=1 m(t) Ω i (t) ( u( t )+1 )v( t )dx ( 1δ ) Ω(t) ( u( t )+1 ) log( u( t )+1 )dx + 1 16π( 1δ ) ( M( t )+| Ω( t ) | ) Ω(t) | v( t ) | 2 dx +( 1δ )log( C Ω(t) )( M( t )+| Ω( t ) | ),

其中 M( t )= Ω(t) u ( t )dx C Ω(t) :=max{ C Ω 1 (t) ,, C Ω m(t) (t) }

结合 M( t ) u 0 L 1 ( 2 ) 和(4.3)式,可以得到,对任意的 t( 0,T )

2 u( t )v( t )dx ( 1δ ) Ω(t) ( u( t )+1 ) log( u( t )+1 )dx + 1 16π( 1δ ) { u 0 L 1 ( 2 ) + 1 s ( v 0 L 1 ( 2 ) + 1 λ u 0 L 1 ( 2 ) ) } Ω(t) | v( t ) | 2 dx +Q, (4.5)

其中常数 Q=log( C Ω(t) ){ u 0 L 1 ( 2 ) + 1 s ( v 0 L 1 ( 2 ) + 1 λ u 0 L 1 ( 2 ) ) }+s u 0 L 1 ( 2 )

进一步选取适当的 s>0 δ>0 ,对任意的 t( 0,T ) ,成立

W( t )= 2 ( 1+u( t ) )log( 1+u( t ) )dx + 1 2 v( t ) L 2 ( 2 ) 2 + λ 2 v( t ) L 2 ( 2 ) 2 2 u( t )v( t )dx δ 2 ( u( t )+1 ) log( u( t )+1 )dx+ λ 2 v( t ) L 2 ( 2 ) 2 Q.

注意到对初始质量 m( u 0 ; 2 ) 的假设,我们可以利用引理4.1建立的能量泛函 W( t ) 的下界估计以及引理3.3得到如下结论。

引理4.2 κ0 χ>1 λ>0 ( u 0 , v 0 )( L 1 ( 2 ) L ( 2 ) )×( L 1 ( 2 ) H 1 ( 2 ) ) ,设初始质量条件满足

M:= 2 u 0 dx<min{ 8π χ , 1 2( κ ( χ 1 ) 2 4 ε * + ( χ 1 ) 2 )max{ C GN 2 ,9 } },

其中 ε ( 0, 1 2 ( χ 1 ) 2 ), C GN 为Gagliardo-Nirenberg常数。假设 ( u,v ) 为模型(2.1)在 2 ×( 0,T ) 上的解,其中 0<T<+ 。则存在依赖于 M v 0 L 1 ( 2 ) T 的常数 R H>0 和依赖于 M 的常数 δ>0 ,使得对任意的 t( 0,T )

δ 2 ( u( t )+1 ) log( u( t )+1 )dx+ κ 2 0 t t v( τ ) L 2 ( 2 ) 2 dτ R,  v( t ) 2 H.

证明 类似(4.5)式的证明,对任意的 t( 0,T ) ,存在常数 Q ,有

1 2 2 | v( t ) | 2 dx W( t )+Q+ χ 16π { u 0 L 1 ( 2 ) + 1 s ( v 0 L 1 ( 2 ) + 1 λ u 0 L 1 ( 2 ) ) } 2 | v( t ) | 2 dx . (4.6)

因为 M< 8π χ ,则

k:= 1 2 χ 16π { u 0 L 1 ( 2 ) + 1 s ( v 0 L 1 ( 2 ) + 1 λ u 0 L 1 ( 2 ) ) }>0,

从而对任意的 t( 0,T ) ,有

v( t ) 2 2 1 k { W( t )+Q } . (4.7)

结合引理3.3和(4.7)式,可以得到对任意的 t( 0,T )

d dt W( t )+κ t v( t ) L 2 ( 2 ) 2 + 2 1 1+u( t ) | u( t ) | 2 dx = 2 u( t )v( t )dxχ 2 u ( t ) | v( t ) | 2 dx+χ 2 u ( t )log( 1+u( t ) )v( t )dx ( 2 χ χ ) 2 4 2 | v | 2 dx + ( χ 1 ) 2 ( κ 2 t v( t )u( t )dx λ 2 u( t )v( t )dx + 2 u 2 ( t )dx ) ( 2 χ χ ) 2 4k { W( t )+Q }+ ( χ 1 ) 2 ( κ 2 t v( t )u( t )dx λ 2 u( t )v( t )dx + 2 u 2 ( t )dx ).

利用Young不等式,对任意 ε >0 ,任意的 t( 0,T ) ,有

d dt W( t )+κ t v( t ) L 2 ( 2 ) 2 + 2 1 1+u( t ) | u( t ) | 2 dx    ε κ ( χ 1 ) 2 t v( t ) L 2 ( 2 ) 2      λ ( χ 1 ) 2 2 u( t )v( t )dx + ( 2 χ χ ) 2 4k { W( t )+Q }      +( κ ( χ 1 ) 2 4 ε + ( χ 1 ) 2 ) u( t ) L 2 ( 2 ) 2 .

对任意的 t( 0,T ) ,存在常数 C 1 ( ε ),  C 2 >0

d dt W( t )+( 1 ε * ( χ 1 ) 2 )κ t v( t ) L 2 ( 2 ) 2 + 2 1 1+u( t ) | u( t ) | 2 dx C 1 ( ε ) u( t ) L 2 ( 2 ) 2 + C 2 { W( t )+Q }, (4.8)

其中 C 1 ( ε ):= κ ( χ 1 ) 2 4 ε * + ( χ 1 ) 2 , C 2 := ( 2 χ χ ) 2 4k +λ ( χ 1 ) 2 .

由引理3.6和 M<min{ 8π χ , 1 2 C 1 ( ε * ) C 3 } 知,存在常数 C 3 , C 4 >0 ,对任意的 t( 0,T ) ,有

d dt W( t )+( 1 ε * ( χ 1 ) 2 )κ t v( t ) L 2 ( 2 ) 2 + C 4 I( u ) C 2 { W( t )+Q }+ C 1 ( ε * ) C 3 M,

其中  C 3 :=max{ C GN 2 ,9 },  C 4 :=1 C 1 ( ε * ) C 3 M

故对每个 χ>1 κ0

d dt W( t )+ κ 2 t v( t ) L 2 ( 2 ) 2 + 1 2 I( u ) C 2 { W( t )+Q }+ C 1 ( ε * ) C 3 M.

再由(3.1)式与引理4.1,对任意的 t( 0,T ) ,可以估计

d dt W( t )+ κ 2 t v( t ) L 2 ( 2 ) 2 = d dt ( W( t )+ κ 2 0 t t v( t ) L 2 ( 2 ) 2 dt +Q )                                        C 2 { W( t )+Q }+ C 1 ( ε * ) C 3 M                                        C 2 { W( t )+ κ 2 0 t t v( t ) L 2 ( 2 ) 2 dt +Q }+ C 1 ( ε * ) C 3 M. (4.9)

根据Gronwall不等式,对任意的 t( 0,T ) ,可以得到

W( t )+ κ 2 0 t t v( t ) L 2 ( 2 ) 2 dt +Q e C 2 t [ W( 0 )+Q+ 0 t C 1 ( ε * ) C 3 Mds ] e C 2 t ( W( 0 )+Q )+ e C 2 t C 1 ( ε * ) C 3 Mt. (4.10)

结合(4.10)式和引理4.1,可知存在依赖于 M v 0 L 1 ( 2 ) T 的常数 R>0 ,使得对任意的 t( 0,T )

δ 2 ( u( t )+1 ) log( u( t )+1 )dx+ κ 2 0 t t v( τ ) L 2 ( 2 ) 2 dτ R .

下面讨论 v( t ) L 2 ( 2 ) 2 的有界性。对任意的 t( 0,T ) ,对任意的 t( 0,T ) ,由(4.7)和(4.10)式,可以得到

v( t ) L 2 ( 2 ) 2 1 k { e C 2 t ( W( 0 )+Q )+ e C 2 t C 1 ( ε * ) C 3 Mt } . (4.11)

故存在依赖于 M v 0 L 1 ( 2 ) T 的常数 H>0 ,使得对任意的 t( 0,T ) ,有

v( t ) L 2 ( 2 ) 2 H.

进一步,通过标准的乘测试函数的过程,现在可以得到关于解 ( u,v ) 更高阶的范数估计。

引理4.3 κ0 χ>0 λ>0 ( u 0 , v 0 )( L 1 ( 2 ) L ( 2 ) )×( L 1 ( 2 ) H 1 ( 2 ) ) ,假设 ( u,v ) 为模型(2.1)在 2 ×( 0,T ) 上的解,其中 0<T<+ 。//

若存在 R δ>0 使得对任意的 t( 0,T ) ,有

δ 2 ( u( t )+1 ) log( u( t )+1 )dx+ κ 2 0 t t v( τ ) L 2 ( 2 ) 2 dτ R, (4.12)

则存在依赖于 T 的常数 G>0 和依赖于 u 0 L 1 ( 2 ) v 0 L 1 ( 2 ) T 的常数 K>0 ,使得对任意的 t[ 0,T ] ,有

u( t ) L 2 ( 2 ) G,  v( t ) L 2 ( 2 ) K.

证明 对于系统(2.1)的局部光滑解 ( u,v ) ,任意的 t( 0,T ) ,对第一个方程乘以 u 并进行分部积分后,结合第二个方程得到

d dt u( t ) L 2 ( 2 ) 2 +2 u( t ) L 2 ( 2 ) 2 =χ 2 u 2 ( t ) ( κ t v( t )+λv( t )u( t ) )dx κχ 2 u 2 ( t ) t v( t )dx λχ 2 u 2 ( t )v( t )dx +χ u( t ) L 3 ( 2 ) 3 . (4.13)

对任意的 t( 0,T ) ,由引理3.5和假设(4.12)可知,任意 ε>0 ,存在常数 C 1 >0 ,有

u( t ) L 3 ( 2 ) ε ( 1+u( t ) )log( 1+u( t ) ) L 1 ( 2 ) 1/3 u( t ) L 2 ( 2 ) 2/3 + C 1 u( t ) L 1 ( 2 ) 1/3 ε u( t ) L 2 ( 2 ) 2/3 + C 1 u( t ) L 1 ( 2 ) 1/3 , (4.14)

再由Hölder不等式和Gagliard-Nirenberg不等式可知,任意的 t( 0,T ) ,存在常数 C 2 >0 ,有

2 u 2 ( t ) t v( t )dx u( t ) L 4 ( 2 ) 2 t v( t ) L 2 ( 2 ) ε u( t ) L 2 ( 2 ) 2 + C 2 t v( t ) L 2 ( 2 ) 2 u( t ) L 2 ( 2 ) 2 . (4.15)

由(4.12)~(4.15)式可知,对任意的 t( 0,T )

d dt u( t ) L 2 ( 2 ) 2 +2 u( t ) L 2 ( 2 ) 2 εχ u( t ) L 2 ( 2 ) 2 +κχ C 2 t v( t ) L 2 ( 2 ) 2 u( t ) L 2 ( 2 ) 2 + C 1 χ u( t ) L 1 ( 2 ) .

取足够小的 ε>0 ,存在常数 C 3 >0 ,对任意的 t( 0,T ) ,有

d dt u( t ) L 2 ( 2 ) 2 + u( t ) L 2 ( 2 ) 2 C 3 χ{ κ t v( t ) L 2 ( 2 ) 2 u( t ) L 2 ( 2 ) 2 +E }

其中 E:= sup t[0,T] ( u( t ) L 1 ( 2 ) + u( t ) L 1 ( 2 ) 2 )

注意到 u( t ) L 1 ( 2 ) = u 0 L 1 ( 2 ) ,并且由Gronwall不等式可知,对任意的 t[ 0,T ] ,有

u( t ) L 2 ( 2 ) 2 u 0 L 2 ( 2 ) 2 e ( C 3 0 t χκ t v( s ) L 2 ( 2 ) 2 ds ) +{ C 3 χE 0 t e ( C 3 s t χκ t v( τ ) L 2 ( 2 ) 2 dτ ) ds }.

从而由假设(4.12)可知,存在依赖于 T 的常数 G>0 ,使得对任意 t[ 0,T ] ,有

u( t ) L 2 ( 2 ) G. (4.16)

进一步,对第二个方程使用分部积分并结合(4.16)式可知,存在常数 C 4 >0 ,对任意的 t[ 0,T ] ,有

κ d dt v( t ) L 2 ( 2 ) 2 =2 2 | v( t ) | 2 dx 2λ 2 v 2 ( t )dx +2 2 u( t )v( t )dx 2 v( t ) L 2 ( 2 ) 2 2λ v( t ) L 2 ( 2 ) 2 +2 v( t ) L 2 ( 2 ) u( t ) L 2 ( 2 ) 2λ v( t ) L 2 ( 2 ) 2 +2 C 4 v( t ) L 2 ( 2 ) ,

即,存在常数 C 5 >0 ,对任意的 t[ 0,T ] ,有

κ d dt v( t ) L 2 ( 2 ) λ v( t ) L 2 ( 2 ) + C 5 .

κ=0 时,引理4.3的结论成立。当 κ>0 时,利用Gronwall不等式,对任意 t[ 0,T ] ,有

v( t ) L 2 ( 2 ) e λ κ t v 0 L 2 ( 2 ) + C 5 λ ( 1 e λ κ t )<

从而存在依赖于 u 0 L 2 ( 2 ) v 0 L 2 ( 2 ) T 的常数 K>0 ,使得对任意 t[ 0,T ] ,有

v( t ) L 2 ( 2 ) K.

最后,我们证明本文的主要结论

定理2.1的证明 根据标准抛物(椭圆)正则性理论,可以由引理4.2和引理4.3得到:对任意的 T>0 ,解 ( u,v ) L ( 2 ×[ 0,T ] ). 因此解在 [ 0,T ] 内不会发生爆破。我们可在 T 时刻后延拓该解,并注意到 u( T ) L 1 ( 2 ) = u 0 L 1 ( 2 ) 保证了先验估计对延拓后的解仍然成立。由此,解可被无限次延拓至任意长时间,这便证明了解的整体存在性。 □

5. 结论

本文在二维全空间的Keller-Segel模型( κ0,χ>1,λ>0 )框架下,基于能量–熵结构与关键积分不等式,建立了小质量条件下经典解整体存在性的统一判据,同时覆盖抛物–抛物与抛物–椭圆两类模型。主要困难在于 χ>1 引致的增益耦合使交叉项符号不定,且与 κ>0 的时间演化项耦合。本文通过构造能量泛函、精确配方与吸收处理耦合项,并联合运用Gagliardo-Nirenberg与Trudinger-Moser型不等式以及对 u L 2 ( 2 ) 的一致先验控制,获得与时间无关的 L 2 与梯度估计,从而实现解的全局延拓。需要指出的是,所得判据为充分条件,其尖锐性仍有进一步提升空间;对初值假设 u 0 L 1 L , v 0 L 1 H 1 且分析目前限于二维全空间与常系数情形。在此基础上,后续工作可致力于提高小质量阈值的尖锐性并刻画临界/次临界质量的长期渐近,同时将本文的能量–熵与插值框架推广至弱初值(仅 L 1 或测度初值)、有界域(Neumann边界)及非常系数模型。

参考文献

[1] Keller, E.F. and Segel, L.A. (1970) Initiation of Slime Mold Aggregation Viewed as an Instability. Journal of Theoretical Biology, 26, 399-415. [Google Scholar] [CrossRef] [PubMed]
[2] Keller, E.F. and Segel, L.A. (1971) Model for Chemotaxis. Journal of Theoretical Biology, 30, 225-234. [Google Scholar] [CrossRef] [PubMed]
[3] Horstmann, D. (2003) From 1970 until Present: The Keller-Segel Model in Chemotaxis and Its Consequences I. Jahresbericht der Deutschen Mathematiker-Vereinigung, 105, 103-165.
[4] Hassan, Z., Shen, W. and Zhang, Y.P. (2023) Global Existence of Classical Solutions of Chemotaxis Systems with Logistic Source and Consumption or Linear Signal Production on .
https://arxiv.org/abs/2310.16001
[5] Goldstein, S. (1948) On Laminar Boundary-Layer Flow Near a Position of Separation. The Quarterly Journal of Mechanics and Applied Mathematics, 1, 43-69. [Google Scholar] [CrossRef
[6] Hillen, T. and Painter, K.J. (2009) A User’s Guide to PDE Models for Chemotaxis. Journal of Mathematical Biology, 58, 183-217. [Google Scholar] [CrossRef] [PubMed]
[7] Biler, P., Boritchev, A. and Brandolese, L. (2022) Large Global Solutions of the Parabolic-Parabolic Keller-Segel System in Higher Dimensions. Journal of Differential Equations, 344, 891-914. [Google Scholar] [CrossRef
[8] Murray, J.D. and Myerscough, M.R. (1991) Pigmentation Pattern Formation on Snakes. Journal of Theoretical Biology, 149, 339-360. [Google Scholar] [CrossRef] [PubMed]
[9] Horstmann, D. (2) From 1970 until Present: The Keller-Segel Model in Chemotaxis and Its Consequences II. Jahresbericht der Deutschen Mathematiker-Vereinigung, 106, 51-69.
[10] Zhao, J. (2023) Global Existence of Large Solutions for the Parabolic-Elliptic Keller-Segel System in Besov Type Spaces. Applied Mathematics Letters, 149, Artilce 108899. [Google Scholar] [CrossRef
[11] Zheng, Z., Wang, J. and Cai, L. (2024) Global Boundedness in a Keller-Segel System with Nonlinear Indirect Signal Consumption Mechanism. Electronic Research Archive, 32, 4796-4808. [Google Scholar] [CrossRef
[12] Verhulst, P.F. (1838) Notice sur la loi que la population suit dans son accroissement. Journal of Mathematical Physics, 10, 113-116.
[13] Nagai, T., Senba, T. and Yoshida, K. (1997) Application of the Trudinger-Moser Inequality to a Parabolic System of Chemotaxis. Funkcialaj Ekvacioj, 40, 411-433.
[14] Mizoguchi, T. and Winkler, M. (2021) Blow-Up in the Two-Dimensional Parabolic Keller-Segel System. Indiana University Mathematics Journal, 70, 2253-2307.
[15] Gajewski, H., Zacharias, K. and Gröger, K. (1998) Global Behaviour of a Reaction‐Diffusion System Modelling Chemotaxis. Mathematische Nachrichten, 195, 77-114. [Google Scholar] [CrossRef
[16] Biler, P. (1998) Local and Global Solvability of Some Parabolic Systems Modeling Chemotaxis. Advances in Mathematical Sciences and Applications, 8, 715-743.
[17] Jäger, W. and Luckhaus, S. (1992) On Explosions of Solutions to a System of Partial Differential Equations Modelling Chemotaxis. Transactions of the American Mathematical Society, 329, 819-824. [Google Scholar] [CrossRef
[18] Nagai, T. (1995) Blow-Up of Radially Symmetric Solutions to a Chemotaxis System. Advances in Mathematical Sciences and Applications, 5, 581-601.
[19] Nagai, T. and Ogawa, T. (2011) Brezis-Merle Inequalities and Application to the Global Existence of the Cauchy Problem of the Keller-Segel System. Communications in Contemporary Mathematics, 13, 795-812. [Google Scholar] [CrossRef
[20] Mizoguchi, N. (2013) Global Existence for the Cauchy Problem of the Parabolic-Parabolic Keller-Segel System on the Plane. Calculus of Variations and Partial Differential Equations, 48, 491-505. [Google Scholar] [CrossRef
[21] Nagai, T., Senba, T. and Suzuki, T. (2000) Chemotactic Collapse in a Parabolic System of Mathematical Biology. Hiroshima Mathematical Journal, 30, 463-497. [Google Scholar] [CrossRef
[22] Moser, J. (1971) A Sharp Form of an Inequality by N. Trudinger. Indiana University Mathematics Journal, 20, 1077-1092.
[23] Trudinger, N. (1967) On Imbeddings into Orlicz Spaces and Some Applications. Indiana University Mathematics Journal, 17, 473-483. [Google Scholar] [CrossRef
[24] Heihoff, F. (2023) Two New Functional Inequalities and Their Application to the Eventual Smoothness of Solutions to a Chemotaxis-Navier-Stokes System with Rotational Flux. SIAM Journal on Mathematical Analysis, 55, 7113-7154. [Google Scholar] [CrossRef