一类强阻尼非线性波动方程解的真空隔离性质
具有非线性阻尼及源项的波动方程解的存在性与爆破
![具有非线性阻尼及源项的波动方程解的存在性与爆破](https://img.taocdn.com/s3/m/a62c5e29b80d6c85ec3a87c24028915f804d8404.png)
【中图分类】O175.29
1 引言
utt+a|ut|m+1ut-Δu=b|u|p+1u为具有非线性阻尼及源项的波动方程.文献[1-2]中,对a=0时波动方程的Dirichlet初边值问题进行了讨论,并对其解在相应空间的存在性与爆破性作了阐述.对方程的Cauchy问题,文献[3-4]就其能量估计进行了深入讨论.文献[5-6]中对型如utt+δut-φ(x)Δu=λ u|u|p-1的非线性波动方程的第一类边值问题解的爆破作了深入讨论.文献[7]中利用能量方法和微分、积分不等式技巧,讨论半线性波动方程具非线性Neumann边界条件的混合问题解的爆破性质.文献[8]中讨论了当a=1、b=1时齐次Dirichlet初边值问题解的存在与爆破问题.本文主要讨论a=1、b=1时,在m和p满足一定条件情况下,如下初边值问题的波动方程解的存在性与爆破:
u(t,x)∈L ,ut(t,x)∈L(0,T;L2(Ω))∩Lm+1(0,T;Lm+1(Ω)).
定理2 假设当N≥3时,当N≤2时,p>1;当1<m<p,且初值足够大,函数q(s)(s∈R1)满足λ Q(s)≤0及λ[q(s)-2Q(s)]≤0,问题(1)的解在有限时间内爆破.其中Q(s)=q(τ)பைடு நூலகம்τ.
具有非线性阻尼及源项的波动方程解的存在性与爆破
李爱萍;杨慧
【摘 要】In this paper,the initial-boundary value problem of wave equation with nonlinear damped and source terms was considered and the existence of the weak solutions for the equation with nonlinear Neumann boundary value was proved by Faedo-Galerkin method. Moreover,the blowing-up conditions of the solution were given.%应用 Faedo-Galerkin方法,证明了一类具有非线性阻尼及源项的波动方程非线性 Neu-mann边界条件初边值问题弱解的存在性,并给出了此问题解的爆破条件。
带有阻尼项的非线性波动方程的精确解
![带有阻尼项的非线性波动方程的精确解](https://img.taocdn.com/s3/m/109caf15964bcf84b9d57bdb.png)
初边 值 问题 , 研 究 了整 体 广 义 解 的存 在 性 和不 存
在性 , 以及 整体 古 典 解 的存 在 性 , 解的 b l o w u p现
象和解 的能 量 衰 减 等.2 0 0 8年 A H MA D 等 在 文献 [ 5 ] 中考 虑 了一个带小 阻尼 的非线性 波动方程
精 确解 , 文献[ 5 ] 化方程 ( 2 ) 为一个 方程 组
m l zm _ 。 =
£ ,
分 方程 的一个 活跃 领域 .正 如人 们所 期待 的那 样 ,
这种 分析 实施 的 函数. 厂 ( U ) 的具 体形式 通 常依 赖 于 实 际情况 或者 依 赖 于对 解 的分 类 的追 求 ( 比如对
带阻 尼 的非线 性 波 动 方程 ( 1 ) 的相 应 解 进 行 了 比
体光滑解的存在唯一性.文献 [ 2— 3 ] 分别利用直
接 相似 约化 法 , 给 出 了具 强 阻 尼项 的 非 线 性 波 动
方 程 的相似 约化 , 获得 了一些 精确 类孤 立波解 .文
献[ 4 ] 考 虑 了一 类 具 有 阻尼 的非 线 性 波 动 方 程 的
本 文 首 先 考 察 了带 阻 尼 的非 线 性 波 动 方 程 ( 2 ) 的行 波解 .对 一般 的 m获 得 了方程 ( 2 ) 行 波解 的隐式 表达 式.讨 论 了带 阻尼 的非 线性 波 动方 程 ( 2 ) 行 波 解 的 极 限. 其 次 借 助 于 分 离 变 量 方 法 卜 将带 阻尼 的非 线 性 波 动 方 程 ( 2 ) 约 化成 非 线性 常微 分方 程组 , 从而 获得 了方 程 ( 2 ) 的一些 显 式精 确整 体光 滑解. 讨论 了这 些解 的极 限 , 并 与不
一类具粘性阻尼项的非线性波动方程的整体吸引子及其维数
![一类具粘性阻尼项的非线性波动方程的整体吸引子及其维数](https://img.taocdn.com/s3/m/dea53341453610661ed9f4f2.png)
一类具粘性阻尼项的非线性波动方程的整体吸引子及其维数本文研究如下的具粘性阻尼项的非线性波动方程初边值问题的解的长时间行为:其中x∈Ω,t∈R<sup>+</sup>,σ(s)=s(?),s≥0,m≥1,Ω是R<sup>N</sup>中具有光滑边界的区域,v是(?)Ω的外法向.g(u)和h(u<sub>t</sub>)是给定的非线性函数,f是自由项.本文分七章:第一章为引言;第二章研究问题(1)-(3)在C(R<sup>+</sup>V<sub>2</sub>)∩C<sup>1</sup>(R<sup>+</sup>;H)中的整体解的存在性和唯一性;第三章研究问题(1)-(3)在相空间X<sub>1</sub>=V<sub>2</sub>×H中整体吸引子的存在性及其维数;第四章研究问题(1)-(3)在C(R<sup>+</sup>;V<sub>2+α</sub>)∩C<sup>1</sup>(R<sup>+</sup>;V<sub>α</sub>)(0<α≤1)空间中解的正则性;第五章研究问题(1)-(3)在相空间X<sub>2</sub>=V<sub>2+α</sub>×V<sub>α</sub>中整体吸引子的存在性及其维数;第六章研究问题(1)-(3)在相空间X<sub>3</sub>=V<sub>3</sub>×V<sub>1</sub>中整体吸引子的存在性及其维数;第七章对抽象条件加以验证并给出具体实例.主要结果如下:定理1假定(H<sub>1</sub>)g:V<sub>2</sub>→V<sub>-2</sub>,其中0<ρ<2,G(s)=(?),1≤m≤(?)(m<∞),(a)<sup>+</sup>=max{0,a},及对任意的,u,v ∈V<sub>2</sub>,||u||<sub>V<sub>2</sub></sub>+||v||<sub>V<sub>2</sub></ sub>≤R,有(H<sub>2</sub>)h=h<sub>1</sub>+h<sub>2</sub>,h<sub>i</sub>:V<sub>1</sub>→V<sub>-1</sub>(i=1,2)且存在常数0<δ<sub>1</sub><1,θ<sub>1</sub>∈(0,1/2),β<sub>1</sub>>0使得(H<sub>3</sub>)f∈V<sub>-1</sub>,(u<sub>0</sub>,u<sub>1</sub>)∈X<sub>1</sub>.则问题(1)-(3)存在唯一解u∈C(R<sup>+</sup>;V<sub>2</sub>)∩C<sup>1</sup>(R<sup>+</sup>;H),且(u,u<sub>t</sub>)在空间X<sub>1</sub>中连续依赖于初值.注1 (H<sub>1</sub>)意味着对任意的η>0,存在常数C<sub>η</sub>及(?)使得注2定理1中的解(u,u<sub>t</sub>)我们用S(t)(u<sub>0</sub>,u<sub>1</sub>)=(u,u<sub>t</sub>)表示.则算子族{S (t)}<sub>t≥0</sub>是空间X<sub>1</sub>中的C<sub>0</sub>-半群.定理2在定理1的假定下,如果存在常数0<δ<sub>2</sub><1/2及σ<sub>1</sub>:0<σ<sub>1</sub><<1使得对任意的v∈V<sub>1</sub>,有(H<sub>4</sub>)(H<sub>5</sub>)f∈V<sub>4σ<sub>1</sub>-1</sub>及对任意的(u,v)∈V<sub>1</sub>,||(u,v)||<sub>X<sub>1</sub></sub>≤R,有则连续半群S(t)(见注2)在X<sub>1</sub>中存在整体吸引子A,A是连通的并有有限的分形维数和Hausdorff维数.定理3在定理1中(H<sub>1</sub>)-(H<sub>2</sub>)成立的条件下,如果(H<sub>6</sub>)映射G(见(4)):V<sub>2</sub>→L<sup>1</sup>且存在常数δ<sub>3</sub>∈(0,1),使得对任意的u∈V<sub>2+α</sub>,v∈V<sub>1+α</sub>,||v||≤R,||u||<sub>V<sub>2</sub></sub>≤R,有(H<sub>7</sub>)f∈V<sub>a-1</sub>,(u<sub>0</sub>,u<sub>1</sub>)∈V<sub>2+α</sub>×V<sub>α</sub>其中0<α≤1.则问题(1)-(3)存在唯一解u∈C(R<sup>+</sup>;V<sub>2+α</sub>)∩C<sup>1</sup>(R<sup>+</sup>;V<sub>α</sub>)且(u,u<sub>t</sub>)在空间X<sub>2</sub>中连续依赖于(u<sub>0</sub>,u<sub>1</sub>).定理4在定理3假定成立的条件下,取0<α<1,如果存在一常数σ<sub>2</sub>:0<σ<sub>2</sub><1-α使得(H<sub>8</sub>)f∈V<sub>-1+α+σ<sub>2</sub></sub>和对任意的(u,v)∈V<sub>2+α</sub>×V<sub>α</sub>,||(u,v)||<sub>X<sub>2</sub></sub>≤R,有则C<sub>0</sub>-半群S(t)(见注2)在X<sub>2</sub>中存在整体吸引子A,A是连通的并有有限的分形维数和Hausdorff维数.定理5在定理3中我们取α=1,m≥2,如果存在δ:0<δ<<1使得任取u∈V<sub>3+δ</sub>,||u||<sub>V<sub>3</sub></sub>≤R,u∈V<sub>1+δ</sub>,||v||<sub>V<sub>1</sub></sub>≤R,都有(H<sub>9</sub>)则注2中定义的C<sub>0</sub>-半群S(t)在X<sub>3</sub>中存在整体吸引子A,A是连通的并有有限的fractal维数和Hausdorff维数.注3由(H<sub>1</sub>)的假定我们可推得m≥2意味着N≤4,特别地,当m=2时N=4.。
保守系统强非线性强迫振动问题的一种解法
![保守系统强非线性强迫振动问题的一种解法](https://img.taocdn.com/s3/m/7790340d4a35eefdc8d376eeaeaad1f346931138.png)
保守系统强非线性强迫振动问题的一种解法
袁镒吾;徐积江
【期刊名称】《力学季刊》
【年(卷),期】2001(22)1
【摘要】本文用新的改进的L-P法求得了一类保守系统强非线性强迫振动问题的一级近似共振周期解,近似解比已有的改进的L-P法的相应解更加准确。
【总页数】4页(P143-146)
【关键词】强非线性;强迫振动;保守系统;共振周期解;L-P法;加权残值法
【作者】袁镒吾;徐积江
【作者单位】中南大学;江西赣州公路分局
【正文语种】中文
【中图分类】O322
【相关文献】
1.一类有阻尼的强非线性自由振动问题的一种解法 [J], 袁镒吾;袁雪辉
2.保守系统中非线性振动问题的数值解法 [J], 陈宜周
3.保守系统弱非线性振动问题的插值摄动解法 [J], 袁镒吾;袁雪辉
4.有阻尼的强非线性Duffing方程强迫振动问题的插值摄动解法 [J], 袁镒吾
5.有正阻尼时强非线性自由振动问题的一种解法 [J], 袁镒吾;曹人清
因版权原因,仅展示原文概要,查看原文内容请购买。
一类具阻尼非线性双曲型方程的初边值问题解的爆破
![一类具阻尼非线性双曲型方程的初边值问题解的爆破](https://img.taocdn.com/s3/m/5e400069ddccda38376baf69.png)
证 明 假设 初边值 问题 ( )~ ( )的解存 在 的最大 时 间是无 限的. 程 ( )的两 边 同乘 以 2 在 上 1 3 方 1 “,
: —
l n
.
2 + √
y日 0 + 矗 0 :() 卢 ()
收 稿 日期 :02— 4—2 21 0 1 基 金 项 目 : 南省 基 础 与 前 沿技 术 研 究 项 目( 0 30 124;13 0 1 13 河 120 4 0 1 12 0 40 9 )
作者简介 : 宋瑞丽(9 8一) , 17 女 河南新野人 , 师, 讲 主要从 事非线性发展方程的研究
非 线性 双 曲型方 程 的三维初 边值 问题
+k V +k + V g V U 1 4 2 4 V ( )=0, t ×( T , ( )∈ , 0, ) () 1
“ :0, V u=0 ( t , ,)∈ a ×( , , 0 ) M , )= ( , ( 0 =U ( , ,)∈ , ( 0 。 ) ,) ) ( t
・
7 6・
河 南工程 学院学报 ( 自然科 学版 )
21 0 2生
使得 当 t 。 , 一 时 日()一 +∞. 中 , 其
2 = Al ±
.
+
2 .
定理 假 设 “ 。∈ n ( , ∈L ( )g o 2 ) , ( )=0 G △u ) ∈ ( )并且 存在 常数 >0使 得 ,( 。
1 整 体 解 不 存 在 定 理
带阻尼项的二阶强迫非线性微分方程解的振动性
![带阻尼项的二阶强迫非线性微分方程解的振动性](https://img.taocdn.com/s3/m/1089a06c1ed9ad51f01df2cf.png)
本 文提 出 了一类 新 的带 阻尼项 的二 阶 强迫非 线
”t ()+q t t 0 () )= (
性 微分 方程
( ( ) ( ( ) k Y( ) ) P tk Y ( ) rt Y t ) ( t ) + ( ) ( t )+
是 振动 的.之后 很 多学 者研 究 了形 式更 为 复杂 的 微 分 方程 的振 动性 质.
关 键 词 :二 阶 微 分 方 程 ;振 动 性 ;变 分 法
中 图分 类 号 :0 7 . 15 1
文献 标志 码 : A
Os i a i n Th o e s f r S c nd. r r No ln a r ur d cl to e r m o e o l O de n i e r Pe t be
W a gl n
.
பைடு நூலகம்
是振 动 的 , 则称 方程 为振 动 的.
P is 于 18 hl … o 9 9年 得 出 当 :
l i m
…
sp u
,
㈤
一
近年 来 , 多 学者 再 次 把 注 意力 投 向了对 非 线 很
÷ ,】 ∞ t s )
Dif r n i lEq to s wih Da pi g f e e ta ua i n t m n
Zh n e g, Zh n p n a gM n a g Li i g
( c o l f ce c , U E S h o o S in e B C A,B in 0 0 4) ej g 1 0 4 i
L 和 R g vh n o 研 究 了非 线 性 方 程 : ooc e k (( ) ( ) + ( )( t )= rt t ) q t_ () 0 厂 的振 动理 论.
带有阻尼项的非线性波动方程的精确解
![带有阻尼项的非线性波动方程的精确解](https://img.taocdn.com/s3/m/34ae159085868762caaedd3383c4bb4cf7ecb7be.png)
带有阻尼项的非线性波动方程的精确解尚亚东【摘要】研究了出现在非线性振动中的一类带阻尼项的非线性波动方程.首先讨论了所论方程的行波解及其极限行为,其次借助于分离变量方法获得了所研究方程的一些显式精确解,讨论了这些解的极限行为.这些解有助于定性或数值分析非线性波动方程解的性态.【期刊名称】《广州大学学报(自然科学版)》【年(卷),期】2013(012)002【总页数】6页(P1-6)【关键词】非线性波动方程;阻尼项;分离变量方法;精确解;整体光滑解【作者】尚亚东【作者单位】广州大学数学与信息科学学院,广东广州510006;广州大学数学与交叉科学广东普通高校重点实验室,广东广州510006【正文语种】中文【中图分类】O175.290 引言对于形如的非线性波动方程的分析,无论从解析角度还是从数值分析看,一直是且还将继续是非线性偏微分方程的一个活跃领域.正如人们所期待的那样,这种分析实施的函数f(u)的具体形式通常依赖于实际情况或者依赖于对解的分类的追求(比如对称分类).特别地,幂函数形式的f,比如f(u)=mum-1,就在许多情况下被考虑.行波形式的解是许多分析中普遍寻求的解.由于实际中粘性耗散的不可避免,经常需要考虑带阻尼项的非线性波动方程[1-4].对于带强粘性阻尼的高阶非线性波动方程,文献[1]考虑了一类具有阻尼项的高阶非线性波动方程的初值问题,证明了一定条件下整体光滑解的存在唯一性.文献[2-3]分别利用直接相似约化法,给出了具强阻尼项的非线性波动方程的相似约化,获得了一些精确类孤立波解.文献[4]考虑了一类具有阻尼的非线性波动方程的初边值问题,研究了整体广义解的存在性和不存在性,以及整体古典解的存在性,解的blow up现象和解的能量衰减等.2008年AHMAD等在文献[5]中考虑了一个带小阻尼的非线性波动方程讨论了所论方程的近似拉格朗日和某种相似约化的不变量.最近苏敬蕊等[6]借助于与部分拉格朗日相关的近似Noether型对称算子,考察了带阻尼的非线性波动方程(2),构造了它的一般形式的近似守恒律.对于带阻尼的非线性波动方程(2)的精确解,文献[5]化方程(2)为一个方程组这个方程组容许有平移对称‘∂t和∂x.这提供了一个能生成波速为c的行波解的组合的Lie点对称X=∂t+c∂x.当 X 有不变量 y=x-ct时,获得了 m=2时的隐式行波解为对于方程(2)在m≠2时的精确行波解以及其他解析形式的显式精确解,就作者所知,还没有文献论及.本文首先考察了带阻尼的非线性波动方程(2)的行波解.对一般的m获得了方程(2)行波解的隐式表达式.讨论了带阻尼的非线性波动方程(2)行波解的极限.其次借助于分离变量方法[7-9]将带阻尼的非线性波动方程(2)约化成非线性常微分方程组,从而获得了方程(2)的一些显式精确整体光滑解.讨论了这些解的极限,并与不带阻尼的非线性波动方程(1)的相应解进行了比较,解释了粘性阻尼的效应.这些解将有助于定性或数值分析带阻尼的非线性波动方程解的性态.1 非线性波动方程的精确行波解考虑非线性波动方程和带阻尼的非线性波动方程(2)的行波解.如所周知,线性波动方程有行波解其中f(ξ)和g(η)为两个任意的二次可微函数.为求方程(5)的行波解,作行波变换ξ=xct,其中c为波的传播速度.这时方程(5)成为方程(6)左右两边关于ξ积分两次得其中k1,k2为任意积分常数.解(9)即为非线性波动方程(5)的隐函数形式的行波解表达式.特别当f(u)=um,m≠1时,可知,非线性波动方程(5)有隐函数形式的行波解这里k1,k2为任意积分常数.如果考虑带阻尼的非线性波动方程(2)的行波解,行波变换使得方程(2)变成方程(11)两边关于ξ积分一次得其中k0为任意积分常数.为了求得方程(2)的行波解,在方程(12)中令k0=0,可得上方程两边同除以u',然后关于ξ积分一次得其中k3为任意积分常数.这正是带阻尼的非线性波动方程(2)的隐函数形式的行波解.明显地,当阻尼项系数ε趋于0时,由带阻尼的非线性波动方程(2)行波解(14)无法得到没有阻尼的非线性波动方程的行波解(10).当m=2时,这里得到的行波解(14)正是文献[5]所得的解.2 非线性波动方程的精确非行波解本节研究带阻尼的非线性波动方程(2)及其对应的无阻尼非线性波动方程的其他形式的显式精确解析解.首先寻求带阻尼的非线性波动方程的和式变量分离解.假设为方程(16)的和式变量分离解,将表达式(17)代入方程(16),则有简单计算整理得到为了获得阻尼非线性波动方程(16)非平凡的分离变量解,要求函数f(x)和g(t)都不恒为常数.方程(19)右端只与自变量x有关,而与自变量t无关.左端为关于自变量t 的一阶线性常微分方程,其系数为x的函数.如果有单变量函数f(x)和g(t)要满足方程(19),必须有(f2(x))xx恒为常数,且f″(x)也恒为常数.分别在不同情况下讨论形如(17)的非平凡解的存在性.情形1 f″(x)=0,但是(f2(x))xx=0.这时方程(19)约化为其解为 g(t)=,这里C1,C2为任意积分常数.由条件f″(x)=0,推知其中a,b是任意待定常数.于是有由条件(f2(x))xx=0,推知a=0.从而f(x)≡常数.所以在此情形下阻尼非线性波动方程(16)有与变量x无关的显式精确解对应的无阻尼非线性波动方程显然有与变量x无关的显式精确解明显地,在阻尼非线性波动方程(16)中让ε趋于0时,方程退化为非线性波动方程(24),但其对应的显式解析解(23)在ε趋于0时并不以非线性波动方程(24)的对应显式解析解(25)为极限.情形2 f″(x)≠0,但是(f2(x))xx=0.这时方程(19)约化为于是有方程(27)左右两端分别为自变量t和x的一元函数.于是f″(x)必须为常数.因此f(x)一定是二次多项式.设其中a,b,c为待定的常数.因此由方程(27),(28),(29)显然有这与方程(19)矛盾.因此在此情形下阻尼非线性波动方程(16)不存在非平凡的和式分离变量解.情形3 f″(x)≠0,但(f2(x))xx≠0.这时不妨设f″(x)=a,a≠0为某个常数.于是有由假设,可知其中b,c为两个任意常数.因此于是,方程(31)左端只与单自变量t有关,而右端却是自变量x的单变量函数.矛盾.说明在此情形下阻尼非线性波动方程(16)也不存在非平凡的和式分离变量解.情形4 f″(x)=0,但(f2(x))xx≠0.这时,方程(19)约化为由条件f″(x)=0,推知其中a,b是任意待定常数.于是有因此,方程(34)变为求得其解为其中C1,C2为任意常数.从而阻尼非线性波动方程(16)有和式变量分离显式精确解这里 C1,C2,a≠0,b 为任意常数.显式精确解析解(39)是阻尼非线性波动方程(16)的整体光滑解,对应于初始值u0(x)显然a=0时,整体光滑解退化为与变量x 无关的精确解(23).对与阻尼非线性波动方程(16)相对应的无阻尼的非线性波动方程(24)做与情形2到情形4相仿的讨论可知,非线性波动方程(24)有和式变量分离显式精确解这里C1,C2,a≠0,b为任意常数.显式精确解析解(40)是非线性波动方程(24)的整体光滑解,对应于初始值u0(x)=ax+b+C2.显然在a=0时,整体光滑解(40)退化为与变量 x无关的精确解(25).在阻尼非线性波动方程(16)中让ε趋于0时,方程退化为非线性波动方程(24),但其对应的显式解析解(39)在ε趋于0时显然并不以非线性波动方程(24)的对应显式解析解(40)为极限.其次,考虑带阻尼的非线性波动方程(2)的乘积形式的变量分离解的显式表达式. 假设方程(2)有乘积形式的变量分离解将表达式(41)代入到方程(2)中,计算得到分离变量得到方程(43)左端为自变量t的单变量函数,而右端却是自变量x的单变量函数,所以只有两端都为常数才能相等.设因此,有非线性常微分方程和为了得到阻尼非线性波动方程(2)的乘积形式分离变量解(41),对于常数λ,需要求解非线性常微分方程(45)和(46).分两种情形来讨论非线性常微分方程(45)和(46)的解.情形1 当λ=0时.这时,非线性常微分方程(45),(46)分别退化为和分别求解线性常微分方程(47)和非线性常微分方程(48),得到精确解和其中 C1,C2,a,b 为任意常数.在此情形下,获得阻尼非线性波动方程(2)的显式精确解析解为对应于初值u0(x)=显式精确解(51)为整体光滑的解析解.进行如同前面一样的分析,可以得到无阻尼的非线性波动方程(24)有相应的乘积形式的变量分离显式解析解其中C1,C2,a,b为任意常数.这个解也是一个整体光滑解析解.同样地,在解(51)中让ε→0,无法其中 n 为待定正整数,而 ai(i=0,1,2,…,m)为待定常数.于是方程(46)左端为一个变量x的nm-2次多项式,而方程(46)右端是变量x的n次多项式.于是必须有解得得到无阻尼的非线性波动方程(24)乘积形式的变量分离显式解析解(52).情形2 当λ≠0时.在此情形下,为了获得阻尼非线性波动方程(2)的乘积形式变量分离解的精确表达式,需要求得非线性常微分方程(45)和(46)的精确解.观察非线性常微分方程(46)的结构,根据多项式函数求导的特点,假设方程(46)有多项式函数形式的解,不妨设当m=2和m=3时,可分别得到n=2和n=1.从而,可假设方程有多项式函数形式解其中a,b,c为待定常数.将表达式(55)代入到方程(54)中,计算得到这里a,c为符号相同的任意常数.因此方程(54),即方程(46)在m=2时有解其中ac>0为任意常数.当m=3时,假设方程有多项式函数形式解其中a,b为待定常数.将表示式(59)代入方程(58),计算得到于是,方程(46)在m=3时有解其中a≠0,b为任意常数.当m=2时,方程(45)对应成为而当m=3,方程(45)对应成为总结前述,可得带阻尼的非线性波动方程(16)的乘积形式的变量分离精确解析解其中ac>0为任意常数,而g(t)满足二阶非线性常微分方程(62).带阻尼的非线性波动方程有乘积形式的变量分离精确解析解其中a≠0,b为任意常数,而g(t)满足二阶非线性常微分方程(63).由以上讨论还可知,无阻尼的非线性波动方程(15)在m=2和m=3时,分别有乘积形式的变量分离显式精确解析解其中ac>0为任意常数,其中a≠0,b为任意常数.解(67)和解(68)都是局部解,在时,这两个解都会出现爆破现象.更一般的,无阻尼的非线性波动方程(15)在m=2和m=3时,分别有乘积形式的变量分离精确解析解其中ac>0为任意常数,k3为任意常数,这里u(x,t)关于t为Weierstrass椭圆函数双周期解.其中a≠0,b,k4为任意常数,这里 u(x,t)关于 t为Weierstrass椭圆函数双周期解.3 结束语非线性振动现象出现于交通科学、声学、材料力学、石油开采等许多科学技术领域,这些现象的数学模型归结为非线性波动方程.由于实际问题中,粘性阻尼的耗散效应不可避免,讨论粘性阻尼对波动的影响机理是非常重要的.反映在数学上,归结为研究阻尼对波动方程解的变化如何影响,阻尼耗散趋于零时,解的极限状态是否与无阻尼时波动方程的解一致.对于出现在非线性振动等应用领域中的非线性波动方程,讨论了一类带阻尼的非线性波动方程和相对应的无阻尼的非线性波动方程精确解析解的存在性.首先获得了这些方程隐函数形式的行波解,接着利用分离变量方法获得了无阻尼非线性波动方程(15)的一些显式精确解析解,这些解既有显式精确整体光滑解,也有在有限时间内发生爆破的显式精确局部解,特别地,获得了关于时间t为双周期函数的Weierstrass椭圆函数解.对于带阻尼的非线性波动方程(2),也找到了一些显式精确解析解,并考察了这些解在ε→0时的极限性态.这些显式精确解对于定性分析非线性波动方程(2)和(15)解的性质以及数值求解将会带来有益帮助.致谢:感谢屈长征教授的热情支持与帮助.参考文献:[1]杨志坚,陈国旺.具有阻尼项的非线性波动方程的初值问题[J].应用数学学报,2000,23(1):45-54.YANG Z J,CHEN G W.Initial value problem for a nonlinear wave equation with damping term[J].Acta Math Appl Sin,2000,23(1):45-54.[2]闫振亚,张鸿庆.具有阻尼项的非线性波动方程的相似约化[J].物理学报,2000,49(11):2113-2117.YAN Z Y,ZHANG H Q.Similarity reductions for a nonlinear wave equation with damping term[J].Acta Phys Sin,2000,49(11):2113-2117.[3] WU H X,FAN T Y.New explicit solutions of the nonlinear wave equations with damping term[J].Appl Math Comput,2007,191(2):457-465.[4] CHEN G W,LU B.The initial boundary value problems for a class of nonlinear wave equations with damping term[J].J Math Anal Appl,2009,351(1):1-15.[5] AHMAD F,KARA A H,BOKHARI A H,et al.On approximate Lagrangians invariants for scaling reductions of a non-linear wave equation with damping[J].Appl Math Comp,2008,206:16-20.[6] SU J R,ZHANG S L,LI J N.Approximate Noether-type symmetries and conservation laws via partial Lagrangians for nonlinear wave equation with damping[J].Commun Theor Phys,2010,53:37-42.[7] EVANS L C.Partial differential equations[M].American Mathematical Society Providence,Rhode Island,Graduate Studies in Mathematics,2002.[8] QU C Z,ZHANG S L,LIU R C.Separation of variables and exact solutions to the quasilinear diffusion equations with nonlinear source [J].Phys D,2000,144:97-123.[9] ESTEVEZ P G,QU C Z,ZHANG S L.Separation of variables of a generalized porous medium equation with nonlinear source[J].J Math Anal Appl,2002,275:44-59.[10]楼森岳,唐晓艳.非线性数学物理方法[M].北京:科学出版社,2006.LOU S Y,TANG X Y.Nonlinear mathematical physics methods[M].Beijing:Science Press,2006.。
一类有正阻尼的强非线性振动问题的解析解
![一类有正阻尼的强非线性振动问题的解析解](https://img.taocdn.com/s3/m/5c2ece651eb91a37f1115c73.png)
吴荣兴,男,宁波大学机械工程与力学学院,博 士研 究生,研 究方 向:非线性振动理论 ( 浙江 宁波 3 5 1) 12 1
第4 期第 6 8页
陈乐平 , 一类有正阻尼 的强非线性振动 问题 的解析解 等:
20 年 l 09 2月
线性解 ,分别为 :
“+8s0+809(3)i+3 。8 =05 (5.)0i13 、8+4 ・1 1 t8 .I+3三 3 ., 1 t . 3+5n 8_ . 。 7i6 130s . _n1  ̄。 . n ・ 2 t 3 3 s
果 的应用 留下 了疑 问。 2 同伦分 析法求解
现在 我们开始 对 ()利用 同伦分析方法 进行求解 ,根据 初始条件 ,可 以选定初 始猜测 解为 , 1
u= ot o Ac s . () 5
根据 ()的性质 ,可 以定义非线性算子 为 , 1
;
[ ’ ,
( 6 )
次和两 次导数 。这 里 当阻尼系数 £≥1 ,方程就变 为强 非线性 和强 阻尼 振动方 程 。 时
袁镒 吾等使用 线化摄 动法 ,权 残法和 插值摄 动法相 结合 的方法 分别求得 了 ()的强非线 性解和 弱非 1
收 稿 日期 : 20 0 9— 0 6— 1 9
作者简介 :陈乐平,男,浙江纺织服装职业技术学院,讲师,在读 博士,研 究方向:机械/ 模具 C AE ( 浙江 宁波 3 5 1) 12 1
摘
要 :阐述利用最新提 出的同伦分析方 法求解有正 阻尼的强非线性振动方程 ,得到 了该方程 的六阶近似解析解 。
计 算结果表明 ,本文得 到的结果与多尺度 法、混合近似法以及数值计算方法得到的结果非常接 近 ,而且 同伦 近似 解提供
阻尼波动方程
![阻尼波动方程](https://img.taocdn.com/s3/m/89b5a3d6112de2bd960590c69ec3d5bbfd0adaea.png)
阻尼波动方程1. 引言阻尼波动方程是描述阻尼介质中波动现象的数学模型。
在物理学和工程学中,阻尼波动方程被广泛应用于描述声波、光波、电磁波等各种波动现象。
本文将详细介绍阻尼波动方程的概念、数学表达以及其在实际应用中的重要性。
2. 阻尼波动方程的概念阻尼波动方程是一种描述阻尼介质中波动传播的偏微分方程。
它考虑了介质中存在的阻尼力对波动传播的影响,使得波动在传播过程中逐渐减弱衰减。
阻尼波动方程的一般形式可以表示为:∂2u ∂t2+2α∂u∂t−c2∂2u∂x2=0其中,u表示波动的振幅,t表示时间,x表示空间坐标,α表示阻尼系数,c表示波速。
3. 阻尼波动方程的解析解阻尼波动方程是一个二阶偏微分方程,可以通过求解特征方程得到解析解。
特征方程的解决方法取决于阻尼系数α的值。
3.1 无阻尼情况当阻尼系数α为零时,阻尼波动方程退化为标准的波动方程:∂2u ∂t2−c2∂2u∂x2=0这是一个齐次线性偏微分方程,可以通过分离变量法或傅里叶变换等方法求解。
3.2 弱阻尼情况当阻尼系数α较小且正值时,阻尼波动方程的解可以近似表示为:u(x,t)=Ae−αt cos(kx−ωt+ϕ)其中,A为振幅,k为波数,ω为角频率,ϕ为相位差。
3.3 强阻尼情况当阻尼系数α较大且正值时,阻尼波动方程的解可以近似表示为:u(x,t)=Ae−αt e−(α2−√α2−ω2)t cos(kx−ωt+ϕ)其中,A为振幅,k为波数,ω为角频率,ϕ为相位差。
4. 阻尼波动方程的应用阻尼波动方程在实际应用中具有重要的意义。
以下是几个典型的应用场景:4.1 声波传播阻尼波动方程可以用于描述声波在介质中的传播过程。
在空气、水、固体等介质中,声波的传播会受到介质的阻尼影响,波动逐渐衰减。
通过求解阻尼波动方程,可以预测声波的传播特性,如衰减程度、频率响应等。
4.2 光波传播阻尼波动方程也可以应用于描述光波在介质中的传播过程。
在光学领域中,介质的吸收和散射会导致光波的衰减。
阻尼波动方程
![阻尼波动方程](https://img.taocdn.com/s3/m/689aeafaba4cf7ec4afe04a1b0717fd5360cb2b9.png)
阻尼波动方程
(原创版)
目录
1.阻尼波动方程的定义与概述
2.阻尼波动方程的求解方法
3.阻尼波动方程的应用领域
正文
阻尼波动方程是物理学中描述阻尼振动系统的偏微分方程。
阻尼振动系统通常由一个振动质点与一个弹性元件以及一个阻尼元件组成。
阻尼波动方程可以用来研究许多实际问题,如弹簧振动、声波传播等。
求解阻尼波动方程的方法有很多,其中最常见的方法是利用拉普拉斯变换和傅里叶变换。
拉普拉斯变换可以将时域上的阻尼波动方程转化为频域上的方程,从而更容易求解。
傅里叶变换则可以将阻尼波动方程从一般坐标变换到特定的坐标系,如傅里叶坐标系,以便于求解。
阻尼波动方程在许多领域都有广泛的应用。
例如,在机械工程中,阻尼波动方程可以用来研究弹簧振动系统的稳定性和振动幅度;在声学中,阻尼波动方程可以用来研究声波在空气中的传播特性;在地球物理学中,阻尼波动方程可以用来研究地震波在地球内部的传播规律。
第1页共1页。
一类非线性波动方程的新的精确解
![一类非线性波动方程的新的精确解](https://img.taocdn.com/s3/m/f5445743767f5acfa1c7cd65.png)
G ( e ) +2 G ( ) +/ z G( ) =0
由方 程 ( 5 ) 易得 , 当 G ( ) ≠0时 ,
收稿 日期 : 2 0 1 3 . 0 1 . 0 7 基金项 目: 国家 自 然科学基金资助项 目( 1 1 1 7 1 0 4 6 , 1 1 0 6 1 0 3 9 )
形为著名 的 L a n d o u — G i n b u r g — H i g g s s 方程 一“ 工 r —m +g U =0 . 主要借助改进的( G / G) 展开法[ 1 0 ] 来求解方程( 1 ) 。首先介绍改进的( G / G) 展开法。
1 预 备 知 识
s i n e 展 开法 … 1, 齐次 平衡 法 [ 2 2, ]D a r b o u r 变 换法 [ 3 3, ]截 断 的 P a i n l e v e 展开法[ , J a c o b i 椭 圆 函数 展 开 法 [ 5, 5 ]t a n h函数
展开方法及 其推 广方法 6, 双 曲正切 函数展 开法_ 7 ] 等 。现在 应用较 为广 泛 的是 王 明亮 提 出的 ( G / G) 展开 法[ 儿 ] , 使用起来简洁 , 高效。
一类带有黏性项的强阻尼耦合波动方程弱解的存在性
![一类带有黏性项的强阻尼耦合波动方程弱解的存在性](https://img.taocdn.com/s3/m/ff4d7029ccbff121dd3683fc.png)
[ 文章 编号] 1 6 7 2 — 2 0 2 7 ( 2 0 1 3 ) 0 3 — 0 0 2 7 - 0 5 [ 中图分 类号]
01 7 4 [ 文献 标识 码] A
0 引 言
本 文考 虑 如下 问题 :
第 1 2卷
第 3期
太 原 师 范 学 院 学 报 ( 自然 科 学 版 ) J OUR NAL O F T AI Y UAN N0R MAI UN I VE R S I TY ( N a t u r a l S c i e n c e E d i t i o n )
J 0 J 0
在本 文 中我们使 用 下列记 号 :
r 。 。 r 。 。
p l 一1 一J g ( s ) d s , , 9 2 —1 一1 ( 5 ) d s , f D —m i n { p  ̄ , P 2 ) , Q —Q×( 0 , £ ) ;
rf
(・,・> 代 表 H ( Q) 与 H ( Q) 空 问 的对 偶对 . <・, ・> 代表 L 的内积.
Hale Waihona Puke C w ( E o , ] , x) 表示其 中的所有 函数 ( f ):[ O , T] 一 x满足在 [ 0 , T] 上 的 连续 映射 £ 一< , ( £ ) ) 对 于
一
△ “ ~a △ “ +l g ( t —r ) △ “ ( r ) d r +e l I m - 1 “ 一f 1 ( “ , ) , ( z , £ ) ∈ Q×( 0 , 丁 )
J 0
rt
一
△ 一a △ +I h ( t —r ) A v ( r ) d r +e I l r _ 一f 2 ( “ , ) , ( L z , £ ) E Q×( o , T )
强阻尼波动方程组解的性质
![强阻尼波动方程组解的性质](https://img.taocdn.com/s3/m/54dc7a08cc17552707220835.png)
Ⅱ ( , 0 ) =/ / , l ( ) ,
∈ n
( , 0 )= 1 ( ) ,
其中, ( u 。 , 。 ) = ( 砩( ) ) , ( u 。 。 ) ∈ ( L ( ) ) .
同文献 [ 2 ] 一样 , F ( , ) 形如 ( 3 ) 式, 被看 作 强
假设 F( u , ) V C l ( u , ) ( , V ) 满 足如下 条件 : C 0 ( 1 U I + l l )≤ F( “ , )≤
C 1 ( 1 M l + l I ) , ( , z , ) ∈R , ( , ) =F ( U , )=
构建 了稳 定集 和不 稳 定集 , 运 用 反证 法 证 明了该 方 程组 的整 体解 存在 性 , 一致 衰减 性 和解 的爆破 . 文献[ 3 ] 研 究 了如下 方程
收 稿 日期 : 2 0 1 0— 0 9—1 9 基金项 目: 国家 自然 科 学 基 金 ( 1 0 9 7 1 1 4 7) 资助项 目
2 0 1 3年 9月
四川师范大学学报 ( 自然科 学版 )
J o u na r l o f S i c h u a n N o r m a l U n i v e r s i t y ( N a t u r a l S c i e n c e )
S e p t . , 2 0 1 3 V o l _ 3 6. N o . 5
设 OCR 是 一个 有 界正则 区域 , 边 界a =F.
r 1 2 , “ 一 △ 一 ∞△u + {u£I m - u£=I u I P _ ,
文献[ 1 ] 考 虑 了如下 方程
” 一
{ ( f , )∈[ 0 , T ) × R “ ,
具变指数非线性项的强阻尼波动方程的爆破
![具变指数非线性项的强阻尼波动方程的爆破](https://img.taocdn.com/s3/m/5a4c4c5af02d2af90242a8956bec0975f565a456.png)
第62卷 第1期吉林大学学报(理学版)V o l .62 N o .12024年1月J o u r n a l o f J i l i nU n i v e r s i t y (S c i e n c eE d i t i o n )J a n 2024d o i :10.13413/j .c n k i .jd x b l x b .2023090具变指数非线性项的强阻尼波动方程的爆破李海霞1,曹春玲2(1.长春师范大学数学学院,长春130032;2.吉林大学数学学院,长春130012)摘要:考虑一类具变指数非线性项的强阻尼波动方程的有限时刻爆破.借助凹方法和适当选取的参数,给出该问题的一个新的爆破准则,并给出爆破时间的上下界估计.结果表明,该爆破准则特别蕴含对任意高的初始能量,该问题存在有限时刻爆破解.关键词:强阻尼;波方程;变指数非线性项;爆破中图分类号:O 175.2 文献标志码:A 文章编号:1671-5489(2024)01-0078-09B l o w -u p t o S t r o n g l y D a m p e d W a v eE qu a t i o nw i t h V a r i a b l e -E x po n e n tN o n l i n e a rT e r m L IH a i x i a 1,C A O C h u n l i n g2(1.S c h o o l o f M a t h e m a t i c s ,C h a n g c h u nN o r m a lU n i v e r s i t y ,C h a n gc h u n 130032,C h i n a ;2.C o l l e g e o f M a t h e m a t i c s ,J i l i nU n i v e r s i t y ,C h a n gc h u n 130012,C h i n a )A b s t r a c t :W e c o n s ide r e dt h ef i n i t e t i m eb l o w -u p t oas t r o ng l y d a m p e d w a v ee q u a t i o n w i t hv a r i a b l e -e x p o n e n t n o n l i n e a r t e r m.W i th t h e h e l p o f c o n c a v em e t h o d a n d a p p r o p ri a t e l y s e l e c t e d p a r a m e t e r s ,w e g a v e an e wb l o w -u p c r i t e r i o n f o r t h i s p r o b l e ma n d e s t i m a t e d t h e u p p e r a n d l o w e r b o u n d s o n t h e b l o w -u p t i m e .T h e r e s u l t s s h o wt h a t t h eb l o w -u p c r i t e r i o n c o n t a i n s s p e c i a l i m p l i c a t i o n s f o r a n y h i g h i n i t i a l e n e r g y ,a n d t h e p r o b l e mh a s a f i n i t e t i m eb l o w -u p s o l u t i o n s .K e y w o r d s :s t r o n g l y d a m p e d ;w a v e e q u a t i o n ;v a r i a b l e -e x p o n e n t n o n l i n e a r t e r m ;b l o w -u p 收稿日期:2023-03-20.第一作者简介:李海霞(1984 ),女,汉族,博士,副教授,从事偏微分方程的研究,E -m a i l :l i h a i x i a 0611@126.c o m.通信作者简介:曹春玲(1971 ),女,汉族,博士,副教授,从事偏微分方程的研究,E -m a i l :c a o c l @jl u .e d u .c n .基金项目:国家自然科学基金(批准号:11626044)和长春师范大学人才引进启动项目(批准号:长师大R C [2016]第008号).0 引 言考虑如下具变指数非线性项的半线性强阻尼波方程的初边值问题:u t t -Δu -Δu t =u p (x )-2u ,x ɪΩ, t >0,u (x ,t )=0,x ɪ∂Ω, t >0,u (x ,0)=u 0(x ), u t (x ,0)=u 1(x ),x ɪΩìîíïïïï,(1)其中Ω⊂ℝn (n ȡ1)是具光滑边界∂Ω的有界区域,p (x )是Ω上的可测函数,u 0ɪH 10(Ω),u 1ɪL 2(Ω).目前,关于初边值问题u t t -Δu -ωΔu t +μu t =u p-2u ,x ɪΩ, t >0,u (x ,t )=0,x ɪ∂Ω, t >0,u (x ,0)=u 0(x ), u t (x ,0)=u 1(x ),x ɪìîíïïïïΩ(2)解的存在与非存在性㊁稳定性与爆破研究已取得了丰富的成果.例如:对非阻尼情形,即当ω=μ=0时,S a t t i n ge r [1]利用稳定集和不稳定集研究了问题(2)局部和整体解的存在性;B a l l [2]得到了该问题在负初始能量时的一个爆破结果;当ω=0,μ>0且适当小时,I k e h a t a [3]给出了问题(2)存在有限时刻爆破解的一个刻画;G a z z o l a 等[4]对问题(2)当参数满足ωȡ0,μ>-ωλ1时进行了系统研究,其中λ1>0是-Δ在H 10(Ω)上的第一特征值,对次临界和临界初始能量的情形,他们证明了整体解的存在与非存在性以及整体解的渐近性.此外,当ω=0,μȡ0时,文献[4]还证明了超临界初始能量情形下问题(2)有限时刻爆破解的存在性,但并未明确当ω>0时,问题(2)是否存在有限时刻爆破解.最近,Y a n g等[5]通过构造合适的泛函并选取恰当参数对上述问题给出了肯定回答.另一方面,具变指数的发展方程近年来也备受关注[6-9].关于具变指数的弱阻尼波动方程的爆破结果,M e s s a o u d i 等[10]考虑了如下波动方程:u t t -Δu +a u t m (x )-2u t =b u p (x )-2u ,(3)其中a 和b 是正常数,m (㊃)和p (㊃)是可测函数.在对问题(3)建立了弱解的存在唯一性之后,文献[10]得到了一个负初始能量时的有限时刻爆破结果.M e s s a o u d i 等[11]将上述爆破结果推广至如下拟线性弱阻尼波动方程上:u t t -d i v (∇u r (x )-2∇u )+a u t m (x )-2u t =b u p (x )-2u ,(4)此外,还将上述问题的爆破条件推广为适当小的正初始能量.对具变指数项的半线性或拟线性强阻尼波动方程的爆破研究目前也有一些结果.例如:A n t o n t s e v [12]研究了拟线性波动方程u t t -d i v (a (x ,t )∇u r (x ,t )-2∇u )-αΔu t =b u p (x ,t )-2u +f (x ,t ),(5)当参数和变指数满足一定条件时,证明了该问题弱解的局部存在性以及负初始能量条件下弱解的有限时刻爆破,但当初始能量非负时,对问题(5)是否存在有限时刻爆破解未给出回答;P a r k [13]研究了问题(5)的一个特例,即问题(1)的爆破性质,利用L e v i n e [14]的凹方法,得到了当初始能量有正上界时问题(1)的一个爆破结果.通过比较文献[4-5]与文献[13]的爆破结果易见,考虑问题(1)在高初始能量条件下的爆破有一定的意义.本文通过将文献[5]中的技巧与L e v i n e 凹方法相结合,给出问题(1)当初始能量可以被 ∇u 0 22+2(u 0,u 1)控制时的一个新的爆破准则.特别地,该爆破准则蕴含问题(1)存在具任意初始能量的有限时刻爆破解.此外,本文还给出爆破时间的上界估计.最后,借助一阶微分不等式技巧,得到爆破时间的一个下界.1 预备知识设p :Ωң[1,ɕ)是一个有界可测函数,并记p 1=e s s i n f x ɪΩp (x ), p 2=e s s s u p x ɪΩp (x ).变指数L e b e s gu e 空间L p (x )(Ω)定义为L p (x )(Ω)=u (x ):u 在Ω上可测,且存在λ>0使得ʏΩλu (x )p (x )d x <{}ɕ.该空间在被赋以范数(L u x e m b u r g 范数)u p (x )=i n f λ>0:ʏΩu (x)λp (x )d x ɤ{}1后成为一个B a n a c h 空间.用 ㊃ q 表示L q (Ω)空间的范数,用(㊃,㊃)表示L 2(Ω)上的内积,在H 10(Ω)中赋以如下内积和范数:<u ,v >=ʏΩ∇u ㊃∇v d x , u 2H 10(Ω)= ∇u 22,并用<㊃,㊃>*表示H -1(Ω)与H 10(Ω)之间的对偶积.此外,用λ1>0表示-Δ在Ω上具有齐次D i r i c h l e t 边界条件的第一特征值.从而97 第1期 李海霞,等:具变指数非线性项的强阻尼波动方程的爆破λ1 u 22ɤ ∇u 22, ∀u ɪH 10(Ω).(6) 引理1[15]设p :Ωң[1,ɕ)是有界可测函数,且满足2ɤp (x )<ɕ, n =1,2;2ɤp 1ɤp (x )ɤp 2<2n n -2,n ȡ3.则嵌入关系H 10(Ω)L p (x )(Ω)是连续且紧的.本文考虑问题(1)在下述意义下的弱解u (x ,t ).在不产生歧义时,u (x ,t )也常被简记为u (t ).弱解的局部存在和唯一性可通过适当修改文献[4,10]中的方法得到.定义1[4,13] 如果u ɪC ([0,T ];H 10(Ω))ɘC 1([0,T ];L 2(Ω))ɘC 2([0,T ];H -1(Ω))满足u t ɪL 2(0,T ;H 10(Ω)),u (0)=u 0,u t (0)=u 1,且对任意的ϕɪH 10(Ω)及几乎所有的t ɪ[0,T ],均有<u t t (t ),ϕ>*+ʏΩ∇u (t )㊃∇ϕd x +ʏΩ∇u t (t )㊃∇ϕd x =(u (t )p (㊃)-2u (t ),ϕ),(7)则称u (x ,t )为问题(1)在[0,T ]上的一个弱解.假设条件:(H 1)设p (x )满足2<p 1ɤp (x )ɤp 2<2*,其中当n =1,2时,2*=ɕ,当n ȡ3时,2*=2n -2n -2;(H 2)对数型H öl d e r 连续性条件:存在A >0及δɪ(0,1),使得对几乎所有满足x -y <δ的x ,y ɪΩ,均有p (x )-p (y )ɤ-Al o g x -y. 定理1[13] 若假设条件(H 1),(H 2)成立,则对任意的u 0ɪH 10(Ω),u 1ɪL 2(Ω),问题(1)存在唯一弱解.记T m a x 为问题(1)弱解u =u (t )的最大存在时间,并定义能量泛函㊁N e h a r i 泛函与不稳定集分别为E (t )=E (u (t ))=12 u t 22+12∇u 22-ʏΩu p (x)p (x )d x , t ɪ[0,T m a x ),(8)I (u )= ∇u 22-ʏΩup (x )d x , u ɪH 10(Ω),(9)N -={u ɪH 10(Ω)I (u )<0}.(10)在式(7)中取ϕ=u t 可得能量恒等式:E (t )+ʏt 0∇u τ 22d τ=E (0), a .e .t ɪ[0,T m a x ).(11)式(11)表明E (t )关于t 在[0,T m a x )上是单调不增的.注1 设T m a x 是u (t )的最大存在时间,若T m a x <ɕ,则必有li m t ңT m a x∇u 2=ɕ.事实上,利用与文献[4]中证明定理3.1类似的讨论和连续延拓定理[16]可知,如果T m a x <ɕ,则l i m t ңT m a x∇u (t ) 22+ u t (t ) 22=+ɕ.(12)于是,由式(8),(11)和S o b o l e v 不等式可得12 u t (t ) 22+12∇u (t ) 22ɤE (0)+ʏΩu p (x )p (x )d x ɤE (0)+1p1ʏΩup (x )d x ɤE (0)+1p1ʏΩ(1+u p 2)d x ɤC 1+C 2 ∇u (t ) p22,(13)其中t ɪ[0,T m a x ),C 1,C 2是不依赖于u 的正常数.结合式(12)和式(13)易见l i m t ңT m a x∇u 2=ɕ.定义2 设u (t )是问题(1)的弱解,T m a x 是它的最大存在时间.若T m a x <ɕ,则称u(t )在有限时刻爆破,并称T m a x 是u 的爆破时间.若T m a x =ɕ,则称u(t )是整体存在的.2 主要结果引理2[14,17]假设ψ(t)是一个正的二次可微函数,且满足08 吉林大学学报(理学版) 第62卷ψᵡ(t )ψ(t )-(1+θ)(ψᶄ(t ))2ȡ0,其中θ>0.如果ψ(0)>0,ψᶄ(0)>0,则当t ңt *ɤt *=ψ(0)θψᶄ(0)时,ψ(t )ңɕ.受文献[5,18]启发,本文证明当初始能量可在上方被 ∇u 0 22+2(u 0,u 1)控制时,问题(1)的解在有限时刻爆破,并给出爆破时间的上界估计.定理2 若变指数p (x )满足假设条件(H 1),(H 2),初值u 0ɪH 10(Ω),u 1ɪL 2(Ω)满足u 0ɪN -,且 ∇u 0 22+2(u 0,u 1)>2p 1(1+λ1)(p1-2)λ1E (0).(14)则问题(1)的解u (t )在有限时刻爆破,且爆破时间T m a x 的上界满足如下估计:T m a x ɤ4[(a 2+(λ-2)2b 0 u 0 22)1/2+a ](λ-2)2b 0,(15)其中λ,a ,b 0是确定的常数.证明:证明过程分为4步.1)证明只要u (t )ɪN -,则函数{t F (t )=ʒ ∇u (t ) 22+2(u ,u t )}在(0,T m a x )上即为严格单调递增的.事实上,将F (t )关于t 求导㊁并结合式(1)及I (u (t ))<0(t ɪ[0,T m a x )),可得F ᶄ(t )=2(∇u ,∇u t )+2 u t 22+2(u ,u t t )=2[(∇u ,∇u t )+ u t 22+(u ,Δu +Δu t +u p (x )-2u )]=2 u t 22- ∇u 22+ʏΩup (x )d ()x =2( u t 22-I (u ))>0.(16)表明F (t )在(0,T m a x )上严格单调递增.2)证明u (t )ɪN -,t ɪ[0,T m a x ).若不然,则由连续性及假设条件I (u 0)<0可知,存在t 0ɪ(0,T m a x ),使得I (u (t ))<0, t ɪ[0,t 0),(17)I (u (t 0))=0.(18)于是,由1)中结论知F (t )= ∇u (t ) 22+2(u ,u t )在区间[0,t 0)上严格递增.将其与式(14)结合,进一步可得∇u (t ) 22+2(u ,u t )> ∇u 0 22+2(u 0,u 1)>2p 1(1+λ1)(p 1-2)λ1E (0), t ɪ(0,t 0).(19)注意到F (t )的连续性和单调性,由式(19)得 ∇u (t 0) 22+2(u (t 0),u t (t 0))>2p 1(1+λ1)(p1-2)λ1E (0).(20)另一方面,由式(8)和式(9)得E (t )ȡ12 u t 22+12 ∇u 22-1p1ʏΩup (x )d x =12 u t 22+p 1-22p 1 ∇u 22+1p1I (u ), t ɪ[0,T m a x ).(21)由式(6),(11),(18),(21)和C a u c h y -S c h w a r z 不等式,可得E (0)ȡE (t 0)ȡ12 u t (t 0) 22+p 1-22p1 ∇u (t 0) 22=12 u t (t 0) 22+p 1-22p 1(1+λ1) ∇u (t 0) 22+(p 1-2)λ12p1(1+λ1) ∇u (t 0) 22ȡ12 u t (t 0) 22+(p 1-2)λ12p 1(1+λ1) u (t 0) 22+(p 1-2)λ12p 1(1+λ1) ∇u (t 0) 22ȡ(p 1-2)λ12p1(1+λ1)[ u t (t 0) 22+ u (t 0) 22+ ∇u (t 0) 22]ȡ(p 1-2)λ12p1(1+λ1)[2(u t (t 0),u (t 0))+ ∇u (t 0) 22],(22)18 第1期 李海霞,等:具变指数非线性项的强阻尼波动方程的爆破与式(20)矛盾.故结论成立.3)通过选取适当的参数并结合L e v i n e 凹方法证明T m a x <ɕ.若不然,假设u 是问题(1)的整体解,则T m a x =ɕ.类似于文献[13]中的证明,对任意T >0,b >0,η>0,定义G (t )=ʏt 0∇u (τ) 22d τ+ u (t ) 22+(T -t ) ∇u 0 22+b (t +η)2, t ɪ[0,T ].(23)则G (t )>0,t ɪ[0,T ],G ᶄ(t )= ∇u (t ) 22+2(u ,u t )- ∇u 0 22+2b (t +η)=2(u ,u t )+2ʏt 0ʏΩ∇u ㊃∇u τdx d τ+2b (t +η),(24)G ᵡ(t )=2(u t,u t)+2(u ,u t t)+2ʏΩ∇u ㊃∇u tdx +2b =2 u t 22+2(u ,Δu +u p (x )-2u )+2b =2[ u t 22-I (u )]+2b .(25)由式(24)得(G ᶄ(t ))2=4(u ,u t )+ʏt 0ʏΩ∇u ㊃∇u τd x d τ+b (t +η[])2=4(u ,u t)2+2(u ,u t)ʏt 0ʏΩ∇u ㊃∇u τd x d τ+2b (t +η)(u ,u t[)+ʏt 0ʏΩ∇u ㊃∇u τd x d ()τ2+2b (t +η)ʏt0ʏΩ∇u ㊃∇u τdx d τ+b 2(t +η)]2.(26)利用C a u c h y -S c h w a r z 不等式知(u ,u t )ɤ u 2 u t 2,(27)ʏt0ʏΩ∇u ㊃∇u τd x d τɤʏt∇u 2∇u τ2d τɤʏt∇u 22d ()τ1/2ʏt∇u τ22d()τ1/2.(28)将式(27),(28)与C a u c h y 不等式相结合,进一步可得(G ᶄ(t ))24ɤ u 22 u t 22+2 u 2 u t 2ʏt∇u 22d ()τ1/2ʏt∇u τ22d()τ1/[2+2b (t +η) u 2u t2+ʏt∇u 22d τʏt∇u τ22dτ+2b (t +η)ʏt∇u 22d ()τ1/2ʏt∇u τ22d()τ1/2+b 2(t +η)]2ɤ u 22u t22+ʏt∇u τ22d()[τ+ʏt∇u 22d τ u t22+ʏt∇u τ22d ()τ+b u 22+ʏt∇u 22d()τ+b (t +η)2ˑu t22+ʏt∇u τ22d()τ+b 2(t +η)]2= u 22+ʏt∇u 22d τ+b (t +η)()2u t22+ʏt∇u τ22dτ+()b .(29)对任意的λɪ(2,p1).注意到式(23),(25),(29)可得G (t )G ᵡ(t )-λ+24(G ᶄ(t ))2ȡG (t )G ᵡ(t )-(λ+2) u t 22+ʏt 0∇u τ22dτ+()[]b =ʒG (t )H (t),(30)其中H (t )=-λ u t 22-2I (u )-(λ+2)ʏt∇u τ 22d τ-b λ.由式(6),(8),(11),(9)知28 吉林大学学报(理学版) 第62卷H (t )ȡ(p 1-λ) u t 22+(p 1-2) ∇u 22-2p 1E (0)+(2p1-λ-2)ʏt∇u τ 22d τ-b λȡ(p 1-λ) u t 22+(p 1-2)-p 1-21+λéëêêùûúú1 ∇u 22+p 1-21+λ1∇u 22-2p 1E (0)-b λȡ(p 1-λ) u t 22+(p 1-2)λ11+λ1( u 22+ ∇u 22)-2p1E (0)-b λ.(31) 下面取λ=p 1-(p 1-2)λ11+λ1ɪ(2,p1).(32)由式(31)及F (t )的单调性知,对任意的b ɪ0,(p 1-2)λ1λ(1+λ1)2(u 0,u 1)+ ∇u 0 22-2p 1(1+λ1)(p 1-2)λ1E (0æèçöø÷æèçùûúú),(33)均有H (t)ȡ(p 1-2)λ11+λ1( u t 22+ u 22+ ∇u 22)-2p1E (0)-b λȡ(p 1-2)λ11+λ12(u ,u t )+ ∇u ()22-2p1E (0)-b λȡ(p 1-2)λ11+λ12(u 0,u 1)+ ∇u 0 22-2p 1(1+λ1)(p 1-2)λ1E (0éëêêùûúú)-b λȡ0.(34)于是,结合式(30)和式(34)可知,对任意满足式(33)的b ,均有G (t )G ᵡ(t )-λ+24(G ᶄ(t ))2ȡ0,t ɪ[0,T ].(35) 选取不依赖于T 的η满足η>m a x0,2 ∇u 0 22-(λ-2)(u 0,u 1)b (λ-2{}),(36)则G (0)= u 0 22+T ∇u 0 22+b η2>0,G ᶄ(0)=2(u 0,u 1)+2b η>4 ∇u 0 22λ-2ȡ0,且对充分大的T 有4G (0)(λ-2)G ᶄ(0)=2( u 0 22+T ∇u 0 22+b η2)(λ-2)[(u 0,u 1)+b η]<T .(37)对G (t )应用引理2可知,存在t *>0满足t *ɤ4G (0)(λ-2)G ᶄ(0)(<T ),(38)使得l i m t ңt -*G (t )=li m t ңt -*ʏt 0∇u (τ) 22d τ+ u (t ) 22+(T -t ) ∇u 0 22+b (t +η)[]2=ɕ.(39)由式(39)易知l i m t ңt -*∇u 22=ɕ.(40)由注1知,式(40)表明u (t )在t *时刻爆破.这与u (t )是问题(1)的整体解矛盾.故T m a x <ɕ.4)给出T m a x 的上界估计.类似于G (t ),定义췍G (t )为췍G (t )=ʏt 0∇u (t ) 22d τ+ u (τ) 22+(T m a x -t ) ∇u 0 22+b (t +η)2, t ɪ[0,T m a x ).利用平行于3)的证明,可得l i m t ңT m a x췍G (t )=ɕ,其中38 第1期 李海霞,等:具变指数非线性项的强阻尼波动方程的爆破T m a x ɤ2( u 0 22+T m a x ∇u 0 22+b η2)(λ-2)[(u 0,u 1)+b η],(41)λ仍由式(32)给出,b 满足式(33).如果仍要求η满足式(36),则式(41)可改写为T m a x ɤT (b ,η)=ʒ2( u 0 22+b η2)(λ-2)[(u 0,u 1)+b η]-2 ∇u 0 22.(42)固定b 满足式(33).记a =2 ∇u 0 22-(λ-2)(u 0,u 1), η0=[a 2+(λ-2)2b u 0 22]1/2+a (λ-2)b.直接验证易知T (b ,η)在η=η0处取到其在η>m a x0,2 ∇u 0 22-(λ-2)(u 0,u 1)b (λ-2{})上的最小值,且m i n ηT (b ,η)=T (b ,η0)=4[(a 2+(λ-2)2b u 0 22)1/2+a ](λ-2)2b.最后,令函数T (b ,η0)关于b 在式(33)上取最小值,得m i n b T (b ,η0)=T (b 0,η0)=4[(a 2+(λ-2)2b 0 u 0 22)1/2+a ](λ-2)2b 0,其中b 0=(p 1-2)λ1λ(1+λ1)2(u 0,u 1)+ ∇u 0 22-2p 1(1+λ1)(p 1-2)λ1E (0éëêêùûúú). 综上,由式(42)可知T m a x ɤ4[(a 2+(λ-2)2b 0 u 0 22)1/2+a](λ-2)2b 0.证毕.注2 定理2表明问题(1)存在具任意高初始能量的有限时刻爆破解.首先选取u 0ɪH 10(Ω),u 1ɪL 2(Ω)满足(u 0,u 1)>0,并令u 0=αu 0,u 1=βu 1,其中α,β是待定正常数.由于p1>2,故存在α1>1,使得I (u 0)=I (αu 0)=α2 ∇u 0 22-ʏΩαp (x )u 0p (x )d x ɤα2 ∇u 0 22-αp 1ʏΩu 0p (x )d x <0, ∀αȡα1.此外,对任意的D >0,存在α2>0,使得 ∇u 0 22= α∇u 0 22>2p 1(1+λ1)(p1-2)λ1D , ∀αȡα2.令α=m a x {α1,α2},对α>0,可选取β>0,使得E (0)=D .注意到(u 0,u 1)>0,显然有 ∇u 0 22+2(u 0,u 1)>2p 1(1+λ1)(p 1-2)λ1D =2p 1(1+λ1)(p1-2)λ1E (0).根据定理2知,问题(1)以(u 0,u 1)为初值的解u (x ,t )在有限时刻爆破,且初始能量满足E (0)=D .下面给出问题(1)爆破时间的一个下界估计.由于爆破时间的下界可以为所考虑的系统提供一个安全(稳定)区间,所以在实际应用中非常重要.定理3 若p (x )满足假设条件(H 1)和(H 2),u (x ,t )是问题(1)的一个在T m a x 时刻爆破的弱解.则T m a x 可从下方估计为T m a x ȡʏɕK (0)d s C 1+C 2s p 2-1,其中K (0)= u 1 22+ ∇u 0 22,C 1,C 2是不依赖于u(x ,t )的正常数.证明:先利用特定泛函满足的一阶微分不等式得到爆破时间的一个下界估计.为避免证明过程冗长,这里只给出n ȡ3时的证明(n =1,2时类似可证).记48 吉林大学学报(理学版) 第62卷K (t )= u t (t ) 22+ ∇u (t ) 22, t ɪ[0,T m a x ).(43)由于u (x ,t )在T m a x 时刻爆破,因此由式(12)知l i m t ңT m a xK (t )=+ɕ.(44)直接计算可得K ᶄ(t )=2[(u t ,u t t )+(∇u ,∇u t )]=2(u t ,u t t -Δu )=2(u t ,Δu t +up (x )-2u )ɤ-2 ∇u t 22+2ʏΩu tu p (x )-1d x ɤ-2 ∇u t 22+2ʏΩu t(1+u p 2-1)d x .(45)对任意的r ɪ1,2n n æèçùûú-2,用S r 表示从H 10(Ω)到L r (Ω)的最佳嵌入常数,即 v r ɤS r ∇v 2, ∀v ɪH 10(Ω).(46)由于p (x )满足假设条件(H 1),直接验证可知2n (p 2-1)n +2<2n n -2,从而H 10(Ω)可以连续嵌入L 2n (p 2-1)/(n +2)(Ω).利用H öl d e r 不等式和带ε的C a u c h y 不等式,可将式(45)中第二项估计如下:2ʏΩut(1+up 2-1)d x ɤ2 u t 2n /(n -2)ʏΩ(1+u p 2-1)2n /(n +2)d ()x (n +2)/(2n )ɤε u t22n /(n -2)+C (ε)ʏΩ(1+u p 2-1)2n /(n +2)d ()x(n +2)/n ɤεS 22n /(n -2) ∇u t 22+C 1(ε)+C 2(ε) u 2(p 2-1)2n (p 2-1)/(n+2)ɤεS 22n /(n -2) ∇u t 22+C 1(ε)+C 2(ε)S 2(p 2-1)2n (p 2-1)/(n +2) ∇u 2(p 2-1)2.(47)选取ε>0适当小使得εS 22n /(n -2)ɤ2,并将式(47)代入(45)可得K ᶄ(t )ɤC 1+C 2K p 2-1(t ),(48)表明对任意的t ɪ[0,T m a x )均有t ȡʏt0K ᶄ(τ)C 1+C 2K p 2-1(τ)d τ.(49)令t ңT m a x 并注意到式(44),得T m a x ȡʏɕK (0)d sC 1+C 2s p 2-1.(50)由于p 2-1>1,故式(50)右端项是有限的.证毕.参考文献[1] S A T T I N G E R D H.O nG l o b a l S o l u t i o n o fN o n l i n e a rH y p e r b o l i cE q u a t i o n s [J ].A r c hR a t i o n a lM e c hA n a l ,1968,30:148-172.[2] B A L LJM.R e m a r k so nB l o w -u p a n d N o n e x i s t e n c eT h e o r e m sf o rN o n l i n e a rE v o l u t i o nE q u a t i o n s [J ].Q u a r tJ M a t hO x f o r dS e r 2,1977,28(112):473-486.[3] I K E HA T A R.S o m eR e m a r k s o n t h eW a v eE q u a t i o n sw i t hN o n l i n e a rD a m p i n g a n dS o u r c eT e r m s [J ].N o n l i n e a r A n a l :T h e o r y ,M e t h o d sA p p l ,1996,27(10):1165-1175.[4] G A Z Z O L A F ,S Q U A S S I N A M.G l o b a lS o l u t i o n sa n d F i n i t e T i m e B l o w u p f o r D a m p e d S e m i l i n e a r W a v e E qu a t i o n s [J ].A n n I n s tH P o i n c a r éA n a lN o nL i n éa i r e ,2006,23(2):185-207.[5] Y A N G YB ,X U RZ .N o n l i n e a r W a v eE q u a t i o n w i t hB o t hS t r o n g l y a n d W e a k l y D a m p e dT e r m s :S u p e r c r i t i c a l I n i t i a l E n e r g y F i n i t eT i m eB l o wu p [J ].C o mm u nP u r eA p p lA n a l ,2019,18(3):1351-1358.[6] A C E R B IE ,M I N G I O N EG.R e g u l a r i t y R e s u l t s f o rS t a t i o n a r y E l e t r o -R h e o l o g i c a lF l u i d s [J ].A r c hR a t i o n M e c h A n a l ,2002,164(3):213-259.[7] A N T O N T S E VSN ,R O D R I G U E S J F .O nS t a t i o n a r y T h e r m o -R h e o l o g i c a l V i s c o u s F l o w s [J ].A n nU n i vF e r r a r a 58 第1期 李海霞,等:具变指数非线性项的强阻尼波动方程的爆破68吉林大学学报(理学版)第62卷S e zⅦS c iM a t h,2006,52(1):19-36.[8] A N T O N T S E VS,S HMA R E VS.C h a p t e rⅠE l l i p t i cE q u a t i o n sw i t h A n i s o t r o p i cN o n l i n e a r i t y a n d N o n s t a n d a r dG r o w t hC o n d i t i o n s[M]//H a n d b o o ko fD i f f e r e n t i a lE q u a t i o n s:S t a t i o n a r y P a r t i a lD i f f e r e n t i a lE q u a t i o n s.N o r t h-H o l l a n d:E l s e v i e r,2006,3:1-100.[9] C H E N Y M,L E V I N ES,R A O M.V a r i a b l eE x p o n e n t,L i n e a rG r o w t hF u n c t i o n a l s i nI m a g eR e s t o r a t i o n[J].S I AMJA p p lM a t h,2006,66(4):1383-1406.[10] M E S S A O U D I SA,T A L A HM E H A A,A l-S MA I LJH.N o n l i n e a rD a m p e dW a v eE q u a t i o n:E x i s t e n c e a n dB l o w-u p[J].C o m p u tM a t hA p p l,2017,74(12):3024-3041.[11] M E S S A O U D I S A,T A L A HM E H A A.B l o w u p i nS o l u t i o n so fa Q u a s i l i n e a r W a v eE q u a t i o n w i t h V a r i a b l e-E x p o n e n tN o n l i n e a r i t i e s[J].M a t h M e t h o d sA p p l S c i,2017,40(18):6976-6986.[12] A N T O N T S E VS.W a v eE q u a t i o nw i t h p(x,t)-L a p l a c i a na n dD a m p i n g T e r m:E x i s t e n c e a n dB l o w-u p[J].D i f f e rE q uA p p l,2011,3(4):503-525.[13] P A R KSH.B l o w-u p o f S o l u t i o n s f o rW a v eE q u a t i o n sw i t hS t r o n g D a m p i n g a n dV a r i a b l e-E x p o n e n tN o n l i n e a r i t y[J].JK o r e a n M a t hS o c,2021,58(3):633-642.[14] L E V I N E H A.S o m eN o n e x i s t e n c e a n d I n s t a b i l i t y T h e o r e m s f o rS o l u t i o n so fF o r m a l l y P a r a b o l i cE q u a t i o no f t h eF o r m P u t=-A u+F(u)[J].A r c hR a t i o n a lM e c hA n a l,1973,51(5):371-386.[15] D I E N I N GL,HA R J U L E H T OP,HÄS TÖP,e t a l.L e b e s g u e a n dS o b o l e vS p a c e sw i t hV a r i a b l eE x p o n e n t s[M].H e i d e l b e r g:S p r i n g e r,2011:1-502.[16] O N O K.O nG l o b a l E x i s t e n c e,A s y m p t o t i cS t a b i l i t y a n dB l o w i n g u p o f S o l u t i o n s f o r S o m eD e g e n e r a t eN o n-l i n e a rW a v eE q u a t i o n so f K i r c h h o f f T y p e w i t ha S t r o n g D i s s i p a t i o n[J].M a t h M e t h o d s A p p lS c i,1997,20(2): 151-177.[17] HA N YZ,G A O WJ,S U NZ,e t a l.U p p e r a n dL o w e rB o u n d so fB l o w-u p T i m e t oaP a r a b o l i cT y p eK i r c h h o f fE q u a t i o nw i t hA r b i t r a r y I n i t i a l E n e r g y[J].C o m p u tM a t hA p p l,2018,76(10):2477-2483.[18] Y A N G H,HA N Y Z.B l o w-u p f o raD a m p e d p-L a p l a c i a nT y p e W a v eE q u a t i o n w i t hL o g a r i t h m i cN o n l i n e a r i t y[J].JD i f f e r e n t i a l E q u a t i o n s,2022,306:569-589.(责任编辑:赵立芹)。
N维空间中一类强阻尼非线性波动方程的解及其性质
![N维空间中一类强阻尼非线性波动方程的解及其性质](https://img.taocdn.com/s3/m/71b3710f0812a21614791711cc7931b764ce7b59.png)
N维空间中一类强阻尼非线性波动方程的解及其性质
廖扬;周晓宇
【期刊名称】《轻工学报》
【年(卷),期】2016(031)006
【摘要】将强阻尼非线性波动方程在三维空间解的性质由三维推广到N(N〉3)维,利用标准的Galerkin方法和Sobolev嵌入定理研究了弱解在该空间下的存在性,运用内积做出了解的耗散估计,并采用Gronwall引理证明了整体吸引子的存在性.
【总页数】5页(P95-99)
【作者】廖扬;周晓宇
【作者单位】河南财经政法大学数学与信息科学学院,河南郑州450002
【正文语种】中文
【中图分类】O176.3
【相关文献】
1.N维空间中一类强阻尼非线性波动方程的解及其性质
2.n维空间中一类耦合系统解的爆破性质
3.强阻尼非线性波动方程的解半群存在吸收集
4.强阻尼非线性波动方程解的渐近性质
5.一类非线性强阻尼扰动发展方程的解
因版权原因,仅展示原文概要,查看原文内容请购买。
具有强阻尼的非线性梁方程的全局吸引子
![具有强阻尼的非线性梁方程的全局吸引子](https://img.taocdn.com/s3/m/0fec34dd8ad63186bceb19e8b8f67c1cfad6ee10.png)
具有强阻尼的非线性梁方程的全局吸引子任永华;张建文【摘要】研究了齐次Dirichlet边界条件下非线性强阻尼梁方程的长时间动力学行为.以经典的算子半群理论和能量的估计技巧为依据,获得了能量的一致衰退估计;通过相空间中新范数的引入,得到了系统的整体解所对应的解半群存在一致有界吸收集;利用半群分解技术,证明了具有强阻尼项的非线性梁方程解的全局吸引子在空间E上的存在性.【期刊名称】《中北大学学报(自然科学版)》【年(卷),期】2014(035)001【总页数】4页(P7-10)【关键词】非线性梁方程;强阻尼;全局吸引子;等价模【作者】任永华;张建文【作者单位】太原理工大学数学学院,山西太原030024;太原理工大学数学学院,山西太原030024【正文语种】中文【中图分类】O1750 引言在非线性发展方程的领域中,系统的长时间动力学行为是由其所对应的半群的吸引子来描述的.一直以来,对系统解的吸引子的研究受到众多科技工作者的高度关注,并涌现出了大量的文献[1-12].当函数,g′(u)有界,g(0)=0,且时,考虑了梯度系统(1)~(3)的渐近行为[11].对于一般抽象函数,由系统(1)~(3)定义的连续半群是点耗散的,并且系统存在整体吸引子[12].目前,弦方程解的吸引子理论已经得到了相当的发展,但对于梁方程的吸引子的研究文献还相对较少.本文在前人研究的基础上,考虑了齐次Dirichlet边界条件下具有强阻尼的梁方程系统在空间E=V×H中全局吸引子的存在性.赋予初始条件和齐次Dirichlet边界条件式中:γ>0,u=u(x,t)是关于变量x和t的Ω×R+ 的实值算子函数.Ω⊂Rn 是具有充分光滑边界Γ的有界开集.1 解的存在唯一性设L2(Ω))是一个自伴正定线性算子,并且算子A的特征值{λi}i∈N满足令H=L2(Ω),V=H20(Ω)和E=V×H,且它们的内积和范数分别为下面考虑系统(1)~(3).为了证明解的存在性,假设函数g(u,v)满足:式中:∀(u,v)∈R×R,Ci为非负常数,0<δ<1.令,则容易将系统(1)~(3)改写为空间E上关于时间的一阶抽象发展方程.系统等价于Hilbert空间E上的初值问题由文献[12]结合上述结论可知,半群eCt的无穷小生成元为C,又由于C是一个扇形算子,则eCt是E中生成的一个解析半群eCt.由参考文献[6],易知N (U)在E上是全局Lipschitz 连续的.再由微分方程的解的存在唯一性理论,有定理1.定理1 对于任意给定的Z0∈E,假设α>0,g(u,v)满足条件(H1)~(H2),那么存在唯一函数Z(t)=Z(t,Z0)∈C(R+,E),使得Z0=Z (0,Z0)且Z(t)满足下面的积分方程Z(t,Z0)关于t和Z0共同连续,根据定理1 解的存在唯一性,对于任意的t≥0,引入一个E上的自治动力系统,则可以定义空间E上的一个连续半群{S(t),t≥0}.其中映射S(t):Z0→Z (t,Z0),且Z(t,Z0)是系统(4)的mild 解.2 全局吸引子的存在性为了得到本文的主要定理,在空间E上定义加权内积和范数如下:其中:易知,函数|·|E等价于E中的通常范数‖·‖E.设,其中设φ=(u,w)T ∈E,则系统(1)~(3)(或(4))可以等价地写成为了研究解的全局吸引子的存在性,首先介绍下面的引理:引理1 对于∀φ=(u,w)T,有证明对于∀φ=(u,w)T,当μ=1-εγ时,有通过简单的计算可得因此,引理得证.接下来讨论空间E上半群{S(t),t≥0}的吸收性质.引理2 对于∀φ=(u,w)T,有或证明∀φ=(u,w)T ∈Ε是系统的解.用φ在E中与问题(6)做内积(.,.)E,得由引理1可知由式(7)可得根据Gronwall不等式,可得到下列在空间(E,|·|E)中的吸收不等式根据上述结论,可直接得出:对应于问题(1)~(3)的半群{S(t),t≥0}存在一致有界吸收集B0.也就是说,吸收集是一致吸收的.因此,当t≥t1(B)时,对于E的任意有界集B,S(t)B⊆B0成立. 引理3 对应于问题(6)的半群{S(t),t≥0}在E中存在有界吸收集.证明设φ=(u,v)T 是问题(6)的解.令φ1(t)和φ2(t)是初值分别为问题(6)中φ1和φ2的两个解,且分别满足下面两个方程:由此可得,Sε(t)=S1(t)+S2(t).类似于引理1和2 的证明,可知成立,因此可得S1(t)是指数衰减的.由于φ2(0)=(0,0)T ∈D(C),用Aφ2在E中与方程(14)作内积,可得类似于引理1的证明,可得另外,通过简单的计算有结合上述计算可得于是,当B0⊂E是一个有界集时,中的有界集.又由于D(A)×紧嵌入E,故是D(A)×中的紧集.结合引理2和引理3,可得下列定理.定理2 假设α>0,g(u,v)满足条件(H1)和(H2),则系统(1)~(3)在空间E中所定义的算子半群{Sε(t),t≥0}存在全局吸引子.参考文献:[1]Temam R.Infinite-Dimensional Dynamical Systems in Mechanics and Physics[M].New York:Springer-Verlag,1988.[2]Zhang J W,Ren Y H,Wu R H,et al.The global attractor of nonlinear thermoelastic coupled Sine-Gordon system[J].Acta Phys.Sin.,2012,61:110404-1-5.[3]Ren Y H.Longtime behavior of a non-autonomous beam equation [J].Scientific Bulletin,Series A:Applied Mathematics and Physics,2013,75:135-146.[4]Zhou S F,Fan X M.Kernel sections for non-autonomous strongly damped wave equations[J].J.Math.Anal.Appl.,2002,275:850-869. [5]Fan X M,Zhou S F.Kernel sections for non-autonomous strongly damped wave equations of non-degenerate Kirchhoff-type[J].Applied Mathematics and Computation,2004,158:253-266.[6]Jiang Y,Xie Y Q.Global attractors for a class nonlinear evolution equation[J].Mathematical Theory and Applications,2010,30:24-28. [7]Wang C M,Wang X,Zhong C K.Existence of global attractors for weak dissipation abstract evolution equations[J].Pure and Applied Mathematics,2012,28:401-411.[8]Meyries M.Global attractors in stronger norms for a class of parabolic systems with nonlinear boundary conditions[J].Nonlinear Analysis,2012,75:2922-2935.[9]Wakasugi Y.Small data global existence for the semilinear wave equation with space-time dependent damping[J].J.Math.Anal.Appl.,2012,393:66-79.[10]Pazy A.Semigroups of Linear Operator and Applications to Partial Differential Equations[M].New York:Springer-Verlag,1983.[11]Webb G F.Existence and asymptotic behavior for a strongly damped nonlinear wave equation[J].Canad.J.Math.,1980,32:334-349. [12]Massatt P.Limiting behavior for strongly damped nonlinear wave equations[J].J.Differential Equations,1983,48:334-349.。
一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程
![一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程](https://img.taocdn.com/s3/m/79f7f771cdbff121dd36a32d7375a417866fc1f7.png)
应用数学MATHEMATICA APPLICATA2022,35(3):511-523一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程杨怡,方钟波(中国海洋大学数学科学学院,山东青岛266100)摘要:本文研究一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程齐次Dirichlet初边值问题解的适定性及定性性质.借助于正则解的适定性并结合稠密性理论导出局部弱解的适定性,且利用修正能量泛函技巧,建立当p<γ时整体适定性.同时,利用反证技巧,证明当p>γ时解的有限时刻爆破现象.关键词:半线性波动方程;基尔霍夫型弱阻尼;对数非线性;整体适定性;爆破中图分类号:O175.29AMS(2010)主题分类:35L20;35A01;35B44文献标识码:A文章编号:1001-9847(2022)03-0511-131.引言我们考虑一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程u tt−∆u−σ(∥∇u∥22)|u t|γ−2u t=|u|p−2u log|u|,(x,t)∈Ω×(0,+∞),(1.1)给出齐次Dirichlet边界条件和初始条件u(x,t)=0,(x,t)∈∂Ω×(0,+∞),(1.2)u(x,0)=u0(x),u t(x,0)=u1(x),x∈Ω,(1.3)其中Ω⊂R N(N≥1)为具有光滑边界∂Ω的有界区域,基尔霍夫型阻尼系数σ,参数p和γ满足(H1)非局部阻尼系数σ∈C1([0,+∞)),对任意s∈[0,+∞),存在正常数b使得σ(s)≥b;;若N=1,2,则2<γ,p<+∞.(H2)参数γ,p满足:若N≥3,则2<γ,p<2N−2N−2我们的模型(1.1)出现于弹性力学理论,量子力学理论等.比如,在一维和二维空间中,具有常数阻尼系数(σ(∥∇u∥2)≡1)的模型(1.1)分别表示服从于粘性效应的均匀弦的横向振动2)=1))的非局部弱阻尼和均匀杆的纵向振动.特别是,当γ=2时,不为常数系数(σ(∥∇u∥22)u t模拟了作用在物体上的摩擦机制,该机制取决于u自身平均值.[1]同时,在量子场项σ(∥∇u∥22理论中,对数非线性项出现在膨胀宇宙学以及超对称场理论中.本文中,我们的目的在于分析在模型(1.1)-(1.3)中基尔霍夫型弱阻尼项σ(∥∇u∥2)|u t|γ−2u t2和对数非线性项|u|p−2u log|u|之间的相互作用与竞争关系对解的适定性及定性性质的影响.为了陈述研究动机,我们回顾此类问题的研究背景.实际上,具有常数系数阻尼项和乘幂型源项的如下半线性波动方程已被许多学者广泛关注且已有许多进展.[2−3]u tt−∆u+g(u t)=|u|p−2u.(1.4)∗收稿日期:2021-06-01基金项目:山东省自然科学基金面上项目(ZR2019MA072);中央高校基本科研基金(201964008)作者简介:杨怡,女,汉族,山东人,研究方向:偏微分方程.通讯作者:方钟波.512应用数学2022比如,Georgiev和Todorova[2]研究了方程(1.4)中g(u t)=a|u t|γ−2u t情形且在齐次Dirichlet边界条件下,证明了问题解的整体存在性及有限时刻爆破现象.之后,文[3]给出了具有任意负初始能量及正初始能量解的爆破现象.关于具有基尔霍夫阻尼系数(σ(∥∇u∥22)=1))的问题的研究方面,大部分集中于基尔霍夫型拟线性波动方程及四阶波动方程解的长时间动力行为,而很少有文献研究爆破现象.比如, Jorge和Narciso[4]在研究如下可扩展梁方程时首次提出了基尔霍夫型阻尼系数u tt+∆2u−κϕ(∥∇u∥22)∆u+σ(∥∇u∥22)g(u t)+f(u)=h(x),其中g(u t)≈|u t|γu t.他们在Dirichlet边界条件和铰接边界条件下建立了强解的适定性并给出了解的长时间动力行为.最近,ZHANG等[5]研究了具有退化非局部非线性阻尼项和乘幂型源项的半线性波动方程Dirichlet初边值问题u tt−∆u+σ(∥∇u∥22)|u t|γ−2u t=|u|p−2u,且利用位势井理论得到了能量的衰减估计和有限时刻解的爆破现象.另一方面,关于具有常数阻尼系数和对数非线性项的半线性波动方程研究也有一些新的进展,其主要难点在于对数非线性源项的单调性和符号无法确定.Cazenave和Haraux[6]首次考虑了具有对数非线性项的Schrˆo dinger方程及Klein-Gordon方程Cauchy问题解的存在性与唯一性.之后,ZHANG等[7]考虑了具有弱阻尼的模型Dirichlet初边值问题并得到了问题解的整体存在性及能量的指数衰减估计值.最近,我们在文[8]中研究了具有对数源和强阻尼的半线性波动方程Dirichlet初边值问题u tt−∆u−∆u t=u log|u|2,(1.5)且利用位势井理论和对数Sobolev不等式,得到了问题的整体可解性及能量衰减和无限爆破结果.文[9]的作者将方程(1.5)推广到具有更一般形式对数非线性项情形.LIAN和XU[10]考虑了具有强弱阻尼和对数源项的非线性波动方程Dirichlet初边值问题u tt−∆u−ω∆u t+µu t=u ln|u|,他们利用压缩映射原理,位势井方法及微分不等式技巧,证明了问题的可解性,能量衰减及解的无限爆破现象.此外,关于具有时滞阻尼的板模型问题的最新进展,我们阅读了文[11].综上所述,关于基尔霍夫型弱阻尼项和对数非线性项竞争的半线性波动方程初边值问题(1.1)-(1.3)的研究尚未得到完善.本文中,我们考虑与文[12]中相同意思下的正则解与弱解且有以下主要难点:1)当p>2时无法利用对数-Sobolev不等式来估计对数所在的项;2)局部弱解的适定性无法直接导出,需考虑更强的解的结果;3)分析两个非线性项之间的竞争关系,即基尔霍夫型非线性弱阻尼项与对数非线性项时遇到难度.为了克服这些困难,启发于文[2,13]的思想,从问题正则解的局部适定性出发,利用稠密性理论和紧性理论导出局部弱解的适定性.并通过修正能量泛函技巧将局部解推广到了整体解.同时利用反证技巧,得到具有负初始能量解的有限时刻爆破现象.本文的剩余部分结构如下:第二节,我们将证明问题(1.1)-(1.3)弱解的局部适定性;第三节中,给出当p<γ时弱解的整体适定性.关于p>γ时具有负初始能量解在有限时刻发生爆破的结论,在第四节中导出.整文中,C及C i(i=1,···)在不同表达式中可能表示不同的正常数.同时,记空间H:= {u∈H10(Ω):∆u∈L2(Ω)},且赋予内积(u,v)H:=(∇u,∇v)+(∆u,∆v),其中(·,·)表示L2(Ω)的内积.此外,H1(Ω)中的模取为∥∇u∥2.2.弱解的局部适定性本节中,我们将先证明问题(1.1)-(1.3)正则解的存在唯一性,之后,通过稠密性理论得到局部弱解的存在唯一性.注意到,得到弱解局部适定性之前,先证明正则解的原因在于:证明基尔霍夫型非局部阻尼项的收敛性时,需用正则解的一些结果.第3期杨怡等:一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程513定理2.1假设(H1),(H2)成立且令(u 0,u 1)∈H 10(Ω)×L 2(Ω),则对某些T >0,问题(1.1)-(1.3)存在唯一弱解并满足正则性u ∈C ([0,T ];H 10(Ω))∩C 1([0,T ];L 2(Ω)),u t ∈L γ(0,T ;L γ(Ω)),u tt ∈L 2(0,T ;H −1(Ω)).证我们先分四个步骤来证明问题(1.1)-(1.3)正则解u ∈L ∞([0,T ];H 10(Ω)),u t ∈L ∞(0,T ;H 10(Ω))∩L γ(0,T ;L γ(Ω)),u tt ∈L ∞(0,T ;L 2(Ω))的存在唯一性结论.第一步逼近问题.令{ωj }∞j =1为H 的一组完备正交基且定义有限维子空间V m :=span {ω1,ω2,···,ωm },m ∈N .我们定义近似解u m (x,t ):=∑mj =1b jm (t )ωj (x ),其中u m (x,t )为如下Cauchy 问题的解:∫Ωu mtt ωj d x +∫Ω∇u m ∇ωj d x +σ(∥∇u m ∥22)∫Ω|u mt |γ−2u mt ωj d x =∫Ω|u m |p −2u m log |u m |ωj d x,(2.1)u m (0)=u 0m =m ∑j =1b jm (0)ωj →u 0,在H 中,(2.2)u mt (0)=u 1m =m ∑j =1b jmt (0)ωj →u 1,在L 2(γ−1)(Ω)∩H 10(Ω)中.(2.3)由ODE 标准理论可知,上述Cauchy 问题(2.1)-(2.3)在区间[0,T m ),T m >0上存在唯一解b jm (t ),我们将通过接下来的先验估计将解延伸到[0,T ].第二步先验估计.第一先验估计对(2.1)两边乘b ′jm (t ),关于j =1,2,···,m 求和并关于时间变量从0到t ≤T m 积分可得¯E m (t )+∫t0σ(∥∇u m (s )∥22)∥u ms (s )∥γγd s =¯E m (0)+∫t 0∫Ω|u m (s )|p −2u m (s )log (|u m (s )|)u ms (s )d x d s,(2.4)其中¯Em (t )=12[∥u mt ∥22+∥∇u m ∥22].且由(2.2)(2.3)可得¯Em (0)<+∞.另一方面,经过简单计算,我们易得log |u (x )|<|u (x )|ααa.e.x ∈Ω,∀α>0.(2.5)结合Young 不等式可以导出∫Ω|u m |p −2u m u mt log |u m |d x ≤1α∫Ω|u m |p +α−1u mt d x,≤14ε1α2γ−1γ∫Ω|u m |(p +α−1)γγ−1d x +ε1γ∫Ω|u mt |γd x.(2.6)根据p 和γ的选择,我们可取α满足:若N ≥3,则0<α<2N(γ−1)(N −2)γ+1−p ;若N =1,2,则0<α<+∞.且应用嵌入H 10(Ω) →L γ(p +α−1)γ−1(Ω),(2.6)可改写为∫Ω|u m |p −2u m u mt log |u m |d x ≤C (p +α−1)γγ−1s (γ−1)4ε1α2γ∥∇u m ∥(p +α−1)γγ−12+ε1γ∥u mt ∥γγ,(2.7)其中C s 为最优嵌入常数.将(2.7)代入(2.4),应用(H1)我们导出12[∥u mt ∥22+∥∇u m ∥22]+(b −ε1γ)∫t 0∥u ms (s )∥γγd s514应用数学2022≤¯E m (0)+C (p +α−1)γγ−1s (γ−1)4ε1α2γ∫t∥∇u m (s )∥(p +α−1)γγ−12d s,≤¯Em (0)+C (p +α−1)γγ−1s(γ−1)4ε1α2γ2β∫t0¯E βm (s )d s,其中β=(p +α−1)γγ−1>1.选取ε1使得ε1γ<b,则由非线性Gronwall 不等式可得:存在正常数L 1使得∥u mt ∥22+∥∇u m ∥22+∫t∥u ms (s )∥γγd s ≤L 1,∀m ∈N ,∀t ∈[0,T ].(2.8)因此,由(2.8)我们得到{u m }在L ∞(0,T ;H 10(Ω))中有界,(2.9){u mt }在L ∞(0,T ;L 2(Ω))中有界,(2.10){u mt }在L γ(0,T ;L γ(Ω))中有界.(2.11)第二先验估计首先我们估计∥u mtt (0)∥22.在(2.1)中,取t =0且ωj =u mtt (0),我们有∥u mtt (0)∥22+σ(∥∇u m (0)∥22)∫Ω|u mt (0)|γ−2u mt (0)u mtt (0)d x −∫Ωu mtt (0)∆u m (0)d x =∫Ω|u m (0)|p −2u m (0)u mtt (0)log |u m (0)|d x.(2.12)利用H¨o lder 不等式,我们可得∥u mtt (0)∥22≤[∥∆u 0m ∥22+σ(∥∇u 0m ∥22)∥u 1m ∥2(γ−1)+∥∇u 0m ∥2]∥u mtt (0)∥2,再由σ的连续性及(2.2),(2.3)可得∥u mtt (0)∥2≤L 2,∀m ∈N ,(2.13)其中L 2为与m 无关的常数.紧接着,将(2.1)两边关于时间t 求导并乘b ′′jm (t ),关于j 求和并关于时间变量从0到t ≤T m 积分可得E m (t )+1γ−1∫t 0σ(∥∇u m (s )∥22)∫Ω|u ms (s )|γ−2|u mss |2d x = E m (0)−2∫t 0σ′(∥∇u m (s )∥22)∫Ω∇u m (s )×∇u ms (s )d x ∫Ω|u mt |γ−2u mt u mss (s )d x d s +∫t 0∫Ω|u m (s )|p −3u m (s )u ms (s )u mss (s )d x d s+(p −1)∫t 0∫Ω|u m (s )|p −3u m (s )u ms (s )u mss (s )log |u m (s )|d x d s,(2.14)其中 E m (t ):=12[∥u mtt ∥22+∥∇u mt ∥22].且由(2.2)和(2.3)可得 E m (0)<+∞.结合(2.9)和(H1),(2.14)可重写为E m (t )+1γ−1b ∫t 0∫Ω|u ms |γ−2|u mss |2d x d s ≤ E m (0)+C 1∫t 0∫Ω∇u m (s )×∇u ms (s )d x3d s +I 1+I 2+I 3,(2.15)其中I 1:=−2∫tσ′(∥∇u m (s )∥22)∫Ω∇u m (s )∇u ms (s )d x∫Ω|u ms |γ−2u ms u mss (s )d x d s,第3期杨怡等:一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程515I 2:=∫t 0∫Ω|u m (s )|p −3u m (s )u ms (s )u mss (s )d x d s,I 3:=(p −1)∫t 0∫Ω|u m (s )|p −3u m (s )u ms (s )u mss (s )log |u m (s )|d x d s.下面,我们估计I 1∼I 3.首先,由H¨o lder 不等式可得∫t 0∫Ω|u ms |γ−2u ms u mss (s )d x d s ≤∫t 0∫Ω|u ms |γ−1|u mss |d x d s≤(∫t 0∥u ms (s )∥γγd s)12(∫t 0∫Ω|u ms (s )|γ−2|u mss (s )|2d x d s )12.(2.16)利用Young 不等式,σ′的连续性及(2.9),我们可以估计I 1如下:I 1≤L 121C 3(∫t 0∫Ω|u ms (s )|γ−2|u mss (s )|2d x d s)12∫t 0∫Ω∇u m (s )∇u ms (s )d x d s ≤L 121C 3ε1∫t 0∫Ω|u ms (s )|γ−2|u mss (s )|2d x d s +L 121C 3ε1∫t 0(∫Ω∇u m (s )∇u ms (s )d x )2d s.(2.17)对I 2应用p −22(p −1)+12(p −1)+12=1的H¨o lder 不等式,嵌入H 10(Ω) →L 2(p −1)(Ω)及(2.9),我们导出I 2≤∫t 0(∫Ω|u m (s )|2(p −1)d x )p −22(p −1)(∫Ω|u ms (s )|2(p −1)d x )12(p −1)(∫Ω|u mss (s )|2d x )12d s≤∫t 0∥∇u m (s )∥p −22∥∇u ms (s )∥2∥u mss (s )∥2d s ≤L p −221∫t∥∇u ms (s )∥2∥u mss (s )∥2d s≤12L p −221∫t∥∇u ms (s )∥22+∥u mss (s )∥22d s.(2.18)接下来,应用(2.5)并选取适当α满足:若N ≥3,则α<2N −2N −2−p ;若N =1,2,则0<α<+∞.我们可类似于(2.18)的过程来导出I 3≤p −12αL p +α−221∫t 0∥∇u ms (s )∥22+∥u mss (s )∥22d s.(2.19)将(2.17)-(2.19)代入(2.15)并整理得E m (t )+(b γ−1−L 121C 3ε1)∫t 0∫Ω|u ms (s )|γ−2|u mss (s )|2d x d s≤12L p −221∫t 0∥∇u ms (s )∥22+∥u mss (s )∥22d s + E m (0)+12(p −1αL p +α−221+L p −221)∫t 0∥∇u ms (s )∥22d s +L 121C 3ε1∫t 0(∫Ω∇u m (s )∇u ms (s )d x )2d s +12(p −1αL p +α−221+L p −221)∫t 0∥u mss (s )∥22d s.选取ε1<b (γ−1)L 121C 3,则存在正常数C 1,C 2使得 E m (t )≤C 1+C 2∫t 0E m d s.由非线性Gronwall 不等式可得,存在正常数L 3使得∥u mtt ∥22+∥∇u mt ∥22≤L 3,∀m ∈N ,(2.20)且知{u mt }在L ∞(0,T ;H 10(Ω))中有界,(2.21){u mtt }在L ∞(0,T ;L 2(Ω))中有界.(2.22)516应用数学2022第三步取极限.结合(2.9)-(2.11)即(2.21)-(2.22),存在函数u 及{u m }∞m =1的子序列(方便起见,仍记为{u m }∞m =1)使得u m W ∗−→u,在L ∞(0,T ;H 10(Ω))中,(2.23)u mt W ∗−→u t ,在L ∞(0,T ;H 10(Ω))中,(2.24)u mt W −→u t ,在L γ(0,T ;L γ(Ω))中,(2.25)u mtt W ∗−→u tt ,在L ∞(0,T ;L 2(Ω))中.(2.26)由(2.23),(2.24)及Aubin-Lions 引理可得u m −→u,在C ([0,T ];H 10(Ω))中.(2.27)因此u m −→u,a.e.(x,t )∈Ω×(0,T ].这表明σ(∥∇u m ∥22)−→σ(∥∇u ∥22),a.e.(x,t )∈Ω×(0,T ],|u m |p −2u m log |u m |−→|u |p −2u log |u |,a.e.(x,t )∈Ω×(0,T ].另一方面,由(2.5),(2.9)及嵌入不等式,我们导出∫Ω |u m |p −2u m log |u m | 2d x ≤∫{x ∈Ω||u m |≤1}|u m |p −2u m log |u m | 2d x +1µ∫{x ∈Ω||u m |>1}|u m |2(p −1+µ)d x,≤[1(p −1)e ]2|Ω|+1µ∥u m ∥2(p −1+µ)2(p −1+µ),≤[1(p −1)e ]2|Ω|+21µB 2(p −1+µ)µ−1∥∇u m ∥2(p −1+µ)2≤C 5,(2.28)其中选择合适的正数µ使其满足:若N ≥3,则0<µ<N2(N −2)+1−p ;若N =1,2,则0<µ<+∞,且对0<x <1,应用不等式|x p −1log x |≤(e (p −1))−1,我们得到|u m |p −2u m log |u m |W ∗−→|u |p −2u log |u |,在L ∞(0,T ;L 2(Ω))中.另外,由H ×(L 2(γ−1)(Ω)∩H 10(Ω))在H 10(Ω)×L 2(Ω)中的稠密性可知(u 0m ,u 1m )→(u 0,u 1),在H 10(Ω))×L 2(Ω)中.在(2.1)-(2.3)中取极限可得u tt −∆u +σ(∥∇u ∥22)g (u t )=|u |p −2u log |u |,在u tt ∈L ∞(0,T ;L 2(Ω))中,u (x,0)=u 0(x ),u t (x,0)=u 1(x ).第四步唯一性.令u 1和u 2是问题(1.1)-(1.3)的解并记z :=u 1−u 2,则由(1.1)我们可知z 满足如下的方程:∫Ωz tt ωj d x +∫Ω∇z ∇ωj d x +∫Ωσ(∥∇u 1∥22)|u 1t |γ−2u 1t ωj d x −∫Ωσ(∥∇u 2∥22)|u 2t |γ−2u 2t ωj d x =∫Ω|u 1|p −2u 1log |u 1|ωj d x −∫Ω|u 2|p −2u 2log |u 2|ωj d x,(2.29)用z t 代替上式中的ωj ,我们得到12d d t (∥z t ∥22+∥∇z ∥22)+I 4=I 5+I 6,(2.30)其中I 4=σ(∥∇u 1∥22)∫Ω(|u 1t |γ−2u 1t −|u 2t |γ−2u 2t )z t d x,I 5=−(σ(∥∇u 1∥22)−σ(∥∇u 2∥22))·∫Ω|u 2t |γ−2u 2t z t d x,I 6=∫Ω(|u 1|p −2u 1log |u 1|−|u 2|p −2u 2log |u 2|)z t d x.下面,我们将估计I 4∼I 6.首先,应用平均值定理,存在θ∈(0,1),我们导出(|u 1t |γ−2u 1t −|u 2t |γ−2u 2t )z t =∫1(θu 1t +(1−θ)u 2t )γ−2z t d θz t第3期杨怡等:一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程517=1γ−1(|θu 1t +(1−θ)u 2t |γ−2(θu 1t +(1−θ)u 2t )) 10z t =1γ−1(12(|u 1t |γ−2+|u 2t |γ−2)(u 1t −u 2t )+12(|u 1t |γ−2−|u 2t |γ−2)(u 1t +u 2t ))z t =12(γ−1)((|u 1t |γ−2+|u 2t |γ−2)|z t |2+(|u 1t |γ−2−|u 2t |γ−2)(|u 1t |2−|u 2t |2))≥12(γ−1)((|u 1t |γ−2+|u 2t |γ−2)|z t |2.结合(H1),我们可以对I 4估计如下:I 4≥b 2(γ−1)∫Ω(|u 1t |γ−2+|u 2t |γ−2)|z 2t |d x.(2.31)然后,利用(H1),H¨o lder 不等式和Young 不等式,我们可估计I 5如下:I 5=−(σ(∥∇u 1∥22)−σ(∥∇u 2∥22))∫Ω|u 2t |γ−2u 2t z t d x≤C ∥∇z ∥2∫Ω|u 2t |z t d x +C ∥∇z ∥2∫Ω|u 2t |γ−1z t d x ≤C ∥u 2t ∥2∥∇z ∥2∥z t ∥2+C ε2∥u 2t ∥γγ∥∇z ∥22+ε2∫Ω(|u 1t |γ−2+|u 2t |γ−2)|z t |2d x.(2.32)对I 6,用H¨o lder 不等式得到I 6≤∥G (u 1)−G (u 2)∥2∥z t ∥2,(2.33)其中G (s )=|s |p −2s log |s |.再由平均值定理及(2.5),存在ξ∈(0,1)使得|G (u 1)−G (u 2)|=|G ′(ξu 1+(1−ξ)u 2)z |≤[1+(p −1)log |ξu 1+(1−ξ)u 2|]|ξu 1+(1−ξ)u 2|p −2|z |≤|ξu 1+(1−ξ)u 2|p −2|z |+(p −1)α3|ξu 1+(1−ξ)u 2|p +α3−2|z |≤|u 1+u 2|p −2|z |+(p −1)α3|u 1+u 2|p +α3−2|z |,(2.34)其中选取合适的α3使其满足:若N ≥3,则0<α3≤p −22(p −1)(2N N −2+2(1−p ));若N =1,2,则0<α3<+∞.对(2.34)右边的两项利用H¨o lder 不等式和Sobolev 嵌入不等式,我们分别导出∫Ω|u 1+u 2|2(p −2)|z |2d x ≤C 1(∫Ω|u 1+u 2|2(p −1)d x )p −2p −1(∫Ω|z |2(p −1)d x)1p −1≤C 1[∥u 1∥2(p −1)2(p −1)+∥u 2∥2(p −1)2(p −1)]p−2p −1∥z ∥22(p −1)≤C 2[∥∇u 1∥2(p −1)2+∥∇u 2∥2(p −1)2]p −2p −1∥∇z ∥22(2.35)∫Ω|u 1+u 2|2(p +α3−2)|z |2d x ≤C 3[∥∇u 1∥p ∗2+∥∇u 2∥p ∗2]p −2p −1∥∇z ∥22,(2.36)其中p ∗=2(p −1)+2α3(p −1)p −2.现在,我们把(2.34)-(2.36)代入到(2.33)并整理得到I 6≤C 4{[∥∇u 1∥2(p −1)2+∥∇u 2∥2(p −1)2]p −2p −1+[∥∇u 1∥p ∗2+∥∇u 2∥p ∗2]p −2p −1}(∥∇z ∥22+∥z t ∥22).(2.37)选取ε2<b2(γ−1),将(2.31),(2.32)和(2.37)代入到(2.30)且结合(2.23)得到d d t(∥∇z ∥22+∥z t ∥22)≤C (1+∥u 2t ∥γγ)(∥∇z ∥22+∥z t ∥22).518应用数学2022对上式从0到t 上积分,应用(2.25)和Gronwall 不等式,可知存在正常数L 4使得∥z t ∥22+∥∇z ∥22≤L 4(∥z 1∥22+∥∇z 0∥22),∀m ∈N ,且∥z t ∥22=∥∇z ∥22=0,即可得唯一性的结论.下面,我们应用稠密性理论,从局部正则解u ∈L ∞([0,T ];H 10(Ω)),u t ∈L ∞(0,T ;H 10(Ω))∩L γ(0,T ;L γ(Ω)),u tt ∈L ∞(0,T ;L 2(Ω))的适定性中导出局部解u ∈C ([0,T ];H 10(Ω))∩C 1([0,T ];L 2(Ω)),u t ∈L γ(0,T ;L γ(Ω)),u tt ∈L 2(0,T ;H −1(Ω))的适定性.由H 在H 10(Ω)中稠密,L 2(γ−1)(Ω)∩H 10(Ω)在L 2(Ω)中稠密及条件(u 0,u 1)∈H 10(Ω)×L 2(Ω)易知,存在{u 0η}⊂H 及{u 1η}⊂(L2(γ−1)(Ω)∩H 10(Ω))使得u 0η→u 0,在H 10(Ω)中,u 1η→u 1,在L 2(Ω)中,当η→+∞.(2.38)且对任意η∈N ,问题(1.1)-(1.3)存在以{u 0η,u 1η}为初值的正则解,且满足u η∈L ∞(0,T ;H 10(Ω)),u ηtt ∈L ∞(0,T ;L 2(Ω)),u ηt ∈L ∞(0,T ;H 10(Ω))∩L γ(0,T ;L γ(Ω)).类似于前述的正则解的存在性证明中第一先验估计的导出过程,并令η2≥η1是两个任取的自然数且记z η:=u η2−u η1则易知∥z ηt ∥22+∥∇z η∥22≤C 1(∥z 1η∥22+∥∇z 0η∥22),∀0≤t <+∞.结合(2.38),我们得到z η(0)=u η1(0)−u η2(0)→0,在H 10(Ω)中,z ηt (0)=u tη1(0)−u tη2(0)→0,在L 2(Ω)中,∥z ηt ∥22+∥∇z η∥22→0,且u η→u,在C 0([0,T ];H 10(Ω))中,(2.39)u ηt →u t ,在C 0([0,T ];L 2(Ω))中.(2.40)因此,以上收敛性并结合(2.28)允许我们对问题(1.1)-(1.3)取极限,且得到弱解满足u tt −∆u +σ(∥∇u ∥22)g (u t )=|u |p −2u log |u |,在u tt ∈L 2(0,T ;H −1(Ω)),u (0)=u 0,u t (0)=u 1,x ∈Ω,此外,关于局部弱解的唯一性需要用正则化方法,且可由Visik-Ladyzenskaya 的标准方法来得到.[14]14−16综上所述,我们得到问题(1.1)-(1.3)存在唯一的局部弱解.注2.1局部弱解的唯一性不可用常见的唯一性证明方法的原因在于:对偶积⟨H −1(Ω),L 2(Ω)⟩没有意义.3.弱解的整体适定性本节中,当p <γ时,结合连续性原理,我们得到与第一节中的局部弱解相同正则性的意思下问题(1.1)-(1.3)整体适定性.我们先给出下面引理,将在证明中起到关键作用.引理3.1假设(H1),(H2)成立,对任意的u ∈H 10(Ω)\{0},存在只依赖于Ω的正常数C >1使得∥u ∥s p ≤C (∥∇u ∥22+∥u ∥pp ),其中2≤s ≤p.证当∥u ∥p ≤1时,由Sobolev 嵌入定理可知∥u ∥s p ≤∥u ∥2p ≤C ∥∇u ∥22.同时,当∥u ∥p >1时,我们有∥u ∥s p ≤∥u ∥pp .由此可知,引理3.1成立.定理3.1假设(H1),(H2)成立,p <γ,令(u 0,u 1)∈H 10(Ω)×L 2(Ω),则问题(1.1)-(1.3)存在唯一整体弱解并满足正则性u ∈C ([0,+∞);H 10(Ω))∩C 1([0,+∞);L 2(Ω)),u t ∈L γ(0,+∞;L γ(Ω)),u tt ∈L 2(0,+∞;H −1(Ω)).第3期杨怡等:一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程519证已证明了弱解的局部存在性定理,因此可证明连续性原理(见文[15]).这表明,解的生命跨度T max =∞或T max 有限且满足lim t →T −max ∥u t ∥22+∥∇u ∥22=+∞.然而,用通常的能量泛函E (t ):=12∥u t ∥22+12∥∇u ∥22−1p ∫Ω|u |p log |u |d x +1p 2∥u ∥pp ,来分析对数源项和基尔霍夫型非线性弱阻尼项之间的相互作用遇到困难,故我们引入如下修正能量泛函:F (t )=12∥u t ∥22+12∥∇u ∥22−1p ∫Ω|u |p log |u |d x +1p 2∥u ∥p p +2pµ∥u ∥p +µp +µ,其中µ满足适当条件且我们只需证明F (t )满足指数形式有界,即通常的能量E (t )得到控制.由简单计算,易知E (t )关于时间变量单调递减,即d d tE (t )=−σ(∥∇u ∥22)∥u t ∥γγ≤0.对F (t )直接求导,结合上式,(H1)和Young 不等式,我们导出F ′(t )=−σ(∥∇u ∥22)∥u t ∥γγ+2(p +µ)pµ∫Ω|u |p +µ−2uu t d x ≤−σ(∥∇u ∥22)∥u t ∥γγ+2(p +µ)pµ∫Ω|u |p +µ−1u t d x ≤−b ∥u t ∥γγ+2(p +µ)pµC ε3∥u ∥p +µp +µ+2(p +µ)pµε3∥u t ∥p +µp +µ,其中ε3>0且C ε3是只依赖于ε3的正常数.注意到p <γ,我们可选µ足够小使得p +µ≤γ.因此,可用嵌入L γ →L p +µ可得F ′(t )≤−b ∥u t ∥γγ+2(p +µ)pµC ε3∥u ∥p +µp +µ+C 1ε3∥u t ∥p +µγ,且F ′(t )≤C 2+2(p +µ)pµC ε3∥u ∥p +µp +µ.(3.1)事实上,当∥u t ∥γγ>1时,我们可选取ε3足够小使得−b ∥u t ∥γγ+C 1ε3∥u t ∥p +µγ≤0且有F ′(t )≤2(p +µ)pµC ε3∥u ∥p +µp +µ;当∥u t ∥γγ≤1时,我们易知F ′(t )≤2(p +µ)pµC ε3∥u ∥p +µp +µ+C 1ε3.由此可知(3.1)显然成立.另一方面,由(2.5)可得如下F (t )的估计:F (t )=12∥u t ∥22+12∥∇u ∥22−1p ∫Ω|u |p log |u |d x +1p 2∥u ∥p p+2pµ∥u ∥p +µp +µ>12∥u t ∥22+12∥∇u ∥22+1p 2∥u ∥p p +1pµ∥u ∥p +µp +µ>1pµ∥u ∥p +µp +µ.(3.2)结合(3.1)和(3.2),可得F ′(t )≤C 2+C 4F (t ),其中C 4=2(p +µ)C ε3pµ且F (t )≤(F (0)+C 2C 4)e C 4t .上式与连续性原理结合,即可得整体解的存在性.事实上,我们只需证明:如果T max <+∞,则lim t →T −max ∥u t ∥22+∥∇u ∥22=+∞就可.利用反证技巧,假设上述结论不成立,即T max <+∞且lim t →T −max ∥u t ∥22+∥∇u ∥22<+∞.则存在一个序列{t n ,n =1,2,···}和一个正常数K 使得当n →+∞时t n →T max 且满足∥u t ∥22+∥∇u ∥22<K,n =1,2,···.由前述证明可知,对每一个n ∈N ,初值为u (x,t n )的问题(1.1)-(1.3)的解在[t n ,t n +T ∗]上存在且唯一,其中正常数T ∗依赖于K 但不依赖于n ∈N .因此,对足够大的n ∈N ,我们可得到T max <t n +T ∗.这与T max 是解的最大存在时间矛盾.定理3.1证毕.4.爆破现象本节中,我们利用反证技巧得到,当p >γ时问题(1.1)-(1.3)具有负初始能量的解在有限时刻发生爆破.引理4.1[11]若满足∫Ω|u |plog |u |d x >0.则存在一个只依赖于Ω的正常数C 使得下式成立∥u ∥p p ≤C [∫Ω|u |plog |u |d x +∥∇u ∥22],∀u ∈L p (Ω).(4.1)520应用数学2022现在,我们陈述有限时刻发生爆破结论.定理4.1假设(H1),(H2)成立,p >γ且令(u 0,u 1)∈H 10(Ω)×L 2(Ω),E (0)<0,则问题(1.1)-(1.3)的解在有限时刻发生爆破,即T max <+∞.证利用反证技巧,假设问题(1.1)-(1.3)解整体存在,即T max =+∞.我们引入辅助函数K (t ):=∥u (t )∥22,H (t ):=−E (t ),∀0≤t ≤T 1,其中正常数T 1将在之后给出.且由E (t )的单调性可知:H ′(t )=−E ′(t )≥0,且H (t )≥H (0)=−E (0)>0,∀0≤t ≤T 1.(4.2)我们记G (t ):=∫Ωσ(∥∇u ∥22)u |u t |γ−2u t d x,∀0≤t ≤T 1.现在,由问题解的整体存在性假设知∥∇u ∥22≤C 1,∥u t ∥22≤C 0,∀0≤t ≤T 1.(4.3)并由条件(H1)中σ(s )的连续性,我们有σ(∥∇u ∥22)≤max 0≤∥∇u ∥22≤C 1σ(∥∇u ∥22):=σ1,∀0≤t ≤T 1.(4.4)结合H¨o lder 不等式和Young 不等式,我们导出|G (t )|= ∫Ωσ(∥∇u ∥22)u |u t |γ−2u t d x =σ(∥∇u ∥22)∫Ω|u ||u t |γ−1d x ≤(σ(∥∇u ∥22)∥u t ∥γγη−γγ−1)γ−1γ(σ(∥∇u ∥22)1γ∥u ∥γη)≤γ−1γη−γγ−1σ(∥∇u ∥22)∥u t ∥γγ+1γηγσ(∥∇u ∥22)∥u ∥γγ≤γ−1γη−γ−1γH ′(t )+1γηγ∥u ∥γγσ1,(4.5)其中正常数η将在之后给出.接下来,对K (t )直接求导得到K ′(t )=2∫Ωuu t d x,K ′′(t )=2dd t ∫Ωuu t d x,且令y (t ):=H 1−ξ+γ2γ(t )+ε4∫Ωuu t d x,其中正常数ε4将在之后给出.由p >γ知,我们可取正常数ξ满足γ(p −γ)p 2<ξ<p −γp<1.(4.6)且对y (t )直接求导,并利用(4.5)式,我们有y ′(t )=1−ξ+γ2γH 1−ξ−γ2γH ′(t )+12K ′′(t )ε4=1−ξ+γ2γH 1−ξ−γ2γH ′(t )+ε4∥u t ∥22−ε4∥∇u ∥22+ε4∫Ω|u |p log |u |d x −ε4σ(∥∇u ∥22)∫Ωu |u t |γ−2u t d x,≥[1−ξ+γ2γH 1−ξ−γ2γ(t )−ε4γ−1γηγ−1γ]H ′(t )+ε4∥u t ∥22−ε4∥∇u ∥22+ε4∫Ω|u |plog |u |d x −ε4ηγγσ1∥u ∥γγ.进一步,引入0<a <1,并将上式改写为y ′(t )≥[1−ξ+γ2γ−ε4γ−1γη−γ−1γH −1+ξ+γ2γ(t )]H 1−ξ−γ2γ(t )H ′(t )+ε4p (1−a )H (t )+ε4p (1−a )+22∥u t ∥22+ε4p (1−a )−22∥∇u ∥22+aε4∫Ω|u |plog |u |d x +(1−a )p ∥u ∥p p −ε41γηγH −ξ(t )H ξ(t )σ1∥u ∥γγ.(4.7)第3期杨怡等:一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程521现在,我们断言∫Ω|u |plog |u |d x >0.事实上,由能量函数E (t )的定义和单调性,我们有1p ∫Ω|u |p log |u |d x =12∥u t ∥22+12∥∇u ∥22+1p2∥u ∥pp −E (t )≥12∥u t ∥22+12∥∇u ∥22+1p2∥u ∥p p −E (0)≥0,(4.8)且引理4.1的条件成立.下面,我们估计(4.7)式右端最后一项.由引理4.1可得∥u ∥γγ≤C 2[(∫Ω|u |plog |u |d x )γp +∥∇u ∥2γp 2],并结合(4.2)和Young 不等式,我们导出H ξ(t )∥u ∥γγ≤(1p ∫Ω|u |plog |u |d x )ξ∥u ∥γγ≤C 3[(∫Ω|u |p log |u |d x )γp +ξ+(∫Ω|u |plog |u |d x )ξ∥∇u ∥2γp2]≤C 4[(∫Ω|u |p log |u |d x )γp +ξ+∥∇u ∥22+(∫Ω|u |p log |u |d x )ξpp −γ].(4.9)由(4.6),知γp<γ+pξp≤1,γp <ξp p −γ≤1,且利用(4.9)可得H ξ(t )∥u ∥γγ≤C 5(∫Ω|u |plog |u |d x +∥∇u ∥22).将上式代入到(4.7)并整理得到y ′(t )≥[1−ξ+γ2γ−ε4γ−1γη−γ−1γH −1+ξ+γ2γ(t )]H 1−ξ−γ2γ(t )H ′(t )+ε4p (1−a )H (t )+ε4p (1−a )+22∥u t ∥22+ε4(p (1−a )−22−1γηγH −ξ(t )σ1C 5)∥∇u ∥22+ε4(a −1γηγH −ξ(t )σ1C 5)∫Ω|u |plog |u |d x +ε4(1−a )p ∥u ∥p p .(4.10)此外,应用(2.5),(4.3),H (t )的定义和嵌入不等式,我们得到H (t )=−12∥u t ∥22−12∥∇u ∥22+1p ∫Ω|u |p log |u |d x −1p2∥u ∥p p≤1p 1δ∥u ∥p +δp +δ≤1p 1δB p +δδ∥∇u ∥p +δ2≤1p 1δB p +δδC p +δ21.且上式与(4.2)和(4.10)结合,我们导出y ′(t )≥ 1−ξ+γ2γ−ε4γ−1γη−γ−1γ(1p 1δB p +δδC p +δ21)−1+ξ+γ2γ H 1−ξ−γ2γ(t )H ′(t )+ε4p (1−a )H (t )+ε4p (1−a )+22∥u t ∥22+ε4(p (1−a )−22−1γηγH −ξ(0)σ1C 5)∥∇u ∥22+ε4(a −1γηγH −ξ(0)σ1C 5)∫Ω|u |plog |u |d x +ε4(1−a )p ∥u ∥p p .(4.11)可选取充分小的正数a 使得p (1−a )−22>0,并可选取η充分小,使得p (1−a )−22−1γηγH −ξ(0)σ1C 5>0,且a −1γηγH −ξ(0)σ1C 5>0.固定η和a 之后,我们现在可选取充分小的ε4使得1−ξ+γ2γ−ε4γ−1γη−γ−1γ(1p 1δB p +δδC p +δ21)−1+ξ+γ2γ>0,522应用数学2022(−E(0))1−ξ+γ2γ−C2ε4C1−12ε4C0≥0,(4.12)且H1−ξ+γ2γ(0)+ε4∫Ωu0u1d x>0.于是,由(4.11)知,存在正常数C6使得y′(t)≥C6[H(t)+∥u t∥22+∥∇u∥22+∥u∥p p]≥0,(4.13)且知y(t)在(0,T max)上单调递增且满足y(t)=H1−ξ+γ2γ(t)+ε4K′(t)≥H 1−ξ+γ2γ(0)+ε4K′(0)>0.再由(4.6)及p>γ知,0<ξ<1且令r:=γ1−ξ+γ2,则根据不等式:|a+b|r≤2r−1(|a|r+|b|r),r≥1,以及Young不等式,我们得到y r(t)≤2r−1(H(t)+ε4∥u(t)∥r2∥u t(t)∥r2)≤C7(H(t)+∥u(t)∥11−ξγ2+∥u t(t)∥22).(4.14)进一步,我们估计(4.14)右端第二项.由(4.6)知,1−ξ>γp,且应用如下不等式:xτ≤(1+1m)(m+x),x≥0,0≤τ≤1,m>0,并取x=∥u(t)∥p2,τ=γ(1−ξ)p<1,m=H(0),得到∥u(t)∥γ1−ξ2≤(1+1H(0))(H(0)+∥u(t)∥p2)≤C8(H(t)+∥u(t)∥pp).(4.15)将(4.15)代入到(4.14)中,我们有y r(t)≤C9(H(t)+∥u(t)∥pp+∥u t(t)∥22),由上式与(4.13),我们得到y′(t)≥C10y r(t),∀0≤t≤T1.(4.16)且对上式关于时间变量从0到t积分,我们有y(t)≥(y(0)1−r−C10(r−1)t)−1r−1,0≤t≤T1.我们可取T1≥T∗=y(0)1−rC10(r−1),且由r:=γ1−ξ+γ2>1和y(0)>0知,在[0,T1]上存在有限时刻T2≤T∗使得y(t)满足limt→T−2y(t)→+∞.(4.17)但是,我们断言:存在正常数使得C11≤y(t)≤C12,∀0≤t≤T1.(4.18)事实上,由(4.3),(4.4),Poincare不等式,H¨o lder不等式和Young不等式和嵌入不等式,我们导出y(t)=H1−ξ+γ2γ(t)+ε4∫Ωuu t d x=(∫tσ(∥∇u∥22)∥uτ∥γγdτ−E(u(0)))1−ξ+γ2γ+ε4∫Ωuu t d x≤(∫tσ1∥uτ∥22dτ−E(u(0)))1−ξ+γ2γ+C2ε4∥∇u∥22+12ε4∥u t∥22≤C11,且y(t)≥(∫tb∥uτ∥γγdτ−E(u(0)))1−ξ+γ2γ−C2ε4∥∇u∥22−12ε4∥u t∥22,≥(−E(0))1−ξ+γ2γ−12ε4C0−C2ε4C1≥C12>0,其中最后一式由(4.12)得到.显然,(4.17)与(4.18)产生矛盾且我们知T max<+∞.定理4.1证毕.参考文献:[1]LANGE H,PERLA M G.Rates of decay of a nonlocal beam equation[J].Differ.Integral Equ.,1997,10(6):1075-1092.第3期杨怡等:一类具有基尔霍夫型弱阻尼和对数非线性项的半线性波动方程523[2]GEORGIEV V,TODOROVA G.Existence of solutions of the wave equation with nonlinear dampingand source terms[J].J.Differential Equations,1994,109:295-308.[3]MESSAOUD S A.Blow up in a nonlinearly damped wave equation[J].Math.Nachr.,2001,231:105-111.[4]JORGE SILVA M A,NARCISO V.Long-time dynamics for a class of extensible beams with nonlocalnonlinear damping[J].Evol.Equ.Control Theory,2017,6(3):437-470.[5]ZHANG H W,LI D H,HU Q Y.Asymptotic stability and blow-up for the wave equation withdegenerate nonlocal nonlinear damping and source terms[J/OL].Appl.Anal.,2020[2021-06-01].DOI:10.1080/00036811.2020.1836354.[6]CAZENAVE T,HARAUX A.´Equations d’´e volution avec non-lin´e arit´e logarithmiqu[J].Ann.Fac.Sci.Toulouse Math.,1980,2(1):21-51.[7]ZHANG H W,LIU G W,HU Q Y.Exponential decay of energy for a logarithmic wave equation[J].J.Partial Diff.Equ.,2015,28(3):269-277.[8]MA L W,FANG Z B.Energy decay estimates and infinite blow-up phenomena for a strongly dampedsemilinear wave equation with logarithmic nonlinear source[J].Math.Method App.Sci.,2018,41(7): 2639-2653.[9]DI H F,SHANG Z F,SONG Z F.Initial boundary value problem for a class of strongly dampedsemilinear wave equations with logarithmic nonlinearity[J].Nonlinear Anal.Real World Appl.,2020, 51:102968.[10]LIAN W,XU R Z.Global well-posedness of nonlinear wave equation with weak and strong damp-ingterms and logarithmic source term[J].Nonlinear Anal.,2020,9:613-632.[11]KAFINI M,MESSAOUDI S.Local existence and blow up of solutions to a logarithmic nonlinearwave equation with delay[J].Appl.Anal.,2020,99(3):530-547.[12]LIU L,WANG M.Global existence and blow-up of solutions for some hyperbolic systems withdamping and source terms[J].Nonlinear Anal.,2006,64(1):69-91.[13]MESSAOUDI S A.Blow up and global existence in a nonlinear viscoelastic wave equation[J].Math.Nachr.,2003,260:58-66.[14]LIONS J L.Quelques Méthodes de R´e solution des Probl´e mes aux Limites Non Lin´eéairess[M].Paris:Dunod,1969.[15]SEGAL I.Nonlinear semigroups[J].Ann.Math.,1963,78:339-364.On a Semilinear Wave Equation with Kirchhoff-type Weak Damping Terms and Logarithmic NonlinearityYANG Yi,FANG Zhongbo(School of Mathematical Sciences,Ocean University of China,Qingdao266100,China) Abstract:Under homogeneous Dirichlet conditions,the well-posedness and qualitative propertiesfor a semilinear wave equation with Kirchhoff-type weak damping terms and logarithmic nonlinearity were considered.By improving the well-posedness for regular solution and density argument,a local existence of weak solutions was proved.Meanwhile,based on modified energy technique and contradiction argument, a global existence with p<γand thefinite time blow-up with p>γwere also established.Key words:Semilinear wave equation;Kirchhoff-type weak damping;Logarithmic nonlinearity; Global existence;Blow-up。
阻尼波动方程
![阻尼波动方程](https://img.taocdn.com/s3/m/3e65b7c2951ea76e58fafab069dc5022aaea4687.png)
阻尼波动方程引言阻尼波动方程是描述阻尼振动的一种数学模型。
在物理学中,许多系统都会受到阻尼的影响,例如弹簧、摆钟、电路等。
通过研究阻尼波动方程,我们可以深入了解这些系统的行为和性质。
阻尼振动简介在自由振动中,一个物体在无外力作用下沿某个轴向做周期性的来回运动。
然而,在实际情况中,摩擦和空气阻力等因素会导致振动能量逐渐减小,最终停止。
这种现象称为阻尼振动。
阻尼振动可以分为三种类型:无阻尼、欠阻尼和过阻尼。
无阻尼振动指的是没有任何外界因素干扰下的理想情况;欠阻尼指的是存在一定的摩擦或空气阻力,但仍能保持周期性运动;过阻尼则是指存在较大的摩擦或空气阻力,导致振荡停止时需要较长时间。
阻尼波动方程在描述阻尼振动时,我们可以使用阻尼波动方程。
阻尼波动方程可以写成如下形式:m * d^2x/dt^2 + c * dx/dt + k * x = 0其中,m是物体的质量,c是阻尼系数,k是弹性系数,x是物体的位移。
该方程包含三个部分:惯性项(m * d2x/dt2),阻尼项(c * dx/dt)和弹性项(k* x)。
这些部分分别描述了物体的惯性、阻尼和弹性特性。
解析解和数值解对于简单的阻尼波动问题,我们可以使用解析方法求得精确解。
然而,对于复杂的系统或非线性情况,通常需要使用数值方法求解。
解析解在某些情况下,可以通过假设特定的位移函数形式,并代入阻尼波动方程中求解。
例如,在无阻尼情况下,可假设位移函数为sin(ωt),然后代入方程进行计算。
通过这种方式可以得到精确解。
数值解当无法找到解析解时,我们可以使用数值方法求解阻尼波动方程。
常用的数值方法包括欧拉法、改进欧拉法和四阶龙格-库塔法等。
这些方法通过将时间连续问题离散化为时间离散的问题,并迭代求解来逼近精确解。
应用领域阻尼波动方程在物理学、工程学和生物学等领域有着广泛的应用。
物理学在物理学中,阻尼振动是研究弹簧、摆钟、声波等现象的基础。
通过解析或数值求解阻尼波动方程,可以预测物体的运动轨迹、频率和振幅等特性。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
(2)
设
u(t )
为问题(1)-(3)的任意解,
若
u0
(x)∈
Bc δ0
,则由
J
(u0
)
≤
E(0)
<
d
(δ
0
)
知 u0
∈ Vδ 0
,即 I (u0 ) <
0
由推论
4
知 u(t)∈Vδ2
,从而 u(t)∈
Bc δ2
即 u(t) 在球
Bδ 2
之外。
定 理 6 令 p , ui (x) (i = 0,1) , e,δi (i = 1,2) , 如 定 理 1 所 定 义 , 若 初 始 能 量 满 足
e<d,
δ1 < δ 2 是方程 d (δ ) = e 的两根,则:
(1) 若 I (u0 ) > 0 或 ∇u0 = 0 则问题(1)-(3)的全部解 u ∈Wδ ;
(2) 若 I (u0 ) < 0 则(1)-(3)的所有解 u 满足 u ∈Vδ 。 定理 3 令 p , ui (x) (i = 0,1) , e,δi (i = 1,2) ,如定理 1 所定义, 0 < E(0) ≤ e ,则任意 δ ∈ (δ1,δ 2 ) 有Wδ ,Vδ 在问题(1)-(3)的流之下是不变的。
u(x,0) = u0 (x) ut (x,0) = u1(x)
x ∈Ω
(5)
u ∂Ω = 0
t ≥0
(6)
其中 Ω ⊂ R n (n = 1,2,3)是有界域。为了得到此问题的整体强解,webb对 f (u) 加了4个条
件,该假设的基本模型方程为 utt − α∆ut − ∆u = − u p−1u 即源项为负,外力与位移方向相
∇u0 ≠ 0 ,由
知 J (u(t0 )) ≠ d(δ 0 ) 。
J (u(t0 )) ≤ E(t0 ) < E(0) < d(δ )
(7)
若 Jδ (u(t0 )) = 0 且 ∇u0 ≠ 0 ,由引理 5 知 J (u(t0 )) ≥ d (δ 0 ) 与(7)式矛盾。从而
u(t )∈Wδ ,δ ∈ (δ1,δ 2 ) 。 (2)令 u(t) 是问题(1)-(3)的一个解,初始能量 E(0) = e ,且 I (u0 ) < 0 。由 I (u0 ) < 0
3. 解的不变集合
定理
1
设
p
满足(H),
u0
(x)
∈
H
1 0
(Ω)
,
u1
(x)
∈
L2
(Ω)
,假设
0
<
e
=
E(0) < d
,
δ1 < δ 2 是方程 d (δ ) = e 的两根,则:
(1) 若 I (u0 ) > 0, ∇u0 = 0 ,则(1)-(3)的所有解 u 满足 u ∈Wδ , δ ∈ (δ1,δ 2 );
(2) 若 I (u0 ) < 0 则(1)-(3)的所有解 u 满足 u ∈Vδ ,δ ∈ (δ1,δ 2 ) 。
证明 (1)令 u(t) 是问题(1)-(3)的一个解,且初始能量 E(0) = e 且 I (u0 ) > 0 或
-2-
∇u0 = 0 , T 为 u(t) 的 存 在 时 间 , 则 有 u0 ∈Wδ 。 事 实 上 , 若 I (u0 ) > 0 ,
知 Jδ0 (u0 ) < 0 且 ∇u0 ≠ 0 ,从而 Jδ (u0 ) 于区间 (δ1,δ 2 )不变号 Jδ (u0 ) < 0,δ ∈ (δ1,δ 2 ), 且 J (u0 ) ≤ E(0) < d(δ )故 u0 ∈Vδ ,δ ∈ (δ1,δ 2 ) 。
下证 u(t)∈Vδ ,δ ∈ (δ1,δ 2 )及 0 < t < T 。 反之 ∃ t0 ∈ (0,T ), δ ∈ (δ1,δ 2 ) 使 u(t0 )∈ ∂Vδ 即 Jδ (u(t0 )) = 0 或 J (u(t0 )) = d (δ )。由 (7)式知 J (u(t0 )) = d (δ )不可能。假设 t0 为第一个使 Jδ (u(t0 )) = 0 的点,则 Jδ (u(t)) < 0,0 ≤ t < t0 由 J (u) ≤ E(0) < d (δ )及引理 2 知
问题(1)-(3)所有解都在球 Bδ2 (或 ∂Bδ2 上)之外,若 u0
x
∈
Bc δ0
。
证明(1)设 u(t) 为问题(1)(- 3)的任意解,若 u0 (x)∈ Bδ0 ,则由 J (u0 ) ≤ E(0) < d (δ 0 )
知 u0 ∈Wδ0 ,即 I (u0 ) > 0 或 ∇u0 = 0 。由推论 4 知 u(t )∈Wδ1 ,从而 u(t )∈ Bδ1 。
0
p +1
定义
-1-
J (u) = 1 ∇u 2 − a u p+1
2
p + 1 p+1
及位势井
I (u ) = ∇u 2 − a u p+1 p +1
{ } W =
u
∈
H
1 0
(Ω)
I
(u)
>
0,
0
<
J
(u)
<
d
U {0}
其中
d
=
inf
J (u) ,其中 u
推论 4 令 p , ui (x) (i = 0,1) , e,δi (i = 1,2) ,如定理 1 所定义, 0 < E(0) ≤ e 则:
(1) 若 I (u0 ) > 0 或 ∇u0 = 0 ,问题(1)-(3)的所有解都属于Wδ1 ;
-3-
(2) 若 I (u0 ) < 0 则问题(1)-(3)的所有解都属于Vδ2 。 证明 令 u(t) 是问题(1)-(3)的任意解,且初始能量 E(0)满足 0 < E(0) ≤ e ,T 为解
从推论 4,引理 3 及 4 可得定理 5。
定理 5 令 p , ui (x) (i = 0,1) , e,δi (i = 1,2) ,如定理 1 所定义, 0 < E(0) ≤ e 则:
( ) (1) 问题(1)-(3)所有解都在球 Bδ1 (或 ∂Bδ1 上)之内,若 u0 x ∈ Bδ0 ;
( ) (2)
反。随着势井理论 [2] 的发展,源项为正的情形也得到了很好的解决 [3] 。并且在位势井理论
基础上随之产生位势井族理论 [4] ,该理论的发展使得该方程又凸显出以往未曾得到的解的
性质.
真空隔离性质由刘亚成教授在文献 [4] 中首先提出,用以解决半线性波动方程的初边值问题。
本文将利用这种理论研究强阻尼非线性波动方程解的真空隔离性质。其中 f (u) 满足条件
0
<
E(0)
≤
e
则
∀δ
∈
(δ1,δ 2
),
Bδ
和
Bc δ
在问题(1)-(3)的流之下是不变的。
( ) 证明
由引理 2 知当 J (u) ≤ E(0) < d (δ ) 时 u ∈ Bδ
Bc δ
⇔ u ∈Wδ (Vδ ) 由定理
3
知
Bδ
,
Bc δ
在问题(1)-(3)的流之下是不变的。
4. 真空隔离现象
证 明 若 u0 (x)∈Wδ ⇒ J (u) < d (δ ), Jδ (u0 ) > 0 或 ∇u0 = 0 。 于 是 有 I (u0 ) > 0 或
∇u0 = 0 。由定理 1 知 u(t)∈Wδ ,δ ∈ (δ1,δ 2 ) ,t ∈ (0,T )。即Wδ 和Vδ 在问题(1)-(3)
的流之下是不变的。
0 < J (u0 ) ≤ E(0) = e < d ⇒ Jδ (u0 ) 于 (δ1,δ 2 ) 不变号,由 Jδ0 (u0 ) > 0 知 Jδ (u0 ) > 0 , δ ∈ (δ1,δ 2 ) 从而 u0 ∈Wδ 。若 ∇u0 = 0 则 u0 ∈Wδ 。
下证 u(t )∈Wδ ,δ ∈ (δ1,δ 2 ) , 0 < t < T 。 反 证 法 , 若 不 然 ∃t0 ,δ 0 使 u(t0 )∈ ∂Wδ0 , 即 J (u(t0 )) = d (δ 0 ), Jδ (u(t0 )) = 0 且
∈
H
1 0
(Ω
),
∇u
≠
0, I (u)
=
0
,
d
的值与文[3]中相同
d
=
1 αC*α
,
α
=
2( p +1)
p −1
,
C*
为
H
1 0
(Ω
)
空间嵌入
L
p+1
(Ω
)
的嵌入常数,即
C*
= sup
u p +1
∇u
。
进一步对于δ ∈ (0,1) ,定义
J
δ
(u
)
=
δ 2
∇u
2
−
a p +1
u
p +1 p +1
2
隔离性质。
关键词:非线性波动方程,位势井族,初边值,强阻尼,真空隔离
中图分类号:O175.29
1. 引言
关于强阻尼非线性波动方程的研究始于1980年,由webb [1] 提出并研究了如下强阻尼非 线性波动方程的初边值问题
utt − α∆ut − ∆u = f (u)