本文中,我们考虑如下具有非线性阻尼项和源项的四阶波动方程
在缺少耗散项,色散项,k -拉普拉斯算子 (β=γ=δ=0), 且取非线性阻尼项 g(ut)=a|ut|m−1ut,非线性源项 f(u)=b|u|p−1u 的情况下,模型 (1.1) 成为如下的非线性双曲方程
如果不考虑 k -拉普拉斯算子和非线性弱阻尼项的影响效应 (γ=0,g(ut)=0),让β=δ=1 且维数 n=1 时,模型 (1.1) 成为如下的四阶色散耗散波动方程
如果不考虑 k -拉普拉斯算子和耗散项的影响效应 (γ=β=0), 且阻尼项为线性阻尼,g(ut)=ut,当n=1 时,模型 (1.1)成为如下四阶非线性色散波方程
当 γ=0,β=δ=1,g(ut)=ut 和 f(u)=|u|p−1u 时, 模型 (1.1)成为如下的四阶耗散色散非线性波方程
在文献[23]中,Messaoudi 考虑了如下的带有非线性阻尼项和非线性源项的拟线性双曲方程
据作者所知,对带有强耗散项 △utt,k -拉普拉斯算子, 非线性弱阻尼项 g(ut) 和非线性源项 f(u)的拟线性波动方程的初边值问题, 迄今还没有任何研究结果. 本文考虑如下具有耗散项、k -拉普拉斯算子、 非线性弱阻尼和源项的四阶波动方程的初边值问题
本文的结构如下: 第2节,我们引入了定理证明所需要的一些数学符号, 基本定义和重要引理; 第3节,得到了方程整体弱解的存在性; 在第4节中, 证明了拥有负初始能量时,方程的解在有限时间内是爆破的.
这一节,我们引入定理证明所需要的数学符号,基本定义和重要引理.
用Lp(Ω)(1≤p<∞) 表示通常的定义在Ω 上的 p 次 Lebesgue 可积的可测函数全体构成的空间,其上的范数定义为 ‖u‖p=(∫Ω|u|pdx)1p. L∞(Ω) 表示定义在 Ω 上的本性有界可测函数全体构成的空间,其上范数定义为 ‖u‖∞=ess⋅supx∈Ω|u(x)|. 在Hilbert空间 L2(Ω) 中定义内积 (u,v)L2=∫Ωuvdx,∀u,v∈L2(Ω). 用 Wm,p(Ω) 表示通常的Sobolev 空间,其上范数定义为 ‖u‖m,p=(∑|α|≤m‖Dαu‖pp)1p. 特别的 p=2 时,Wm,p(Ω) 简记为 Hm(Ω), 其上的内积和范数分别为 (u,v)Hm=∑|α|≤m(Dαu,Dαv)L2 和 ‖u‖Hm=(∑|α|≤m‖Dαu‖22)12. △=∂2∂x21+⋯+∂2∂x2n, ∇=(∂∂x1,⋯,∂∂xn).
C∞0(Ω) 是由在 Ω 中具有紧支集的全体 C∞ 函数组成的; Wm,p0(Ω) 可以看成是 C∞0(Ω) 在 Wm,p(Ω) 中的闭包; Wm,p0(Ω) 是 Sobolev 空间 Wm,p(Ω) 的一个完备的子空间.
为了获得本文的主要结果,我们首先引入如下能量泛函
简单起见,我们在 [0,T) 区间上给出了问题 (1.9) 弱解的定义,其中 T 可以理解为无穷大或解存在的时间区间的上限.
定义2.1 称函数 u(x,t) 是初边值问题 (1.9) 的一个弱解,若 u∈L∞(0,T;W1,k0(Ω)), ut∈L∞(0,T;H10(Ω))∩Lm(QT) 并且满足如下的条件
(i)~ 对 ∀ v∈W1,∞(0,T;W1,k0(Ω))∩Lm(QT),有
(ii)~ u(x,0)=u0(x) 在 W1,k0(Ω) 中,ut(x,0)=u1 在 H10(Ω)∩Lm(Ω) 中.
引理2.1 假设 p,m,k>2,a,b>0,若 u(x,t) 是初边值问题 (1.9) 的一个弱解, 则 E(t) 是一个非增的函数,即
证 在 (1.9) 式两边同时乘以 ut,并在 Ω 上积分, 通过分部积分法和类似于文献 [24] 的推导,我们可证明对方程的任意正则解, (2.3),(2.4) 式是成立的. 采用稠密性方法,可将以上结论推广到对任意的弱解也是成立的. 关于稠密性方法的相关理论知识,可见参考文献 [25, 26].
接下来,我们定义如下的一些泛函和集合 I(t)=I(u(t))=‖∇u‖kk+‖∇u‖22−b‖u‖pp, J(t)=J(u(t))=1k‖∇u‖kk+12‖∇u‖22−bp‖u‖pp, E(t)=E(u(t))=12‖ut‖22+12‖∇ut‖22+J(t) 和 W={u∈W1,k0(Ω)|I(t)>0}∪{0}.
注2.1 泛函 J(u(t)) 和 I(u(t)) 在空间 W1,k0(Ω) 中是连续的, 关于其连续性的严格证明可以直接从它们的定义中推导得到.
引理2.2 \label{le2.2} 假设a,b>0,p>k>2,并且满足
证 因为 I(u0)>0,由 I(u(t)) 的连续性,存在时间 t1 使得 I(u(t))≥0 对所有的 t∈[0,t1] 成立. 这样,我们可看到
我们将文献 [27] 第二章中的引理 6.1 稍微改动,即得到下面的引理.
引理2.3 \label{le2.2} 假设 m,k>2,并假设函数 u 满足 u∈L∞(0,∞;W1,k0(Ω)),ut∈L∞(0,∞; H10(Ω))∩Lm(Q∞),u(x,0)=u0, ut(x,0)=u1. 进一步假设 utt−△u−△ut−△utt=g, g∈L2(0,T;H−1(Ω))+Lm′(QT), 则对 a.e. t∈[0,∞),接下来的不等式成立
注2.2 引理2.3的证明过程完全类似于文献 [27] 第二章引理 6.1的证明.
接着,令 {wj} 为空间 W1,k0(Ω)∩Lm(Ω) 中的一组基函数,通过 Galerkin 方法,我们构造初边值问题 (1.9) 的近似解 ul=ul(x,t),它具有如下的形式
关于非线性常微分方程组初值问题的一般结果保证了问题 (2.13) 的解在区间 [0,tl] 上存在. 在第3节中得到的一系列先验估计表明了问题 (2.13) 解的时间区间可以扩展到 [0,∞). 更进一步,通过先验估计,也可证明近似解 ul(x,t) 极限逼近值即为初边值问题 (1.9) 的整体弱解.
本节我们证明初边值问题 (1.9) 解的整体存在性.
定理3.1 假设 a,b>0,m,k>2,p>k 且 (2.5) 式成立. 假设 (u0,u1)∈W×H10(Ω) 并满足 (2.6) 式, 则问题 (1.9) 存在整体弱解 u∈L∞(0,∞;W1,k0(Ω)), ut∈L∞(0,∞;H10(Ω))∩Lm(Q∞), 并且对所有 t≥0,有 u∈W.
证 在 (2.13) 式两边同乘以 djl(t)′,j=1,⋯,l,并对 j 求和,得
利用 Sobolev 嵌入不等式,可推得
于是由 (3.3)-(3.6) 式,我们有
因此,我们可得到 ‖ult‖2H1, ‖∇ul‖22, ‖∇ul‖kk, ∫t0‖∇ult‖22dτ, ∫t0‖ult‖mmdτ 是有界的且都可以被一个与 l 和 t 无关的常数所控制. 因而推出 ul∈L∞(0,∞;W1,k0(Ω)), ult∈L∞(0,∞;H10(Ω))∩Lm(Q∞). 注意到 (a|ult|m−2ult,ult)=a‖ult‖mm, (b|ul|p−2ul,ul)=b‖ul‖pp 和 (‖∇ul‖k−2∇ul,∇ul)=‖∇ul‖kk, 所以可推出 a|ult|m−2ult∈Lm′(Q∞), b|ul|p−2u∈L∞(0,∞;Lp′(Ω)), ‖∇ul‖k−2∇ul∈L∞(0,∞;Lk′(Ω)). 这样,通过以上的讨论,由致密性原理,可知存在 {ul} 的一个子序列(在这仍用 {ul} 来表示其子序列) 和函数 u(x,t),ξ,η,γ,使得 ul(x,t)⟶u(x,t) 在 L∞(0,∞;W1,k0(Ω)) 弱∗收敛, l⟶∞, ul(x,t)⟶u(x,t) 在 Q∞=Ω×(0,∞)上几乎处处收敛, l⟶∞, ∇ul(x,t)⟶∇u(x,t) 在 Q∞=Ω×(0,∞)上几乎处处收敛, l⟶∞, ult(x,t)⟶ut(x,t) 在 L∞(0,∞;H10(Ω)) 弱∗收敛, l⟶∞, ult(x,t)⟶ut(x,t) 在 Lm(Q∞) 弱收敛, l⟶∞, a|ult|m−2ult⟶ξ 在 Lm′(Q∞) 弱收敛, l⟶∞, b|ul|p−2ul⟶η 在 L∞(0,∞;Lp′(Ω)) 弱∗收敛, l⟶∞, b|ul|p−2ul⟶η 在 Q∞=Ω×(0,∞)上几乎处处收敛, l⟶∞, |∇ul|k−2∇ul⟶γ 在 L∞(0,∞;Lk′(Ω)) 弱∗收敛, l⟶∞, |∇ul|k−2∇ul⟶γ 在 Q∞=Ω×(0,∞)上几乎处处收敛, l⟶∞, 根据文献 [27] 第一章引理 1.3,推出 b|ul|p−2ul⟶b|u|p−2u=η, |∇ul|k−2∇ul⟶|∇u|k−2∇u=γ. 考虑到基函数 {wj} 在函数空间 W1,k0(Ω)∩Lm(Ω) 的稠密性,所以我们可选择函数 v∈W1,∞(0,∞;W1,k0(Ω))∩Lm(Q∞),它具有形式 v(x,t)=∞∑j=1kj(t)wj, 其中 kj(t) 是给定的已知函数. 在 (2.13) 式两边乘以 kj(t),j=1,2,⋯, 并对 j 求和,我们可得
接下来,我们需要证明的是
在上式中,当 l⟶∞时 取极限得到
本节在初始能量 E(0)<0 的条件下,我们证明了问题 (1.9) 在有限的时间内爆破解的存在性.
定理4.1 假设 a,b>0,p,m,k>2,m<k<p 且 (2.5) 式成立. 若初始能量
证 通过与引理 2.1 证明中类似的讨论,我们可得到 E′(t)+‖∇ut‖22=−a‖ut‖mm≤0, 令 H(t)=−E(t),则有
利用 Sobolev 嵌入不等式和 (4.2) 式,可得到估计
注4.1 由爆破时间的估计式 (4.35),我们可以看出 F(0) 的取值越大, 那么问题 (1.9) 的解爆破发生越快.
注4.2 文献 [28, 29] 中,在低初始能量 (E(0)<d) 的情况下, 已经研究和讨论了具有强阻尼项的非线性双曲方程解的适定性. 文献 [5] 中,Xu 在临界的初始能量 (E(0)=d) 情况下研究了半线性的双曲方程的初边值问题, 证明了方程的解在有限的时间内是爆破的. 而本文所研究的问题 (1.9), 不管是低初始能量 (E(0)<d) 还是临界的初始能量 (E(0)=d), 都没有类似于文献 [5, 28, 29] 中的结果. 在初始能量 (E(0)<d) 或 (E(0)=d) 的情况下问题 (1.9) 的适定性的研究,我们将在其他文章中讨论.