Processing math: 3%

数学物理学报, 2025, 45(3): 726-747

一类半线性波动方程弱耦合系统解的破裂

冯振东,1, 郭飞,2,3,*, 李岳群,2

1丽水职业技术学院智能制造学院 浙江丽水 323000

2南京师范大学数学科学学院 南京 210023

3南京师范大学大规模复杂系统数值模拟教育部重点实验室 南京 210023

Breakdown of Solutions to a Weakly Coupled System of Semilinear Wave Equations

Feng Zhendong,1, Guo Fei,2,3,*, Li Yuequn,2

1School of Intelligent Manufacturing, Lishui Vocational and Technical College, Zhejiang Lishui 323000

2School of Mathematical Sciences, Nanjing Normal University, Nanjing 210023

3Key Laboratory of NSLSCS (NNU), Ministry of Education, Nanjing 210023

通讯作者: E-mail: guof@njnu.edu.cn

收稿日期: 2024-01-25   修回日期: 2024-10-15  

基金资助: 国家自然科学基金(11731007)
江苏省高校优势学科建设工程
江苏省自然科学基金(BK20221320)
和江苏省研究生科研与实践创新计划项目(KYCX24_1786)

Received: 2024-01-25   Revised: 2024-10-15  

Fund supported: NSFC(11731007)
Priority Academic Program Development of Jiangsu Higher Education Institutions, the NSF of Jiangsu Province(BK20221320)
Postgraduate Research and Practice Innovation Program of Jiangsu Province(KYCX24_1786)

作者简介 About authors

E-mail:492024816@qq.com;

yqli1214@163.com

摘要

该文关注了一类带有尺度不变阻尼项, 质量项和一般的非线性记忆项的半线性波动方程弱耦合系统的 Cauchy 问题. 首先, 在指数 p,q 和系数 μ1,μ2 的适当假设下, 利用 Banach 不动点定理建立了该问题的局部解. 这里, pq 为非线性记忆项的指数, μ1μ2 分别为阻尼项和质量项的系数. Palmieri 对相应线性齐次方程解的衰减估计在证明局部解过程中起了重要作用. 之后, 采用迭代技术并结合试验函数方法, 得到了能量解的破裂.

关键词: 半线性波动方程; 弱耦合系统; 破裂; 试验函数; 迭代技术

Abstract

This paper addresses the Cauchy problem for a weakly coupled system of semilinear wave equations with scale-invariant dampings, mass, and general nonlinear memory terms. Firstly, a local (in time) existence result for this problem is established using Banach's fixed point theorem, subject to suitable assumptions on the exponents p,q and coefficients μ1,μ2.Here, p and q represent the powers of the nonlinear memory terms, while μ1 and μ2 denote the coefficients of the dampings and mass terms, respectively. It is noteworthy that Palmieri's decay estimates for the solution to the corresponding linear homogeneous equation play a crucial role in proving the local well-posedness result. Subsequently, employing an iteration technique in conjunction with the test function method, we obtain a blowup result for energy solutions.

Keywords: semilinear wave equation; weakly coupled system; blowup; test function method; iteration skill

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

本文引用格式

冯振东, 郭飞, 李岳群. 一类半线性波动方程弱耦合系统解的破裂[J]. 数学物理学报, 2025, 45(3): 726-747

Feng Zhendong, Guo Fei, Li Yuequn. Breakdown of Solutions to a Weakly Coupled System of Semilinear Wave Equations[J]. Acta Mathematica Scientia, 2025, 45(3): 726-747

1 引言

本文研究如下带有尺度不变阻尼项, 质量项和一般的非线性记忆项的半线性波动方程弱耦合系统的 Cauchy 问题

{uttΔu+μ11+tut+μ22(1+t)2u=g1|v|p,t>0,vttΔv+μ11+tvt+μ22(1+t)2v=g2|u|q,t>0,(u,ut,v,vt)(0,x)=(u0,u1,v0,v1)(x),xRn,
(1.1)

其中 u=u(t,x)v=v(t,x) 是未知函数, t 是时间变量,

xRn, 系数 μ1>0μ2, 指数 p,q>1, (时间) 卷积非线性项是

(g_k\ast |w|^r)(t,x)\triangleq \int _{0}^{t} g_k(t-\eta)|w(\eta,x)|^r {\rm d}\eta, w=v, u,

这里的时变记忆核 (或称松弛函数) g_k=g_k(t)>0, k=1,2.

(1.1) 式中的方程称为尺度不变的, 是因为其相应的线性方程在所谓的双曲尺度变换 (hyperbolic scaling)

\begin{equation*} \widetilde{u}(t,x)=u(\lambda (1+t)-1,\lambda x),\quad \lambda >0 \end{equation*}之下是不变的, 见文献 [23].

方程

\begin{equation*}\label{wuli1} m\dot{V}(t)=-\int_{-\infty}^{t}k(t-\tau)V(\tau){\rm d}\tau + F_{R}(t) \end{equation*}

在物理学中用来描述大颗粒在溶剂中的运动. 这里 m 是质量, k 是记忆项, F_{R}(t) 是随机力. 由 V(t) 表征的稳态解可视为一中心高斯过程, 它的自相关函数 r(t) 满足确定性时滞微分方程

m\dot{r}(t)=-\int_{0}^{t}k(t-\tau)r(\tau){\rm d}\tau, \quad r(0)=I.

方程中的右端项就是特殊的非线性记忆项, 更多细节参见文献 [18].

方程

\begin{equation}\label{Euler} u_{tt}-\Delta u=\frac1{\Gamma(1-\gamma)} \int_0^t (t-s)^{-\gamma} |u(s)|^p {\rm d}s, \end{equation}
(1.2)

其中 0<\gamma<1,\ p>1, \Gamma 是 Euler-Gamma 函数, 可视为经典的半线性波动方程

\begin{equation}\label{1.3} u_{tt}-\Delta u=|u|^p \end{equation}
(1.3)

的一个近似, 这是因为在分布意义下成立

\lim_{\gamma\rightarrow1} \frac1{\Gamma(1-\gamma)} s^{-\gamma}_+ = \delta (s).

研究方程 (1.3) 的 Cauchy 问题的解的破裂与小初值整体解的存在性问题称为 Strauss 猜想.具体来说, 若 1< p\leqslant p_S(n), 则具非负初值的解在有限时刻内破裂, 否则, 当 p>p_S(n) 时, 具小初值的解整体存在, 其中临界指数 p=p_S(n) 是 Strauss 指数: 当 n=1 时, p_S(1)=\infty; 当 n\geqslant2 时则是如下方程之正根

\begin{equation}\label{gamma S} \gamma_S(n,r)\triangleq(n-1)r^2 - (n+1) r -2 =0. \end{equation}
(1.4)

更多细节可参见文献 [7-9,12,13,17,27-29,32-34].

回顾有关无记忆项的半线性波方程的一些结果. 对尺度不变的阻尼型波方程

\begin{equation} \left\{ \begin{array}{ll} u_{tt}-\Delta u+\frac{\mu_1}{1+t}u_t=|u|^p, &\quad t>0, \\ (u,u_t)(0,x)=(u_0,u_1)(x), &\quad x\in \mathbb{R}^n, \end{array} \right. \end{equation}
(1.5)

其中 p>1, \mu_1>0, 当 1< p\leqslant p_F(n), \mu_1>11< p\leqslant p_F(n+\mu_1-1), 0<\mu_1\leqslant1 时, Wakasugi[31] 得出了解的破裂结果, 这里 p_F(n)=1+\frac{2}{n} 是 Fujita 指数. 对方程 (1.5) 解的适定性和渐近行为也已有广泛研究[5,10,15].

若考虑质量项, (1.5) 式化为

\begin{equation} \left\{ \begin{array}{ll} u_{tt}-\Delta u+\frac{\mu_1}{1+t}u_t+\frac{{\mu_2}^2}{(1+t)^{2}}u=|u|^p, &\quad t>0,\\ (u,u_t)(0,x)=(u_0,u_1)(x), &\quad x\in \mathbb{R}^n, \end{array} \right. \end{equation}
(1.6)

其中 p>1, \mu_1, \mu_2 为非负常数. 记

\begin{equation}\label{delta def} \delta=(\mu_1-1)^2-4\mu_2^2, \end{equation}
(1.7)

它用来刻画阻尼项和质量项之间的相互作用. 一方面, 当 \delta \geqslant(n+1)^2 时, 阻尼起主导作用, 临界指数为 Fujita 型的 p_F\big(n+{\mu_1-1-\sqrt{\delta}}{2}\big)[21],此时, 方程性状表现地像抛物型方程. 另一方面, 当 \delta \geqslant0, 1< p< p_S(n+\mu_1)0\leqslant \delta< n^2, p=p_S (n+\mu_1) 时, 文献 [25] 证明了 (1.6) 式的能量解将在有限时刻内破裂, 并得到了生命跨度的上界估计. 因此, 至少在考虑解的破裂现象时, 此时 (1.6) 式性状表现地类似波动方程.

下面考察两个带有非线性幂次项的半线性波动方程的弱耦合系统.

\begin{equation} \left\{ \begin{array}{ll} u_{tt}-\Delta u=|v|^p, & t>0,\\ v_{tt}-\Delta v=|u|^q, & t>0, \\ (u,u_t,v,v_t)(0,x)=(u_0,u_1,v_0,v_1)(x), & x\in \mathbb{R}^n, \end{array} \right. \end{equation}
(1.8)

其中 p,\,q>1, 在 p-q 平面上的临界曲线由

\label{1.8 blow up} F(p,q)\triangleq \max\left\{{p+2+q^{-1}}{pq-1},{q+2+p^{-1}}{pq-1} \right \}-{n-1}{2}

决定. 大致来说, 当 F(p,q)\geqslant 0 时, (1.8) 式的局部解将在有限时刻内破裂; 而当 F(p,q)< 0 时, 其小初值整体解将存在且唯一. 应当指出弱耦合系统的研究并不是对单个半线性方程的相应结果的简单推广. 事实上, 易知

1< p\leqslant p_S(n) \ \ \Leftrightarrow\ \ {1+p^{-1}}{p-1}\geqslant {n-1}{2}.

注意到不等式

\label{1.83} \max\bigg \{{p+2+q^{-1}}{pq-1},{q+2+p^{-1}}{pq-1} \bigg \}\leqslant \max\bigg \{{1+p^{-1}}{p-1},{1+q^{-1}}{q-1} \bigg \}

中, 等号当且仅当 p=q 时成立. 因而, 对 p\neq q 有可能出现

F(p,q)\triangleq \max\bigg \{{p+2+q^{-1}}{pq-1},{q+2+p^{-1}}{pq-1} \bigg \}\geqslant {n-1}{2},

(这意味着在 p-q 平面内 (p,q) 属于系统 (1.8) 的解的破裂区域), 但 pq 中可能有一个比 Strauss 指数 p_S(n) 大, 即单个方程 (1.3) 的解全局存在. 在文献 [1,11,14,26] 中可查阅更多细节.

若考虑到线性阻尼项的影响, 则 (1.8) 式化为

\begin{equation}\label{1.9} \left\{ \begin{array}{ll} u_{tt}-\Delta u+u_t=|v|^p, & t>0,\\ v_{tt}-\Delta v+v_t=|u|^p, & t>0,\\ (u,u_t,v,v_t)(0,x)=(u_0,u_1,v_0,v_1)(x), & x\in \mathbb{R}^n, \end{array} \right. \end{equation}
(1.9)

其中 p,q>1. 系统 (1.9) 的临界曲线为

\max \bigg\{{p+1}{pq-1}, {q+1}{pq-1}\bigg\}={n}{2}.

Sun 和 Wang[30] 分别研究了当 n=1,3 时 (1.9) 式的解的整体存在性和 n\geqslant 1 的整体解的不存在性. Narazaki[19] 将 Sun 和 Wang 的整体解的结果推广到 n=1,2,3 并改进了 n=3 时的时间衰减估计. 另外,Nishihara[20] 对任意空间维数在加权能量空间中给出了 (1.9) 式的解的整体存在性和破裂结果.

现在, 考虑具有非线性卷积记忆项的半线性波动方程或方程组的相关结果. Chen 与 Palmieri[3] 对带有特殊的记忆项 (1.2) 的问题

\begin{equation}\label{1.61} \left\{ \begin{array}{ll} u_{tt}-\Delta u=g \ast |u|^p, & t>0,\\ (u,u_t)(0,x)=(u_0,u_1)(x), \ & x\in \mathbb{R}^n, \end{array} \right. \end{equation}
(1.10)

确定了一个广义指数 p_{0,S}(n,\gamma), 它是下述方程的最大正根

(n-1)p^2-\big( n+1+2(1-\gamma) \big) p-2=0.

他们进一步指出, 在次临界情况 (即 1 < p < p_{0,S}(n,\gamma), n\geqslant2p > 1, n=1) 与临界情况下 (即 p=p_{0,S}(n,\gamma), n\geqslant 2), 问题 (1.10) 的能量解在有限时间内破裂.

Chen[2] 研究了一类具有一般非线性记忆项的弱耦合系统

\begin{equation}\label{1.10} \left\{ \begin{array}{ll} u_{tt}-\Delta u=g_1 \ast |v|^p, & t>0,\\ v_{tt}-\Delta v=g_2 \ast |u|^p, & t>0,\\ (u,u_t,v,v_t)(0,x)=(u_0,u_1,v_0,v_1)(x), & x\in \mathbb{R}^n, \end{array} \right. \end{equation}
(1.11)

其中 p,\,q>1. 若记忆核 g_k (k=1,2) 满足两种不同的增长条件 (见文献 [2,定理 2.2]), 他证明了方程 (1.11) 的能量解 (u,v) 必在有限时间内破裂. 值得注意的是, 系统 (1.8), (1.9) 和 (1.11) 中波传播的速度相同. Li 和 Wang[16] 研究了一组具有不同波速的耦合波动方程的部分同步系统.

本文在假设 p,q>1 和常数 \mu_1>0, \mu_2\geqslant0 使得 \delta > 0 情形下研究 Cauchy 问题 (1.1), 另外, 还假定 g_k(t)\in C^1([0,\infty)), k=1,2.

为了研究解的破裂现象, 需要建立解的局部存在性结果 (定理 2.1). 我们将利用 Banach 不动点定理来证明定理 2.1. 为此, 需要用到 Palmieri 对相应线性齐次方程解的衰减估计来构造工作空间 (见引理 3.1-3.2). 然而, 由于参数 \mu_1\mu_2 的取值变化, 导致相互作用量 \delta 的范围不确定, 故需要根据 \delta 的值进行分类分析.

研究半线性波方程解的破裂的一种有效方法是试验函数方法并结合 Kato 型微分不等式, 参见文献 [28,32]. 然而, 由于记忆项的影响, 得到破裂函数的精确下界估计具有挑战性 (见文献 [32,(2.7) 式]). 受文献 [2] 和 [25] 的启发, 我们将采用试验函数结合迭代方法 (见 Palmieri 与 Takamura[24], Chen 与 Reissig[4] 以及其参考文献) 来研究 (1.1) 式解的破裂. 由于方程组 (1.1) 较 (1.11) 式来说更为复杂, 在 Chen[2] 中所选取的试验函数不能直接运用, 本文选取函数

\begin{equation*} \Phi (t,x):=\lambda(t)\varphi(x) \end{equation*}

作为替代, 这里

\lambda(t):=(1+t)^{{\mu_1+1}{2}}K_{{\sqrt{\delta}}{2}}(t+1),

其中的 K_{{\sqrt{\delta}}{2}}(t+1) 是修正的第二类 Bessel 函数 (t\geqslant 0);

\begin{equation*} \varphi := \begin{cases} \int_{\mathbb{S}^{n-1}}{\rm e}^{x\cdot \omega}{\rm d} \omega, &n\geqslant 2,\\ {\rm e}^{x}+{\rm e}^{-x},& n=1, \end{cases} \end{equation*}

\mathbb{S}^{n-1}n-1 维单位球面. 此外, 由于非线性记忆项的一般性, 文献 [28] 中所建立的 Kato 引理对 (1.1) 式不再适用. 受 Chen[2] 的启发, 我们引入函数 G_k(t)(k=1,2) 并采用迭代过程进行证明.

本文安排如下. 在第 2 节, 给出本文的主要结果, 包括局部适定性 (定理 2.1) 和解的破裂结果 (定理 2.2). 在对线性方程结果做初步介绍之后, 第 3 节将给出局部适定性结果的证明. 最后在第四节, 将通过试验函数方法结合迭代技术来证明主要结果: 解的破裂 (定理 2.2).

记号说明. 本文中 f\lesssim (\text{或}\gtrsim)\ g 表示存在常数 C>0 使得 f\leqslant (\text{或}\geqslant)\ Cg, g\lesssim f \lesssim g 记为 f\sim g. 记自然数集为 N_0=N\cup \{0\}. 此外, B_R=B_R(0)\mathbb{R}^n 中以 R 为半径的球. 函数空间均定义在 \mathbb{R}^n 上, 在不引起混淆的情况下省略 \mathbb{R}^n.

2 主要结果

先给出系统 (1.1) 的局部适定性结果

定理 2.1 局部适定性 假设 (u_0,u_1,v_0,v_1) \in (H^1 \times L^2)^2, 其紧支集含在球 B_R(0) 中;

n=1,2 \text { 时 } 1<p, q<\infty, n \geqslant 3 \text { 时 } 1<p, q \leqslant \frac{n}{n-2} \text {; }
(2.1)

\mu_1 >0, \mu_2\geqslant 0 使得

\delta =(\mu_1-1)^2-4\mu_2^2>0;

非负时变记忆核 g_k(t) \in L_{loc}^1([0,\infty)) (k=1,2). 则存在 T\in (0,\infty] 和系统 (1.1) 的唯一温和解 (mild solution)

(u,v) \in \left(C([0,T),H^1)\cap C^1([0,T),L^2)\right)^2,

使得对任意的 t \in [0,T), \text{supp} u(t,\cdot),\ \text{supp} v(t,\cdot) \subset B_{R+t}.

解的破裂结果与能量解有关, 其定义如下

定义 2.1(u_0,u_1,v_0,v_1) \in (H^1 \times L^2)\times (H^1 \times L^2), (u,v) 是弱耦合系统 (1.1) 在 [0,T) 上的能量解, 如果

\begin{equation} \begin{split} &u\in C([0,T),H^1)\cap C^1([0,T),L^2) \ \mbox{且} g_2(t)\ast |u|^q \in L_{loc}^1([0,\infty)\times \mathbb{R}^n);\\ &v\in C([0,T),H^1)\cap C^1([0,T),L^2)\ \mbox{且} g_1(t)\ast |v|^p \in L_{loc}^1([0,\infty)\times \mathbb{R}^n), \end{split} \end{equation}
(2.2)

且对任意的试验函数 \phi,\psi \in C_0^{\infty}([0,\infty)\times \mathbb{R}^n) 和任意的 t\in (0,T), 成立

\begin{equation}\label{energy solution u} \begin{split} & \int u_t(t,x)\phi(t,x){\rm d}x-\int u_t(0,x)\phi(0,x){\rm d}x\\ & +\int_0^t\int (-u_s(s,x)\phi_s(s,x)+\nabla u(s,x)\cdot\nabla \phi(s,x)){\rm d}x {\rm d}s\\ & +\int_0^t\int \bigg (\frac{\mu_1}{1+s}u_s(s,x)+\frac{{\mu_2}^2}{(1+s)^{2}}u(s,x)\bigg )\phi (s,x){\rm d}x {\rm d}s\\ &=\int_0^t\int (g_1\ast |v|^p)(s,x)\phi(s,x){\rm d}x{\rm d}s \end{split} \end{equation}
(2.3)

\begin{equation}\label{energy solution v} \begin{split} & \int v_t(t,x)\psi(t,x){\rm d}x-\int v_t(0,x)\psi(0,x){\rm d}x\\ & +\int_0^t\int (-v_s(s,x)\psi_s(s,x)+\nabla v(s,x) \cdot \nabla \psi(s,x)){\rm d}x {\rm d}s\\ & +\int_0^t\int \bigg (\frac{\mu_1}{1+s}v_s(s,x)+\frac{{\mu_2}^2}{(1+s)^{2}}v(s,x)\bigg )\psi (s,x){\rm d}x {\rm d}s\\ &=\int_0^t\int (g_2\ast |u|^q)(s,x)\psi(s,x){\rm d}x{\rm d}s. \end{split} \end{equation}
(2.4)

基于局部适定性结果和能量解的定义, 解的破裂结果是

定理 2.2 (解的破裂) 假设 n=1, 21< p,\ q<\infty, n\geqslant31< p,\ q\leqslant\frac{n}{n-2}; \mu_1 >0, \mu_2\geqslant 0 使得 \delta> 0; 非负时变记忆核 g_k(t)\in L^1_{loc}([0,\infty))\cap C^1([0,\infty)), \ k=1,2; p,\ q 满足

\begin{equation}\label{blowup2} \max\bigg \{{p+2+q^{-1}}{pq-1}-{\mu_1}{2},{q+2+p^{-1}}{pq-1}-{\mu_1}{2} \bigg \}> {n-1}{2}; \end{equation}
(2.5)

(u_0,u_1,v_0,v_1) \in \big(H^1 \times L^2 \big)^2 非负且其紧支集含在 B_R(0) 中, u_1,v_1 不恒为零且

\begin{equation}\label{condition1} u_0(x)\geqslant0 \mbox{且} u_1(x)+{\mu_1-1-\sqrt{\delta}}{2}u_0(x)\geqslant0, \end{equation}
(2.6)
\begin{equation} v_0(x)\geqslant0 \mbox{且} v_1(x)+{\mu_1-1-\sqrt{\delta}}{2}v_0(x)\geqslant0, \end{equation}
(2.7)

则弱耦合系统 (1.1) 在 [0,T) 上的局部能量解 (u,v) 满足

\text{supp} u,\ \text{supp} v \subset \{(t,x)\in [0,T)\times \mathbb{R}^n | |x|\leqslant t+R\},

且必在有限时刻内破裂, 这里 T 是解 (u,v) 的生命跨度.

例 2.1 考虑如下带有多项式衰减的非线性记忆项的弱耦合系统

\begin{equation}\label{ex1.1} \left\{ \begin{array}{ll} u_{tt}-\Delta u+\frac{\mu_1}{1+t}u_t+\frac{{\mu_2}^2}{(1+t)^{2}}u=\int _{0}^{t} {|v(\eta,x)|^p}{(1+t-\eta)^{\gamma_1}} {\rm d}\eta, & t>0,\\ v_{tt}-\Delta v+\frac{\mu_1}{1+t}v_t+\frac{{\mu_2}^2}{(1+t)^{2}}v=\int _{0}^{t} {|u(\eta,x)|^q}{(1+t-\eta)^{\gamma_2}} {\rm d}\eta, &t>0,\\ (u,u_t,v,v_t)(0,x)=(u_0,u_1,v_0,v_1)(x), & x \in \mathbb{R}^n. \end{array} \right. \end{equation}
(2.8)

在定理 2.2 的假设下, 可知弱耦合系统 (2.8) 的解 (u,v) 在有限时间内破裂.

3 定理 2.1 的证明

3.1 预备工作

线性齐次方程

\begin{equation} \label{L W EQ} \left\{ \begin{array}{ll} w_{tt}-\Delta w+\frac{\mu_1}{1+t}w_t+\frac{{\mu_2}^2}{(1+t)^{2}}w=0, & t>0,\\ (w,w_t)(s,x)=(w_0,w_1)(x), & x \in \mathbb{R}^n \end{array} \right. \end{equation}
(3.1)

的解 w=w(t,x) 可表成

w(t,x)=E_0(t,s,x)\ast_{(x)} w_0(x)+E_1(t,s,x)\ast_{(x)} w_1(x),

这里 \mu_1 >0, \mu_2\geqslant 0, E_0(t,s,x)E_1(t,s,x) 是在分布意义下初值分别为 (w_0,w_1)=(\delta_0,0)(0,\delta_0) 的解, \delta_0 是变量 x 的 Dirac 分布函数, \ast_{(x)} 表示对变量 x 的卷积. 另外, 要求初值属于具有 L^1 正则性的经典能量空间, 即

(w_0,w_1) \in (H^1 \cap L^1)\times (L^2\cap L^1).

对任意的 k \in [0,1], 令

D^k(\mathbb{R}^n):=(H^k\cap L^1)\times (L^2\cap L^1),\ D(\mathbb{R}^n):= D^1 (\mathbb{R}^n).

下面两个衰减估计在局部适定性 (定理 2.1) 的证明中起了重要作用. 若初始时间为 s=0, 则有

引理 3.1\mu_1 >0, \mu_2\geqslant 0 使得 \delta =(\mu_1-1)^2-4\mu_2^2>0, (w_0,w_1) \in D(\mathbb{R}^n), 则对任意的 k \in [0,1], (3.1) 式的能量解 w 满足衰减估计

\begin{array}{l} \|w\|_{\dot{H}^{k}} \lesssim\left\|\left(w_{0}, w_{1}\right)\right\|_{D^{k}}\\ \times\left\{\begin{array}{ll} (1+t)^{-k-\frac{n+\mu_{1}}{2}+\frac{1+\sqrt{\delta}}{2}}, & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}>k, \ \ \ \ (3.2a)\\ (1+t)^{-\frac{\mu_{1}}{2}} \ell(t), & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}=k, \ \ \ \ (3.2b)\\ (1+t)^{-\frac{\mu_{1}}{2}}, & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}< k, \ \ \ \ (3.2c) \end{array}\right. \end{array}

这里 \ell(t)=1+(\ln(1+t))^{{1}{2}}. 此外, 若在 (3.2)诸式中取 k=1, 则 \|w_t(t,\cdot)\|_{L^2}\|\nabla w(t,\cdot)\|_{L^2} 满足同样的衰减估计.

若初始时间为 s\geqslant 0, 且第一个初值为零, 则有

引理 3.2\mu_1 >0, \mu_2\geqslant 0 使得 \delta =(\mu_1-1)^2-4\mu_1^2>0, w_0=0, w_1\in L^2\cap L^1, 则对 t\geqslant s 以及 k \in [0,1], (3.1) 式的能量解 w 满足

\begin{array}{l} \|w\|_{\dot{H}^{k}} \lesssim\left(\left\|w_{1}\right\|_{L^{1}}+(1+s)^{\frac{n}{2}}\left\|w_{1}\right\|_{L^{2}}\right)(1+s)^{\frac{1+\mu_{1}}{2}-\frac{\sqrt{\delta}}{2}}\\ \times\left\{\begin{array}{ll} (1+t)^{-k-\frac{n+\mu_{1}}{2}+\frac{1+\sqrt{\delta}}{2}}, & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}>k, \ \ \ (3.3a)\\ (1+t)^{-\frac{\mu_{1}}{2} \tilde{\ell}(t, s),} & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}=k, \ \ \ (3.3b)\\ (1+t)^{-\frac{\mu_{1}}{2}}(1+s)^{-k-\frac{n}{2}+\frac{1+\sqrt{\delta}}{2}}, & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}< k, \ \ \ (3.3c) \end{array}\right. \end{array}

其中 \tilde{\ell}(t,s)=1+(\ln ({1+t}{1+s}))^{{1}{2}}. 此外, 若在 (3.3) 诸式中取 k=1, 则 \|w_t(t,\cdot)\|_{L^2}\|\nabla w(t,\cdot)\|_{L^2} 满足同样的衰减估计.

引理 3.3 (Gagliardo-Nirenberg 不等式 [6])p,q,r (1\leqslant p,q,r\leqslant \infty)\sigma \in [0,1] 满足

\begin{equation} {1}{p}=\sigma ({1}{r}-{1}{n})+(1-\sigma){1}{q} \end{equation}
(3.4)

(除 p=\inftyr=n, n\geqslant 2 外). 则对任意的 h \in C_0^1, 存在常数 C=C(p,q,r,n)>0, 使得

\begin{equation} \|h\|_{L^p}\leqslant C \|h\|_{L^q}^{1-\sigma}\|\nabla h\|_{L^r}^{\sigma}. \end{equation}
(3.5)

现在证明定理 2.1. 定义

\begin{equation}\label{space def}\begin{split} X(T):=\bigg\{&(u,v)\in \left(C([T],H^1)\cap C^1([T],L^2)\right)^2 \\ &\text{且} \text{supp}\ u(t,\cdot), \text{supp} v(t,\cdot) \subset B_{R+t},\ \forall\ t \in [T]\bigg\}, \end{split} \end{equation}
(3.6)

并赋以范数

\|(u,v)\|_{X(T)}\triangleq \underset{t\in[T]}\max \big(M[u](t)+M[v](t)\big),

其中

M[w](t)\triangleq \|w(t,\cdot)\|_{L^2(\mathbb{R}^n)}+\|w_t(t,\cdot)\|_{L^2(\mathbb{R}^n)}+\|\nabla w(t,\cdot) \|_{L^2(\mathbb{R}^n)}, w=u,v.

容易验证 (X(T),\|\cdot\|_{X(T)}) 是 Banach 空间. 对 T,K>0, 定义

\begin{equation}\label{X(T,K) def} X(T,K):= \left\{(u,v) \in X(T) \bigg| \|(u,v)\|_{X(T)}\leqslant K\right\}. \end{equation}
(3.7)

由Duhamel 原理, 可定义算子

N: (u,v) \in X(T)\longrightarrow N(u,v)=(u^l+u^n,v^l+v^n),

这里

u^l=E_0(t,s,x)\ast u_0(x)+E_1(t,s,x)\ast u_1(x),
u^n=\int _0^t E_1(t-\tau,0,x)\ast (g_{1} \ast |v|^p)(\tau,x) {\rm d}\tau;
v^l=E_0(t,s,x)\ast v_0(x)+E_1(t,s,x)\ast v_1(x),
v^n=\int _0^t E_1(t-\tau,0,x)\ast (g_{2} \ast |u|^q)(\tau,x) {\rm d}\tau,

Nu= u^l + u^n 是如下问题的唯一解

\begin{equation}\label{Nu} \left\{ \begin{array}{ll} u_{tt}-\Delta u+\frac{\mu_1}{1+t}u_{t}+\frac{{\mu_2}^2}{(1+t)^{2}}u=g_1\ast|v|^p, & t>0,\\ (u,u_{t})(0,x)=(u_0,u_1)(x), & x \in \mathbb{R}^n, \end{array} \right. \end{equation}
(3.8)

Nv= v^l + v^n 是如下问题的唯一解

\begin{equation}\label{Nv} \left\{ \begin{array}{ll} v_{tt}-\Delta v+\frac{\mu_1}{1+t}v_{t}+\frac{{\mu_2}^2}{(1+t)^{2}}v=g_2\ast|u|^q, & t>0,\\ (v,v_{t})(0,x)=(v_0,v_1)(x), & x \in \mathbb{R}^n. \end{array} \right. \end{equation}
(3.9)

我们的目标是通过选取合适的 TK, 使映射 N 是从 X(T,K) 到其自身的压缩映射. 这意味着问题 (1.1) 的局部解为算子 N 的不动点. 为此, 需建立如下不等式

\begin{equation}\label{well-def} \|N(u,v)\|_{X(T)}\leqslant C_{0,T}+C_{T}(\|(u,v)\|_{X(T)}^p+\|(u,v)\|_{X(T)}^q), \end{equation}
(3.10)
\begin{equation}\label{contractive} \begin{split} \|N(u,v)-N(\overline{u},\overline{v})\|_{X(T)}\leqslant& C_{T}^{'} \|(u,v)-(\overline{u},\overline{v})\|_{X(T)} \\ &\times \underset{r=p,q}{\sum} (\|(u,v)\|_{X(T)}^{r-1}+\|(\overline{u},\overline{v})\|_{X(T)}^{r-1}), \end{split} \end{equation}
(3.11)

这里 C_{0,T} 仅与初值 u_0,u_1,v_0,v_1 的范数大小有关且当 T\to 0 是有界的. 另外, C_{T}C_{T}^{'} 均趋于 0 (T\to 0). 事实上, 因 T\rightarrow 0C_T \rightarrow 0, 故对于充分大的 K, 可选取充分小的 T 使得 (3.10) 式右端的两项均小于 {K}{2}. 因此, NX(T,K) 映射到自身. 而对上述充分小的 T, 易从 (3.11) 式推出 N 的压缩性. 因此, 在 X(T) 中, (1.1) 式局部解的存在性与唯一性就由 Banach 不动点定理得到.

考虑到 Nu = u^l + u^n Nv = v^l +v^n, 若要建立 (3.10) 与 (3.11) 式需要对 u^l,\ v^l, u^nv^n 做估计. 简便起见, 下面只给出 u^lu^n 的证明, v^lv^n 的估计可类似得出.

3.2 (3.10) 式的证明

在引理 3.1 中分别取 k=0,1, 可得

\begin{array}{l} \|w(t, \cdot)\|_{L^{2}} \leqslant C\left\|\left(w_{0}, w_{1}\right)\right\|_{D^{0}}\\ \times\left\{\begin{array}{ll} (1+t)^{-\frac{n+\mu_{1}}{2}+\frac{1+\sqrt{\delta}}{2}}, & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}>0, \ \ \ \ \ (3.12a)\\ (1+t)^{-\frac{\mu_{1}}{2}} \ell(t), & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}=0, \ \ \ \ \ (3.12b)\\ (1+t)^{-\frac{\mu_{1}}{2}}, & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}< 0 \ \ \ \ \ \ \ (3.12c) \end{array}\right. \end{array}

\begin{array}{l} \left\|w_{t}(t, \cdot)\right\|_{L^{2}},\|\nabla w(t, \cdot)\|_{L^{2}} \leqslant C\left\|\left(w_{0}, w_{1}\right)\right\|_{D}\\ \times\left\{\begin{array}{ll} (1+t)^{-1-\frac{n+\mu_{1}}{2}+\frac{1+\sqrt{\delta}}{2}}, & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}>1, \ \ \ \ \ (3.13a)\\ (1+t)^{-\frac{\mu_{1}}{2}} \ell(t), & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}=1, \ \ \ \ \ (3.13b) \\ (1+t)^{-\frac{\mu_{1}}{2}}, & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}<1, \ \ \ \ \ (3.13c) \end{array}\right. \end{array}

这里 \ell(t)=1+(\ln(1+t))^{{1}{2}}.

对应于 \delta 的不同取值, 下面的证明将分五种情形, 如下图所示.

情形 1 {1+\sqrt{\delta}-n}{2}>1>0, 即 \delta >(n+1)^2>(n-1)^2.

由 (3.12a) 与 (3.13a) 式得

\begin{equation}\label{u11} \|u^l\|_{L^2}\leqslant C (1+t)^{-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}}\|(u_0,u_1)\|_{D^0}, \end{equation}
(3.14)
\begin{equation}\label{ut ux11} \begin{split} \|u_t^l\|_{L^2},\|\nabla u^l\|_{L^2}&\leqslant C (1+t)^{-1-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}} \|(u_0,u_1)\|_{D} \\ &\leqslant C (1+t)^{-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}} \|(u_0,u_1)\|_{D}, \end{split} \end{equation}
(3.15)

\begin{equation}\label{u0 11} M[u^l](t)\leqslant C(1+t)^{-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}} \left( \|(u_0,u_1)\|_{D^0}+\|(u_0,u_1)\|_{D}\right). \end{equation}
(3.16)

从而

\begin{equation}\label{(u0,v0) 11} \|(u^l,v^l)\|_{X(T)} \leqslant C (1+t)^{-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}} \|(u_0,u_1,v_0,v_1)\|_{D\times D} \leqslant C_{0,T}. \end{equation}
(3.17)

容易看出 C_{0,T} 仅依赖于初值的范数且当 T\rightarrow 0 时有界. 类似地可以处理

情形 2 0<{1+\sqrt{\delta}-n}{2}<1, 即 (n-1)^2<\delta <(n+1)^2 (由 (3.12a) 和 (3.13c) 式可得);

情形 3 {1+\sqrt{\delta}-n}{2}<0<1, 即 \delta <(n-1)^2 <(n+1)^2 (由 (3.12c) 和 (3.13c) 式可得);

情形 4 0={1+\sqrt{\delta}-n}{2}<1, 即 (n-1)^2=\delta <(n+1)^2 (由 (3.12b) 和 (3.13c) 式可得);

情形 5 0<{1+\sqrt{\delta}-n}{2}=1, 即 (n-1)^2<\delta =(n+1)^2 (由 (3.12a) 和 (3.13b) 式可得).

因此

\begin{equation}\label{(u0,v0)1'} \begin{split} \|(u^l,v^l)\|_{X(T)}\leqslant C_{0,T}. \end{split} \end{equation}
(3.18)

接下来估计 u^n. 对任意的 \eta \in [0,T), 由 Gagliardo-Nirenberg 不等式 (引理 (3.3)) 可得

\begin{equation}\label{GN} \begin{split} \|w(\eta,\cdot)\|^r_{L^{2r}}\leqslant C\|w(\eta,\cdot)\|_{L^2}^{(1-{n}{2}(1-{1}{r}))r} \|\nabla w(\eta,\cdot)\|_{L^2}^{{n}{2}(1-{1}{r})r}\leqslant C(M[w](\eta))^r, \end{split} \end{equation}
(3.19)

这里 n=1,2r>1, n\geqslant31< r\leqslant \frac{n}{n-2}.

在引理 3.2 中分别取 s=0, k=0,1, 可得

\begin{array}{l} \|w(t, \cdot)\|_{L^{2}} \leqslant C\left(\left\|w_{1}(t, \cdot)\right\|_{L^{1}}+\left\|w_{1}(t, \cdot)\right\|_{L^{2}}\right)\\ \times\left\{\begin{array}{ll} (1+t)^{-\frac{n+\mu_{1}}{2}+\frac{1+\sqrt{\delta}}{2}}, & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}>0, \ \ \ \ \ \ (3.20a)\\ (1+t)^{-\frac{\mu_{1}}{2} \tilde{\ell}(t),} & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}=0, \ \ \ \ \ \ (3.20b)\\ (1+t)^{-\frac{\mu_{1}}{2},} & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}<0 \ \ \ \ \ \ \ \ (3.20c) \end{array}\right. \end{array}

\begin{aligned}\left\|w_{t}(t, \cdot)\right\|_{L^{2}},\|\nabla w(t, \cdot)\|_{L^{2}} \leqslant & C\left(\left\|w_{1}(t, \cdot)\right\|_{L^{1}}+\left\|w_{1}(t, \cdot)\right\|_{L^{2}}\right) \\& \times\left\{\begin{array}{ll}(1+t)^{-1-\frac{n+\mu_{1}}{2}+\frac{1+\sqrt{\delta}}{2}}, & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}>1,(3.21 \mathrm{a}) \\(1+t)^{-\frac{\mu_{1}}{2} \tilde{\ell}(t),} & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}=1,(3.21 \mathrm{~b}) \\(1+t)^{-\frac{\mu_{1}}{2},} & \text { 若 } \frac{1+\sqrt{\delta}-n}{2}<1,(3.21 \mathrm{c})\end{array}\right.\end{aligned}

这里 \tilde{\ell}(t)=1+(\ln(1+t))^{{1}{2}}.

情形 1 {1+\sqrt{\delta}-n}{2}>1>0, 即 \delta >(n+1)^2>(n-1)^2.

注意到 (3.6) 式中要求函数有紧支集, 利用 \rm H\ddot{o}lder 不等式, 可得

\begin{equation} \|g_{1} \ast |v|^p\|_{L^1}\leqslant (R+\tau)^{{n}{2}} \|g_{1} \ast |v|^p\|_{L^2}. \end{equation}
(3.22)

由 Gagliardo-Nirenberg 不等式得到

\begin{equation}\label{fact1} \begin{split} \| |v(\eta,\cdot)|^p\|_{L^2}&=\|v(\eta,\cdot)\|^p_{L^{2p}}{\leqslant} C\|v(\eta,\cdot)\|_{L^2}^{(1-{n}{2}(1-{1}{p}))p} \|\nabla v(\eta,\cdot)\|_{L^2}^{{n}{2}(1-{1}{p})p}\\ &\leqslant C(M[v](\eta))^p \leqslant C\|(u,v)\|_{X(T)}^p, \end{split} \end{equation}
(3.23)

这里当 n=1, 21< p<\infty, 当 n\geqslant 3 1< p\leqslant\frac{n}{n-2}. 再由 Duhamel 原理以及 (3.20a) 式可得

\begin{equation}\label{ur11} \begin{split} \|u^n\|_{L^2}& = \|\int _0^t E_1(t-\tau,0,x)\ast_{(x)} (g_{1} \ast |v|^p)(\tau,x) {\rm d}\tau\|_{L^2}\\ &\leqslant C\int _0^t (1+t-\tau)^{-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}} (\|g_{1} \ast |v|^p\|_{L^1}+\|g_{1} \ast |v|^p\|_{L^2}) {\rm d}\tau\\ & \leqslant C\int _0^t (1+t-\tau)^{-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}} (1+(R+\tau)^{{n}{2}})\|g_{1} \ast |v|^p\|_{L^2} {\rm d}\tau\\ &\leqslant C\int _0^t (1+t-\tau)^{-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}}(R+\tau)^{{n}{2}} \\ &\quad\times \int _{0}^{\tau} g_1(\tau-\eta)\| |v(\eta,\cdot)|^p\|_{L^2}{\rm d}\eta {\rm d}\tau\\ &\leqslant C\int _0^t (1+t)^{{1+\sqrt{\delta}+n}{2}} \|(u,v)\|_{X(T)}^p {\rm d}\tau \\ &\leqslant CT(1+T)^{{1+\sqrt{\delta}+n}{2}}\|(u,v)\|_{X(T)}^p, \end{split} \end{equation}
(3.24)

这里

\begin{equation}\label{fact2} \bigg| \int _{0}^{\tau} g_1(\tau-\eta) {\rm d}\eta\bigg|=\int _{0}^{\tau} g_1(\eta){\rm d}\eta<\int _{0}^{T} g_1(\eta){\rm d}\eta<\infty \end{equation}
(3.25)

是因为 g_1(t) \in L_{loc}^1([0,\infty)).

\nabla^j\partial_t^m u^n=\nabla^j\partial_t^m \int _0^t E_1(t-\tau,0,x)\ast (g_{1} \ast |v|^p)(\tau,x) {\rm d}\tau, 这里 j+m=1, j,m\in N\bigcup \{0\}. 易见 j=0,m=1 时, \nabla^j\partial_t^m u^n=u_t^n, 而 j=1,m=0 时, \nabla^j\partial_t^m u^n=\nabla u^n.

类似于 (3.24) 式, 由 (3.21a), (3.22) 和 (3.23) 式可得

\begin{equation}\label{urt urx11} \begin{split} \|\nabla^j\partial_t^m u^n\|_{L^2}&=\|\nabla^j\partial_t^m\int _0^t E_1(t-\tau,0,x)\ast_{(x)} (g_{1} \ast |v|^p)(\tau,x) {\rm d}\tau\|_{L^2}\\ &\leqslant C \int_0^t(1+t-\tau)^{-1-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}}(1+(R+\tau)^{{n}{2}})\|g_{1} \ast |v|^p\|_{L^2} {\rm d}\tau\\ &\leqslant C \int _0^t(1+t-\tau)^{-1-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}}\\ &\quad \cdot (R+\tau)^{{n}{2}}\bigg| \int _{0}^{\tau} g_1(\tau-\eta)\|v(\eta,\cdot)|^p\|_{L^2} {\rm d}\eta \bigg| {\rm d}\tau\\ &\leqslant C (1+t)^{{-1+\sqrt{\delta}}{2}}t\|(u,v)\|^{p}_{X(T)} \leqslant CT(1+T)^{{-1+\sqrt{\delta}}{2}}\|(u,v)\|^{p}_{X(T)}, \end{split} \end{equation}
(3.26)

上式中用到的

\begin{equation*} \begin{split} (1+t-\tau)^{-1-{n+\mu_1}{2}+{1+\sqrt{\delta}}{2}} <(1+t-\tau)^{{1+\sqrt{\delta}-n}{2}-1}<(1+t)^{{1+\sqrt{\delta}-n}{2}-1}=(1+t)^{{-1-n+\sqrt{\delta}}{2}} \end{split} \end{equation*} 是因为 \delta >(n+1)^2$ 的缘故.

结合 (3.24) 和 (3.26) 式, 得到

\begin{equation}\label{(ur,vr) 11} \begin{split} \|(u^n,v^n)\|_{X(T)}&\leqslant CT(1+T)^{{1+\sqrt{\delta}+n}{2}}(\|(u,v)\|^{p}_{X(T)}+\|(u,v)\|^{q}_{X(T)}).\\ \end{split} \end{equation}
(3.27)

可利用相同的方式来处理情形 2-5 中的 u^n. 因此

\left\|\left(u^{n}, v^{n}\right)\right\|_{X(T)} \leqslant C_{T}\left(\|(u, v)\|_{X(T)}^{p}+\|(u, v)\|_{X(T)}^{q}\right),
(3.28)

其中 C_{T}\rightarrow 0 (T\rightarrow 0).

至此, 不等式 (3.10) 可以由 (3.17) 和 (3.28) 式推出.

3.3 (3.11)式的证明

对任意的 (u,v),\ (\overline{u},\overline{v})\in X(T), 易得

\begin{equation}\label{Nuv-Nu'v'} \|N(u,v)- N(\overline{u},\overline{v})\|_{X(T)} =\left\|\int _0^t E_1(t-\tau,0,x)\ast_{(x)} (H_p(\tau,x),H_q(\tau,x)) {\rm d}\tau \right\|_{X(T)}, \end{equation}
(3.29)

其中

\begin{equation} H_p(\tau,x)=g_{1}\ast (|v|^p-|\overline{v}|^p)(\tau,x)=\int _{0}^{\tau} g_1(\tau-\eta)(|v(\eta,x)|^p-|\overline{v}(\eta,x)|^p) {\rm d}\eta \end{equation}
(3.30)

\begin{equation}\label{Hq} H_q(\tau,x)=g_{2}\ast (|u|^q-|\overline{u}|^q)(\tau,x)=\int _{0}^{\tau} g_2(\tau-\eta)(|u(\eta,x)|^q-|\overline{u}(\eta,x)|^q){\rm d}\eta. \end{equation}
(3.31)

注意到

| |u|^r-|\overline{u}|^r|\leqslant C|u-\overline{u}|(|u|^{r-1}+|\overline{u}|^{r-1}),

\begin{equation}\label{wr-w'r} \begin{split} \| |u|^r-|\overline{u}|^r\|_{L^2} &\leqslant C\|u-\overline{u}\|_{L^{2r}}\cdot \| |u|^{r-1}+|\overline{u}|^{r-1}\|_{L^{2r'}}\\ &\leqslant C\|u-\overline{u}\|_{L^{2r}}\cdot (\| |u|^{r-1}\|_{L^{2r'}}+\| |\overline{u}|^{r-1}\|_{L^{2r'}})\\ &\leqslant C M[u-\overline{u}](t)((M[u])^{r-1}+(M[\overline{u}])^{r-1}), \end{split} \end{equation}
(3.32)

这里当 n=1,2r>1, 而当 n\geqslant31< r\leqslant \frac{n}{n-2}, r'r 的共轭指数, 且用到了对任意的 \eta \in [T] 成立

\begin{array}{c}\|u-\bar{u}\|_{L^{2 r}} \leqslant M[u-\bar{u}](\eta), \\\left\||u|^{r-1}\right\|_{L^{2 r^{\prime}}}=\|u\|_{L^{2 r}}^{r-1} \leqslant C(M[u](\eta))^{r-1}, \\\left\||\bar{u}|^{r-1}\right\|_{L^{2 r^{\prime}}}=\|\bar{u}\|_{L^{2 r}}^{r-1} \leqslant C(M[\bar{u}](\eta))^{r-1}.\end{array}

由 (3.31) 式可得

\begin{equation} \begin{aligned} \| H_q(\tau,\cdot)\|_{L^2} &= \| g_2 \ast (|u|^q-|\overline{u}|^q) (\tau,\dot)\|_{L^2}\\ & \leqslant \int_0^{\tau} g_2(\tau-\eta) \||u|^q-|\overline{u}|^q \|_{L^2} {\rm d}\eta\\ &\leqslant C M[u-\overline{u}](\tau)((M[u])^{r-1}+(M[\overline{u}])^{r-1}). \end{aligned} \end{equation}
(3.33)

因此, 从 (3.29) 式得

\begin{equation} \begin{aligned} & \|N(u,v) - N(\overline{u}, \overline{v})\|_{L^2}\\ &\leqslant \int_0^t \| E_1(t-\tau,0,x)\ast_{(x)} (H_p(\tau,x),H_q(\tau,x))\|_{L^2} {\rm d}\tau\\ &\leqslant C\int (1+t-\tau)^{-\frac{n+\mu_1}2 + \frac{1+\sqrt{\delta}}2} \|(H_p,H_q)\|_{L^2}\\ &\leqslant C'_T \|(u,v) - (\overline{u}, \overline{v})\|_{X(T)} \sum_{r=p,q} (\|(u,v)\|_{X(T)}^{r-1} + \|(\overline{u}, \overline{v})\|_{X(T)}^{r-1}), \end{aligned} \end{equation}
(3.34)

其中 C'_T\rightarrow 0 (T\rightarrow 0). \|\nabla^j\partial_t^m(u^n,v^n)\|_{L^2} (j+m=1j,m\in N_0) 的估计可用类似方法得到. 至此, (3.11) 式证毕.

4 定理 2.2 的证明

4.1 准备工作

在构造试验函数之前, 先回顾第二类 \zeta 阶修正 Bessel 函数

K_{\zeta}(t)=\int_0^{\infty}\exp(-t \cosh z)\cosh (\zeta z){\rm d}z, \zeta \in \mathbb{R},

它是如下方程的解

\left( t^2 {{\rm d}^2}{{\rm d}t^2}+t{\rm d}{{\rm d}t}-(t^2+\zeta ^2) \right) K_{\zeta}(t)=0, t>0.

以下是 K_{\zeta}(t) (\zeta 是实参数) 一些有用的性质, 更多细节可参考文献 [25].

\bullet 极限性质

\begin{equation}\label{limit behavior} K_{\zeta}(t)=\sqrt{{\pi}{2t}}{\rm e}^{-t}[1+O(t^{-1})], t\rightarrow \infty. \end{equation}
(4.1)

\bullet 导数恒等式

\begin{equation}\label{deriative behavior} {\rm d}{{\rm d}t}K_{\zeta}(t)=-K_{\zeta+1}(t)+{\zeta}{t}K_{\zeta}(t). \end{equation}
(4.2)

\lambda(t):=(1+t)^{{\mu_1+1}{2}}K_{{\sqrt{\delta}}{2}}(t+1), t\geqslant 0.

它满足

\left( {{\rm d}^2}{{\rm d}t^2}-{\mu_1}{1+t}{\rm d}{{\rm d}t}+{\mu_1+\mu_2^2}{(1+t)^2}-1 \right)\lambda(t)=0, t>0.

仿照 Yordanov 和 Zhang[32], 引入函数

\begin{equation} \varphi:= \begin{cases} \int_{\mathbb{S}^{n-1}}{\rm e}^{x\cdot \omega}{\rm d} \omega,&n\geqslant 2,\\ {\rm e}^{x}+{\rm e}^{-x},&n=1, \end{cases} \end{equation}
(4.3)

这里 \mathbb{S}^{n-1}n-1 维单位球面. 函数 \varphi 满足

\Delta \varphi (x)=\varphi (x), x \in \mathbb{R}^n

\varphi (x)\thicksim C_n|x|^{-{n-1}{2}}{\rm e}^{|x|}, |x|\rightarrow \infty.

利用文献 [25,引理 2.1] 的证明方法, 注意到系统 (1.1) 中非线性项的非负性, 选取

\Phi (t,x) = \lambda(t) \varphi(x)

作为试验函数去推导 \int_{\mathbb{R}^n} |u(t,x)|^q{\rm d}x\int_{\mathbb{R}^n} |v(t,x)|^p{\rm d}x 的下界,可得

引理 4.1 假定u_0,u_1,v_0,v_1 非负且其紧支集含在球 B_R(0) 中, 满足

\begin{equation} u_0(x)\geqslant 0 \mbox{且} u_1(x)+{\mu_1-1-\sqrt{\delta}}{2}u_0(x)\geqslant0, \end{equation}
(4.4)
\begin{equation} v_0(x)\geqslant 0 \mbox{且} v_1(x)+{\mu_1-1-\sqrt{\delta}}{2}v_0(x)\geqslant0, \end{equation}
(4.5)

则 (1.1) 式的局部解 (u,v) 满足

\text{supp} u,\ \text{supp} v \subset \{(t,x)\in [0,T)\times \mathbb{R}^n | |x|\leqslant t+R\}.

同时存在一个充分大且不依赖于 u_0,u_1,v_0,v_1T_0, 使得对任意的 t>T_0, 当 p,q>1 时, 成立

\begin{equation} \int_{\mathbb{R}^n} |u(t,x)|^q {\rm d}x\geqslant C_1(1+t)^{n-1-{n+\mu_1-1}{2}q}, \end{equation}
(4.6)
\begin{equation} \int_{\mathbb{R}^n} |v(t,x)|^p {\rm d}x\geqslant \widetilde {C_1}(1+t)^{n-1-{n+\mu_1-1}{2}p}. \end{equation}
(4.7)

4.2 迭代分析

为了研究解的破裂, 我们将研究如下与解相关的泛函的演化情况

\begin{equation}\label{U V def} U(t)\triangleq \int_{\mathbb{R}^n}u(t,x){\rm d}x, V(t)\triangleq \int_{\mathbb{R}^n}v(t,x){\rm d}x. \end{equation}
(4.8)

将证明: 当变量 t 趋于某个有限时间, 这两个泛函趋于无穷,从而解 (u,v) 必在有限时刻内破裂. 以下证明分为四个步骤进行.

首先, 推导迭代框架. 在 (3.2) 和 (2.4) 式中选取试验函数 \phi \psi 满足

\phi =\psi \equiv 1 \mbox{于} \big\{ (s,x)\in [t]\times \mathbb{R}^n:|x|\leqslant s+R\big\},

使得

\begin{equation}\label{U u} \begin{split} &\int_{\mathbb{R}^n}u_t(t,x){\rm d}x-\int_{R^n}u_t(0,x){\rm d}x\\ &+\int_0^t\int_{\mathbb{R}^n}\bigg (\frac{\mu_1}{1+s}u_s(s,x)+\frac{{\mu_2}^2}{(1+s)^{2}}u(s,x)\bigg ){\rm d}x {\rm d}s\\ =&\int_0^t\int_{\mathbb{R}^n}\int_0^s g_1(s-\tau)|v(\tau,x)|^p{\rm d}\tau {\rm d}x{\rm d}s \end{split} \end{equation}
(4.9)

\begin{equation}\label{V v} \begin{split} &\int_{\mathbb{R}^n}v_t(t,x){\rm d}x-\int_{R^n}v_t(0,x){\rm d}x\\ &+\int_0^t\int_{\mathbb{R}^n}\bigg (\frac{\mu_1}{1+s}v_s(s,x)+\frac{{\mu_2}^2}{(1+s)^{2}}v(s,x)\bigg ){\rm d}x {\rm d}s\\ =&\int_0^t\int_{\mathbb{R}^n}\int_0^s g_2(s-\tau)|u(\tau,x)|^q{\rm d}\tau {\rm d}x{\rm d}s, \end{split} \end{equation}
(4.10)

此即

\begin{equation}\label{U} \begin{split} & U^{'}(t)- U^{'}(0)+\int_0^t \frac{\mu_1}{1+s}U^{'}(s){\rm d}s+\int_0^t\frac{{\mu_2}^2}{(1+s)^{2}}U(s){\rm d}s\\ & =\int_0^t\int_{\mathbb{R}^n}\int_0^s g_1(s-\tau)|v(\tau,x)|^p{\rm d}\tau {\rm d}x{\rm d}s \end{split} \end{equation}
(4.11)

\begin{equation}\label{V} \begin{split} & V^{'}(t)- V^{'}(0)+\int_0^t\int_{\mathbb{R}^n} \frac{\mu_1}{1+s}V^{'}(s){\rm d}s+\int_0^t\frac{{\mu_2}^2}{(1+s)^{2}}V(s){\rm d}s\\ & =\int_0^t\int_{\mathbb{R}^n}\int_0^s g_2(s-\tau)|u(\tau,x)|^q{\rm d}\tau {\rm d}x{\rm d}s. \end{split} \end{equation}
(4.12)

t 求导得

\begin{equation}\label{U''} \begin{split} U^{''}(t)+\frac{\mu_1}{1+t}U^{'}(t)+\frac{{\mu_2}^2}{(1+t)^{2}}U(t)=\int_{\mathbb{R}^n}\int_0^t g_1(t-\tau)|v(\tau,x)|^p{\rm d}\tau {\rm d}x, \end{split} \end{equation}
(4.13)
\begin{equation}\label{V''} \begin{split} V^{''}(t)+\frac{\mu_1}{1+t}V^{'}(t)+\frac{{\mu_2}^2}{(1+t)^{2}}V(t)=\int_{\mathbb{R}^n}\int_0^t g_2(t-\tau)|u(\tau,x)|^q{\rm d}\tau {\rm d}x. \end{split} \end{equation}
(4.14)

二次方程

r^2-(\mu_1-1)r+\mu_2^2=0

\delta=(\mu_1-1)^2-4\mu_2^2> 0 时有两个互异实根

\begin{equation} r_{1,2}=\frac{\mu_{1}-1\mp\sqrt{\delta}}{2}. \end{equation}
(4.15)

显然有

\begin{equation} \left. \begin{aligned} &\quad \mu_{1}>1 &\Rightarrow &\quad r_{1,2}>0\\ &0\leqslant\mu_{1}< 1&\Rightarrow &-1< r_{1,2}< 0\\ &\quad \mu_{1}=1&\Rightarrow &\quad r_{1,2}=0 \\ \end{aligned} \right\} \Rightarrow r_{1,2}+1>0. \end{equation}
(4.16)

从而 (4.13) 式可改写为

\begin{equation}\label{RE U''} \begin{split} & \bigg(U^{'}(t)+\frac{r_1}{1+t}U(t)\bigg)^{'}+\frac{{r_2}+1}{1+t}\bigg(U^{'}(t)+\frac{r_1}{1+t}U(t)\bigg)\\ &=\int_{R^n}\int_0^t g_1(t-\tau)|v(\tau,x)|^p{\rm d}\tau {\rm d}x. \end{split} \end{equation}
(4.17)

两边同乘 (1+t)^{r_2+1}, 在 [t] 上积分, 利用 (2.6) 式得

\begin{equation} U^{'}(t) +\frac{r_1}{1+t}U(t) \geqslant(1+t)^{-r_2-1}\int_0^t(1+s)^{r_2+1} \int_{\mathbb{R}^n}\int_0^s g_1(s-\tau)|v(\tau,x)|^p{\rm d}\tau {\rm d}x{\rm d}s. \end{equation}
(4.18)

两边再次同乘 (1+t)^{r_2+1} 并在 [t] 上积分, 则有

\begin{equation} \begin{split} & (1+t)^{r_1}U(t)+U(0)\\ &\geqslant \int_0^t(1+\tau)^{r_1-r_2-1}\int_0^{\tau}(1+s)^{r_2+1}\int_{\mathbb{R}^n}\int_0^s g_1(s-\eta)|v(\eta,x)|^p {\rm d}\eta {\rm d}x{\rm d}s {\rm d}\tau, \end{split} \end{equation}
(4.19)

这里由 u_0(x)\geqslant 0 可知

U(0)=\int_{\mathbb{R}^n}u(0,x){\rm d}x=\int_{\mathbb{R}^n}u_0dx 是非负的, 因此

\begin{equation}\label{frame} \begin{split} U(t)&\geqslant (1+t)^{-r_1}\\ & \cdot \int_0^t(1+\tau)^{r_1-r_2-1}\int_0^{\tau}(1+s)^{r_2+1}\int_0^s g_1(s-\eta)\int_{\mathbb{R}^n} |v(\eta,x)|^p {\rm d}x {\rm d}\eta {\rm d}s {\rm d}\tau. \end{split} \end{equation}
(4.20)

注意到 r_1-r_2-1<0, 将

\begin{equation} \int_{\mathbb{R}^n}|v(\eta,x)|^p {\rm d}x\geqslant C_0|V(\eta)|^p \big(\int_{B_{R+\eta}}dx \big)^{-(p-1)}=C_0(\eta +R)^{-n(p-1)}|V(\eta)|^p \end{equation}
(4.21)

代入到 (4.20) 式, 得

\begin{equation}\label{Di Tui U} U(t)\geqslant C_0(1+t)^{-r_2-1}\int_0^t\int_0^{\tau}(1+s)^{r_2+1-n(p-1)}\int_0^s g_1(s-\eta)|V(\eta)|^p {\rm d}\eta {\rm d}s {\rm d}\tau. \end{equation}
(4.22)

类似地, 可得

\begin{equation}\label{Di Tui V} \begin{split} &V(t)\geqslant C_0(1+t)^{-r_2-1}\int_0^t\int_0^{\tau}(1+s)^{r_2+1-n(q-1)}\int_0^s g_2(s-\eta)|U(\eta)|^q {\rm d}\eta {\rm d}s {\rm d}\tau. \end{split} \end{equation}
(4.23)

(4.22) 和 (4.23) 式两个不等式将作为迭代框架在后续证明中使用.

其次, 推导 U(t)V(t) 的第一个下界估计. 对 g_k\in C^1([0,\infty)), k=1,2, 引入

\begin{equation} G_k(t):=\int_0^t g_k(\eta){\rm d}\eta. \end{equation}
(4.24)

它满足

\bulletG_k^{'}(t)=g_k(t)>0,\ G_k(0)=0.

\bulletG_k(t) 是严格单调递增函数且 G_k(t)\geqslant 0.

\bullet 对任意的 t\geqslant t_0>0, \int_0^tg_k(t-\eta)\eta^{\alpha}{\rm d}\eta= \int_0^tg_k(\eta)(t-\eta)^{\alpha}{\rm d}\eta

\begin{equation}\label{G(t)} \begin{split} =[G_k(\eta)(t-\eta)^{\alpha}]\big|_{\eta=0}^{\eta=t}+\alpha \int_0^t G_k(\eta)(t-\eta)^{\alpha -1}{\rm d} \eta \quad\quad\quad\quad\quad\quad\,\,\\ \geqslant\, \begin{cases} G_k(t_0),& \alpha=0,\\ \alpha \int_{t_0}^{t} G_k(\eta)(t-\eta)^{\alpha -1}{\rm d}\eta, &\alpha >0,\,\,\\ \end{cases} \geqslant G_k(t_0)(t-t_0)^{\alpha},\quad\quad\, \end{split} \end{equation}
(4.25)

这里 \alpha \geqslant 0.

由 (4.20) 式, 引理 4.1 以及 (4.25) 式, 当t\geqslant \max\{T_0, t_0\} 时, 可得

\begin{aligned}U(t) \geqslant & \widetilde{C_{1}}(1+t)^{-r_{1}} \\& \times \int_{0}^{t}(1+\tau)^{r_{1}-r_{2}-1} \int_{0}^{\tau}(1+s)^{r_{2}+1} \int_{0}^{s} g_{1}(s-\eta)(1+\eta)^{n-1-\frac{n+\mu_{1}-1}{2} p} \mathrm{~d} \eta \mathrm{~d} s \mathrm{~d} \tau \\\geqslant & \widetilde{C_{1}}(1+t)^{-r_{2}-1-\frac{n+\mu_{1}-1}{2} p} \int_{0}^{t} \int_{0}^{\tau}(1+s)^{r_{2}+1} \mathrm{~d} s \mathrm{~d} \tau \int_{0}^{s} g_{1}(s-\eta) \eta^{n-1} \mathrm{~d} \eta \\\geqslant & \widetilde{C_{1}}(1+t)^{-r_{2}-1-\frac{n+\mu_{1}-1}{2} p} G_{1}\left(t_{0}\right) \int_{t_{0}}^{t} \int_{t_{0}}^{\tau}\left(s-t_{0}\right)^{r_{2}+1}\left(s-t_{0}\right)^{n-1} \mathrm{~d} s \mathrm{~d} \tau \\\geqslant & \widetilde{C_{1}} G_{1}\left(t_{0}\right) \\\left(n+r_{2}+2\right)\left(n+r_{2}+1\right) & (1+t)^{-r_{2}-1-\frac{n+\mu_{1}-1}{2} p}\left(t-t_{0}\right)^{r_{2}+n+2}\end{aligned}
(4.26)

类似地, 当 t\geqslant \max\{T_0, t_0\} 成立

\begin{equation} V(t)\geqslant {C_1G_2(t_0)}{(n+r_2+2)(n+r_2+1)} (1+t)^{-r_2-1-{n+\mu_1-1}{2}q} (t-t_0)^{r_2+n+2}. \end{equation}
(4.27)

下面, 开始迭代过程. 假设

\begin{equation}\label{U first item} U(t)\geqslant Q_j (1+t)^{-\theta_{j}} (t-L_{j}t_0)^{\sigma_j}, \text{对} t\geqslant L_{j}t_0 \end{equation}
(4.28)

\begin{equation}\label{V first item} V(t)\geqslant \widetilde{Q}_{j} (1+t)^{-\widetilde{\theta}_{j}} (t-L_{j}t_0)^{\widetilde{\sigma}_{j}}, \text{对} t\geqslant L_{j}t_0, \end{equation}
(4.29)

这里 \{Q_j\}_{j\geqslant 1},\{\widetilde{Q}_j\}_{j\geqslant 1},\{\theta_{j}\}_{j\geqslant 1},\{\widetilde{\theta}_{j}\}_{j\geqslant 1},\{\sigma_j\}_{j\geqslant 1},\{\widetilde{\sigma}_j\}_{j\geqslant 1}是非负实数列 (在后续归纳步骤中待定). 当 j=1 时, (4.26)-(4.27) 式给出

\begin{equation} \begin{array}{ll} Q_1={\widetilde{C_1}G_1(t_0)}{(n+r_2+2)(n+r_2+1)}, &\widetilde{Q}_{1}={C_1G_2(t_0)}{(n+r_2+2)(n+r_2+1)},\\ \theta_{1}=r_2+1+{(n+\mu_1-1)}{2}p, &\widetilde{\theta}_{1}=r_2+1+{(n+\mu_1-1)}{2}q,\\ \sigma_1=r_2+n+2, & \widetilde{\sigma}_{1}=r_2+n+2. \end{array} \end{equation}
(4.30)

受 Chen[2] 的启发, 构造 \{L_j\}_{j\geqslant1} 如下

\begin{equation}\label{Lj} L_j:=\prod_{k=0}^{j}l_k, \ \forall j\geqslant 1, \ \text{其中}\ l_k:=1+(pq)^{-{k-1}{2}}>1, \forall k\geqslant 1. \end{equation}
(4.31)

易知 L_j 是单调递增的, \prod_{k=0}^{\infty}l_k 是收敛的 (由比式判别法可证), \lim\limits_{k\to \infty}(\ln l_{k+1})/(\ln l_k)=(pq)^{-{1}{2}} < 1; 进一步, 利用数学归纳法可知 L\triangleq\prod_{k=0}^{\infty}l_k>1L_j \in [L].

将 (4.29) 式代入 (4.22) 式, 得

\begin{equation}\label{Uj+1 2} \begin{split} U(t)&\geqslant C_0(1+t)^{-r_2-1}\int_0^t \int_0^{\tau}(1+s)^{r_2+1-n(p-1)}\\ &\quad\int_0^s g_1(s-\eta)\big|\widetilde{Q}_{j} (1+\eta)^{-\widetilde{\theta}_{j}} (\eta-L_{j}t_0)^{\widetilde{\sigma}_{j}}\big|^p {\rm d}\eta {\rm d}s {\rm d}\tau\\ &\geqslant C_0\widetilde{Q}_{j}^p(1+t)^{-r_2-1-n(p-1)-p \widetilde{\theta}_{j}}\\ & \times \int_{L_j t_0}^t \int_{L_j t_0}^{\tau} (s-L_{j+1}t_0)^{r_2+1} \int_{L_j t_0}^s g_1(s-\eta)(\eta-L_{j}t_0)^{p\widetilde{\sigma}_{j}}{\rm d}\eta {\rm d}s {\rm d}\tau. \end{split} \end{equation}
(4.32)

换元并分部积分, 对任意的 s\geqslant L_{j+1}t_0, 可得

\begin{equation} \begin{split} \int_{L_j t_0}^s g_1(s-\eta)(\eta-L_{j}t_0)^{p\widetilde{\sigma}_{j}}{\rm d}\eta &=\int_{0}^{s-L_jt_0}g_1(\eta)(s-L_jt_0-\eta)^{p\widetilde{\sigma}_{j}}{\rm d}\eta\\ &=p\widetilde{\sigma}_{j} \int_{0}^{s-L_jt_0}G_1(\eta)(s-L_jt_0-\eta)^{p\widetilde{\sigma}_{j}-1}{\rm d}\eta\\ &\geqslant p \widetilde{\sigma}_{j}\int_{L_jt_0(l_{j+1}-1)}^{s-L_jt_0}G_1(\eta)(s-L_jt_0-\eta)^{p\widetilde{\sigma}_{j}-1}{\rm d}\eta\\ &\geqslant G_1\big(L_jt_0(l_{j+1}-1)\big)\cdot(s-L_jt_0 l_{j+1})^{p\widetilde{\sigma}_{j}}, \end{split} \end{equation}
(4.33)

这里用到了

s\geqslant L_{j+1}t_0\Rightarrow s-L_{j}t_0\geqslant L_jt_0(l_{j+1}-1)=L_jt_0(pq)^{-{j}{2}}>0.

注意到 g_1(t) \in C^1([0, \infty)), L>1, L_j \in [L] (j\geqslant 1), 可知存在正整数 j_0 使得当 j\geqslant j_0 时, 成立

0< L_jt_0(l_{j+1}-1)=L_jt_0(pq)^{-{j}{2}}\ll1.

因此

\begin{equation} \begin{split} & G_1(L_jt_0(l_{j+1}-1))\\ &=G_1(0)+g_1(\xi)L_jt_0(l_{j+1}-1)\\ &=G_1(0)+g_1(0)L_jt_0(pq)^{-{j}{2}}+g_1^{'}(\zeta)O\bigg( \big( (pq)^{-{j}{2}}L_jt_0 \big)^2\bigg)\\ &= (pq)^{-j}L_jt_0\bigg((pq)^{{j}{2}} g_1(0)+g_1^{'}(\zeta)O(1) \bigg)\geqslant C_2(pq)^{-j},\ j\geqslant \max\{j_0,j_1\}, \end{split} \end{equation}
(4.34)

式中 \xi \in (0,L_jt_0(pq)^{-{j}{2}}), \zeta \in (0, \xi) , 选取的正整数 j_1 保证了 (pq)^{{j}{2}} g_1(0)+g_1^{'}(\zeta)O(1)\geqslant C >0 (由于 g_1(0)>0). 联合 (4.32)-(4.34) 式, 对 t \geqslant L_{j+1}t_0, j \geqslant \max\{j_0,j_1\}, 得到

\begin{align} U(t)&\geqslant C_0\widetilde{Q}_{j}^p(1+t)^{-r_2-1-n(p-1)-p\widetilde{\theta}_{j}}\nonumber\\ &\quad\times \int_{L_j t_0}^t \int_{L_j t_0}^{\tau} (s-L_{j+1}t_0)^{r_2+1+p\widetilde{\sigma}_{j}}C_3 (pq)^{-j}{\rm d}s{\rm d}\tau \\ &\geqslant {C_0C_2(pq)^{-j}\widetilde{Q}_{j}^p}{(r_2+3+p\widetilde{\sigma}_{j})(r_2+2+p\widetilde{\sigma}_{j})} (1+t)^{-r_2-1-n(p-1)-p\widetilde{\theta}_{j}}(t-L_{j+1}t_0)^{r_2+3+p\widetilde{\sigma}_{j}}\nonumber. \end{align}
(4.35)

t \geqslant L_{j+1}t_0, j \geqslant \max\{j_0,\widetilde{j_1}\}, 类似可得

\begin{equation} V(t)\geqslant {\widetilde{C_0}\widetilde{C_2}(pq)^{-j}Q_{j}^q}{(r_2+3+ q\sigma_{j})(r_2+2+q\sigma_{j})}(1+t)^{-r_2-1-n(q-1)-q\theta_j} \times (t-L_{j+1}t_0)^{r_2+3+q\sigma_{j}}, \end{equation}
(4.36)

其中选取 \widetilde{j_1} 使得 (pq)^{{j}{2}} g_2(0)+g_2^{'}(0)O(1)\geqslant C>0 (由于 g_2(0)>0).从而, 令

\begin{equation}\label{j+1 2} \begin{split} &Q_{j+1}:={C_0C_2(pq)^{-j}\widetilde{Q_j}^p}{(r_2+3+p\widetilde{\sigma}_{j}) (r_2+2+p\widetilde{\sigma}_{j})}, \\ &\widetilde{Q}_{j+1}:={\widetilde{C_0}\widetilde{C_2}(pq)^{-j}Q_j^q}{(r_2+3+q\sigma_{j}) (r_2+2+q\sigma_{j})},\\ &\theta_{j+1}:=r_2+1+n(p-1)+p\widetilde{\theta}_j,\quad\quad\widetilde{\theta}_{j+1}:=r_2+1+n(q-1)+q\theta_{j},\\ &\sigma_{j+1}:=p\widetilde{\sigma}_j+r_2+3, \quad\quad\quad\quad\quad\quad\,\,\,\,\widetilde{\sigma}_{j+1}:=q\sigma_j+r_2+3, \end{split} \end{equation}
(4.37)

则迭代过程 (4.28) 与 (4.29) 式是合理的.

最后, 证明解 (u,v) 必在有限时刻内破裂. 对任意奇数j\geqslant\max\{j_0, j_1, \widetilde{j_1}\}, 有

\begin{equation}\label{theta j} \begin{split} &\theta_{j}= \left({A_1}{pq-1}+\theta_1\right) (pq)^{{j-1}{2}}-{A_1}{pq-1}, A_1=r_2+1-n+(r_2+1)p+npq,\\ &\widetilde{\theta}_{j}= \left({\widetilde{A}_1}{pq-1}+\widetilde{\theta}_1\right) (pq)^{{j-1}{2}} -{\widetilde{A}_1}{pq-1}, \widetilde{A}_1=r_2+1-n+(r_2+1)q+npq \end{split} \end{equation}
(4.38)

\begin{equation}\label{sigma j} \begin{split} &\sigma_{j}=\left({B_1}{pq-1}+\sigma_1\right) (pq)^{{j-1}{2}}-{B_1}{pq-1}, B_1=r_2+3+(r_2+3)p,\\ &\widetilde{\sigma}_{j}= \left({\widetilde{B}_1}{pq-1}+\widetilde{\sigma}_1\right) (pq)^{{j-1}{2}} -{\widetilde{B}_1}{pq-1}, \widetilde{B}_1=r_2+3+(r_2+3)q. \end{split} \end{equation}
(4.39)

同样地, 对于任意偶数 j\geqslant\max\{j_0, j_1, \widetilde{j_1}\}, 则有

\begin{aligned} \sigma_{j}&=p\widetilde{\sigma}_{j-1}+r_2+3 =p \left[\left( {\widetilde{B}_1}{pq-1}+\widetilde{\sigma}_1\right)(pq)^{{j-2}{2}} -{\widetilde{B}_1}{pq-1}\right] +r_2+3\\ &={1}{q} \left({\widetilde{B}_1}{pq-1}+\widetilde{\sigma}_1\right) (pq)^{{j}{2}} -{p\widetilde{B}_1}{pq-1}+r_2+3\\ &={1}{q} \left({\widetilde{B}_1}{pq-1}+\widetilde{\sigma}_1\right) (pq)^{{j}{2}} +{-p(r_2+3)-(r_2+3)}{pq-1} \end{aligned}

\widetilde{\sigma}_{j}={1}{p} \left({B_1}{pq-1}+\sigma_1 \right) (pq)^{{j}{2}} +{-q(r_2+3)-(r_2+3)}{pq-1}.

从而对任意的正整数 j\geqslant \max\{j_0,j_1,\widetilde{j_1}\}, 存在与 j 无关的正常数 B_1^{'}, \widetilde{B_1^{'}} 使得

\sigma_{j}< B_1^{'}(pq)^{\frac{j}{2}}, \widetilde{\sigma}_{j}<\widetilde{B_1^{'}}(pq)^{\frac{j}{2}}.

由 (4.37) 式可知

\begin{equation} \begin{split} Q_{j}& \geqslant {C_0C_2(pq)^{-(j-1)}}{B_{1}^{'2}(pq)^{j}}\widetilde{Q}_{j-1}^p \geqslant {C_0\widetilde{C_0}^pC_2\widetilde{C_2}^p(pq)^{3p+1}}{B_{1}^{'2}\widetilde{B_{1}^{'}}^{2p}} (pq)^{-2j(p+1)}Q_{j-2}^{pq} \\ &\triangleq E_2(pq)^{-2j(p+1)}Q_{j-2}^{pq} \end{split} \end{equation}
(4.40)

以及

\begin{equation} \widetilde{Q}_{j}\geqslant {C_0^q\widetilde{C_0}C_2^q\widetilde{C_2}(pq)^{3p+1}}{B_{1}^{'2p}\widetilde{B_{1}^{'}}^{2}} (pq)^{-2j(q+1)}\widetilde{Q}_{j-2}^{pq} \triangleq \widetilde{E}_2(pq)^{-2j(q+1)}\widetilde{Q}_{j-2}^{pq}. \end{equation}
(4.41)

故对任意奇数 j\geqslant \max\{j_0,j_1,\widetilde{j_1}\}, 成立

\begin{equation} \begin{split} \ln Q_j&\geqslant (pq)^{{j-1}{2}} \ln Q_1-2(p+1) \ln (pq)\sum_{k=0}^{{j-3}{2}}((j-2k)(pq)^k) + \ln E_2 \sum_{k=0}^{{j-3}{2}}(pq)^k\\ &=(pq)^{{j-1}{2}} \bigg(\ln Q_1-{2(p+1)(3pq-1) \ln(pq)}{(pq-1)^2}+{\ln E_2}{pq-1}\bigg) \\ &\quad+{2(p+1)(2pq+j(pq-1)) \ln (pq)}{(pq-1)^2}-{\ln E_2}{pq-1} \end{split} \end{equation}
(4.42)

\begin{equation} \begin{split} \ln \widetilde{Q}_j&\geqslant (pq)^{{j-1}{2}} \bigg(\ln \widetilde{Q}_1-{2(q+1)(3pq-1) \ln (pq)}{(pq-1)^2}+{\ln \widetilde{E}_2}{pq-1}\bigg) \\ &\quad +{2(q+1)(2pq+j(pq-1))\ln (pq)}{(pq-1)^2}-{\ln \widetilde{E}_2}{pq-1}. \end{split} \end{equation}
(4.43)

若取 j_2, \widetilde{j_2} 为分别满足

{2(p+1)(2pq+j(pq-1)) \ln (pq)}{(pq-1)^2}-{\ln E_2}{pq-1}>0,
{2(q+1)(2pq+j(pq-1))\ln (pq)}{(pq-1)^2}-{\ln \widetilde{E}_2}{pq-1}>0,

的最小正整数, 则对任意的奇数

j\geqslant \max\{j_0,j_1,\widetilde{j_1},j_2,\widetilde{j_2}\}, 成立

\begin{equation}\label{logQj} \begin{split} \ln Q_j\geqslant (pq)^{{j-1}{2}} \ln \big(Q_1(pq)^{-{2(p+1)(3pq-1)}{(pq-1)^2}}E_2^{{1}{pq-1}}\big) \triangleq (pq)^{{j-1}{2}}\ln E_3 \end{split} \end{equation}
(4.44)

\begin{equation} \begin{split} \ln \widetilde{Q}_j\geqslant (pq)^{{j-1}{2}} \ln \big(\widetilde{Q}_1(pq)^{-{2(q+1)(3pq-1)}{(pq-1)^2}}\widetilde{E}_2^{{1}{pq-1}}\big) \triangleq (pq)^{{j-1}{2}} \ln \widetilde{E_3}, \end{split} \end{equation}
(4.45)

其中 E_3,\ \widetilde{E_3} 是与 j 无关的正常数.

t\geqslant \max\{T_0,1,2Lt_0\} 时, 易得

\begin{equation} \ln (1+t)\leqslant \ln (2t),\ \ln (t-Lt_0)\geqslant \ln({t}{2}), \end{equation}
(4.46)

从而对任意的奇数

j\geqslant \max\{j_0,j_1,\widetilde{j_1},j_2,\widetilde{j_2}\}, 注意到

(4.38), (4.39) 及 (4.44) 式, 考虑到 L>1L_j \in [L], (见 (4.31) 式), 从 (4.28) 式中可推出: 当 t\geqslant\max\{T_0,1,2Lt_0\} 时,

\begin{equation}\label{U2 lower} \begin{split} U(t)&\geqslant \exp\big((pq)^{{j-1}{2}}\ln E_3\big)\\ &\quad \times(1+t)^{-({A_1}{pq-1}+\theta_1)(pq)^{{j-1}{2}}+{A_1}{pq-1}}(t-t_0)^{({B_1}{pq-1}+\sigma_1)(pq)^{{j-1}{2}}-{B_1}{pq-1}}\\ &= \exp\bigg((pq)^{{j-1}{2}}\ln \big(E_3(1+t)^{-({A_1}{pq-1}+\theta_1)}(t-t_0)^{({B_1}{pq-1}+\sigma_1)}\big)\bigg)\\ & \quad\times (1+t)^{{A_1}{pq-1}} (t-t_0)^{-{B_1}{pq-1}}\\ &\geqslant \exp\bigg((pq)^{{j-1}{2}}J_2(t)\bigg)(1+t)^{{A_1}{pq-1}}(t-t_0)^{-{B_1}{pq-1}}, \end{split} \end{equation}
(4.47)

其中

\begin{equation} \begin{split} J_2(t)&\triangleq \ln\big(E_3 2^{-({B_1}{pq-1}+\sigma_1)-({A_1}{pq-1}+\theta_1)} t^{({B_1}{pq-1}+\sigma_1)-({A_1}{pq-1}+\theta_1)}\big)\\ &=\ln \big(E_3 2^{-{2(r_2+2)(p+1)}{pq-1}-({(n+\mu_1-1)}{2}p+2r_2+2n+3)} t^{-{n+\mu_1-1}{2}p+1+{2(p+1)}{pq-1}} \big). \end{split} \end{equation}
(4.48)

类似地, 可得到: 当 t\geqslant\max\{T_0,1,2Lt_0\} 时,

\begin{equation}\label{V2 lower} \begin{split} V(t)\geqslant \exp\bigg((pq)^{{j-1}{2}}\widetilde{J_2}(t) \bigg)(1+t)^{{\widetilde{A_1}}{pq-1}}(t-t_0)^{-{\widetilde{B_1}}{pq-1}}, \end{split} \end{equation}
(4.49)

其中

\begin{equation} \begin{split} \widetilde{J_2}(t)&\triangleq \ln\big(\widetilde{E}_3 2^{-({\widetilde{B}_1}{pq-1}+\widetilde{\sigma}_1)-({\widetilde{A}_1}{pq-1}+\widetilde{\theta}_1)} t^{({\widetilde{B}_1}{pq-1}+\widetilde{\sigma}_1)-({\widetilde{A}_1}{pq-1}+\widetilde{\theta}_1)}\big)\\ &=\ln \big(\widetilde{E}_3 2^{-{2(r_2+2)(q+1)}{pq-1}-({(n+\mu_1-1)}{2}q+2 r_2+2n+3)} t^{-{n+\mu_1-1}{2}q+1+{2(q+1)}{pq-1}} \big). \end{split} \end{equation}
(4.50)

对任意的奇数

j\geqslant \max\{j_0,j_1,\widetilde{j_1},j_2,\widetilde{j_2}\}, 条件 (2.5) 保证了函数 J_2(t)\widetilde{J_2}(t)t 的幂次是正的, 从而存在 t_1,\ \widetilde{t_1} 使得

\begin{equation}\label{t1 range} t_1^{-{n+\mu_1-1}{2}p+1+{2(p+1)}{pq-1}}>\bigg(E_3 2^{-{2(r_2+2)(p+1)}{pq-1}-({(n+\mu_1-1)}{2}p+2r_2+2n+3)}\bigg)^{-1} \end{equation}
(4.51)

\begin{equation}\label{t1' range} \widetilde{t_1}^{-{n+\mu_1-1}{2}q+1+{2(q+1)}{pq-1}}>\bigg(\widetilde{E}_3 2^{-{2(r_2+2)(q+1)}{pq-1}-({(n+\mu_1-1)}{2}q+2 r_2+2n+3)}\bigg)^{-1}. \end{equation}
(4.52)

由此, 若取 t\geqslant \max\{T_0,1,2Lt_0\}t\geqslant \min \{t_1,\widetilde{t_1}\}, 则当 j\to \infty 时, (4.47) 和 (4.59) 式中的泛函 U(t),\ V(t) 的下界必在有限时刻内破裂, 这就完成了定理 2.2 的证明.

参考文献

Agemi R, Kurokawa Y, Takamura H.

Critical curve for p-q systems of nonlinear wave equations in three space dimensions

J Differential Equations, 2000, 167: 87-133

[本文引用: 1]

Chen W H.

Interplay effects on blow-up of weakly coupled systems for semilinear wave equations with general nonlinear memory terms

Nonlinear Anal, 2021, 202: 112160

[本文引用: 6]

Chen W H, Palmieri A. Blow-up result for a semilinear wave equation with a nonlinear memory term.// Cicognani M, Del Santo D, Parmeggiani A, Reissig M. Anomalies in Partial Differential Equations. Cham: Springer, 2021, 43: 77-97

[本文引用: 1]

Chen W H, Reissig M.

Blow-up of solutions to Nakao's problem via an iteration argument

J Differential Equations, 2021, 275: 733-756

[本文引用: 1]

D'Abbicco M.

The threshold of effective damping for semilinear wave equations

Math Methods Appl Sci, 2015, 38(6): 1032-1045

[本文引用: 1]

Friedman A. Partial Differential Equations. New York: Krieger, 1976

[本文引用: 1]

Georgiev V, Lindblad H, Sogge C D.

Weighted Strichartz estimates and global existence for semilinear wave equations

Amer J Math, 1997, 119(6): 1291-1319

[本文引用: 1]

Glassey R T.

Existence in the large for u_{tt}-\Delta u = F (u) in two space dimensions

Math Z, 1981, 178: 233-261

Glassey R T.

Finite-time blow-up for solutions of nonlinear wave equations

Math Z, 1981, 177(3): 323-340

[本文引用: 1]

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(3/4): 1017-1040

[本文引用: 1]

Ikeda M, Sobajima M, Wakasa K.

Blow-up phenomena of semilinear wave equations and their weakly coupled systems

J Differential Equations, 2019, 267(9): 5165-5201

[本文引用: 1]

Jiao H, Zhou Z.

An elementary proof of the blow-up for semilinear wave equation in high space dimensions

J Differential Equations, 2003, 189(2): 355-365

[本文引用: 1]

John F.

Blow-up of solutions of nonlinear wave equations in three space dimensions

Manuscripta Math, 1979, 28: 235-268

[本文引用: 1]

Kurokawa Y, Takamura H, Wakasa K.

The blow-up and lifespan of solutions to systems of semilinear wave equation with critical exponents in high dimensions

Differential Integral Equations, 2012, 25(3/4): 363-382

[本文引用: 1]

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(9): 5377-5394

[本文引用: 1]

Li Y C, Wang C M.

On a partially synchronizable system for a coupled system of wave equations in one dimension

Commun Anal Mech, 2023, 15(3): 470-493

[本文引用: 1]

Lindblad H, Sogge C D.

Long-time existence for small amplitude semilinear wave equations

Amer J Math, 1996, 118(5): 1047-1135

[本文引用: 1]

Martin H.

Mathematical analysis of some iterative methods for the reconstruction of memory kernels

Electron Trans, 2021, 54: 483-498

[本文引用: 1]

Narazaki T.

Global solutions to the Cauchy problem for the weakly coupled system of damped wave equations

Dynamical systems, differential equations and applications, 7th AIMS Conference, 2009, special: 592-601

[本文引用: 1]

Nishihara K, Wakasugi Y.

Critical exponent for the Cauchy problem to the weakly coupled damped wave system

Nonlinear Anal, 2014, 108(5): 249-259

[本文引用: 1]

Nunes do Nasci\acute{m}eto W, Palmieri A, Reissig M.

Semi-linear wave models with power non-linearity and scale-invariant time-dependent mass and dissipation

Math Nachr, 2017, 290(11/12): 1779-1805

[本文引用: 1]

Palmieri A.

Global existence of solutions for semi-linear wave equation with scale-invariant damping and mass in exponentially weighted spaces

J Math Anal Appl, 2018, 461(2): 1215-1240

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(2/3): 1176-1220

[本文引用: 1]

Palmieri A, Takamura H.

Blow-up for a weakly coupled system of semilinear damped wave equations in the scattering case with power nonlinearities

Nonlinear Anal, 2019, 187(2): 467-492

[本文引用: 1]

Palmieri A, Tu Z.

Lifespan of semilinear wave equation with scale invariant dissipation and mass and sub-Strauss power nonlinearity

J Math Anal Appl, 2019, 470: 447-469

[本文引用: 4]

Said-Houari B, Kirane M.

Global nonexistence results for a class of hyperbolic systems

Nonlinear Anal, 2011, 74(17): 6130-6143

[本文引用: 1]

Schaeffer J.

The equation u_{tt}-\Delta u = |u|^p for the critical value of p.

Proc Roy Soc Edinburgh, 1985, 101(1/2): 31-44

[本文引用: 1]

Sideris T C.

Nonexistence of global solutions to semilinear wave equations in high dimensions

J Differential Equations, 1984, 52(3): 378-406

[本文引用: 2]

Strauss W A.

Nonlinear scattering theory at low energy

J Functional Analysis, 1981, 41(1): 110-133

[本文引用: 1]

Sun F, Wang M.

Existence and nonexistence of global solutions for a nonlinear hyperbolic system with damping

Nonlinear Anal, 2007, 66(12): 2889-2910

[本文引用: 1]

Wakasugi Y. Critical exponent for the semilinear wave equation with scale invariant damping. Fourier Analysis: Pseudo-differential Operators, Time-Frequency Analysis and Partial Differential Equations. Springer International Publishing, 2014: 375-390

[本文引用: 1]

Yordanov B T, Zhang Q S.

Finite time blow up for critical wave equations in high dimensions

J Funct Anal, 2006, 231(2): 361-374

[本文引用: 4]

Zhou Y.

Cauchy problem for semilinear wave equations in four space dimensions with small initial data

J Partial Differential Equations, 1995, 8(2): 135-144

Zhou Y.

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

Chin Ann Math, 2007, 28(2): 205-212

[本文引用: 1]

/