一类非线性反应扩散耦合系统的整体解和爆破
Global Solution and Blow-Up for a Class of Nonlinear Reaction Diffusion Coupled Systems
DOI: 10.12677/AAM.2021.101022, PDF, HTML, XML, 下载: 707  浏览: 874 
作者: 樊佳幸:西南交通大学数学学院,四川 成都
关键词: 非线性抛物方程组Galerkin方法整体解爆破Nonlinear Parabolic Equations Galerkin Method Global Solution Blow up
摘要: 本文研究一类带有扩散项的非线性抛物方程组的初边值问题,借助Sobolev嵌入定理,Gagliardo-Niren- berg不等式,利用Galerkin方法,建立了弱解整体存在和爆破的充分条件,结合伯努利不等式得到了弱解整体存在时的有界性。
Abstract: In this paper, the initial boundary value problem for a class of Nonlinear Parabolic Equations with diffusion term is studied. By means of Sobolev embedding theorem, Gagliardo Nirenberg inequality and Galerkin method, sufficient conditions for the global existence and blow up of weak solutions are established. The boundedness of global existence of weak solutions is obtained by combining Bernoulli inequality.
文章引用:樊佳幸. 一类非线性反应扩散耦合系统的整体解和爆破[J]. 应用数学进展, 2021, 10(1): 189-202. https://doi.org/10.12677/AAM.2021.101022

1. 引言

本文研究如下带有扩散项的非线性抛物方程组的初边值问题

{ u t ( Ω | u | 2 d x ) γ Δ u = | u | p 1 u + α ( u v ) , ( x , t ) Ω × [ 0 , T ] , v t ( Ω | u | 2 d x ) γ Δ v = | v | p 1 v + α ( v u ) , ( x , t ) Ω × [ 0 , T ] , u ( x , 0 ) = u 0 ( x ) , x Ω , v ( x , 0 ) = v 0 ( x ) , x Ω , u ( x , t ) = v ( x , t ) = 0 , ( x , t ) Ω × [ 0 , T ] , (1.1)

其中, ( Ω | u | 2 d x ) γ Δ u , ( Ω | v | 2 d x ) γ Δ v 为扩散项,可知,其为非线性的, Ω R n 为边界充分光滑的有界

区域。常数 α , γ R p > 1

(1.1)可用来描述热传导现象、半导体中的电子与空穴流、液体在多孔介质中的运动规律等问题。当 γ = 0 时,(1.1)是经典的抛物型方程组,与之相关的抛物方程组整体解的存在和爆破已被众多学者进行了广泛的研究。当 γ 0 时,(1.1)即为反应扩散耦合系统,与之相关的抛物方程组的研究也有一些成果。例如,在文献 [1] 中,Carlos等考虑了如下反应扩散耦合系统

{ u t a ( l ( u ) ) Δ u + f ( u v ) = α ( u v ) , ( x , t ) Ω × [ 0 , T ] , v t a ( l ( v ) ) Δ v f ( u v ) = α ( v u ) , ( x , t ) Ω × [ 0 , T ] , u ( x , 0 ) = u 0 ( x ) , x Ω , v ( x , 0 ) = v 0 ( x ) , x Ω , u ( x , t ) = v ( x , t ) = 0 , ( x , t ) Ω × [ 0 , T ] ,

其中, Ω R n 为边界充分光滑的有界区域,常数 α > 0 f : R R a : R R 是利普希茨连续函数,且 a ( s ) m > 0 ( s R ) l : L 2 ( Ω ) R 是线性连续函数。采用Galerkin方法和Aubin紧性定理,研究了弱解的整体存在性、唯一性及指数衰减。

Rui等在文献 [2] 中考虑了有界域上带有更一般的非局部扩散项的非线性反应扩散耦合系统,其扩散项作用于两个线性形式 l 1 , l 2 上,即

{ u t a 1 ( l 1 ( u ) , l 2 ( u ) ) Δ u + λ 1 | u | p 2 u = f 1 , ( x , t ) Ω × [ 0 , T ] , v t a 2 ( l 1 ( v ) , l 2 ( v ) ) Δ v + λ 2 | v | p 2 v = f 2 , ( x , t ) Ω × [ 0 , T ] , u ( x , 0 ) = u 0 ( x ) , x Ω , v ( x , 0 ) = v 0 ( x ) , x Ω , u ( x , t ) = v ( x , t ) = 0 , ( x , t ) Ω × [ 0 , T ] ,

其中, Ω R n 为边界充分光滑的有界区域,常数 λ 1 , λ 2 > 0 p > 1 f 1 , f 2 L 2 ( 0 , T ; L 2 ( Ω ) ) l i : L 2 ( Ω ) R , i = 1 , 2 是线性连续函数,即存在一个函数 g i L 2 ( Ω ) ,使得对所有的 u L 2 ( Ω ) ,都有 l i ( u ) = l g i ( u ) = Ω g i ( x ) u ( x ) d x , i = 1 , 2 成立。 0 < m a i ( s ) M ( s , r R , i = 1 , 2 ) ,且 a i : R 2 R , i = 1 , 2 是整体利普希茨连续函数,同样采用Galerkin方法和Aubin紧性定理,研究了弱解的整体存在性、唯一性以及指数衰减。

事实上,在以往大多是文献中,均假设扩散项系数的有界的,即存在正常数 m , M ,使得

0 < m a ( s ) M < ( s , r R ) ,此时问题总是非退化的。但在刻画扰动传播的有限性等问题的反应扩散方程组时总会出现退化情形,因此有必要研究具有退化性的反应扩散方程组。本文旨在研究一种扩散项

( Ω | u | 2 d x ) γ Δ u , ( Ω | v | 2 d x ) γ Δ v ( γ R ) 的情形,此时扩散项系数 ( Ω | u | 2 d x ) γ , ( Ω | v | 2 d x ) γ ( γ R ) 可能

是零或无穷大。受文献 [3] 的启发,结合Sobolev嵌入定理、Gagliardo-Nirenberg不等式和伯努利不等式,利用Galerkin方法,建立了弱解的全局存在性、有界性及爆破的充分条件。

首先给出弱解的定义如下:

定义1 u , v L ( 0 , T ; H 1 ( Ω ) ) u t v t L 2 ( 0 , T ; L 2 ( Ω ) ) ,对任意的 ω H 0 1 ( Ω ) , t [ 0 , T ] ,都有等式

Ω u t ω + ( Ω | u | 2 d x ) u ω | u | p 1 u ω α ( u v ) ω d x = 0

Ω v t ω + ( Ω | v | 2 d x ) v ω | u | p 1 v ω α ( v u ) ω d x = 0

成立,且

u ( x , 0 ) = u 0 ( x ) , v ( x , 0 ) = v 0 ( x ) , x Ω ,

则称 ( u , v ) 是问题(1.1)的一个弱解。

现给出本文主要结论如下:

定理1 假设 u 0 ( x ) , v 0 ( x ) H 0 1 ( Ω ) α > 0 γ > p 1 2 ,当 n > 2 时, 1 < p n + 2 n 2 ;当 n 2 时, 1 < p <

则问题(1.1)存在一个整体弱解 ( u , v ) ,满足

u , v L ( 0 , T ; L 2 ( Ω ) ) , u , v L ( 0 , T ; L 2 ( Ω ) ) , u t , v t L 2 ( 0 , T ; L 2 ( Ω ) ) ,

n < 2 ,则由初值函数 ( u 0 ( x ) , v 0 ( x ) ) 决定的整体弱解 ( u , v ) 是唯一的。

定理2 假设定理1的条件成立,进一步假设 γ 2 α C 5 1 ,且

u 0 L 2 2 + v 0 L 2 2 C 4 + 2 C 5 γ C 5 ( γ + 1 ) 2 α ,

则有

u L 2 2 + v L 2 2 u 0 L 2 2 + v 0 L 2 2 ,

其中, C 4 = 2 γ P + 1 2 γ + 2 ( γ + 1 p + 1 ) p + 1 2 γ p + 1 ( C 2 ( p + 1 ) ( 2 γ + 2 ) 2 γ p + 1 + C 3 ( p + 1 ) ( 2 γ + 2 ) 2 γ p + 1 ) , C 2 , C 3 , C 5 为嵌入定理中不等式的系数。

定理3 假设 u 0 ( x ) , v 0 ( x ) H 0 1 ( Ω ) α < 0 0 < γ < p 1 2 1 < p < n + 4 n ,则存在一个正常数

T 0 = n n p + 4 γ + 4 c ( p 1 ) ( 4 γ 2 n γ + 4 ) ( u 0 L 2 2 + v 0 L 2 2 ) ( 1 p ) ( 2 γ n γ + 2 ) n n p + 4 γ + 4 , c R + ,

使得当 0 t T 0 时,问题(1.1)存在一个局部弱解 ( u , v ) ,满足

u , v L ( 0 , T ; L 2 ( Ω ) ) , u , v L ( 0 , T ; L 2 ( Ω ) ) ; u t , v t L 2 ( 0 , T ; L 2 ( Ω ) ) .

注1 由 1 < p < n + 4 n 可推出当 n > 2 时, 1 < p n + 2 n 2 ;当 n 2 时, 1 < p <

注2 因为 p < n + 4 n ,故有 n p n 4 4 < 0 < γ ,所以集合 ( n p n 4 4 , p 1 2 ) 不是空集。

定理4 假设定理3的条件成立,且

1 2 γ + 2 u 0 L 2 2 γ + 2 + 1 2 γ + 2 v 0 L 2 2 γ + 2 1 p + 1 u 0 L p + 1 p + 1 1 p + 1 v 0 L p + 1 p + 1 α 2 u 0 v 0 L 2 2 < 0 ,

则问题(1.1)满足初值函数的弱解不是全局的,即在有限时间

2 c ( p 1 ) ( u 0 L 2 2 + v 0 L 2 2 ) p 1 2 < T

内爆破。其中, c = 2 p 4 γ 2 p + 1 ( 2 | Ω | ) 1 p 2

全文将 L p ( Ω ) H s ( Ω ) 分别简记为 L p 和, H s ,记 ( u , v ) = Ω u v d x ,c表示不依赖于未知函数的正常数,在不同地方不一定相同。

2. 准备工作

本节主要给出结论证明过程中需要的三个引理。

引理1 [4] (Sobolev嵌入定理)对任意函数 u 0 ( x ) H 0 1 ( Ω ) ,有

u L p + 1 C 1 u L 2 ,

其中,当 n > 2 时, 0 p n + 2 n 2 ;当 n 2 时, 0 p < ,常数 C 1 仅依赖于 n , p , Ω 。当 n < 2 时, H 0 1 ( Ω )

中的函数的连续的,且

s u p | u | C 2 u L 2 ,

常数 C 2 仅依赖于 n , Ω

引理2 [4] (Gagliardo-Nirenberg不等式)对任意函数 u ( x ) W 0 1 , p ( Ω ) p 1 r 1 ,有如下不等式成立,

u L q C 3 u L p θ u L r 1 θ ,

其中

θ = ( 1 r 1 q ) ( 1 n 1 p + 1 r ) 1 ,

1) 当 p n = 1 时, r q

2) 当 n > 1 p < n 时,若 r n p n p ,则 r q n p n p ;若 r n p n p ,则 n p n p q r

3) 当 p = n > 1 时, r q

4) 当 p > n > 1 时, r q

常数 C 3 仅依赖于 n , p , q , r

引理3 [5] (伯努利不等式) 若 0 x > 1 α R ,则有

( 1 + x ) α > 1 + α x ,当 α < 0 α > 1 时;

( 1 + x ) α < 1 + α x ,当 0 < α < 1 时。

3. 当 γ > p 1 2 时,解的全局性和有界性

证明定理1:利用Galerkin方法。取空间 H 0 1 ( Ω ) 的一组基 { ω i } i = 1 ,且 ω i Δ 算子的狄利克雷问题的特征函数:

{ Δ ω i = λ ω i ω i | Ω = 0 (3.1)

将(3.1)两边同时与 ω i 做内积,可知特征值 λ i 0 , i N * ,将 ω i 标准化,使得 ω i L 2 = 1 。接下来,寻找(1.1)具有以下形式的解

u m ( x , t ) = i = 1 m θ i m ( t ) ω i ( x ) ,

v m ( x , t ) = i = 1 m ϕ i m ( t ) ω i ( x ) ,

其中 θ i m ( t ) , ϕ i m ( t ) C 1 ( [ 0 , T ] ) ,且由以下常微分方程决定

( u m t , ω j ) + ( ( Ω | u m | 2 d x ) γ u m , ω j ) = ( | u m | p 1 u m , ω j ) + ( α ( u m v m ) , ω j ) , 1 j m , (3.2)

( v m t , ω j ) + ( ( Ω | v m | 2 d x ) γ v m , ω j ) = ( | v m | p 1 v m , ω j ) + ( α ( v m u m ) , ω j ) , 1 j m , (3.3)

且初值条件满足

u m ( 0 ) = u 0 m u 0 m = i = 1 m α i m ( t ) ω i ( x ) u 0 H 0 1 ( Ω ) 中, m (3.4)

v m ( 0 ) = v 0 m v 0 m = i = 1 m β i m ( t ) ω i ( x ) v 0 H 0 1 ( Ω ) 中, m (3.5)

因此,由常微分方程中Picard迭代法,(3.2)~(3.5)在某一区间 [ 0 , t m ] , 0 < t m < T 中存在一个局部解 ( u m ( t ) , v m ( t ) ) 。下面证明,对任意的 T > 0 ,这样的解可以通过先验估计延拓到整个区间 [ 0 , T ] 上。

在(3.2)的两边同时乘以 θ i m ,对i从1到m求和,得到

0 t u m t L 2 2 d τ + 1 2 0 t u m L 2 2 γ d d τ u m L 2 2 d τ = 1 p + 1 u m L p + 1 p + 1 1 p + 1 u 0 m L p + 1 p + 1 + α 2 u m L 2 2 α 2 u 0 L 2 2 + α 0 t Ω v m u m t d x d τ , (3.6)

f ( t ) = u m L 2 2 ,则(3.6)左边第二项

1 2 0 t u m L 2 2 γ d d τ u m L 2 2 d τ = 1 2 0 t f ( τ ) γ d f ( τ ) = 1 2 γ + 2 f ( t ) γ + 1 1 2 γ + 2 f ( 0 ) γ + 1 ,

换回原函数名,即

1 2 0 t u m L 2 2 γ d d τ u m L 2 2 d τ = 1 2 γ + 2 u m L 2 2 γ + 2 1 2 γ + 2 u 0 m L 2 2 γ + 2 , (3.7)

将(3.7)代入(3.6)整理可得

0 t u m t L 2 2 d τ + 1 2 γ + 2 u m L 2 2 γ + 2 1 p + 1 u m L p + 1 p + 1 α 2 u m L 2 2 = 1 2 γ + 2 u 0 m L 2 2 γ + 2 1 p + 1 u 0 m L p + 1 p + 1 α 2 u 0 m L 2 2 α 0 t Ω v m u m t d x d τ , (3.8)

同理,在(3.3)的两边同时乘以 ϕ i m ,可以得到

0 t v m t L 2 2 d τ + 1 2 γ + 2 v m L 2 2 γ + 2 1 p + 1 v m L p + 1 p + 1 α 2 v m L 2 2 = 1 2 γ + 2 v 0 m L 2 2 γ + 2 1 p + 1 v 0 m L p + 1 p + 1 α 2 v 0 m L 2 2 α 0 t Ω u m v m t d x d τ , (3.9)

将(3.8)和(3.9)左右两边相加有

0 t u m t L 2 2 d τ + 0 t v m t L 2 2 d τ + 1 2 γ + 2 u m L 2 2 γ + 2 + 1 2 γ + 2 v m L 2 2 γ + 2 1 p + 1 u m L p + 1 p + 1 1 p + 1 v m L p + 1 p + 1 α 2 u m L 2 2 α 2 v m L 2 2 = c α 0 t Ω ( v m u m t + u m v m t ) d x d τ = c α 0 t Ω ( v m u m ) t d x d τ = c ( α Ω v m u m d x α Ω v 0 m u 0 m d x )

整理化简即

0 t u m t L 2 2 d τ + 0 t v m t L 2 2 d τ + 1 2 γ + 2 u m L 2 2 γ + 2 + 1 2 γ + 2 v m L 2 2 γ + 2 = 1 p + 1 u m L p + 1 p + 1 + 1 p + 1 v m L p + 1 p + 1 + α 2 u m L 2 2 + α 2 v m L 2 2 α Ω v m u m d x + c = 1 p + 1 u m L p + 1 p + 1 + 1 p + 1 v m L p + 1 p + 1 + α 2 Ω | u m v m | 2 d x + c 1 p + 1 u m L p + 1 p + 1 + 1 p + 1 v m L p + 1 p + 1 + α 2 u m L 2 2 + α 2 v m L 2 2 + c , (3.10)

由于 γ > p 1 2 ,当 n > 2 时, 1 < p 2 n 2 ;当 n 2 时, 1 < p < ,利用引理1和Young不等式,(3.10)

可化为

0 t u m t L 2 2 d τ + 0 t v m t L 2 2 d τ + 1 2 γ + 2 u m L 2 2 γ + 2 + 1 2 γ + 2 v m L 2 2 γ + 2 c + c ( u m L 2 2 γ + 2 ) p + 1 2 γ + 2 + c ( v m L 2 2 γ + 2 ) p + 1 2 γ + 2 + α 2 ( u m L 2 2 γ + 2 ) 2 2 γ + 2 + α 2 ( v m L 2 2 γ + 2 ) 2 2 γ + 2 c + 1 8 γ + 8 u m L 2 2 γ + 2 + 1 8 γ + 8 v m L 2 2 γ + 2 + 1 8 γ + 8 u m L 2 2 γ + 2 + 1 8 γ + 8 v m L 2 2 γ + 2 ,

0 t u m t L 2 2 d τ + 0 t v m t L 2 2 d τ + 1 4 γ + 4 u m L 2 2 γ + 2 + 1 4 γ + 4 v m L 2 2 γ + 2 c

由此可以得到

0 t u m t L 2 2 d τ c , 0 t v m t L 2 2 d τ c , (3.11)

u m L 2 2 γ + 2 c , v m L 2 2 γ + 2 c , (3.12)

由(3.12)和Poincaré不等式,

u m L 2 2 c u m L 2 2 c , (3.13)

v m L 2 2 c v m L 2 2 c . (3.14)

因此,由先验估计(3.11)~(3.14)可知,存在一个函数 u , v { u m } , { v m } 的子序列,不妨仍记 { u m } , { v m } ,使得

u m u , v m v 弱*收敛于 L ( 0 , T ; L 2 ( Ω ) ) 中, (3.15)

u m u , v m v 弱*收敛于 L ( 0 , T ; L 2 ( Ω ) ) 中, (3.16)

u m t u t , v m t v t 弱收敛于 L 2 ( 0 , T ; L 2 ( Ω ) ) 中。 (3.17)

一方面,由Aubin紧性定理, W = { Φ | Φ L ( 0 , T ; H 0 1 ( Ω ) ) , Φ t L 2 ( 0 , T ; L 2 ( Ω ) ) } 紧嵌入到 L ( 0 , T ; H 0 1 ( Ω ) ) 中,根据(3.15)~(3.17)可以得到

u m u L ( 0 , T ; H 0 1 ( Ω ) ) 中, (3.18)

于是,存在 { u m } 的子序列,不妨仍记为 { u m } ,使得

u m u 几乎处处收敛于 Ω × [ 0 , T ]

因为 s | s | p 1 s 是一个连续函数,所以

| u m | p 1 u m | u | p 1 u 几乎处处收敛于 Ω × [ 0 , T ]

u m L ( 0 , T ; H 0 1 ( Ω ) ) 中的一致有界性和引理1, | u m | p 1 u m L p + 1 p ( 0 , T ; L p + 1 p ( Ω ) ) 上一致有界,则

| u m | p 1 u m | u | p 1 u 弱收敛于 L p + 1 p ( 0 , T ; L p + 1 p ( Ω ) )

另一方面,由(3.18)可知

u m L 2 u L 2 L [ 0 , T ]

为了讨论方便,记 a ( u ) = ( Ω | u | 2 d x ) γ ,由于a是连续的,可以得到

a ( u m ) a ( u ) L [ 0 , T ]

固定(3.2)中的j,由于 p + 1 p < 2 ,且每一项均在 L p + 1 p [ 0 , T ] 中弱收敛,同理可讨论序列 { v m } 的收敛性。

于是令 m ,可以得到

Ω ( u t ω j + a ( u ) u ω j | u | p 1 u ω j α ( u v ) ω j ) d x = 0 , j = 1 , 2 , .

根据基 { ω i } i = 1 的稠密性,

Ω ( u t ω + a ( u ) u ω | u | p 1 u ω α ( u v ) ω ) d x = 0 , ω H 0 1 ( Ω ) .

由(3.4),(3.16),(3.17)知,在 L 2 ( Ω ) 中, u m ( 0 ) u ( 0 ) u m ( 0 ) u 0 ,根据极限的唯一性, u ( 0 ) = u 0 。以类似的方式讨论有 v ( 0 ) = v 0 。因此 ( u , v ) 是问题(1.1)的整体弱解。

唯一性的证明如下:

( u 1 , v 1 ) ( u 2 , v 2 ) 是(1.1)的两个满足初值条件的整体弱解,且 ( u 1 u 2 ) ( x , 0 ) = 0 ( v 1 v 2 ) ( x , 0 ) = 0 ,则有

d d t u 1 ( Ω | u 1 | 2 d x ) γ Δ u 1 = | u 1 | p 1 u 1 + α ( u 1 v 1 ) , (3.19)

d d t v 1 ( Ω | v 1 | 2 d x ) γ Δ v 1 = | v 1 | p 1 v 1 + α ( v 1 u 1 ) , (3.20)

d d t u 2 ( Ω | u 2 | 2 d x ) γ Δ u 2 = | u 2 | p 1 u 2 + α ( u 2 v 2 ) , (3.21)

d d t v 2 ( Ω | v 2 | 2 d x ) γ Δ v 2 = | v 2 | p 1 v 2 + α ( v 2 u 2 ) , (3.22)

进一步的,由(3.19)~(3.22)有

d d t ( u 1 u 2 , h 1 ) + ( Ω | u 1 | 2 d x ) γ Ω u 1 h 1 d x ( Ω | u 2 | 2 d x ) γ Ω u 2 h 1 d x = ( | u 1 | p 1 u 1 | u 2 | p 1 u 2 , h 1 ) + α ( ( u 1 v 1 ) ( u 2 v 2 ) , h 1 ) ,

d d t ( v 1 v 2 , h 2 ) + ( Ω | v 1 | 2 d x ) γ Ω v 1 h 2 d x ( Ω | v 2 | 2 d x ) γ Ω v 2 h 2 d x = ( | v 1 | p 1 v 1 | v 2 | p 1 v 2 , h 2 ) + α ( ( v 1 u 1 ) ( v 2 u 2 ) , h 2 ) ,

接下来,取 h 1 = u 1 u 2 , h 2 = v 1 v 2 可得

1 2 d d t u 1 u 2 L 2 2 + ( Ω | u 1 | 2 d x ) γ ( u 1 , ( u 1 u 2 ) ) ( Ω | u 2 | 2 d x ) γ ( u 2 , ( u 1 u 2 ) ) = ( | u 1 | p 1 u 1 | u 2 | p 1 u 2 , u 1 u 2 ) + α ( u 1 v 1 , u 1 u 2 ) α ( u 2 v 2 , u 1 u 2 ) c Ω sup x Ω ( | u 1 | p 1 , | u 2 | p 1 ) ( u 1 u 2 ) 2 d x + α u 1 u 2 L 2 2 α ( v 1 v 2 , u 1 u 2 ) c Ω ( | u 1 | p 1 + | u 2 | p 1 ) ( u 1 u 2 ) 2 d x + α u 1 u 2 L 2 2 + α 2 v 1 v 2 L 2 2 + α 2 u 1 u 2 L 2 2 c Ω ( u 1 u 2 ) 2 d x sup x Ω ( | u 1 | p 1 + | u 2 | p 1 ) + 3 α 2 u 1 u 2 L 2 2 + α 2 v 1 v 2 L 2 2 . (3.23)

同理,

1 2 d d t v 1 v 2 L 2 2 + ( Ω | v 1 | 2 d x ) γ ( v 1 , ( v 1 v 2 ) ) ( Ω | v 2 | 2 d x ) γ ( v 2 , ( v 1 v 2 ) ) c Ω ( v 1 v 2 ) 2 d x sup x Ω ( | v 1 | p 1 + | v 2 | p 1 ) + 3 α 2 v 1 v 2 L 2 2 + α 2 u 1 u 2 L 2 2 . (3.24)

而(3.23)左边第二项和第三项利用Hölder不等式可化简为

( ( Ω | u 1 | 2 d x ) γ u 1 ( Ω | u 2 | 2 d x ) γ u 2 , u 1 u 2 ) = u 1 L 2 2 γ + 2 + u 2 L 2 2 γ + 2 u 1 L 2 2 γ Ω u 1 u 2 d x u 2 L 2 2 γ Ω u 1 u 2 d x u 1 L 2 2 γ + 2 + u 2 L 2 2 γ + 2 u 1 L 2 2 γ + 1 u 2 L 2 u 2 L 2 2 γ + 1 u 2 L 2 = ( u 1 L 2 2 γ + 1 u 2 L 2 2 γ + 1 ) ( u 1 L 2 u 2 L 2 ) > 0.

同理,(3.24)左边第二项和第三项

( ( Ω | v 1 | 2 d x ) γ v 1 ( Ω | v 2 | 2 d x ) γ v 2 , v 1 v 2 ) > 0.

于是,将(3.23)和(3.24)左右两边分别相加,由引理1可以得到

1 2 ( d d t u 1 u 2 L 2 2 + v 1 v 2 L 2 2 ) c Ω ( u 1 u 2 ) 2 d x sup x Ω ( | u 1 | p 1 + | u 2 | p 1 ) + c Ω ( v 1 v 2 ) 2 d x sup x Ω ( | v 1 | p 1 + | v 2 | p 1 ) + 2 ( u 1 u 2 L 2 2 + v 1 v 2 L 2 2 ) c ( u 1 u 2 L 2 2 + v 1 v 2 L 2 2 ) .

所以,根据Growall不等式,

u 1 u 2 L 2 2 + v 1 v 2 L 2 2 = 0.

证明定理2:将问题(1.1)中两个方程分别与 u , v 做内积可得

1 2 d d t u L 2 2 + u L 2 2 γ + 2 = u L p + 1 p + 1 + α u L 2 2 α Ω u v d x , (3.25)

1 2 d d t v L 2 2 + v L 2 2 γ + 2 = v L p + 1 p + 1 + α v L 2 2 α Ω u v d x , (3.26)

将(3.25)与(3.26)左右两边相加有

1 2 d d t u L 2 2 + 1 2 d d t v L 2 2 + u L 2 2 γ + 2 + v L 2 2 γ + 2 = u L p + 1 p + 1 + v L p + 1 p + 1 + α u v L 2 2 , (3.27)

因为 γ > p 1 2 ,且由引理1和Young不等式得到

1 2 d d t u L 2 2 + 1 2 d d t v L 2 2 + u L 2 2 γ + 2 + v L 2 2 γ + 2 C 2 P + 1 ( u L 2 2 γ + 2 ) p + 1 2 γ + 2 + C 3 P + 1 ( v L 2 2 γ + 2 ) p + 1 2 γ + 2 + α u v L 2 2 1 2 u L 2 2 γ + 2 + 1 2 v L 2 2 γ + 2 + α u v L 2 2 + C 4 ,

d d t ( u L 2 2 + v L 2 2 ) + u L 2 2 γ + 2 + v L 2 2 γ + 2 2 α u v L 2 2 + C 4 2 α ( u L 2 2 + v L 2 2 ) + C 4 . (3.28)

其中, C 4 = 2 γ p + 1 2 γ + 2 ( γ + 1 p + 1 ) p + 1 2 γ p + 1 ( C 2 ( p + 1 ) ( 2 γ + 2 ) 2 γ p + 1 + C 3 ( p + 1 ) ( 2 γ + 2 ) 2 γ p + 1 ) C 2 , C 3 , C 5 为嵌入定理中不等式的系数。由

引理3可知

u L 2 2 γ + 2 = ( 1 + ( u L 2 2 1 ) ) γ + 1 > 1 + ( γ + 1 ) ( u L 2 2 1 ) , (3.29)

v L 2 2 γ + 2 = ( 1 + ( v L 2 2 1 ) ) γ + 1 > 1 + ( γ + 1 ) ( v L 2 2 1 ) , (3.30)

将(3.29)与(3.30)代入(3.28)中,即 y = u L 2 2 + v L 2 2 满足微分不等式

y + ( C 5 ( γ + 1 ) 2 α ) y C 4 + 2 C 5 γ , (3.31)

在(3.31)两边同时乘以 e ( c 1 ( γ + 1 ) 2 α ) t ,则有

e ( C 5 ( γ + 1 ) 2 α ) t y + ( C 5 ( γ + 1 ) 2 α ) e ( C 5 ( γ + 1 ) 2 α ) t y ( C 4 + 2 C 5 γ ) e ( C 5 ( γ + 1 ) 2 α ) t , (3.32)

( e ( C 5 ( γ + 1 ) 2 α ) t y ) ( C 4 + 2 C 5 γ ) e ( C 5 ( γ + 1 ) 2 α ) t ,

两边同时对t积分,可以得到

e ( C 5 ( γ + 1 ) 2 α ) t y y ( 0 ) C 4 + 2 C 5 γ C 5 ( γ + 1 ) 2 α ( e ( C 5 ( γ + 1 ) 2 α ) t 1 ) ,

整理化简则有

y ( y ( 0 ) C 4 + 2 C 5 γ C 5 ( γ + 1 ) 2 α ) e ( C 5 ( γ + 1 ) 2 α ) t + C 4 + 2 C 5 γ C 5 ( γ + 1 ) 2 α .

由于 γ 2 α C 5 1 ,故 e ( C 5 ( γ + 1 ) 2 α ) t 1 ,又因为 C 4 + 2 C 5 γ C 5 ( γ + 1 ) 2 α y ( 0 ) ,所以

y ( y ( 0 ) C 4 + 2 C 5 γ C 5 ( γ + 1 ) 2 α ) e ( C 5 ( γ + 1 ) 2 α ) t + C 4 + 2 C 5 γ C 5 ( γ + 1 ) 2 α y ( 0 ) .

换回原函数名,即结论成立。

注3 由于 p > 1 ,故 γ + 1 > p + 1 2 > 1 ,又 u L 2 2 0 ,所以 u L 2 2 1 1 ,当 u L 2 2 1 > 1 时,由引理3,

(3.29)显然成立;当 u L 2 2 1 = 1 时,代入(3.29)两边验证可知其也成立,当 u L 2 2 1 = 0 时,(3.29)等号成立。(3.30)同理可得。

4. 当 0 < γ < p 1 2 时,解的局部存在性和爆破

证明定理3:依旧采用Galerkin方法, { ω i } i = 1 { u 0 m } , { u m } 与定理1中所述相同。同理,因为 γ > 0 ,由常微分方程中Picard迭代法,(3.2)~(3.5)在某一区间 [ 0 , t m ] , 0 < t m < T 中存在一个局部解 ( u m ( t ) , v m ( t ) )

在(3.2),(3.3)的两边分别同时乘以 θ i m , ϕ i m ,并对对i从1到m求和,得到

1 2 d d t u m L 2 2 + u m L 2 2 γ + 2 = u m L p + 1 p + 1 + α u m L 2 2 α Ω u m v m d x , (4.1)

1 2 d d t v m L 2 2 + v m L 2 2 γ + 2 = v m L p + 1 p + 1 + α v m L 2 2 α Ω u m v m d x , (4.2)

由于 α < 0 ,于是将(4.1),(4.2)左右两边分别相加可以得到

1 2 d d t u m L 2 2 + 1 2 d d t v m L 2 2 + u m L 2 2 γ + 2 + v m L 2 2 γ + 2 = u m L p + 1 p + 1 + v m L p + 1 p + 1 + α u m v m L 2 2 u m L p + 1 p + 1 + v m L p + 1 p + 1 . (4.3)

因为 1 < p < n + 4 n ,所以 1 2 1 n < 1 p + 1 < 1 2 ,于是根据引理2可以得到

| Ω u m p + 1 d x | C 3 ( m L 2 2 γ + 2 ) δ u m L 2 δ , (4.4)

由于 γ > n p n 4 4 ,所以

δ = ( p + 1 ) θ 2 γ + 2 = n p n 4 γ + 4 < 1 ,

又因为 γ < p 1 2 ,故有

δ = ( 1 θ ) ( p + 1 ) = n n p + 2 p + 2 2 > 0 ,

其中, θ 是由 1 p + 1 = ( 1 2 1 n ) θ + 1 2 ( 1 θ ) 唯一确定的[0, 1]中的数,于是根据Young不等式,

| Ω u m p + 1 d x | 1 2 u m L 2 2 γ + 2 + c 2 u m L 2 2 δ , (4.5)

其中 δ = δ 2 ( 1 δ ) > 1 。同理

| Ω v m p + 1 d x | 1 2 v m L 2 2 γ + 2 + c 2 v m L 2 2 δ . (4.6)

将(4.5),(4.6)代入(4.3)中则有

d d t u m L 2 2 + d d t v m L 2 2 c u m L 2 2 δ + c v m L 2 2 δ c ( u m L 2 2 + v m L 2 2 ) δ , (4.7)

于是,微分不等式(4.7)的解 u m L 2 2 + v m L 2 2 可由如下初值问题的解来表示,

d y d t = c y δ , y ( 0 ) = y 0 = u 0 m L 2 2 + v 0 m L 2 2 , (4.8)

如果

t < t = 1 c ( δ 1 ) y 0 δ 1 , (4.9)

则有(4.8)的解是有限的,因此,在区间 0 t T 0 = t 2 内,可以得到

u m ( t ) L 2 2 + v m ( t ) L 2 2 c , (4.10)

于是根据(3.10),(4.4),(4.10)有

u m ( t ) L ( 0 , T 0 ; L 2 ( Ω ) ) c , v m ( t ) L ( 0 , T 0 ; L 2 ( Ω ) ) c , (4.11)

u m t ( t ) L 2 ( 0 , T 0 ; L 2 ( Ω ) ) c , v m t ( t ) L 2 ( 0 , T 0 ; L 2 ( Ω ) ) c , (4.12)

u m ( t ) L ( 0 , T 0 ; L 2 ( Ω ) ) c , v m ( t ) L ( 0 , T 0 ; L 2 ( Ω ) ) c . (4.13)

由(4.11)~(4.13)和Aubin紧性定理,重复定理1的证明过程,可以得到 ( u , v ) 是问题(1.1)在区间 0 t T 0 上的弱解。

下面证明问题(1.1)满足初值函数的弱解不是全局的,即在有限时间爆破。

证明定理4:将问题(1.1)中两个方程分别与 u t , v t 做内积可得

u t L 2 2 + 1 2 γ + 2 d d t u L 2 2 γ + 2 = 1 p + 1 d d t u L p + 1 p + 1 + α 2 d d t u L 2 2 α Ω v u t d x , (4.14)

v t L 2 2 + 1 2 γ + 2 d d t v L 2 2 γ + 2 = 1 p + 1 d d t v L p + 1 p + 1 + α 2 d d t v L 2 2 α Ω u v t d x , (4.15)

将(4.14)与(4.15)左右两边分别相加有

u t L 2 2 + v t L 2 2 + 1 2 γ + 2 d d t u L 2 2 γ + 2 + 1 2 γ + 2 d d t v L 2 2 γ + 2 = 1 p + 1 d d t u L p + 1 p + 1 + 1 p + 1 d d t v L p + 1 p + 1 + α 2 d d t u L 2 2 + α 2 d d t v L 2 2 α d d t Ω u v d x ,

u t L 2 2 + v t L 2 2 + 1 2 γ + 2 d d t u L 2 2 γ + 2 + 1 2 γ + 2 d d t v L 2 2 γ + 2 = 1 p + 1 d d t u L p + 1 p + 1 + 1 p + 1 d d t v L p + 1 p + 1 + α 2 d d t u v L 2 2 ,

于是

d d t ( 1 2 γ + 2 u L 2 2 γ + 2 + 1 2 γ + 2 v L 2 2 γ + 2 1 p + 1 u L p + 1 p + 1 1 p + 1 v L p + 1 p + 1 α 2 u v L 2 2 ) < 0.

两边关于t积分得到

u L 2 2 γ + 2 v L 2 2 γ + 2 2 γ + 2 p + 1 u L p + 1 p + 1 2 γ + 2 p + 1 v L p + 1 p + 1 α ( γ + 1 ) u v L 2 2

u 0 L 2 2 γ + 2 v 0 L 2 2 γ + 2 + 2 γ + 2 p + 1 u 0 L p + 1 p + 1 + 2 γ + 2 p + 1 v 0 L p + 1 p + 1 + α ( γ + 1 ) u 0 v 0 L 2 2 . (4.16)

F ( t ) = u L 2 2 + v L 2 2 ,

E ( t ) = 1 2 γ + 2 u L 2 2 γ + 2 + 1 2 γ + 2 v L 2 2 γ + 2 1 p + 1 u L p + 1 p + 1 1 p + 1 v L p + 1 p + 1 α 2 u v L 2 2 ,

则由(4.16)可以得到

F ( t ) = 2 ( u , u t ) + 2 ( v , v t ) = 2 ( u , ( Ω | u | 2 d x ) γ Δ u + | u | p 1 u + α ( u v ) ) + 2 ( v , ( Ω | v | 2 d x ) γ Δ v + | v | p 1 v + α ( u v ) ) = 2 ( u L 2 2 γ + 2 v L 2 2 γ + 2 + u L p + 1 p + 1 + v L p + 1 p + 1 + α u v L 2 2 ) ( 4 γ + 4 ) E ( 0 ) + 2 p 4 γ 2 p + 1 u L p + 1 p + 1 + 2 p 4 γ 2 p + 1 v L p + 1 p + 1 2 α γ u v L 2 2 , (4.17)

由Hölder不等式有

u L p + 1 p + 1 | Ω | 1 p 2 ( u L 2 2 ) p + 1 2 . (4.18)

因为 α < 0 0 < γ < p 1 2 E ( 0 ) < 0 ,将(4.18)代入(4.17),根据不等式

| a + b | n 2 n 1 ( a n + b n )

可以得到

F ( t ) c F ( t ) p + 1 2 ,

F ( t ) p 1 2 1 F ( 0 ) p 1 2 p 1 2 c t ,

其中, c = 2 p 4 γ 2 p + 1 ( 2 | Ω | ) 1 p 2 。因此,当 t 2 c ( p 1 ) ( u 0 L 2 2 + v 0 L 2 2 ) p 1 2 < T 时, ( u L 2 2 + v L 2 2 ) p 1 2

故结论成立。

参考文献

[1] Raposo, C.A., et al. (2008) Solution and Asymptotic Behaviour for a Nonlocal Coupled System of Reaction-Diffusion. Acta Applicandae Mathematicae, 102, 37-56.
https://doi.org/10.1007/s10440-008-9207-5
[2] Rui, M.P. (2016) A Reaction-Diffusion Model for the Non-Local Coupled System: Existence, Uniqueness, Long-Time Behaviour and Localization Properties of Solutions. IMA Journal Applied Mathematics, 81, 344-364.
https://doi.org/10.1093/imamat/hxv041
[3] Tsutsumi, M. (1972) Existence and Nonexistence of Global Solutions for Nonlinear Parabolic Equations. Publicational of the Research Institute for Mathematical Sciences, 8, 211-229.
https://doi.org/10.2977/prims/1195193108
[4] Ladyzhenskaya, O.A., Solonnikov, V.A. and Uralceva, N.N. (1968) Linear and Quasi-Linear Equations of Parabolic Type. American Mathematical Society.
https://doi.org/10.1090/mmono/023
[5] 葛健芽. Bernoulli不等式的几个标记[J]. 浙江师范大学学报(自然科学版), 2003, 26(2): 119-122.