本文主要讨论在强拓扑空间
罗双利,张江卫,李亚楠
长沙理工大学数学与统计学院,湖南 长沙
收稿日期:2018年10月1日;录用日期:2018年10月15日;发布日期:2018年10月22日
本文主要讨论在强拓扑空间 H 0 1 ( Ω ) ∩ H 2 ( Ω ) × L μ 2 ( ℝ + ; H 0 1 ( Ω ) ∩ H 2 ( Ω ) ) 具有衰退记忆的非经典扩散方程解的存在性、唯一性和解对初值的连续依赖性。值得注意的是非线性项满足任意阶多项式指数增长条件。
关键词 :衰退记忆,非经典扩散方程,强拓扑空间,适定性
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
在本文中,我们研究如下具有衰退记忆的非经典扩散方程在强拓扑空间的适定性:
{ u t − Δ u t − Δ u − ∫ 0 ∞ k ( s ) Δ u ( t − s ) d s + f ( u ) = g ( x , t ) ∈ Ω × ( 0 , ∞ ) u ( x , t ) | ∂ Ω = 0 x ∈ ∂ Ω , t ∈ ℝ u ( x , 0 ) = u 0 ( x , 0 ) x ∈ Ω . (1.1)
其中 Ω ⊂ ℝ 3 是具有适当光滑边界的有界域。外力项 g ( x ) ∈ L 2 ( Ω ) 。对于非线性项, f ∈ C 1 ( ℝ , ℝ ) , f ( 0 ) = 0 且满足
γ 1 | s | p − β 1 ≤ f ( s ) s ≤ γ 2 | s | p + β 2 , s ∈ ℝ , p ≥ 2 (1.2)
f ′ ( s ) ≥ − l (1.3)
其中 γ i , β i ( i = 1 , 2 ) ,l为正常数。令
F ( s ) = ∫ 0 s f ( y ) d y
对于记忆核 μ ( s ) (参见文献 [
令 k ( s ) = ∫ s ∞ μ ( s ) d s , μ ( s ) = − k ′ ( s ) , μ ∈ C 1 ( ℝ + ) ∩ L 1 ( ℝ + ) , μ ( s ) ≥ 0 , μ ′ ( s ) ≤ 0 ,对于所有的 s ∈ ℝ + ,
m 0 = ∫ 0 ∞ μ ( s ) d s < ∞ (1.4)
μ ′ ( s ) + δ μ ( s ) ≤ 0 (1.5)
对于所有的 s ∈ ℝ + , δ 为正常数。
我们引入表示位移历史变量,定义为
η t = η t ( x , s ) = ∫ 0 s u ( x , t − r ) d r , s ∈ ℝ + (1.6)
η s t ( x , s ) = u ( x , t − s ) , η t t ( x , s ) = u ( x , t ) − η s t ( x , s ) (1.7)
那么方程(1.1)可转化为
{ u t − Δ u t − Δ u − ∫ 0 ∞ μ ( s ) Δ η t ( s ) d s + f ( u ) = g , η t t = − η s t + u , (1.8)
初边值条件为
该方程作为通常的非经典扩散方程的推广情形来源并被运用于流体力学、固体力学和热传导理论领域 [
u t ( t ) − Δ u ( t ) = f ( u ) + g (1.10)
u t ( t ) − Δ u ( t ) − Δ u t ( t ) = f ( u ) + g (1.11)
在反应扩散过程中如果我们考虑传播媒质的黏弹性,就需要在方程(1.11)中增加衰退记忆项,产生的模型即为我们将要研究的具有衰退记忆的非经典扩散方程。
u t ( t ) − Δ u ( t ) − Δ u t ( t ) − ∫ 0 ∞ k ( s ) Δ u ( t − s ) d s + f ( u ( t ) ) = g ( t ) (1.12)
据我们了解,对具有衰退记忆的非经典扩散方程解的长时间行为的研究主要集中在弱拓扑空间上。而在强拓扑空间 H 0 1 ( Ω ) ∩ H 2 ( Ω ) × L μ 2 ( ℝ + ; H 0 1 ( Ω ) ∩ H 2 ( Ω ) ) 解的动力学行为考虑的人特别少。这里必须指出的是首先对非线性项假设仅为临界指数增长而不是任意阶指数,其次对于外力项假设具有较强的正则性,即 g ( x ) ∈ L 2 ( Ω ) 所以本文主要讨论在强拓扑空间下,非经典扩散方程强解的适定性。
u t ( t ) − Δ u t ( t ) − Δ u ( t ) − ∫ 0 ∞ k ( s ) Δ u ( t − s ) d s + f ( u ( t ) ) = g (1.13)
为了方便,常用的一些记号和函数空间进行约定。在本文中用C记任意正常数,在不同的行中C所表示的正常数可能不同,甚至同一行中也可表示不同的正常数。记 | · | 2 为 L 2 ( Ω ) 的模,其内积记为 ( · , · ) , ‖ · ‖ 0 = | ∇ • | 2 为 H 0 1 ( Ω ) 的范数, ‖ · ‖ 1 = | Δ · | 2 为 H 0 1 ( Ω ) ∩ H 2 ( Ω ) 的范数。令X为Hilbert空间其内积与范数分别记为
〈 φ , ψ 〉 μ , ν r = ∫ 0 ∞ μ ( s ) 〈 φ ( s ) , ψ ( s ) 〉 ν r d s
‖ φ ‖ μ , ν r 2 = ∫ 0 ∞ μ ( s ) ‖ φ ( s ) ‖ ν r 2 d s
我们定义一族Hilbert空间: M r = ν r × L μ 2 ( ℝ + ; ν r ) ,并且赋予范数
‖ z ‖ M r 2 = ‖ ( u , η t ) ‖ M r 2 = 1 2 ( ‖ u ‖ 0 2 + α ‖ u ‖ 1 2 + ‖ η t ‖ μ , ν r 2 )
引理1.1:根据文献 [
〈 η t , η s t 〉 μ , ν ≥ δ 2 ‖ η t ‖ μ , ν 2 (1.14)
成立。
若 z 0 ∈ M 2 ,显然 z 0 ∈ M 1 ,于是由 [
引理1.2:设 g ∈ L 2 ( Ω ) ,f满足(1.2)和(1.3)。初始值 z 0 ∈ M 2 ,对任意 T > 0 ,系统(1.8)存在唯一解 z = ( u , η t ) ,且满足以下估计:存在仅依赖于假设(1.2)~(1.5)中参数与 | g | 2 的正常数 ρ 0 和 γ ,对任意的 z 0 ∈ M 2 及 t ∈ [ 0 , T ] ,有
‖ z ‖ M 1 2 ≤ ρ 0 + Q ( ‖ z 0 ‖ M 2 ) e − r t (1.15)
其中Q为严格单增的正函数。
定义2.1:对任意的 T > 0 ,记 I = [ 0 , T ] ,设 Ω ⊂ ℝ 3 是有界光滑区域, f ( u ) 满足(1.2)和(1.3), z 0 = ( u 0 , η 0 ) ∈ M 2 , z = ( u , η t ) 是方程(1.1)的整体强解,如果z满足方程(1.1)并且
u ∈ C ( I ; ν 1 ) , u , u t ∈ L 2 ( I ; ν 2 ) , η t ∈ C ( I ; L μ 2 ( ℝ + ; ν 2 ) ) ;
η t t + η s t ∈ L ∞ ( I ; L μ 2 ( ℝ + ; ν 1 ) ) ∩ L 2 ( I ; L μ 2 ( ℝ + ; ν 2 ) ) 。
并且对于任意的 ω ( x ) ∈ ν 0 , φ ( x , t ) ∈ L μ 0 ( ℝ + ; ν 0 ) 有
( u t , ω ) + ( α Δ u t , ω ) + 〈 Δ η t , ω 〉 μ , ν 0 + ( Δ u , ω ) + ( f ( u ) , ω ) = ( g , ω )
〈 Δ η t t + Δ η s t , φ 〉 μ , ν 0 = ( u , φ )
对于 t ∈ I 几乎处处成立。
引理2.1:(强解的存在性)对于任意 T > 0 , I = [ 0 , T ] ,对任意 t ∈ I , z 0 = ( u 0 , η 0 ) ∈ M 2 ,问题(1.8)有唯一强解
z = ( u ( x , t ) , η t )
满足
证明:设 ω j ( x ) 是特征方程
− Δ ω j = λ j ω j , ω j | ∂ Ω = 0
对应于特征值 λ j 的特征函数,因此 { ω j } j = 1 ∞ 组成 L 2 ( Ω ) 的标准正交基,即
当 Ω ∈ C 2 时,有 ω j ( x ) ∈ ν 2 ,且 { ω j } j = 1 ∞ 是 H 0 1 ( Ω ) 的正交基(但非标准正交基), ν 2 即为 { ω j } j = 1 ∞ 在 H 2 ( Ω ) 中的闭线性扩张。因此 { ω j } j = 1 ∞ 也是 ν 2 的正交基。
接下来我们将构造 L μ 2 ( ℝ + ; ν 2 ) 的一组基 { ς i } i = 1 ∞ 。假设 { l k } k = 1 ∞ 是 L μ 2 ( ℝ + ) 的一组正交基,则选择 l k ω j , k , j = 1 , 2 , ⋯ , ∞ 作为 L μ 2 ( ℝ + ; ν 2 ) 的正交基 { ς i } i = 1 ∞ 。
设 E n 1 , E n 2 分别为由 ω 1 , ω 2 , ⋯ , ω n 和 ς 1 , ς 2 , ⋯ , ς n 生成的线性子空间,即
E n 1 = s p a n { ω 1 , ω 2 , ⋯ , ω n } , E n 2 = s p a n { ς 1 , ς 2 , ⋯ , ς n }
记 P n 为 ν 2 到 E n 1 的正交投影, Q n 为 L μ 2 ( ℝ + ; ν 2 ) 到 E n 2 的正交投影,即对任意的 z = ( u , η t ) ∈ M 1 有
u n ( t ) = P n u = ∑ j = 1 n a j ( t ) ω j , u t n ( t ) = P n u t = ∑ j = 1 n b j ( t ) ω j , η n t ( s ) = Q n η t = ∑ j = 1 n c j ( t ) ς j (s)
这里
a j ( t ) = 〈 u , ω j 〉 ν 2 , b j ( t ) = 〈 u t , ω j 〉 ν 2 = d d t a j ( t ) , b j ( t ) = 〈 η t ( s ) , ς j ( s ) 〉 μ , ε 1 。
定义 z n = ( u n , η n t ) 是方程(1.1)的解,形式如下:
u n ( t ) = ∑ j = 1 n a j ( t ) ω j 和 η n t ( s ) = ∑ j = 1 n c j ( t ) ς j ( s )
由Galerkin方法可知,它满足如下的非线性常微分方程组:
{ u n t − Δ u n t − Δ u n − ∫ 0 ∞ μ ( s ) Δ η n t ( s ) d s + f ( u n ) = g n , η n t t = − η n s t + u n , (2.1)
方程(2.1)两边用 − Δ u n 作用,并在 Ω 上对x积分得
1 2 d d t ‖ z n ‖ M 2 2 + ‖ u n ‖ 1 2 + 〈 η n s t , η n t 〉 μ , ν 2 = + 〈 f ( u n ) , Δ u n 〉 − 〈 g n , Δ u n 〉 (2.2)
根据引理1.1有
〈 η n s t , η n t 〉 μ , ν 2 ≥ δ 2 ‖ η n t ‖ μ , ν 2 2 (2.3)
利用Cauchy不等式,有
− 〈 g n , Δ u n 〉 ≤ 1 2 | g n | 2 2 + 1 2 ‖ u n ‖ 1 2
由(1.3)式,我们有
〈 f ( u n ) , Δ u n 〉 = ∫ Ω f ( u n ) Δ u n = − ∫ Ω f ′ ( u n ) | ∇ u n | 2 ≤ l ‖ u n ‖ 0 2
所以,我们得到
d d t ‖ z n ‖ M 2 2 + ‖ u n ‖ 1 2 + δ ‖ η n t ‖ μ , ν 2 2 ≤ 2 l ‖ u n ‖ 0 2 + | g n | 2 2 (2.4)
在 ( 0 , t ) 上对t积分,有
‖ z n ( t ) ‖ M 2 2 + ∫ 0 t ‖ u n ( s ) ‖ 1 2 d s + δ ∫ 0 t ‖ η n t ( s ) ‖ μ , ν 2 2 d s ≤ ∫ 0 T ( 2 l ‖ u n ( s ) ‖ 0 2 + | g n | 2 2 ) d s + ‖ z n ( 0 ) ‖ M 2 2 ≤ C ( | g | 2 2 + 1 ) + ‖ z n ( 0 ) ‖ M 2 2 (2.5)
因此,可以得到
Δ u n 在 L 2 ( 0 , T ; L 2 ( Ω ) ) 中关于 n , t 是一致有界的,
Δ η n t 在 L 2 ( 0 , T ; L μ 2 ( ℝ + , L 2 ( Ω ) ) ) 中关于 n , t 是一致有界的。
因此,通过Alaoglu弱紧性定理和自反性,存在 { z n = ( u n , η n t ) } n = 1 ∞ 的子列 { z n j = ( u n j , η n j t ) } j = 1 ∞ 使得
Δ u n j ⇀ Δ u 在 L 2 ( 0 , T ; L 2 ( Ω ) ) 是弱收敛,
Δ η n t ⇀ Δ η t 在 L 2 ( 0 , T ; L μ 2 ( ℝ + , L 2 ( Ω ) ) ) 是弱收敛。
为了方便起见,在不至于引起误解的情况下,在接下来的讨论中我们将 { z n = ( u n , η n t ) } n = 1 ∞ 的子列仍然用 { z n = ( u n , η n t ) } n = 1 ∞ 表示。
下面将证明 f ( u n ) 在 L 2 ( Ω ) 中是有界的。由 n = 3 知,Sobolev嵌入 D ( A ) ⊂ L ∞ ( Ω ) ,我们有
| u n | 2 ( p − 1 ) ≤ C λ ‖ u n ‖ 1 ≤ C λ M
| u n | 2 ( p − 1 ) 2 ( p − 1 ) ≤ ( C λ M ) 2 ( p − 1 )
∫ 0 T ∫ Ω | u n ( t ) | 2 ( p − 1 ) d t ≤ ( C λ M ) 2 ( p − 1 ) T (2.6)
其中 C λ 是嵌入常数。
因此,
u n ∈ L 2 ( p − 1 ) ( 0 , T ; L 2 ( p − 1 ) (Ω))
由(1.2)和(2.6),使得
| f ( u n ) | 2 2 ≤ C ( | u n | 2 ( p − 1 ) + 1 ) (2.7)
故
∫ 0 T ∫ Ω | f ( u n ) | 2 d t ≤ C [ ∫ 0 T ∫ Ω | u n | 2 ( p − 1 ) d x d t + m e s Ω T ] ≤ C
f ( u n ) ∈ L 2 ( 0 , T ; L 2 ( Ω ) ) = L 2 (ΩT)
f ( u n ) 在 L 2 ( Ω T ) 关于 n , t 是有界的,
f ( u n ) ⇀ χ 在 L 2 ( Ω T ) 是弱收敛的,这里 Ω T = Ω × [ 0 , T ] 。
由于对任意的 t ∈ [ 0 , T ] , u n → u 在 L 2 ( 0 , T ; D ( A ) ) ,由紧性定理, u n → u ( n → ∞ ) 在 L 2 ( 0 , T ; L 2 ( Ω ) ) 是强收敛的。因此 u n → u ( n → ∞ ) 在 L 2 ( Ω ) 几乎处处收敛,利用f的连续性,有:
并有 f ( u n ) 在 [ 0 , T ] 上关于n一致有界,由Lebesgue逐项积分定理,对任意 φ ( x , t ) ∈ C 0 ( [ 0 , T ] ; L 2 ( Ω ) ) 有
lim n → ∞ ∫ Ω × [ 0 , T ] f ( u n ) φ d x d t = ∫ Ω × [ 0 , T ] f ( u ) φ d x d t
则 f ( u n ) ⇀ f ( u ) 在 L 2 ( Ω T ) 中收敛。由弱极限的唯一性知, χ = f ( u ) 。
接下来估计 u n t 。方程(2.1)两边用 − Δ u n t 作用,得
1 2 d d t ‖ u n ‖ 1 2 + ‖ u n t ‖ 0 2 + ‖ u n t ‖ 1 2 = 〈 Δ η n t , − Δ u n t 〉 μ , ν 2 + 〈 f ( u n ) , Δ u n t 〉 − 〈 g n , Δ u n t 〉 (2.8)
用Holder和Young不等式,有
〈 Δ η n t , − Δ u n t 〉 μ , ν 2 ≤ m 0 ‖ η n t ‖ μ , ν 2 2 + 1 4 ‖ u n t ‖ 1 2
和
〈 f ( u n ) , Δ u n t 〉 ≤ | f ( u n ) | 2 2 + 1 4 ‖ u n t ‖ 1 2
− 〈 g n , Δ u n t 〉 ≤ | g n | 2 2 + 1 4 ‖ u n t ‖ 1 2
综上
两边在 ( 0 , t ) 上对t积分,得
‖ u n ( t ) ‖ 1 2 + 2 ∫ 0 t ‖ u n t ( s ) ‖ 0 2 d s + 1 2 ∫ 0 t ‖ u n t ( s ) ‖ 1 2 d s ≤ C ∫ 0 t ( ‖ η n t ( s ) ‖ μ , ν 2 2 + | f ( u n ( s ) ) | 2 2 + | g n | 2 2 ) d s + ‖ u n ( 0 ) ‖ 1 2
u n t 在 L 2 ( 0 , T ; H 0 1 ( Ω ) ) 中有界,在前面已经证明了 − Δ η n t 在 L 2 ( 0 , T ; L μ 2 ( ℝ + ; L 2 ( Ω ) ) ) 有界, f ( u n ) 在 L 2 ( Ω T ) 有界。因此我们有:
u n t 在 L 2 ( 0 , T ; D ( A ) ) 中一致有界,
Δ u n t 在 L ∞ ( 0 , T ; L 2 ( Ω ) ) 中一致有界。
由Alaoglu弱紧性定理和自反性,知
u n t ⇀ u t 在 L 2 ( 0 , T ; D ( A ) ) 是弱收敛的, Δ u n t ⇀ Δ u t 在 L 2 ( 0 , T ; L 2 ( Ω ) ) 是弱收敛。
为了证明 z = ( u , η t ) 是该方程的一个强解,我们还需要验证初值条件 u ( x , 0 ) = u 0 ( x , 0 ) 于 H 0 1 ( Ω ) ∩ H 2 ( Ω ) × L μ 2 ( ℝ + ; H 0 1 ( Ω ) ∩ H 2 ( Ω ) ) 成立。但基于上述估计,可由标准的方法得到(参见文献 [
设 z 1 = ( u 1 , η 1 t ) , z 2 = ( u 2 , η 2 t ) 是方程(1.1)的两个解,且初值 z 10 , z 20 ∈ M 2 ,令
{ ω t − Δ ω t − Δ ω − ∫ 0 ∞ μ ( s ) Δ ς t ( s ) d s + f ( u 1 ) − f ( u 2 ) = 0 , ς t t = − ς s t + ω , (3.1)
方程(3.1)两边用 − Δ ω 作用,有
1 2 d d t ‖ z ′ ‖ M 2 2 = − ‖ ω ‖ 1 2 − 〈 ς s t , ς t 〉 μ , ν 2 + 〈 f ( u 1 ) − f ( u 2 ) , Δ ω 〉 (3.2)
这里 z ′ = ‖ ω ‖ 0 2 + ‖ ω ‖ 1 2 + ‖ ς t ‖ μ , ν 2 2 .
用(1.14),有:
− 〈 ς s t , ς t 〉 μ , ν 2 ≤ − δ 2 ‖ ς t ‖ μ , ν 2 2
由Holder和Young不等式和(2.7)式,得
〈 f ( u 1 ) − f ( u 2 ) , Δ ω 〉 = ∫ Ω ( f ( u 1 ) − f ( u 2 ) ) Δ ω ≤ C ∫ Ω ( | u 1 | p − 2 + | u 2 | p − 2 ) | ω | | Δ ω | + C ∫ Ω | ω | | Δ ω | ≤ C | u n | 2 ( p − 1 ) p − 2 | Δ ω | 2 2 + C | u 2 | 2 ( p − 1 ) p − 2 | Δ ω | 2 2 + C | Δ ω | 2 2 ≤ C | Δ ω | 2 2
综上,我们可以得到
1 2 d d t ‖ z ′ ‖ Μ 2 2 + 1 2 ‖ ω ‖ 1 2 + δ 2 ‖ ς t ‖ μ , ν 2 2 ≤ C | Δ ω | 2 2
d d t ‖ z ′ ‖ M 2 2 ≤ C ℜ ‖ z ′ ‖ M 2 2
其中 C ℜ 是与初值及常数 ρ 有关而与t无关的常数。由Gronwall引理得
‖ z 1 − z 2 ‖ M 2 2 ≤ e C ℜ T ‖ z 1 ( 0 ) − z 2 ( 0 ) ‖ M 2 2
当且仅当 ω ( 0 ) = 0 ,等号成立。所以证明解的唯一性,和对初值的连续依赖性。
本篇论文,我要感谢我的导师谢永钦教授和我的同门师姐们。从论文选题、论文的撰写、修改到最终完稿,他们都给予了我很大的帮助,再次感谢导师为我所做的一切。
罗双利,张江卫,李亚楠. 带记忆项非经典扩散方程整体强解的存在性Existence of Global Strong Solutions for Nonclassical Diffusion Equations with Memory[J]. 应用数学进展, 2018, 07(10): 1233-1240. https://doi.org/10.12677/AAM.2018.710143