Processing math: 33%

数学物理学报, 2020, 40(5): 1319-1332 doi:

论文

常微分不等式及其在半线性波动方程的应用

黄守军,, 孟希望,

Improved Ordinary Differential Inequality and Its Application to Semilinear Wave Equations

Huang Shoujun,, Meng Xiwang,

通讯作者: 黄守军, E-mail: sjhuang@ahnu.edu.cn

收稿日期: 2019-03-29  

基金资助: 国家自然科学基金.  11301006
安徽省自然科学基金.  1408085MA01

Received: 2019-03-29  

Fund supported: the NSFC.  11301006
the NSF of Anhui Province.  1408085MA01

作者简介 About authors

孟希望,E-mail:3572861950@qq.com , E-mail:3572861950@qq.com

摘要

该文首先得到两类变系数的常微分不等式的爆破结果,可视为文献[3,定理3.1]的推广.其次,作为改进的常微分不等式的一个应用,考虑具有尺度不变阻尼项的半线性波动方程的柯西问题,给予初值合理假设,得到当μ>11<p<1+2n时解的生命跨度的上界估计.该结果的证明方法主要来自于文献[11].

关键词: 常微分不等式 ; 半线性波动方程 ; 生命跨度

Abstract

In this paper, we first derive some blow-up results for two ordinary differential inequalities with variable coefficients, which are the generalizations of Theorem 3.1 in Li and Zhou[3]. Second, as an application of the improved ordinary differential inequality, we consider the Cauchy problem for the semilinear wave equation with scale-invariant damping and deduce the upper bound of the lifespan for the case μ>1 and 1<p<1+2n under some suitable assumptions for the initial data. The method for the latter result is due to Lai and Zhou[11].

Keywords: Ordinary differential inequality ; Semilinear wave equation ; Lifespan

PDF (388KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

黄守军, 孟希望. 常微分不等式及其在半线性波动方程的应用. 数学物理学报[J], 2020, 40(5): 1319-1332 doi:

Huang Shoujun, Meng Xiwang. Improved Ordinary Differential Inequality and Its Application to Semilinear Wave Equations. Acta Mathematica Scientia[J], 2020, 40(5): 1319-1332 doi:

1 引言

本文讨论如下具有尺度不变耗散和幂次方非线性项的半线性波动方程的柯西问题

{vττΔyv+μ1+τvτ=|v|p,τ>0,yRn,v(0,y)=εv0(y),vτ(0,y)=εv1(y),yRn,
(1.1)

其中p>1, μ>0是常数, ε>0是衡量初值大小的参数.该模型常称为尺度不变的,因为相应的线性波动方程在尺度变换˜v=v(λ(1+τ)1,λy),λ>0下是不变的.方程(1.1)可用来描述具有摩擦和热传导的波的传播,在磁流体力学、粘弹性力学等领域也有很多应用[1-2].

柯西问题(1.1)的一个基本问题是确定临界值pc(n),使得当p>pc(n)时, (1.1)存在小初值的整体解,而当1<ppc(n)时,解发生爆破.据我们所知,该问题至今仍然没有得到完全解决.下面我们给出一些有关的研究工作.

柯西问题(1.1)已有大量研究结果,参见文献[4-6]及其参考文献.容易看出,该方程的解强烈地依赖于参数μ的值.特别地,当0<μ<1时(此时阻尼项称为non-effective),解的渐近行为与相应的齐次波动方程的解的性质有关[4, 7-8]. Wirth[13-14]给出相应的线性波动方程阻尼的分类. Lai, Takamura和Wakasa[4]猜测方程解的行为依赖于μ的临界值为μ=1.文献[4]的主要工作为:对于pF(n)p<p0(n+2μ),其中μ(0,μ0(n)),在初值作合理假设后,柯西问题(1.1)解的生命跨度的上界为ε2p(p1)γ(n+2μ;p),其中pF(n):=1+2n是Fujita指标, μ0(n):=n2+n+2n+2, p0(n+2μ)是方程γ(n+2μ;p)=0的正根,这里γ(n+2μ;p)由二次多项式γ(n;p):=2+(n+1)p(n1)p2来定义.

后来, Ikeda和Sobajima[5]重新考虑柯西问题(1.1),并研究超Fujita指标情形时解的爆破现象.文献[5]的结果可总结为:对于pF(n)pp0(n+μ),这里p0(n)是Strauss指标(即方程γ(n;p)=0的正根),如果具有紧支撑的初值是光滑的,且为非负,那么解将在有限时间内发生爆破,解的生命跨度的上界可以得到.可以看出,当解发生爆破时, (p,μ)的范围比文献[4]中对应的范围要稍微大一些,详细结果参见文献[5].

另一方面,对于模型(1.1),利用Zhang[9]的测试函数方法, Wakasugi[15]得到1<ppF(n) (其中μ>1),以及1<ppF(n+μ1) (其中0<μ<1)时的爆破结果.生命跨度的上界由下式给出

Cεp12n(p1),其中1<p<pF(n),μ1;
(1.2)

Cεp12(n+μ1)(p1)其中1<p<pF(n+μ1),0<μ<1,
(1.3)

这里C是一个与ε无关的正常数.

现在我们给出如下阻尼波动方程有关的结果

{vττΔyv+vτ=|v|p,τ>0,yRn,v(0,y)=εv0(y),vτ(0,y)=εv1(y),yRn.
(1.4)

当空间维数n=1,2, p属于(1,pF(n)]时, Li和Zhou[3]得到上述方程解的生命跨度的上界sharp估计.对于整数的p以及n1时,解的生命跨度的下界估计由Li[16]给出,从而对于n=1,2情形解的生命跨度的sharpness得到证明.临界指标p=pF(n)由Todorova和Yordanov[17]发现. Nishihara[18]给出三维情形解的生命跨度的上界的sharp估计.对于n4的高维情形,生命跨度的sharpness最近才被Lai和Zhou[11]所证明.文献[11]中证明的关键之处在于利用经典热核的半群性质得到一个微分不等式,而该微分不等式的爆破结果已由文献[3]给出.这方面的其他相关工作,可以参见文献[11-12]以及其参考文献.

现在我们回到柯西问题(1.1).正如Takamura[19]所指出, Kato引理在研究半线性波动方程时起着非常重要的作用.同时,我们发现Li和Zhou[3]中的常微分不等式的有关结果在Lai和Zhou[11]证明(1.4)解的生命跨度的sharpness时是一个重要的工具.受到Takamura[19]和Lai和Zhou[11]工作的启发,本文我们将Li和Zhou[3]中有关常微分不等式的结果推广到变系数的情形,然后利用所得结果去研究变系数的半线性波动方程柯西问题(1.1)解的爆破现象.

本文主要考虑μ>1的情形.此时, Wirth[13-14]指出阻尼项在相应的线性波动方程中是有效的(effective).因此,类似于Lai和Zhou[11]的方法,选取热核作为测试函数去研究柯西问题(1.1)的解的爆破将是可行的.关于柯西问题(1.1)的主要结果如下定理.

定理1.1 设μ>11<p<pF(n)=1+2n.令初值满足v0(y),v1(y)H1(Rn)×L2(Rn)supp(v0,v1){yRn||y|μμ+1}.此外,假设

v0(y)0,Rn[(μ1)v0(y)+v1(y)]dy>0.
(1.5)

v(τ,y)是柯西问题(1.1)在[0,T(ε))上的解,那么我们可以得到如下关于生命跨度的估计

T(ε)Cεp12n(p1),
(1.6)

这里C是一个与ε无关的正常数.

注1.1 尽管初值的假设(1.5)与文献[15]中的不同,但生命跨度的上界估计(1.6)与(1.2)是一致的.

注1.2 对于μ=1情形,我们相信本文改进的常微分不等式的有关结果仍然适用.然而,利用本文选取热核作为测试函数的方法去得到柯西问题(1.1)当0<μ<1时解的爆破结果,看起来不是非常显然.

注1.3 对于情形p=pF(n)=1+2n,本文的方法将失效,因为(3.28)式右端的最后一项将变为正数,相应的后续讨论将无法进行.这个情况与Lai和Zhou[11]不同,文献[11]研究了临界情形p=pF(n).之所以方法失效,推测起来原因是此时的柯西问题(1.1)中的方程是变系数的,而不是常系数的.

本文结构如下.下一节我们给出两个改进的常微分不等式的爆破结果.在第3节,利用该推广的常微分不等式的爆破结果研究具有尺度不变的半线性波动方程,进而得到柯西问题(3.5)的爆破结果.

2 改进的常微分不等式

本节我们将建立两个常微分不等式,它们改进了Li和Zhou[3]相关的结果.文献[3]的结果在研究常系数的半线性波动方程时非常有用[3, 11].下面的讨论方法来自于Li和Zhou[3].

定理2.1 设I(t)满足

a(t)I(t)+b(t)I(t)C0I1+α(t)(1+t)β
(2.1)

I(0)>0,I(0)0,
(2.2)

这里α>0,β0, C0>0都是常数, a(t)b(t)是光滑的正函数.假设存在常数m12,C1>0C2>0,使得

a(t)0,C1<a(t),0<b(t)C2a(t)m,t0.
(2.3)

那么,当0β1, I=I(t)必在有限时间内爆破.而且,如果I(0)=ε,其中ε>0是小参数,那么I=I(t)的生命跨度T(ε)具有如下的上界估计

T(ε){exp(Aεα),β=1,Bεα1β,0β<1,
(2.4)

这里AB是与ε无关的正常数.

 类似于文献[3]中定理3.1的证明,我们仅需要证明I(t)必在有限时间内爆破,从而生命跨度估计(2.4)可由文献[3]的讨论容易得到.

类似文献[3],不失一般性,我们假设I(0)>0.下面我们仅考虑情形β=0,而剩余情形β=10<β<1可作类似讨论.对于现在的情形, (2.1)式可写为

a(t)I(t)+b(t)I(t)C0I1+α(t).
(2.5)

我们声称I(t)必在有限时间内发生爆破.为此,我们只需要证明:对充分小的η>0,记I1(t)是如下柯西问题的解

{I1(t)=ηa(t)mI1(t)1+α2,t=0:I1=I1(0),
(2.6)

这里

0<I1(0)<I(0),
(2.7)

在解的存在范围内满足

I(t)>I1(t).
(2.8)

事实上,对充分小的η>0,成立

I1(0)=ηa(0)mI1+α21(0)<I(0)
(2.9)

I1(t)=η2a2m(t)(1+α2)I1+α1mηa(t)m1I1+α21(t)a(t).

于是,注意到I1(t)I1(0)以及(2.3)式中的假设,对充分小的η>0,成立

a(t)I1(t)+b(t)I1(t)=(η2a(t)2m1(1+α2)+η(b(t)ma(t))a(t)m1Iα21(t))I1+α1(t)(η2C12m1(1+α2)+ηb(t)a(t)mIα21(t))I1+α1(t)(η2C12m1(1+α2)+ηC2Iα21(0))I1+α1(t)C0I1+α1(t).
(2.10)

由(2.5), (2.10), (2.7)和(2.9)式以及文献[3]中的引理3.1,我们得到

I(t)>I1(t),t0.

因此,利用(2.7)式,我们得到(2.8)式.容易看出I1(t)必在有限时间内爆破,从而上述论断成立,也就是说I(t)必在有限时间内爆破.证毕.

注2.1 如果在(2.1)式中令a(t)b(t)1,那么上述结果可以回到Li和Zhou[3]的结果.

注2.2 如果在(2.3)式中取m=12,那么a(t)具有正下界的假设可以去掉,见(2.10)式的第一个等式.

注2.3 公式(2.3)中的假设a(t)0是为了保证文献[3]中的Lemma 3.2成立,而该引理对于β=10<β<1情形的证明具有非常重要的作用,见文献[3]中的定理3.1的讨论.

推论2.1 设I(t)满足

(1+t)γ1I(t)+(1+t)γ2I(t)C0I1+α(t)(1+t)β
(2.11)

I(0)>0,I(0)0,
(2.12)

这里γ10,γ20,α>0,β0C0>0是常数.假设

γ1γ2β.
(2.13)

那么,当0β<1时, I=I(t)必在有限时间内爆破.此外,如果I(0)=ε,其中ε>0是小参数,那么I=I(t)的生命跨度T(ε)具有如下的上界估计

T(ε)Cεα1β,
(2.14)

这里C是与ε无关的正常数.

 尽管证明方法类似于文献[3]的定理3.1,但是需要指出证明是非平凡的.此处我们仅仅考虑情形0<β<1,因为情形β=0的证明更简单.

J(t)=I((1+t)1/(1β)1).
(2.15)

于是,由(2.11)–(2.12)式,我们得到

(1β)(1+t)(γ1β)/(1β)J(t)+(1+t)γ2/(1β)J(t)β(1+t)γ111βJ(t)C01βJ1+α,
(2.16)

以及

J(0)>0,J(0)0.
(2.17)

由文献[3,推论3.1],在(2.11)式中我们仍然可以得到I(t)>0.于是,利用如下关系I(τ)=(1β)(1+t)β1βJ(t),这里τ=(1+t)11β1,我们可以看出J(t)>0.从而,公式(2.16)化为

(1β)(1+t)(γ1β)/(1β)J(t)+(1+t)γ2/(1β)J(t)C01βJ1+α.
(2.18)

对(2.18)和(2.17)式,再次利用文献[3]中推论3.1,我们可以得到

J(t)>0,t>0.
(2.19)

a(t)=(1β)(1+t)(γ1β)/(1β),b(t)=(1+t)γ2/(1β),

并在柯西问题(2.6)中假设m=1.利用关键条件(2.13),我们得到

a(t)1β>0.
(2.20)

注意到

b(t)a(t)=(1+t)(γ2γ1+β)/(1β)1β,
(2.21)

并利用(2.20)和(2.21)式,当γ2γ1+β0时, (2.3)式中的不等式必定成立.而条件γ2γ1+β0可由假设条件(2.13)来保证.因此,类似于定理的证明,我们得到J(t)必定在有限时间内爆破,也就是说, I(t)必在有限时间内爆破.最后,我们得到生命跨度估计(2.14).证毕.

注2.4 注意到条件(2.13),我们可以发现上述结果不可能通过令γ1=γ2=0回到文献[3]中定理3.1.

3 在半线性波动方程的应用

本节我们将应用上述改进的微分不等式的有关结果,去讨论具有尺度不变阻尼项的柯西问题(1.1)解的爆破,并证明定理1.1.本文我们仅考虑μ>1的情形.类似于文献[11]的方法,我们将利用经典的热核作为测试函数来讨论柯西问题(1.1).

n维热核为

E(t,x)=1(4πt)n2e|x|24t.
(3.1)

热核E(t,x)具有半群性质.下面的结果来自于文献[11].

引理3.1 热核满足

RnE(t,x)dx=1
(3.2)

E(t,x)E(s,x)=E(t+s,x),
(3.3)

这里表示关于x的卷积.

定理1.1的证明  为证明该定理,引入

t=12μ((1+τ)21),y=μ1+μx.
(3.4)

在新坐标(t,x)下,柯西问题(1.1)可以写成如下的形式

{utΔxu=μ1+μ|u|p2μt+1μ(1+μ)utt,t>0,xRn,u(0,y)=εf(x):=εv0(μ1+μx),xRn,ut(0,y)=εg(x):=εμv1(μ1+μx),xRn.
(3.5)

利用Duhamel原理,柯西问题(3.5)的解可写为

u(t,x)=εE(t)f+μ1+μt0E(tτ)|u(τ)|pdτt0E(tτ)2μτ+1μ(1+μ)uττdτ.
(3.6)

由分部积分、初值以及引理3.1,经计算可知

t0E(tτ)2μτ+1μ(1+μ)uττdτ=1μ(1+μ)(RnE(0,xy)(2μt+1)ut(t,y)dyRnE(t,xy)ut(0,y)dyRnt0τuτ(E(tτ,xy)(2μτ+1))dτdy)=1μ(1+μ)((2μt+1)ut(t,x)εE(t)gRnt0τuτ(E(tτ,xy)(2μτ+1))dτdy).
(3.7)

将(3.7)式带入(3.6)式,我们得到

u+2μt+1μ(1+μ)ut1μ(1+μ)Rnt0τuτ(E(tτ,xy)(2μτ+1))dτdy=εE(t)f+εμ(1+μ)E(t)g+μ1+μt0E(tτ)|u(τ)|pdτ.
(3.8)

注意到

Rnt0τuE(tτ,xy)dτdy=Rnu(t,y)E(0,xy)dyRnu(0,y)E(t,xy)dyRnt0u(τ,y)τE(tτ,xy)dτdy=u(t,x)εE(t)fRnt0uτE(tτ,xy)dτdy,

(3.8)式可改写为

μ(μ1)u+(2μt+1)utRnt0(2μτ+1)τuτE(tτ,xy)dτdy+2μRnt0u(τ,y)τE(tτ,xy)dτdy=μ(μ1)εE(t)f+εE(t)g+μ2t0E(tτ)|u|pdτ.
(3.9)

在(3.9)式两端同时乘以 (4\pi(1+t))^{\frac{n}{2}}E(t+1) ,并在 {{\Bbb R}} ^n 作积分,我们得到

\begin{eqnarray} &&\mu(\mu-1)\int_{{{\Bbb R}} ^n}e^{-\frac{|x|^2}{4(t+1)}}u(t, x){\rm d}x+(2\mu t+1) \int_{{{\Bbb R}} ^n}e^{-\frac{|x|^2}{4(t+1)}}u_t(t, x){\rm d}x\\ &&+2\mu(4\pi(t+1))^{\frac{n}{2}}\int_{{{\Bbb R}} ^n}\int_0^t\int_{{{\Bbb R}} ^n}E(t+1, x)\partial_{\tau} E(t-\tau, x-y)u(\tau, y){\rm d}y{\rm d}\tau {\rm d}x\\ &&-(4\pi(t+1))^{\frac{n}{2}}\int_{{{\Bbb R}} ^n}\int_0^t\int_{{{\Bbb R}} ^n}(2\mu\tau+1)E(t+1, x)\partial_{\tau} E(t-\tau, x-y)\partial_{\tau}u(\tau, y){\rm d}y{\rm d}\tau {\rm d}x\\ & = &\varepsilon (4\pi(t+1))^{\frac{n}{2}}\int_{{{\Bbb R}} ^n}E(2t+1, x)[\mu(\mu-1)f(x)+g(x)]{\rm d}x\\ &&+\mu^2 (4\pi(t+1))^{\frac{n}{2}}\int_0^t\int_{{{\Bbb R}} ^n}E(2t+1-\tau, x)|u|^p(\tau, x){\rm d}x {\rm d}\tau. \end{eqnarray}
(3.10)

引入

\begin{equation} G(t) = \int_{{{\Bbb R}} ^n}e^{-\frac{|x|^2}{4(t+1)}}u(t, x){\rm d}x, \end{equation}
(3.11)

\begin{equation} F(t) = \left(\int_{{{\Bbb R}} ^n}e^{-\frac{|x|^2}{4(t+1)}}|u(t, x)|^p{\rm d}x\right)^{\frac{1}{p}}(t+1)^{\frac{n(p-1)}{2p}}, \end{equation}
(3.12)

\begin{equation} A(t) = \int_{{{\Bbb R}} ^n}e^{-\frac{|x|^2}{4(t+1)}}u_t(t, x){\rm d}x, \end{equation}
(3.13)

\begin{eqnarray} B(t)& = &2\mu(4\pi(t+1))^{\frac{n}{2}}\int_{{{\Bbb R}} ^n}\int_0^t\int_{{{\Bbb R}} ^n}E(t+1, x)\partial_{\tau} E(t-\tau, x-y)u(\tau, y){\rm d}y{\rm d}\tau {\rm d}x\\&& -(4\pi(t+1))^{\frac{n}{2}}\int_{{{\Bbb R}} ^n}\int_0^t\int_{{{\Bbb R}} ^n}(2\mu\tau+1)E(t+1, x)\\ &&\times \partial_{\tau} E(t-\tau, x-y)\partial_{\tau}u(\tau, y){\rm d}y{\rm d}\tau {\rm d}x \end{eqnarray}
(3.14)

\begin{equation} D(t) = (4\pi(t+1))^{\frac{n}{2}}\int_0^t\int_{{{\Bbb R}} ^n}E(2t+1-\tau, x)|u|^p(\tau, x){\rm d}x {\rm d}\tau. \end{equation}
(3.15)

于是,类似于文献[11],我们得到

\begin{equation} G(t) \lesssim F(t), \quad A(t)\lesssim G'(t)+F(t)(1+t)^{-1} \end{equation}
(3.16)

\begin{equation} D(t)\gtrsim \int_0^t\frac{F^p(\tau)}{(1+\tau)^{\frac{n(p-1)}{2}}}{\rm d}\tau. \end{equation}
(3.17)

结合(3.10)–(3.17)式得到

\begin{eqnarray} &&\mu(\mu-1)G(t)+(2\mu t+1)G'(t)+\frac{2\mu t+1}{t+1}F(t)+B(t)\\&\gtrsim &\varepsilon \left(\frac{t+1}{2t+1}\right)^{\frac{n}{2}}\int_{{{\Bbb R}} ^n}e^{-\frac{|x|^2}{4(2t+1)}}\left[\mu(\mu-1)f(x)+g(x)\right]{\rm d}x+\mu^2 D(t). \end{eqnarray}
(3.18)

我们现在估计 B(t) .

\begin{equation} B(t) = B_1(t)+B_2(t), \end{equation}
(3.19)

这里

\begin{eqnarray} B_1(t)& = &-\left(4\pi(t+1)\right)^{\frac{n}{2}}\left\{(2\mu t+1)\int_{{{\Bbb R}} ^n}2\pi n \left(4\pi(t+1)\right)^{-\frac{n}{2}-1}e^{-\frac{|x|^2}{4(t+1)}}u(t, x){\rm d}x\right.\\ &&-(2\mu t+1)\int_{{{\Bbb R}} ^n} \left(4\pi(t+1)\right)^{-\frac{n}{2}}\frac{|x|^2}{4(t+1)^2}e^{-\frac{|x|^2}{4(t+1)}}u(t, x){\rm d}x\\ &&-\int_{{{\Bbb R}} ^n}2\pi n \left(4\pi(2t+1)\right)^{-\frac{n}{2}-1}e^{-\frac{|x|^2}{4(2t+1)}}u(0, x){\rm d}x\\ &&+\int_{{{\Bbb R}} ^n} \left(4\pi(2t+1)\right)^{-\frac{n}{2}}\frac{|x|^2}{4(2t+1)^2}e^{-\frac{|x|^2}{4(2t+1)}}u(0, x){\rm d}x\\ &&\left.-4\mu\int_0^t\int_{{{\Bbb R}} ^n}E_{\tau}(2t+1-\tau, x)u(\tau, x){\rm d}x{\rm d}\tau\right\}\\ &\equiv& B_{11}(t)+B_{12}(t)+B_{13}(t)+B_{14}(t)+B_{15}(t)+B_{16}(t), \end{eqnarray}
(3.20)

具体写出来为

B_{11} = -\left(4\pi(t+1)\right)^{\frac{n}{2}}(2\mu t+1)\int_{{{\Bbb R}} ^n}2\pi n \left(4\pi(t+1)\right)^{-\frac{n}{2}-1}e^{-\frac{|x|^2}{4(t+1)}}u(t, x){\rm d}x,

B_{12} = \left(4\pi(t+1)\right)^{\frac{n}{2}}(2\mu t+1)\int_{{{\Bbb R}} ^n} \left(4\pi(t+1)\right)^{-\frac{n}{2}}\frac{|x|^2}{4(t+1)^2}e^{-\frac{|x|^2}{4(t+1)}}u(t, x){\rm d}x,

B_{13} = \left(4\pi(t+1)\right)^{\frac{n}{2}}\int_{{{\Bbb R}} ^n}2\pi n \left(4\pi(2t+1)\right)^{-\frac{n}{2}-1}e^{-\frac{|x|^2}{4(2t+1)}}u(0, x){\rm d}x,

B_{14} = \left(\frac{t+1}{2t+1}\right)^{\frac{n}{2}}\int_{{{\Bbb R}} ^n} \frac{|x|^2}{4(2t+1)^2}e^{-\frac{|x|^2}{4(2t+1)}}u(0, x){\rm d}x,

B_{15} = 2\mu\left(4\pi(t+1)\right)^{\frac{n}{2}}\int_0^t\int_{{{\Bbb R}} ^n}2\pi n\left(4\pi(2t+1-\tau)\right)^{-\frac{n}{2}-1} e^{-\frac{|x|^2}{4(2t+1-\tau)}}u(\tau, x){\rm d}x{\rm d}\tau

B_{16} = -2\mu\left(4\pi(t+1)\right)^{\frac{n}{2}}\int_0^t\int_{{{\Bbb R}} ^n}\left(4\pi(2t+1-\tau)\right)^{-\frac{n}{2}} \frac{|x|^2}{4(2t+1-\tau)^2}e^{-\frac{|x|^2}{4(2t+1-\tau)}}u(\tau, x){\rm d}x{\rm d}\tau,

同时

\begin{eqnarray} B_2(t)& = &\left(4\pi(t+1)\right)^{\frac{n}{2}}\int_0^t(2\mu \tau+1)\int_{{{\Bbb R}} ^n}\partial_{\tau\tau} E(2t+1-\tau)u(\tau, x){\rm d}x{\rm d}\tau{}\\ & = &\left(4\pi(t+1)\right)^{\frac{n}{2}}\int_0^t(2\mu\tau+1)\int_{{{\Bbb R}} ^n}\left\{8\pi^2 n (1+\frac{n}{2})\left(4\pi(2t+1-\tau)\right)^{-\frac{n}{2}-2}e^{-\frac{|x|^2}{4(2t+1-\tau)}}u(\tau, x)\right.{}\\ &&-2n\pi\left(4\pi(2t+1-\tau)\right)^{-\frac{n}{2}-1}\frac{|x|^2}{4(2t+1-\tau)^2}e^{-\frac{|x|^2}{4(2t+1-\tau)}}u(\tau, x){}\\ && +\left(4\pi(2t+1-\tau)\right)^{-\frac{n}{2}}\frac{|x|^4}{16(2t+1-\tau)^4}e^{-\frac{|x|^2}{4(2t+1-\tau)}}u(\tau, x){}\\ &&\left.-(2+\frac{n}{2})\frac14 \left(4\pi(2t+1-\tau)\right)^{-\frac{n}{2}}\frac{|x|^2}{4(2t+1-\tau)^3}e^{-\frac{|x|^2} {4(2t+1-\tau)}}u(\tau, x)\right\}{\rm d}x{\rm d}\tau{}\\ &\equiv &B_{21}(t)+B_{22}(t)+B_{23}(t)+B_{24}(t). \end{eqnarray}
(3.21)

由(3.16)式,直接计算可得

\begin{equation} B_{11}(t), \;B_{12}(t)\lesssim \frac{2\mu t+1}{t+1}F(t), \end{equation}
(3.22)

\begin{equation} B_{13}(t) = \varepsilon\frac{n}{2}\frac{1}{2t+1}\left(\frac{t+1}{2t+1}\right)^{\frac{n}{2}} \int_{{{\Bbb R}} ^n}e^{-\frac{|x|^2}{4(2t+1)}}f(x){\rm d}x \end{equation}
(3.23)

\begin{equation} B_{14}(t) = -\varepsilon\frac{1}{2t+1}\left(\frac{t+1}{2t+1}\right)^{\frac{n}{2}} \int_{{{\Bbb R}} ^n}e^{-\frac{|x|^2}{4(2t+1)}}\frac{|x|^2}{4(2t+1)^2}f(x){\rm d}x. \end{equation}
(3.24)

类似于文献[11],由Young不等式,我们得到

\begin{equation} B_{15}(t), \;B_{16}(t)\lesssim \frac{\mu}{6}D(t)+C_0(t+1)^{\frac{n(p-1)-2}{2(p-1)}}, \end{equation}
(3.25)

这里 C_0 是依赖于 \pi, n, p \mu 的正常数.此外,由Hölder不等式可得

\begin{eqnarray} B_{21}(t)&\lesssim& \left(4\pi(t+1)\right)^{\frac{n}{2}}(t+1)^{-2} \int_0^t(2\mu\tau+1)\int_{{{\Bbb R}} ^n}E(2t+1-\tau, x)u(\tau, x){\rm d}x{\rm d}\tau\\ &\lesssim & \left(4\pi(t+1)\right)^{\frac{n}{2}}(t+1)^{-2} \left(\int_0^t\int_{{{\Bbb R}} ^n}E(2t+1-\tau, x)|u|^p{\rm d}x{\rm d}\tau\right)^{\frac{1}{p}}\\ &&\cdot\left(\int_0^t\int_{{{\Bbb R}} ^n}(2\mu \tau+1)^{1-\frac{1}{p}}E(2t+1-\tau, x){\rm d}x{\rm d}\tau\right)^{1-\frac{1}{p}}\\ &\lesssim& \left(4\pi(t+1)\right)^{\frac{n}{2}} \left(\int_0^t\int_{{{\Bbb R}} ^n}E(2t+1-\tau, x)|u|^p{\rm d}x{\rm d}\tau\right)^{\frac{1}{p}}\cdot (2\mu t+1)(t+1)^{1-\frac{1}{p}-2}\\ &\lesssim& \frac{\mu}{6}D(t)+C_0(t+1)^{1+\frac{n}{2}-\frac{p}{p-1}}\\ &\lesssim&\frac{\mu}{6}D(t)+C_0(t+1)^{\frac{n(p-1)-2}{2(p-1)}}, \end{eqnarray}
(3.26)

上式我们利用了Young不等式,这里 C_0 是仅依赖于 \pi, n, p \mu 的正常数.类似地,我们得到

\begin{equation} B_{22}(t), \;B_{23}(t), \;B_{24}(t)\lesssim \frac{\mu}{6}D(t)+C_0(t+1)^{\frac{n(p-1)-2}{2(p-1)}}. \end{equation}
(3.27)

结合(3.19)–(3.27)式,由(3.18)式可以推得

\begin{eqnarray} &&\mu(\mu-1)G(t)+(2\mu t+1)G'(t)+\frac{2\mu t+1}{t+1}F(t)\\ &\gtrsim& \mu(\mu-1)D(t)+\varepsilon\left(\frac{t+1}{2t+1}\right)^{\frac{n}{2}}\int_{{{\Bbb R}} ^n} e^{-\frac{|x|^2}{4(2t+1)}}\frac{|x|^2}{4(2t+1)^2}f(x){\rm d}x\\ &&+\varepsilon \left(\frac{t+1}{2t+1}\right)^{\frac{n}{2}}\int_{{{\Bbb R}} ^n} e^{-\frac{|x|^2}{4(2t+1)}}[\mu(\mu-1)f(x)+g(x)]{\rm d}x\\ &&-\varepsilon \frac{n}{2(2t+1)}\left(\frac{t+1}{2t+1}\right)^{\frac{n}{2}}\int_{{{\Bbb R}} ^n} e^{-\frac{|x|^2}{4(2t+1)}}f(x){\rm d}x-C_0(t+1)^{\frac{n(p-1)-2}{2(p-1)}}, \end{eqnarray}
(3.28)

其中 C_0 是仅依赖于 \pi, n, p \mu 的正常数.由(3.18)式,我们将(3.28)式改写为

\begin{equation} \mu(\mu-1)G(t)+(2\mu t+1)G'(t)+\frac{2\mu t+1}{t+1}F(t)\gtrsim \mu(\mu-1)\int_0^t\frac{F^p(\tau)} {(1+\tau)^{\frac{n(p-1)}{2}}}{\rm d}\tau+\varepsilon H(t)-I(t), \end{equation}
(3.29)

这里

\begin{eqnarray} H(t)& = &\left(\frac{t+1}{2t+1}\right)^{\frac{n}{2}}\int_{{{\Bbb R}} ^n} e^{-\frac{|x|^2}{4(2t+1)}}\frac{|x|^2}{4(2t+1)^2}f(x){\rm d}x\\ &&+ \left(\frac{t+1}{2t+1}\right)^{\frac{n}{2}}\int_{{{\Bbb R}} ^n} e^{-\frac{|x|^2}{4(2t+1)}}[\mu(\mu-1)f(x)+g(x)]{\rm d}x\\ &&- \frac{n}{2(2t+1)}\left(\frac{t+1}{2t+1}\right)^{\frac{n}{2}}\int_{{{\Bbb R}} ^n} e^{-\frac{|x|^2}{4(2t+1)}}f(x){\rm d}x, \end{eqnarray}
(3.30)

\begin{equation} I(t) = C_0(t+1)^{\frac{n(p-1)-2}{2(p-1)}}. \end{equation}
(3.31)

在(3.29)式两端同乘以 (2\mu t+1)^{\frac{\mu-3}{2}} ,所得方程在 [0, t] 上积分可得

\begin{eqnarray} &&(2\mu t+1)^{\frac{\mu-1}{2}}G(t)+\int_0^t\frac{(2\mu \tau+1)^{\frac{\mu-1}{2}}}{1+\tau}F(\tau){\rm d}\tau\\ &\gtrsim & G(0)+\mu(\mu-1)\int_0^t(2\mu s+1)^{\frac{\mu-3}{2}}\int_0^s\frac{F^p(\tau)}{(1+\tau)^{\frac{n(p-1)}{2}}}{\rm d}\tau {\rm d}s\\ &&+\int_0^t(2\mu \tau+1)^{\frac{\mu-3}{2}}[\varepsilon H(\tau)-I( \tau)]{\rm d}\tau. \end{eqnarray}
(3.32)

注意到(1.5)式,我们有

G(0) = \varepsilon\int_{{{\Bbb R}} ^n}e^{-\frac{|x|^2}{4}}f(x){\rm d}x\geq0.

于是由(3.16)式,我们可以从(3.32)式推导得到

\begin{eqnarray} &&(2\mu t+1)^{\frac{\mu-1}{2}}F(t)+\int_0^t\frac{(2\mu \tau+1)^{\frac{\mu-1}{2}}}{1+\tau}F(\tau){\rm d}\tau\\ &\gtrsim & \int_0^t\left[(2\mu t+1)^{\frac{\mu-1}{2}}-(2\mu \tau+1)^{\frac{\mu-1}{2}}\right]\frac{F^p(\tau)}{(1+\tau)^{\frac{n(p-1)}{2}}}{\rm d}\tau+J(t), \end{eqnarray}
(3.33)

其中

\begin{equation} J(t) = \int_0^t(2\mu \tau+1)^{\frac{\mu-3}{2}}[\varepsilon H(\tau)-I( \tau)]{\rm d}\tau. \end{equation}
(3.34)

由假设(1.5),我们得到

\begin{equation} \lim\limits_{t\rightarrow\infty}H(t) = H_0\equiv 2^{-\frac{n}{2}}\int_{{{\Bbb R}} ^n}\left[\mu(\mu-1)f(x)+g(x)\right]>0, \end{equation}
(3.35)

也就是说,存在时间 t_1 > 0 使得

H(t)>\frac{H_0}{2}>0, \;\;\forall\; t>t_1.

由于我们仅考虑次Fujita指标情形,即 1 < p < 1+\frac{2}{n} ,于是 \frac{n(p-1)-2}{2(p-1)} < 0 .从而存在时间 t_2 > 0 使得

I(t)\leq\frac{H_0}{4}\varepsilon, \;\;\forall\;t>t_2.

t_3 = \max\{t_1, t_2\} , { } M_H = \max_{[0, \infty)}|H(t)| < \infty 以及 { } M_I = \max_{[0, \infty)}I(t) > 0 .于是

\begin{eqnarray} J(t)& = &\int_0^{t_3}(2\mu\tau+1)^{\frac{\mu-3}{2}}\left[\varepsilon H(\tau)-I(\tau)\right]{\rm d}\tau +\int_{t_3}^t(2\mu\tau+1)^{\frac{\mu-3}{2}}\left[\varepsilon H(\tau)-I(\tau)\right]{\rm d}\tau\\ &\gtrsim&\int_0^{t_3}(2\mu\tau+1)^{\frac{\mu-3}{2}}\left[-\varepsilon M_H-M_I\right]+\frac{H_0}{4}\varepsilon\int_{t_3}^t(2\mu\tau+1)^{\frac{\mu-3}{2}}{\rm d}\tau\\ & = &\frac{H_0\varepsilon}{8\mu(\mu-1)}(2\mu t+1)^{\frac{\mu-1}{2}}+\bigg[\frac{H_0\varepsilon}{8\mu(\mu-1)}(2\mu t+1)^{\frac{\mu-1}{2}}{}\\ &&-\frac{1}{\mu(\mu-1)}\left(\frac{H_0}{4}\varepsilon+M_H+M_I\right) (2\mu t_3+1)^{\frac{\mu-1}{2}}\bigg]+\frac{1}{\mu(\mu-1)}(\varepsilon M_H+M_I).{\qquad} \end{eqnarray}
(3.36)

由于 \mu > 1 ,如果取 t 充分大,即 t > t_4(> t_3), 那么(3.36)式中的最后两项是正的.故

\begin{equation} J(t)\gtrsim \varepsilon (2\mu t+1)^{\frac{\mu-1}{2}}, \quad\forall\;t>t_4. \end{equation}
(3.37)

利用(3.37)式,由(3.33)式可得

\begin{eqnarray} &&\frac{\rm d}{{\rm d}t}\left[(1+t)\int_0^t\frac{(2\mu \tau+1)^{\frac{\mu-1}{2}}}{1+\tau}F(\tau){\rm d}\tau\right] {}\\ &\gtrsim & \varepsilon (2\mu t+1)^{\frac{\mu-1}{2}}+\int_0^t\left[(2\mu t+1)^{\frac{\mu-1}{2}}-(2\mu \tau+1)^{\frac{\mu-1}{2}}\right]\frac{F^p(\tau)}{(1+\tau)^{\frac{n(p-1)}{2}}}{\rm d}\tau, \;\forall \;t>t_4.{\qquad} \end{eqnarray}
(3.38)

\alpha(t) = (1+t)\int_0^t\frac{(2\mu \tau+1)^{\frac{\mu-1}{2}}}{1+\tau}F(\tau){\rm d}\tau,

\beta(t) = \int_{t_0}^t\int_0^s\left[(2\mu s+1)^{\frac{\mu-1}{2}}-(2\mu \tau+1)^{\frac{\mu-1}{2}}\right]\frac{F^p(\tau)}{(1+\tau)^{\frac{n(p-1)}{2}}}{\rm d}\tau {\rm d}s,

其中 t_0 > t_4 待定.将(3.38)式在 [t_4, t] 上作积分可得

\begin{equation} \alpha(t)-\alpha(t_4)\gtrsim \varepsilon \int_{t_4}^t(2\mu\tau+1)^{\frac{\mu-1}{2}}{\rm d}\tau +\beta(t)-\beta(t_4), \end{equation}
(3.39)

其中 \alpha(t_4)\geq0 , \beta(t_4)\leq0 .于是,由(3.39)式可知,当 t 充分大时,即 t > t_0 > t_4 ,有

\begin{eqnarray} \alpha(t)&\gtrsim& \frac{\varepsilon}{2}(2\mu t+1)^{\frac{\mu+1}{2}} +\left[ \frac{\varepsilon}{2}(2\mu t+1)^{\frac{\mu+1}{2}}- \varepsilon(2\mu t_4+1)^{\frac{\mu+1}{2}}\right]+\beta(t)\\ &\gtrsim& \varepsilon(2\mu t+1)^{\frac{\mu+1}{2}}+\beta(t). \end{eqnarray}
(3.40)

此外,直接计算可得

\begin{equation} (2\mu t+1)\beta''(t)-\mu(\mu-1)\beta'(t) = \mu(\mu-1) \int_0^t\frac{(2\mu \tau+1)^{\frac{\mu-1}{2}}}{(1+\tau)^{\frac{n(p-1)}{2}}}F^p(\tau){\rm d}\tau. \end{equation}
(3.41)

由Hölder不等式,我们发现

\begin{eqnarray} \alpha(t)\lesssim (1+t)\left(\int_0^t(2\mu\tau+1)^{\frac{\mu-1}{2}}\frac{F^p(\tau)}{(1+\tau)^{\frac{n(p-1)}{2}}}{\rm d}\tau \right)^{\frac{1}{p}}\left(\int_0^t\frac{(2\mu\tau+1)^{\frac{\mu-1}{2}}}{(1+\tau)^{p'-\frac{n(p-1)}{2p}p'}}{\rm d}\tau\right)^{\frac{1}{p'}}, \end{eqnarray}
(3.42)

这里 p' = \frac{p}{p-1}. 很容易看出

\begin{eqnarray} \int_0^t\frac{(2\mu\tau+1)^{\frac{\mu-1}{2}}}{(1+\tau)^{p'-\frac{n(p-1)}{2p}p'}}{\rm d}\tau & = &\int_0^t(2\mu\tau+1)^{\frac{\mu-1}{2}}(1+\tau)^{\frac{n}{2}-\frac{p}{p-1}}{\rm d}\tau\\ &\lesssim&\int_0^t(1+\tau)^{\frac{\mu-1}{2}+\frac{n}{2}-\frac{p}{p-1}}{\rm d}\tau\lesssim (1+t)^A, \end{eqnarray}
(3.43)

这里

\begin{equation} A = \left\{\begin{array}{ll}\Upsilon+1, &\mbox{if}\;\;\Upsilon+1>0, \\ 1, &\mbox{if}\;\;\Upsilon+1 = 0, \\0, &\mbox{if}\;\;\Upsilon+1<0, \end{array}\right. \end{equation}
(3.44)

其中

\begin{equation} \Upsilon = \frac{\mu-1}{2}+\frac{n}{2}-\frac{p}{p-1}. \end{equation}
(3.45)

将(3.43)式代入到(3.42)式可得

\alpha(t)\lesssim (1+t)^{1+\frac{A}{p'}}\left(\int_0^t(2\mu\tau+1)^{\frac{\mu-1}{2}}\frac{F^p(\tau)}{(1+\tau)^{\frac{n(p-1)}{2}}}{\rm d}\tau \right)^{\frac{1}{p}},

\begin{eqnarray} \int_0^t\frac{(2\mu\tau+1)^{\frac{\mu-1}{2}}}{(1+\tau)^{\frac{n(p-1)}{2}}}F^p(\tau){\rm d}\tau &\gtrsim & \alpha^p(t)(1+t)^{-p-\frac{p}{p'}A}\\ &\gtrsim & (1+t)^{-p-(p-1)A}\left[\varepsilon(2\mu t+1)^{\frac{\mu+1}{2}}+\beta(t)\right]^p, \end{eqnarray}
(3.46)

这里我们利用了(3.40)式.将(3.46)式代入(3.41)式可以得到

\begin{equation} (2\mu t+1)\beta''(t)-\mu(\mu-1)\beta'(t)\gtrsim \mu(\mu-1) (1+t)^{-p-(p-1)A}\left[\varepsilon(2\mu t+1)^{\frac{\mu+1}{2}}+\beta(t)\right]^p. \end{equation}
(3.47)

\beta(t) = (2\mu t+1)^{\frac{\mu+1}{2}}(\gamma(t)-\varepsilon) .于是, (3.47)式可写为

\begin{equation} (2\mu t+1)\gamma''(t)+\mu(\mu+3)\gamma'(t)\gtrsim \mu(\mu-1) \frac{\gamma^p(t)}{(1+t)^{p+(p-1)\left(A-\frac{\mu+1}{2}\right)}}. \end{equation}
(3.48)

而且,容易看出

\begin{equation} \gamma(t_0) = \varepsilon>0, \quad \gamma'(t_0)>0. \end{equation}
(3.49)

在(3.48)式两端同时乘以 (2\mu t+1)^m ,其中 m 待定,于是我们得到

\begin{equation} (2\mu t+1)^{m+1}\gamma''(t)+(2\mu t+1)^m\gamma'(t)\gtrsim \frac{\gamma^p(t)}{(1+t)^{p-m+(p-1)\left(A-\frac{\mu+1}{2}\right)}}\equiv\frac{\gamma^{1+(p-1)}(t)}{(1+t)^{\beta}}. \end{equation}
(3.50)

为了应用推论2.1,我们需要保证

0\leq\beta = p-m+(p-1)\left(A-\frac{\mu+1}{2}\right)<1,

等价地

\begin{equation} 0\leq\beta = (p-1)[A-(\Upsilon+1)]+\frac{n}{2}(p-1)-m<1. \end{equation}
(3.51)

由于在本文中我们仅考虑情形: 1 < p < 1+\frac{2}{n} ,这意味着 0 < \frac{n}{2}(p-1) < 1 ,于是我们可得

m = (p-1)[A-(\Upsilon+1)].

由(3.44)式,我们观察到 m\geq0 ,从而 \beta = \frac{n}{2}(p-1) .进而,取 \alpha = p-1 并通过(3.4)式回到原来的坐标 (\tau, y) ,我们就可以运用推论2.1证明定理1.1的结果.

参考文献

Almenar P , J'odar L , Matin J A .

Mixed problems for the time-dependent telegraph equation: continuous numerical solutions with a priori error bounds

Math Comput Modelling, 1997, 25, 31- 44

URL     [本文引用: 1]

Kato M , Sakuraba M .

Global existence and blow-up for semilinear damped wave equations in three space dimensions

Nonlinear Analysis, 2019, 182, 209- 225

URL     [本文引用: 1]

Li T T , Zhou Y .

Breakdown of solutions to \Box u+u_t=|u|^{1+\alpha}

Discrete Contin Dyn Syst, 1995, 1, 503- 520

DOI:10.3934/dcds.1995.1.503      [本文引用: 21]

Lai N A , Takamura H , Wakasa K .

Blow-up for semilinear wave equations with the scale invariant damping and super-Fujita exponent

J Differential Equations, 2017, 263, 5377- 5394

DOI:10.1016/j.jde.2017.06.017      [本文引用: 5]

Ikeda M , Sobajima M .

Life-span of solutions to semilinear wave equation with time-dependent critical damping for specially localized initial data

Math Ann, 2018, 372, 1017- 1040

DOI:10.1007/s00208-018-1664-1      [本文引用: 3]

Palmieri A , Reissig M .

A competition between Fujita and Strauss type exponents for blow-up of semi-linear wave equations with scale-invariant damping and mass

J Differential Equations, 2019, 266, 1176- 1220

DOI:10.1016/j.jde.2018.07.061      [本文引用: 1]

D'Abbicco M .

The threshold of effective damping for semilinear wave equation

Math Methods Appl Sci, 2015, 38, 1032- 1045

DOI:10.1002/mma.3126      [本文引用: 1]

D'Abbicco M , Lucente S .

A modified test function method for damped wave equation

Adv Nonlinear Stud, 2013, 13, 863- 889

URL     [本文引用: 1]

Zhang Q S .

A blow-up result for a nonlinear wave equation with damping: the critical case

C R Acad Sci Paris Ser I Math, 2001, 333, 109- 114

DOI:10.1016/S0764-4442(01)01999-1      [本文引用: 1]

Zhou Y .

Blow up of solutions to semilinear wave equations with critical exponent in high dimensions

Chin Ann Math Ser B, 2007, 28, 205- 212

DOI:10.1007/s11401-005-0205-x     

Lai N A , Zhou Y .

The sharp lifespan extimate for semilinear damped wave equation with Fujita critical power in higher dimensions

J Math Pures Appl, 2019, 123, 229- 243

DOI:10.1016/j.matpur.2018.04.009      [本文引用: 15]

Yordanov B T , Zhang Q S .

Finite time blow up for critical wave equations in high dimensions

Journal of Functional Analysis, 2006, 231, 361- 374

URL     [本文引用: 1]

Wirth J .

Wave equations with time-dependent dissipation I. Non-effective dissipation

J Differential Equations, 2006, 222, 487- 514

DOI:10.1016/j.jde.2005.07.019      [本文引用: 2]

Wirth J .

Wave equation with time-dependent dissipation Ⅱ. Effective dissipation

J Differential Equations, 2007, 232, 74- 103

DOI:10.1016/j.jde.2006.06.004      [本文引用: 2]

Wakasugi Y. Critical Exponent for the Semilinear Wave Equation with Scale Invariant Damping//Ruzhansky M, Turunen V, et al. Fourier Analysis. Boston: Birkhäuser, 2014: 375-390

[本文引用: 2]

Li Y C .

Classical solutions for fully nonlinear wave equuations with dissipaction (in Chinese)

Chin Ann Math Ser A, 1996, 17, 451- 466

[本文引用: 1]

Todorova G , Yordanov B .

Critical exponent for a nonlinear wave equation with damping

J Differ Equ, 2001, 174, 464- 489

DOI:10.1006/jdeq.2000.3933      [本文引用: 1]

Nishihara K .

L_p-L_q estimates of solutions to the damped wave quation in 3-dimensional space and their application

Math Z, 2003, 244, 631- 649

DOI:10.1007/s00209-003-0516-0      [本文引用: 1]

Takamura H .

Improved Kato's lemma on ordinary differential inequality and its application to semilinear wave equations

Nonlinear Analysis, 2015, 125, 227- 240

DOI:10.1016/j.na.2015.05.024      [本文引用: 2]

/