一类带不连续值的二维趋化系统解的全局存在性
Global Existence of a Two-Dimension Chemotaxis System with Discontinuous Data
DOI: 10.12677/PM.2023.134107, PDF, HTML, XML, 下载: 166  浏览: 266 
作者: 钟 华, 彭红云:广东工业大学数学与统计学学院,广东 广州
关键词: 趋化渐近稳定不连续初值有效粘性通量Chemotaxis Asymptotic Stability Discontinous Initial Data Effective Viscous Flux
摘要: 本文研究了描述肿瘤血管生成的二维趋化系统解的适定性和大时间行为。证明了当时间趋于无穷时,该系统的解收敛到一个常平衡态。同以往的结果相比,我们研究了不连续初值解的适定性。利用有效粘性通量和加权能量估计,得到了解的一系列先验估计,通过这些先验估计可以得到弱解的全局存在性。以往的研究都是关于连续初值的研究,得到的都是经典解。而实际情况中不连续的情形更加普遍,而本文的结果表明,在不连续初始值的情况下,我们能够得到二维趋化系统的全局弱解。
Abstract: This paper is concerned with the well-posedness and large-time behavior of a two-dimensional PDEODE hybrid chemotaxis system describing the initiation of tumor angiogenesis. We prove the solution converge to the constant equilibrium when the time tends to infinity. In contrast to the ex-isting related results, where continuous initial data is imposed, we are able to prove the asymptotic stability for discontinuous initial data with large oscillations. The key ingredient in our proof is the use of the so-called “effective viscous flux”, which enables us to obtain the desired energy estimates and regularity. The previous results are all about the continuous initial data and the classical solutions are obtained. However, discontinuity is more common in practice, and the results of this paper show that we can obtain the global weak solution of the two-dimensional chemotaxis system under the discontinuous initial data.
文章引用:钟华, 彭红云. 一类带不连续值的二维趋化系统解的全局存在性[J]. 理论数学, 2023, 13(4): 1018-1032. https://doi.org/10.12677/PM.2023.134107

1. 前言和主要的结果

1.1. 背景

本文研究以下趋化模型

{ u t = Δ u ( ξ u ln c ) , c t = μ u c + ϵ Δ u . (1)

该模型描述了在血管形成期间,血管内皮细胞生长因子与血管内皮细胞之间的关系 [1] [2] 。其中 u ( x , t ) 表示血管内皮细胞的密度, c ( x , t ) 表示血管内皮细胞生长因子的浓度。参数 ξ > 0 代表趋化系数, μ 表示趋化血管内皮细胞生长因子的消耗率, ϵ 表示生长因子的扩散率。

模型(1)中,第一个方程含有一个对数敏感函数 ln c ,为了消去对数奇性,可以通过如下Cole-Hopf变换 [3] [4]

v = 1 μ ( ln c ) = 1 μ c c , (2)

和伸缩变换 t ˜ = ϵ μ t x ˜ = ϵ μ v ˜ = ϵ μ v ,将(1.1)转换为如下系统(为了方便起见,还是用 t , v , x 分别代替 t ˜ , v ˜ , x ˜ )

{ u t ( u v ) = Δ u , v t ( ϵ | v | 2 + u ) = ϵ Δ v , ( u , v ) ( x , 0 ) = ( u 0 , v 0 ) = ( u 0 , v 0 ) ( x ) , (3)

无穷远状态为

( u , v ) ( ± , t ) = ( u ¯ , 0 ) , (4)

其中 u ¯ > 0

对于模型(1)已经有了大量的研究,比如在一维空间中该系统行波解稳定性的研究 [5] [6] [7] [8] ,关于解的全局存在性以及大时间行为的研究 [9] - [14] 。对于高维无界区域的情形,研究者在文献 [1] 中给出了该模型在二维空间非线性行波解的稳定性,在文献 [15] [16] 中得到了经典解的爆破准则,在文献 [4] 中得到了该系统的全局适定性以及解的大时间行为。对于高维有界区域 Ω d ( d 3 ) ,在Neuman边界条件下,研究者在文献 [16] 中得到了具有指数衰减的全局解。在文献 [11] 中,研究者在二维空间研究了模型(1)的边界层问题。关于模(1)的型其它研究请参考文献 [17] [18] 。

以上的研究大都是关于连续初值的研究,在连续初始值的情况下,可以要求初始值足够光滑,这时得到的都是经典解。而在实际应用中不连续的情形更加普遍(比如分片光滑的初值就是不连续的),因此研究带不连续初值的趋化模型是一个重要而有意义的题。本文将研究二维趋化系统在不连续初值下解的存在性问题,进而考虑解的大时间行为,这有助于研究该模型在不连续初值下解的性态。

首先我们先介绍弱解的定义:

定义1.1我们说 ( u , v ) 是系统(3)~(4)的弱解,如果 ( u , v ) 满足,对于所有的测试函数 ϕ C 0 ( 2 × [ 0 , ) ) 满足以下关系

2 u 0 ϕ 0 + 0 2 ( u ϕ t u ϕ ) d x d t = 0 2 u v ϕ d t

2 v 0 j ϕ 0 d x + 0 2 ( v j ϕ t u ϕ x j + ϵ ϕ x j | v j | 2 ϵ ϕ x j x j v j ) d x d t = 0 ,

其中 i , j = 1 , 2 ϕ 0 ( x ) = ϕ ( x , 0 )

1.2. 主要结果

定理1.2如果 4 < p 0 < 如果初始值满足

u 0 1 L 2 ( 2 ) , v 0 L 2 ( 2 ) L p 0 ( 2 ) , v 0 0 , v 0 = 0 , (5)

= ( 2 , 1 ) 是旋度算子。然后对于任意的 M > 0 v 0 L p 0 ( 2 ) M ,那么存在 η > 0 ,使得

u 0 1 L 2 ( 2 ) 2 + v 0 L 2 ( 2 ) 2 = θ 0 η ,

那么Cauchy问题(3)-(4)有一个全局的弱解 ( u , v ) ,使得

{ u 1 L ( [ 0 , ) ; L 2 ( 2 ) ) C ( ( 0 , ) ; C ( 2 ) ) , v , u L 2 ( [ 0 , ) ; L 2 ( 2 ) ) v L ( [ 0 , ) ; L 2 ( 2 ) L p 0 ( ℝ 2 )

2 < p 1 p 2 ,则

u ( x , t ) 1 L p 1 ( 2 ) 0 , v ( x , t ) L p 2 ( 2 ) 0 , t .

2. 预备知识

下面介绍一些概念:

H k ( 2 ) 代表 2 上指标为k的Sobolev空间,其范数定义为

f H k ( 2 ) : = ( j = 0 k 2 | x j f | 2 d x ) 1 2 ,

旋度算子

= ( 2 , 1 ) , (7)

θ 0 = u 0 1 L 2 2 + v 0 L 2 2 . (8)

其中 θ 0 < 1 ,下面我们将简写 L 2 。首先介绍以下两个引理 [5] [16] 。

引理2.1令 s > d 2 ( u 0 1 , v 0 ) H s ( d ) 。那么存在一个最大时间 T * = ( u 0 1 H s , v 0 H s ) > 0 使得Cauchy问题(3)~(4),存在一个唯一解使得

( u 1 , v ) L ( ( 0 , T * ] , H s ( d ) ) .

引理2.2令 s > d 2 + 1 ,且 ( u 0 1 , v 0 ) H s ( d ) ( u , v ) 是引理2.1中得到的解,如果

0 T v L q 2 q q d d t < ( d < q < ) , (9)

那么解可以延拓到 T * 之外。

作变量替换 u ˜ = u 1 ,系统(3)可转化为

{ u ˜ t Δ u ˜ = ( u ˜ v ) + v , v t u ˜ = ϵ ( | v | 2 + Δ v ) , ( u ˜ , v ) ( x , 0 ) = ( u 0 1 , v 0 ) ( x ) . (10)

定义有效粘性通量 F 为:

F = u ˜ + ( u ˜ + 1 ) v . (11)

那么(10)中的第一个方程可以写成

F = u ˜ t . (12)

由此可以得到关于F的估计:

引理2.3 [19] 如果 ( u ˜ , v ) 是(10)的光滑解,那么存在常数C使得

F L p C ( u ˜ t L p ) + u ˜ v L p ( p > 1 ) .

引理2.4 [20] 令 1 q r 0 < a 1 使得

1 p = a ( 1 q 1 2 ) + ( 1 a ) 1 r ,

那么对于 u W 1 , p ( 2 ) L r ( 2 ) ,然后存在一个常数只依赖于 q , r , n 的C,使得

u L p C u L q ( 2 ) a u L r ( 2 ) 1 a . (13)

3. 渐进解的先验估计

首先对初值进行磨光

u ˜ 0 δ = j δ u ¯ 0 , v 0 δ = j δ v 0 .

考虑以下渐进系统

{ u ˜ t δ Δ u ˜ δ = ( u ˜ δ v δ ) + v δ v t δ u ˜ δ = ϵ ( ( | v | δ ) 2 + Δ v δ ) , (14)

其中初始值为

( u ˜ 0 δ , v 0 δ ) H 3 ( ℝ 2 )

u ˜ 0 δ 2 + v 0 δ 2 u ˜ 0 2 + v 0 2 = θ 0 .

为方便起见,下面用 ( u ˜ , v ) 代替 ( u ˜ δ , v δ ) 。令 T > 0 ( u ˜ , v ) 是(14)在 2 × [ 0 , T ] 上的光滑解,令 σ ( t ) = min { 1 , t } ,定义

{ A 1 ( T ) = sup t [ 0 , T ] ( u ˜ 2 + v 2 ) + 0 T u ˜ 2 d t + ϵ 0 T v 2 d t , A 2 ( T ) = sup t [ 0 , T ] ( σ u ˜ 2 + σ 2 u ˜ t 2 + σ 2 v t 2 + σ ϵ v 2 ) + 0 T σ 2 u ˜ t 2 d t + ϵ 0 T σ 2 v 2 d t + 0 T σ v t 2 d t + 0 T σ u ˜ t 2 d t , A 3 ( T ) = sup t [ 0 , T ] v L 4 4 + 0 T v L 4 4 d t . (15)

下面我们将通过先验假设的方法来得到解的先验估计,我们首先假设

A 1 ( T ) 2 θ 0 , A 2 ( T ) 2 θ 0 1 2 , A 3 ( T ) 2 θ 0 η 0 , sup t [ 0 , T ] v L p 0 6 M , (16)

其中

η 0 = p 0 4 2 ( p 0 2 ) ( 0 , 1 2 ) . (17)

下面我们从基本能量估计开始。

引理3.1如果定理1.2中的条件满足, ( u ˜ , v ) 是满足初值条件并且满足(16)的光滑解,则有下面不等式成立

u ˜ 2 + v 2 + 0 T u ˜ 2 d t + ϵ 0 T v 2 d t 3 2 θ 0 , (18)

证明 用 u ˜ 乘以(10)式中的第一个方程,用 σ v 乘以(10)中的第二个方程,把两式所得结果相加,然后在 2 上作积分,则有然后用类似文献 [19] 引理3.1的方法,容易得到所需的结果。

接下来我们将给出解的一阶能量估计,得到解的一阶能量估计时需要利用加权能量的技巧。

引理3.2如果定理1.2中的条件满足, ( u ˜ , v ) 是满足初值条件并且满足(16)的光滑解,则有

σ u ˜ 2 + σ u ˜ t 2 + σ ϵ v t 2 + σ ϵ v 2 + v t 2 d t + 0 T σ u ˜ t 2 d t + ϵ 0 T σ 2 v t 2 d t + 1 4 0 T σ 2 u ˜ t 2 d t θ 0 1 2 , (19)

其中 σ = σ ( t ) = min { 1 , t }

证明 第一步:用 σ u ˜ t σ v t ,分别乘以(14)中第一个方程,第二个方程,然后在 2 × ( 0 , T ] 作积分,可得到

1 2 σ u ˜ 2 + 0 T σ u ˜ t 2 d t = 1 2 0 σ ( T ) u ˜ 2 d t 0 T σ 2 u ˜ v u ˜ t d x d t 0 T σ 2 v u ˜ t d x d t (20)

1 2 σ ϵ v 2 + 0 T σ v t 2 d t = ϵ 0 T σ 2 ( v t | v | 2 ) d x d t + ϵ 0 σ ( T ) v 2 d t + 0 T σ 2 v t u ˜ d x d t . (21)

其中上式用了 1 2 0 T σ t u ˜ 2 d t = 1 2 0 σ ( T ) u ˜ 2 d t 1 2 0 T σ t v 2 d t = 1 2 0 σ ( T ) v 2 d t

首先处理(20)式右边的第二项,根据G-N不等式有

0 T σ 2 u ˜ v u ˜ t d x d t λ 0 T σ 2 u ˜ t 2 d t + C 0 T u ˜ v 2 d t λ 0 T σ 2 u ˜ t 2 d t + C 0 T u ˜ L 4 2 v L 4 2 d t λ 0 T σ 2 u ˜ t 2 d t + C 0 T u ˜ 2 d t + C 0 T u ˜ 2 v L 4 4 d t λ 0 T σ 2 u ˜ t 2 d t + C θ 0 . (22)

对于(21)中右边最后一项,通过计算可得

0 T 2 v u ˜ t d x d t = 0 T ( σ 2 v u ˜ d x ) t d t + 0 σ ( T ) 2 v u ˜ d x d t + 0 T σ 2 v t u ˜ d x d t = : I 1 + I 2 + I 3 , (23)

由Cauchy不等式,(15),(16),及 0 σ 1 ,则有

I 1 = σ 2 v u ˜ d x σ 4 u ˜ 2 + σ v 2 σ 4 u ˜ 2 + 2 θ 0 , (24)

I 2 = 0 σ ( T ) 2 v u ˜ d x d t 1 2 0 σ ( T ) u ˜ 2 d t + 1 2 0 σ ( T ) v 2 d t 1 2 0 T u ˜ 2 d t + 1 2 sup t [ 0 , T ] v 2 2 θ 0 . (25)

I 3 = 0 T σ 2 v t u ˜ d x d t = 0 T σ 2 | u ˜ | 2 d x d t + ϵ 0 T σ 2 | v | 2 u ˜ d x d t + ϵ 0 T σ 2 Δ v u ˜ d x d t ,

对于 I 3 中的第二项,由Cauchy不等式以及(16)有

ϵ 0 T σ 2 | v | 2 u ˜ d x d t C 0 T σ 2 ϵ | v | | v | | u ˜ | d x d t C θ 0 1 2 + η 0 .

I 3 中的第三项,由(14)中的第一个等式有

I 3 = 0 T σ 2 v t u ˜ d x d t C θ 1 2 + η 0 . (26)

由(20),(22),(23),(24),(25),(26)得

1 4 σ u ˜ 2 + 0 T σ u ˜ t 2 d t λ 0 T σ 2 u ˜ t 2 d t + C θ 0 1 2 + η 0 . (27)

并且

0 T σ 2 u v t d x d t 1 8 0 T σ 2 v t 2 d x d t + C 0 T σ 2 | u ˜ | 2 d x d t 1 8 0 T σ 2 v t 2 d x d t + C θ 0 1 2 + η 0 , (28)

ϵ 0 T 2 σ | v | 2 v t d x d t ϵ 0 T 2 σ | v | | σ | | v | | v t | d x d t ϵ σ v L ( 1 4 0 T 2 σ | v t | 2 d x d t + 0 T 2 | v | 2 d x d t ) ϵ 4 0 T σ v t 2 d t + C θ 0 1 2 + η 0 , (29)

则有

σ u ˜ 2 + σ ϵ v 2 + 0 T σ v t 2 d t + 0 T σ u ˜ t 2 d t λ 0 T σ 2 u ˜ t 2 d t + C θ 0 1 2 + η 0 . (30)

步骤2. 估计 0 T σ 2 u ˜ t 2 d t ,对时间求导得

{ u ˜ t t Δ u ˜ t = ( u ˜ v ) t + v t v t t u ˜ t = 2 ϵ ( | v | 2 ) t + ϵ Δ v t . (31)

所得方程第一个方程乘以 σ 2 u ˜ t ,第二个方程乘以 σ 2 v t ,把所得结果相加并在作积分可得

σ 2 2 u ˜ t 2 + σ 2 2 v t + 0 T σ 2 u ˜ t 2 d t + ϵ 0 T σ 2 v t 2 d t 0 σ ( T ) σ u ˜ t 2 d t + 0 σ ( T ) σ v t 2 d t + ϵ 0 σ ( T ) σ 2 2 v t ( | v | 2 ) d x d t + 0 T σ 2 2 | ( u ˜ v ) t | | u ˜ t | d x d t 0 T σ u ˜ t 2 d t + 0 σ ( T ) σ v t 2 d t + 0 T σ 2 2 | u ˜ v t u ˜ t | d x d t + 0 T σ 2 2 | u ˜ t v u ˜ t | d x d t + ϵ | 0 σ ( T ) σ 2 2 v t ( | v | 2 ) t d x d t | . (32)

因此有

σ 2 2 u ˜ t 2 + σ 2 2 v t 2 + 0 T σ 2 u ˜ t 2 d t + ϵ 0 T σ 2 v t 2 d t 1 4 θ 0 1 2 + 0 σ ( T ) σ v t 2 d t + 0 T σ u ˜ t 2 d t + 0 T σ 2 2 | u ˜ v t u ˜ t | d x d t + 0 T 2 σ 2 u ˜ t v u ˜ t d x d t . (34)

下面估计(34)右边第四项,由

0 T σ 2 2 | u ˜ v t u ˜ t | d x d t 1 4 0 T σ 2 2 | u ˜ t | 2 d x d t + 0 T σ 2 2 | u ˜ v t | 2 d x d t 1 4 0 T σ 2 2 | u ˜ t | 2 d x d t + ( σ u ˜ L ) 2 0 T 2 σ | v t | 2 d t d x 1 4 0 T σ 2 2 | u ˜ t | 2 d x d t + C θ 0 1 2 + η 0 . (35)

下面估计(34)中右边最后一项,由(13)和(16)可得

0 T 2 u ˜ t v u t d x d t 1 4 0 T σ 2 u ˜ t 2 d t + 0 T σ 2 u ˜ t v 4 d t 1 4 0 T σ 2 u ˜ t 2 d t + 0 T σ 2 u ˜ t L 4 2 v L 4 2 d t 1 4 0 T σ 2 u ˜ t 2 d t + C 0 T σ 2 u ˜ t u ˜ t v L 4 2 d t 1 2 0 T σ 2 u ˜ t 2 d t + C 0 T σ 2 u ˜ t 2 v L 4 4 d t 1 2 0 T σ 2 u ˜ t 2 d t + C θ 0 η 0 0 T σ 2 u ˜ t 2 d t . (36)

根据 0 σ 1 θ 0 1 ,可得

1 2 σ 2 u ˜ t 2 + 1 2 σ 2 v t 2 + 1 4 0 T σ 2 u ˜ t 2 d t + ϵ 0 T σ 2 v t 2 d t C λ 0 T σ 2 u ˜ t 2 d t + 0 T σ v t 2 d t + C θ 0 1 2 + η 0 + 1 4 θ 0 1 2 ,

选取 ( C + 4 ) λ 1 8 ,可得

σ u ˜ 2 + σ ϵ v 2 + σ 2 u ˜ t 2 + σ 2 v t 2 + 100 0 T σ v t 2 d t + 0 T σ u ˜ t 2 d t + 1 4 0 T σ 2 u ˜ t d t + ϵ 0 T σ 2 v t 2 d t C θ 0 1 2 + η 0 + 1 2 θ 0 1 2 . (37)

选取 θ 0 ,使得 C θ 0 η 0 1 2 ,则有

σ u ˜ 2 + σ ϵ v 2 + σ 2 u ˜ t 2 + σ 2 v t 2 + 100 0 T σ v t 2 d t + 0 T σ u ˜ t 2 d t + 1 4 0 T σ 2 u ˜ t d t + ϵ 0 T σ 2 v t 2 d t θ 0 1 2 .

因此完成引理3.2的证明。

接下来,我们将给出 u ˜ ( x , t ) 的估计,从而得到 u ˜ ( x , t ) 的大时间行为。

引理3.3 [19] 如果定理1.2中的条件满足, ( u ˜ , v ) 是满足初值条件并且满足(16)的光滑解,则有

σ 4 u ˜ L 4 C θ 0 1 2 + η 0 ,

t ( σ ( T ) , T ] ,则

1 4 u ˜ ( x , t ) 1 4 . (38)

证明:引理3.3的证明可参考文献 [19] 。

以下引理是关于v的估计。

引理3.4如果定理1.2中的条件满足, ( u ˜ , v ) 是满足初值条件并且满足(16)的光滑解,那么

sup t [ 0 , T ] 2 v 4 d x + 0 T 2 v 4 d x d t θ 0 η 0 ,

其中 η 0 = p 0 4 2 ( p 0 2 )

证明:由 v t = u ˜ + ϵ Δ v + ϵ | v | 2 ,有

Δ v ϵ | v | 2 + ( u ˜ + 1 ) v = F , (39)

两边同时乘以 | v | 2 v ,然后在 2 作积分,则

1 4 ( 2 v 4 d x ) t + 2 ( u ˜ + 1 ) v 4 d x = ϵ 2 | v | 2 v | v | 2 d x + ϵ 2 Δ v | v | 2 v d x + 2 F | v | 2 v d x .

[ σ ( T ) , T ] 作积分,可得

1 4 2 v 4 d x + 3 4 σ ( T ) T 2 v 4 d x d t sup t [ 0 , σ ( T ) ] ( 1 4 2 v 4 d x ) + σ ( T ) T | F | | v | 3 d x d t + C ϵ σ ( T ) T v L 2 2 | v | 2 | v | d t d x 3 ϵ σ ( T ) T 2 | v | 2 | v | 2 d x d t sup t [ 0 , σ ( T ) ] ( 2 v 4 d x ) + C ϵ σ ( T ) T 2 ( | v | 4 + | v | 2 ) d t d x + 1 4 σ ( T ) T 2 v 4 d x d t + C σ ( T ) T | F | 4 d x d t sup t [ 0 , σ ( T ) ] ( 1 4 2 v 4 d x ) + C θ 0 1 2 + η 0 + 1 4 σ ( T ) T 2 v 4 d x d t + C σ ( T ) T 2 | F | 4 d x d t . (40)

由G-N不等式(13)有

σ ( T ) T F L 4 4 d t C σ ( T ) T F 2 F 2 d t C ( sup σ ( T ) t T F 2 ) σ ( T ) T F 2 d t . (41)

由于

F u ˜ + u ˜ v + v u ˜ + u ˜ v + v u ˜ + u ˜ L 4 v L 4 + v u ˜ + C u ˜ 1 2 u ˜ 1 2 v L 4 + v

u ˜ + C u ˜ v L 4 2 + v u ˜ + C θ 0 1 4 + η 0 2 + C θ 0 1 2 u ˜ + C θ 0 1 4 + η 0 2 , (42)

σ ( T ) t T σ = 1 ,则

sup σ ( T ) t T F 2 4 σ u ˜ 2 + C ( θ 0 1 4 + η 0 ) 2 4 θ 0 1 2 + C θ 0 1 2 + η 0 C θ 0 1 2 . (43)

那么由(13)有

σ ( T ) T F 2 d t C σ ( T ) T σ u ˜ t 2 d t + C σ ( T ) T σ u ˜ L 4 2 v L 4 2 d t C θ 0 1 2 + C σ ( T ) T σ u ˜ L 4 2 v L 4 2 d t . (44)

0 T σ 2 u ˜ L 4 4 d t C θ 0 η 0 和(16)可得

C σ ( T ) T σ u ˜ L 4 2 v L 4 2 d t C σ ( T ) T σ 2 u ˜ L 4 4 d t + C σ ( T ) T v L 4 4 d t C θ 0 η 0 , (45)

因此

σ ( T ) T F 2 d t C θ 0 η 0 , (46)

由(41),(44),(46)得

σ ( T ) T F L 4 4 d t C θ 0 1 2 + η 0 .

那么

sup t [ σ ( T ) , T ] 2 v 4 d x + σ ( T ) T v 4 d x d t C ( M ) θ 0 4 η 0 + C θ 0 1 2 + η 0 . (47)

0 t σ ( T ) ,由插值不等式,(3.3)可得

v L p v 2 ( p 0 p ) p ( p 0 2 ) v L p 0 p 0 ( p 2 ) p ( p 0 2 ) C θ 0 p 0 p p ( p 0 2 ) M p 0 ( p 2 ) p ( p 0 2 )

v L 4 C θ 0 p 0 4 4 ( p 0 2 ) M p 0 ( p 2 ) 4 ( p 0 2 ) C ( M ) θ 0 η 0 , (48)

因此

sup t [ 0 , σ ( T ) ] 2 v 4 d x 2 sup t [ 0 , σ ( T ) ] 2 v 4 d x C ( M ) θ 0 4 η 0 . (49)

由(47),(48),(49)取 C ( M ) θ 0 3 η 0 1 2 C θ 0 1 2 1 2 ,则

sup t [ 0 , T ] 2 v 4 d x + 0 T 2 v 4 d t d x C ( M ) θ 0 4 η 0 + 1 2 θ 0 1 2 + η 0 θ 0 η 0 .

引理3.5如果定理1.2中的条件满足, ( u ˜ , v ) 是满足初值条件并且满足(16)的光滑解,那么

sup t [ 0 , T ] v L p 0 3 M . (50)

证明:用 σ σ | v | p 0 2 | v | 乘以(3.27),然后在 上作积分,则有

1 p 0 ( 2 σ p 0 | v | p 0 d x ) t + 2 ( u ˜ + 1 ) ( σ | v | ) p 0 d x = ϵ 2 σ p 0 | v | p 0 2 v | v | 2 d x + ϵ 2 σ p 0 | v | p 0 2 v Δ v d x + 2 σ p 0 | F | | v | p 0 2 d x + 2 σ p 0 1 | v | p 0 d x . (51)

将上式在 [ 0 , T ] 上积分有

2 σ p 0 | v | p 0 d x C ϵ 0 T 2 σ p 0 | v | p 0 | v | d x d t + C ϵ 0 T 2 σ p 0 | v | 2 | v | p 0 2 d x d t + 0 T 2 σ p 0 | F | | v | p 0 1 d x d t + 0 T 2 σ p 0 1 | v | p 0 d x d t ϵ 0 T σ p 0 2 2 v L p 0 2 2 | v | 2 | v | d x d t + C ϵ 0 T ( σ p 0 2 2 v L p 0 2 ) 2 | v | 2 d x d t + 0 T ( v L p 0 4 ) σ p 0 4 2 2 | v | d x d t + 0 T | σ F | | v | p 0 1 d x d t ϵ ( 0 T 2 v 4 + | v | 2 d x d t ) + C θ 0 1 2 + η 0 + C θ 0 η 0 + 0 T | σ F | | σ v | p 0 1 d x d t . (52)

由不等式

a b ϵ a q + ( ϵ q ) r q r 1 b r , a , b 0 , ϵ > 0 , q , r > 0 , 1 q + 1 r = 1 ,

0 T 2 | σ F | | σ v 0 | p 0 1 d x d t 1 2 0 T 2 | σ v | p 0 d x d t + 1 p 0 ( 2 p 0 2 p 0 ) p 0 1 0 T | σ F | p 0 d x d t 1 2 0 T σ p 0 4 2 v L p 0 4 2 | v | 4 d t d x + 1 p 0 ( 2 p 0 2 p 0 ) p 0 1 0 T 2 | σ F | p 0 d x d t ,

因此

2 | v | p 0 d x C θ 0 η 0 + ( 2 p 0 2 p 0 ) p 0 1 0 T 2 | σ F | p 0 d x d t .

由G-N不等式(13)可得

0 T σ F L p 0 p 0 d t C 0 T σ F 2 σ F p 0 2 d t C sup t [ 0 , T ] σ F 2 σ F p 0 4 0 T σ F 2 d t . (53)

sup t [ 0 , σ ( T ) ] σ 1 2 F C σ 1 2 u ˜ + C σ 1 2 θ 0 1 4 C θ 0 1 4

0 T σ F 2 d t C θ 0 η 0 ,

可得

sup t [ 0 , T ] σ F 2 sup t [ 0 , T ] σ F 2 C θ 0 1 2 .

由(10),(13)可得

u ˜ L 4 F u ˜ v v L 4 F L 4 + u ˜ v L 4 + v L 4 C ( F 1 2 F 1 2 + u ˜ L v L 4 + v L 4 ) ( F 1 2 F 1 2 + u ˜ L 4 1 2 u ˜ L 4 1 2 v L 4 + v L 4 ) .

利用上式和(13)可得

下面对 I 4 做估计,由(13)和Hölder不等式可得

F C u ˜ t + C u ˜ v C u ˜ t + C u ˜ L 4 v L 4 .

由(15),(16),(43)及以上不等式,可得

I 4 C 0 T σ ( u ˜ 2 + θ 0 1 2 + η 0 ) σ ( u ˜ t 2 + u ˜ L 4 v L 4 2 ) d t σ ( A 2 ( T ) + θ 0 1 2 + η 0 ) 0 T σ ( u t 2 + u ˜ L 4 2 v L 4 2 ) d t C θ 0 1 2 0 T σ u ˜ t 2 d t + C θ 0 1 2 0 T σ 2 u ˜ L 4 4 d t + C θ 0 1 2 A 3 ( T ) C θ 0 1 2 A 2 ( T ) + C θ 0 1 2 0 T σ 2 u ˜ L 4 4 d t + C θ 0 1 2 A 3 ( T ) C θ 0 + C θ 0 1 2 0 T σ 2 u ˜ L 4 4 d t + C θ 0 1 2 + η 0 C θ 0 1 2 + C θ 0 1 2 0 T σ 2 u ˜ L 4 4 d t .

以下估计 I 5 ,由(13),(16),(17)可得

I 5 = C 0 T σ 2 u ˜ L 4 2 u ˜ L 4 2 d t 1 4 0 T σ 2 u ˜ L 4 4 d t + C A 3 2 ( T ) 0 T 0 T σ 2 u ˜ L 4 4 d t 1 4 0 T σ 2 u ˜ L 4 4 d t + C A 3 ( T ) 2 0 T σ 2 u ˜ 2 u ˜ 2 d t 1 4 0 T σ 2 u ˜ L 4 4 d t + C θ 0 1 2 + 3 η 0 1 4 0 T σ 2 u ˜ L 4 4 d t + C θ 0 η 0 .

θ 0 使得 C θ 0 1 2 1 4 ,那么

0 T σ 2 u ˜ L 4 4 d t C θ 0 η 0 + C θ 0 1 2 C θ 0 η 0 .

0 T σ F 2 d t C 0 T σ u ˜ t 2 d t + C 0 T σ u ˜ L 4 2 v L 4 2 d t 0 T σ u ˜ t 2 d t + C 0 T ( σ 2 u ˜ L 4 4 + v L 4 4 ) d t C θ 0 1 2 + C θ 0 η 0 C θ 0 η 0 .

由于

sup t [ 0 , T ] σ F C θ 0 η 0 2 + C σ u ˜ t + C σ 1 2 u ˜ C σ ( θ 0 η 0 2 + θ 0 1 4 ) C θ 0 η 0 2 ,

那么

sup t [ 0 , T ] σ F p 0 4 C θ 0 η 0 ( p 0 4 ) 2 .

把以上三个不等式代入(53)可得

0 T σ F L p 0 p 0 d t C θ 0 η 0 ( p 0 2 ) + 1 2 = C θ 0 p 0 2 4 , (54)

因此取 θ 0 足够下,当 C θ 0 η 0 + C θ 0 p 0 2 4 ( 3 M ) p 0 时,

sup t [ 0 , T ] v L p 0 3 M .

因此完成了引理3.5的证明。

根据引理3.1~3.5,可得

A 1 ( T ) 3 2 θ 0 , A 2 ( T ) θ 0 1 2 , A 3 ( T ) θ 0 η 0 , sup t [ 0 , T ] v L p 0 3 M . (55)

4. 主定理证明

下面利用爆破准则得到全局解

命题4.1假设 ( u ˜ 0 δ , v 0 δ ) H 3 ( 2 ) ,那么系统(16)有一个唯一解,使得 ( u ˜ δ , v δ ) L ( [ 0 , ) , H 3 )

证明:通过引理2.1,则存在 T * > 0 ,使得满足初始条件的系统(3.1)在 2 × [ 0 , T * ] 上的Cauchy问题有唯一解。首先,从(3.3)可知

A 1 ( 0 ) θ 0 , A 2 ( 0 ) < 2 θ 0 1 2 , v 0 δ L p 0 M .

C ( M ) θ 0 3 η 0 1

A 3 ( 0 ) = v 0 δ L 4 4 C ( M ) θ 0 4 η 0 θ 0 η 0 . (56)

因此,存在一个 T 1 > 0 ,使得 T 1 T * ,(3.3)成立。

反证,如果 ,由爆破准则,则

0 T * v δ L 4 4 d t < .

根据引理3.1-3.5,当 T = T * 时,上式成立,所以存在 T ** > T * ,使得

( u δ 1 , v δ ) L ( [ 0 , T ** ] , H 3 ) ,

和当 T = T ** 时(16)成立,与最大存在时刻矛盾。所以根据连续性技巧,可证命题4.1。

从(55)得到下面关于 δ 的一致估计

{ sup t [ 0 , T ] ( u ˜ δ 2 + v δ 2 ) + ϵ v δ 2 + 0 T u ˜ δ 2 d t + ϵ 0 T v δ 2 d t 3 θ 0 2 , sup t [ 0 , T ] ( σ u ˜ δ 2 + σ 2 u ˜ t δ 2 + σ 2 v t δ 2 + σ ϵ v δ 2 ) + 0 T σ 2 u ˜ t δ 2 d t + ϵ 0 T σ 2 v δ 2 d t + 0 T σ v t δ 2 d t + 0 T σ u ˜ t δ 2 d t θ 0 1 2 , sup t [ 0 , T ] u ˜ δ L 4 4 + 0 T v δ L 4 4 d t θ 0 η 0 , sup t [ 0 , T ] v δ L p 0 3 M . (57)

由(13),(57),可得

σ u δ 1 L 4 2 C u δ 1 σ u δ C , σ 2 u δ L 4 2 C . (58)

所以由(4.2),(4.3),则有

{ u σ 1 L ( [ 0 , ) , L 2 ( 2 ) ) , ( v δ , u δ ) L 2 ( [ 0 , ) , L 2 ( 2 ) ) , u δ 1 L ( ( 0 , ) , W 1 , 4 ( 2 ) ) , u t δ , v t δ L 2 ( ( 0 , ) , H 1 ( 2 ) ) , v L ( [ 0 , ) ; L ( 2 ) L p 0 ( 2 ) ) L 4 ( [ 0 , ) ; L 4 ( 2 ) ) . (59)

W 1 , 4 ( 2 ) ↪↪ C ( 2 ) L 2 ( 2 ) ↪↪ H 1 ( 2 ) ,因此Aubin-Lions-Simon引理,可以抽取一个子列 ( u δ , v δ ) ,当 δ 0 ,则

{ v δ ( , t ) v stronglyin C ( [ 0 , ) , H 1 ( 2 ) ) , u δ ( , t ) u ( , t ) Stronglyin C ( [ 0 , ) , C ( 2 ) ) , u δ ( , t ) u ( , t ) weaklyin L 2 ( [ 0 , ) , L ( 2 ) ) , v δ ( , t ) v ( , t ) weaklyin L 2 ( [ 0 , ) , L ( 2 ) ) . (60)

那么

1 u ˜ 2 d t C ,

1 | ( u ˜ 2 ) | t d t 2 1 u ˜ u ˜ t d t ( 1 u ˜ 2 d t ) 1 2 ( 1 σ 2 u ˜ t 2 d t ) 1 2 C ,

那么

u ˜ 0 , t , (61)

据此可得

u 1 L 0 , t . (62)

根据插值不等式,当 2 < p 1 时,可得

u 1 L p 1 0 , t .

类似地

1 | ( v 2 ) t | d t C ,

那么

v L p 1 0 , t .

至此完成定理1.2的证明。

参考文献

[1] Corrias, L., Perthame, B. and Zaag, H. (2003) A Chemotaxis Model Motivated by Angiogenesis. Comptes Rendus de l’Académie des Sciences Paris, 336, 141-146.
https://doi.org/10.1016/S1631-073X(02)00008-0
[2] Corrias, L., Perthame, B. and Zaag, H. (2004) Global Solutions of Some Chemotaxis and Angiogenesis Systems in High Space Dimensions. Milan Journal of Mathematics, 72, 1-28.
https://doi.org/10.1007/s00032-003-0026-x
[3] Levine, H.A., Sleeman, B.D. and Nilsen-Hamilton, M. (1997) A System of Reaction Diffusion Equations Arising the Theory of Reinforced Random Walks. SIAM Journal on Applied Mathematics, 47, 683-730.
https://doi.org/10.1137/S0036139995291106
[4] Othmer, H.G. and Stevens, A. (1997) Aggregation, Blowup, and Collapse: The ABCs of Taxis in Reinforced Random Walks. SIAM Journal on Applied Mathematics, 57, 1044-1081.
https://doi.org/10.1137/S0036139995288976
[5] Jin, H.Y., Li, J.Y. and Wang, Z.A. (2013) Asymptotic Stability of Traveling Waves of a Chemotaxis Model with Singular Sensitivity. Journal of Differential Equations, 255, 193-219.
https://doi.org/10.1016/j.jde.2013.04.002
[6] Mei, M., Peng, H.Y. and Wang, Z.A. (2015) Asymptotic Profile of a Parabolic-Hyperbolic System with Boundary Effect Arising from Tumor Angiogenesis. Journal of Differential Equations, 259, 5168-5191.
https://doi.org/10.1016/j.jde.2015.06.022
[7] Li, T. and Wang, Z.A. (2009) Nonlinear Stability of Traveling Waves to a Hyperbolic-Parabolic System Modeling Chemotaxis. SIAM Journal on Applied Mathematics, 70, 1522-1541.
https://doi.org/10.1137/09075161X
[8] Peng, H.Y. and Wang, Z.A. (2020) On a Parabolic-Hyperbolic Chemotaxis System with Discontinuous Data: Well-Posedness, Stability and Regularity. Journal of Differential Equations, 268, 4374-4415.
https://doi.org/10.1016/j.jde.2019.10.025
[9] Guo, J., Xiao, J.X., Zhao, H.J. and Zhu, C.J. (2009) Global Solu-tions to Hyperbolic-Parabolic Coupled System with Large Initial Data. Acta Mathematica Scientia. Series B. English Edition, 29, 629-641.
https://doi.org/10.1016/S0252-9602(09)60059-X
[10] Li, D., Pan, R.H. and Zhao, K. (2015) Quantitative Decay of a One-Dimensional Hybrid Chemotaxis Model with Large Data. Nonlinearity, 28, 2181-2210.
https://doi.org/10.1088/0951-7715/28/7/2181
[11] Li, T., Li, T. and Zhao, K. (2011) On a Hyperbolic-Parabolic System Modeling Chemotaxis. Mathematical Methods in the Applied Sciences, 21, 1631-1650.
https://doi.org/10.1142/S0218202511005519
[12] Li, D., Pan, R.H. and Zhao, K. (2012) Global Dynamics of a Hyperbolic-Parabolic Model Arising from Chemotaxis. SIAM Journal on Applied Mathematics, 72, 417-443.
https://doi.org/10.1137/110829453
[13] Li, T. and Wang, Z.A. (2010) Nonlinear Stability of Large Amplitude Viscous Shock Waves of a Generalized Hyperbolic-Parabolic System Arising in Chemotaxis. Mathematical Models and Methods in Applied Sciences, 20, 1967-1998.
https://doi.org/10.1142/S0218202510004830
[14] Martinez, V., Wang, Z.A. and Zhao, K. (2018) Asymptotic and Viscous Stability of Large Amplitude Solutions of a Hyperbolic System Arising from Biology. Indiana University Mathematics Journal, 67, 1383-1424.
https://doi.org/10.1512/iumj.2018.67.7394
[15] Fan, J. and Zhao, K. (2012) Blow up Criterion for a Hyperbol-ic-Parabolic System Arising from Chemotaxis. Journal of Mathematical Analysis and Applications, 394, 687-695.
https://doi.org/10.1016/j.jmaa.2012.05.036
[16] Wang, Z.A. and Hou, Q.Q. (2019) Convergence of Boundary Layers for the Keller-Segel System with Singular Sensitivity in the Half-Plane. Journal de Mathématiques Pures et Appliquées, 130, 251-287.
https://doi.org/10.1016/j.matpur.2019.01.008
[17] Jin, H. and Zou, F. (2023) Nonlinear Stability of Traveling Waves to a Parabolic-Hyperbolic System Modeling Chemotaxis with Periodic Perturbations. Journal of Differential Equations, 352, 23-66.
https://doi.org/10.1016/j.jde.2022.12.033
[18] Jin, H. and Xu, K. (2023) Boundedness of a Chemotaxis-Convection Model Describing Tumor-Induced Angiogenesis. Acta Mathematica Scientia. Series B. English Edition, 43, 156-168.
https://doi.org/10.1007/s10473-023-0110-y
[19] Peng, H.Y., Wang, Z.A. and Zhu, C.J. (2022) Global Weak Solutions and Asymptotics of a Singular PDE-ODE Chemotaxis System with Discontinuous Data. Science China Mathematics, 65, 269-290.
https://doi.org/10.1007/s11425-019-1754-0
[20] Friedman, A. (1969) Partial Differential Equation. Academic Press, New York.