扩散方程是物理学中描述物质和能量随着时间变化在空间中传播的现象, 其在材料科学、地质勘探中有着广泛的应用. 近些年关于非线性扩散方程的研究主要体现在求其解的全局存在性、解的爆破性和爆破时间的上下界, 见文献 [1],[2],[3],[4]. 分数阶 $p$-Laplace 方程作为整数阶的推广在近些年被很多学者广泛关注, 其理论研究方向多集中于研究爆破解的最大存在时间问题, 见文献 [5],[6],[7],[8],[9],[10],[11],[12],[13],[14],[15],[16],[17]. 在分数阶偏微分方程中引入对数非线性源项也是这些年研究的热点方向, 在描述扩散的过程中引入对数非线性源项可以反映介质的特殊扩散特性, 在生物数学中引入对数非线性源项可以更直接的描述生物种群的增长和衰减, 近些年的研究集中于对其解的全局存在性、 唯一性以及爆破性的探讨, 见文献 [18],[19],[20],[21],[22],[23],[24],[25],[26],[27],[28],[29],[30],[31],[32],[33].

非线性伪抛物型方程是在伪抛物型方程的基础上引入非线性源项的一类偏微分方程, 其研究领域包括非线性扩散, 非线性波的传播以及描述多孔介质中液体的流动. 近年来, 对非线性伪抛物型方程的研究也取得了许多新的进展. 在文献 [33] 中, 作者温兰, 杨晗研究了带有对数非线性源项的 $p$-Laplace 伪抛物方程解的存在性和爆破

作者应用了 Galerkin 近似解的方法证明了方程 (1.1) 局部解的存在性, 然后应用势阱理论得到了初值在稳定集中时方程解的存在性, 之后通过构造辅助函数, 应用 Hölder 不等式和 Young 不等式等方法证明了解在有限时刻的爆破.

分数阶 $p$-Laplace 方程是对整数阶方程更深一层的探讨, 两类方程都使用 Galerkin 近似的方法探究初始能量在适当条件下解的全局存在性. 在探究解的爆破等问题的时候都会采用构造辅助函数并结合微分不等式的方法得到结果, 然而分数阶 $p$-Laplace 方程由于其分数阶导数的存在, 在理论分析上比整数阶更为复杂. 对于此类方程的研究, Liao M 等在文献 [5] 中研究了在任意初始能量下方程 (1.2) 的全局存在性, 解的唯一性和在有限时间内的爆破性

其中函数 $f$ 满足

(a) $f\in C^{1}$ 并且 $f(0)=f^{'}(0)=0$;

(b) 当 $u>0$ 时, 函数 $f(u)$ 单调递增并且是凸的, 当 $u<0$ 时是凹的;

(c) $(q+1)F(u)\le uf(u)$ 并且 $\left | uf(u) \right |\le r\left | F(u) \right |$, 其中 $F(u)=\int_{0}^{u}f(s)\mathrm{d}s$ 并且

其中

(d) $u(uf^{'}(u)-(p-1)f(u))\ge 0$, 等号仅在 $u=0$ 时成立.

在文中, 作者运用了改进势阱理论方法, 设阱的深度为 $d=\inf_{u\in N}J(u)$, 通过构造 Galerkin 近似解求出了方程 (1.2) 在次临界 $J(u_{0})<d$、临界 $J(u_{0})=d$ 和超临界 $J(u_{0})>d$ 状态下解的全局存在性、唯一性和在有限时间内的爆破性. 作者在整数阶的基础上更进一步地研究了分数阶 $p$-Laplace 方程解的性质, 两者在数学性质和物理实际应用领域中都有着重要意义, 但分数阶 $p$-Laplace 方程的研究和应用则更加具有挑战性.

Kirchhoff 型方程源自 Kirchhoff 在 1883 年对弦振动方程的推广, 与经典的 $p$-Laplace 方程不同, Kirchhoff 型方程能够描述更多现实世界中非线性和非均匀的物理现象, 由于其特殊的性质引来了众多学者对其的研究. 在文献 [18] 中, 作者 Zeng F 等研究了下列含有对数非线性源项的分数阶 $p$-Laplace Kirchhoff 型扩散方程的初边值问题. 作者在方程中加入 Kirchhoff 型函数 $M(\cdot )$ 可以更深入的研究各类物理系统中非线性扩散现象

作者利用 Galerkin 近似的方法和势阱理论, 并结合 Sobolev 空间相关知识, 证明了方程 (1.3) 在次临界和临界状态下解的整体存在性和有限时间内的爆破性, 讨论了弱解的增长率, 研究了基态解的性质.

近些年, 在偏微分方程中引入对数非线性源项也是众多学者研究的热门方向, 尤其是在 Kirchhoff 型方程中加入对数非线性源项, 不仅在解的稳定性和复杂分析的研究中更具挑战性, 而且为物理学研究领域提供了更为准确的数学工具, 这使得方程在处理复杂系统的建模中具有更广泛的应用. 在文献 [19] 中, Zhao Q 等研究了下列含有对数非线性源项的 Kirchhoff 型方程

这篇论文中引入对数非线性源项 $\left | u \right |^{q-1}u\text{log}\left | u \right |$ 可以更好的描述复杂的物理现象. 作者通过构造 Galerkin 近似解的方法研究了方程 (1.4) 的解在任意初始能量条件下的全局存在性以及解的爆破性, 得到了全局解的指数衰减率和爆破解的最大存在时间. 同时研究了方程基态解的性质, 并建立了基态解与全局解之间的收敛关系.

本文参考方程 (1.3) 和 (1.4), 在其研究的问题上进一步延伸, 考虑了下列分数阶 $p$-Laplace 扩散方程的初边值问题

其中 $s\in (0,1)$, $2< p< N/s$, $1\le \theta < p_{s}^{*}/p$, $p_{s}^{*}=\frac{Np}{N-sp}$, $\Omega \subset R^{N}(N\ge 1)$ 为一个具有 Lipschitz 光滑边界的有界区域, $u_{0}\ne0$ 为 $\Omega$ 上的初始函数, 我们令 $M(t)=a+bt^{\theta -1}$, 其中 $a=0$, $b=1$, $t=\left [ u \right ]_{s,p}^p$, $(-\Delta) ^s_p$ 为分数阶 $p$-Laplace 算子, 并且满足

本文的具体结构如下：在第二节中, 我们介绍了一些符号以及引理；在第三和第四节中, 我们证明了在次临界 $J(u_{0} )<h$, 和临界 $J(u_{0} )=h$ 状态下解的全局存在性和在有限时间内的爆破性.

为了阐述本文的结果, 我们给出了如下定义

定义1.1 如果函数 $u(x,t)$ 在 $\Omega \times [0,T)$ 上满足 $u\in L^\infty (0,T;W_{0}^{s,p}(\Omega ) )$, $u_t\in L^2(0,T;$ $H^{s}_{0}(\Omega ))$, $u(x,0)=u_0\in W_{0}^{s,p}(\Omega )$, 且对于任意的 $v\in W_{0}^{s,p}(\Omega )$ 在 $t\in(0,T)$ 上都有

其中

则称 $u$ 为方程 (1.5) 在 $\Omega \times [0,T)$ 上的一个弱解.

定义1.2 设 $u$ 为方程 (1.5) 的一个弱解, 若 $u$ 的最大存在时间 $T< \infty$, 并且满足

则称 $u$ 在有限时间内爆破.

在本节我们给出了定理的证明过程中所用到的符号以及引理. 在本文中 $\left \|\cdot \right \| _p$ 来表示 Lebesgue 空间 $L^p(\Omega )$ 的空间范数, 其中 $p\in (0,+\infty )$, 用 $(\cdot,\cdot)$ 来表示 $L^2(\Omega )$ 的内积, $\left [ u \right ]_{s,p}$ 为 $(s,p)\text{-Gagliardo}$ 半范数

定义

其中 $\varrho \in (0,p_{s}^{\ast }-\theta p )$, $D_{\ast }$ 和 $C_{\ast }$ 为 $W_{0}^{s,p}(\Omega )\hookrightarrow L^{\theta p+\varrho } (\Omega )$ 以及 $W_{0}^{s,p}(\Omega )\hookrightarrow L^{\theta p} (\Omega )$ 的最佳嵌入常数. 对于 $u(x)\in W_{0}^{s,p}(\Omega )$, 定义的泛函如下所示

根据 (2.1), (2.2) 式可以得到

定义

其中

下面将给出证明过程中所需要的引理.

引理2.1 ^[18] 如果 $\varrho \in (0,p_{s}^{\ast }-\theta p )$, 则有 $\text{ln}t \le \frac{1}{e\varrho }t ^{\varrho }\text{, }$ 其中 $t\in [1,+\infty )$.

引理2.2 设 $u\in W_{0}^{s,p}(\Omega )\setminus\left \{ 0 \right \}$, $\lambda \in (0,+\infty )$, 则

(I) $\lim\limits_{\lambda \to 0^+}J (\lambda u)=0$, $\lim\limits_{\lambda \to +\infty } J(\lambda u)=-\infty$;

(II) 存在唯一的 $\lambda^{\ast} \in (0,+\infty )$, 使得 $\frac{\mathrm{d} }{\mathrm{d} \lambda } J(\lambda u)\mid _{\lambda =\lambda ^{\ast } }=0\text{, } \nonumber$ 并且 $J(\lambda u)$ 在 $(0,\lambda ^{\ast } )$ 上递增, 在 $(\lambda ^{\ast },+\infty )$ 上递减, 在 $\lambda=\lambda ^{\ast }$ 处取得最大值;

(III) 对于任意的 $0< \lambda < \lambda ^{\ast }$ 有 $I(\lambda u)> 0$, 对于任意的 $\lambda > \lambda ^{\ast }$ 有 $I(\lambda u)< 0$, 并且 $I(\lambda ^{\ast }u)=0$.

证 对于任意的 $u\in W_{0}^{s,p} (\Omega )$, 有

其中 $\lambda>0$, 则 (I) 成立.

通过简单的计算可以得到

根据 (2.7) 式, 当 $0<\lambda<\lambda^{\ast }$ 时有 $\frac{\mathrm{d} }{\mathrm{d} \lambda } J(\lambda u(x))>0\text{, } \nonumber$ 当 $\lambda ^{\ast } < \lambda <+\infty$ 时有 $\frac{\mathrm{d} }{\mathrm{d} \lambda } J(\lambda u(x))< 0\text{, } \nonumber$

因此 (II) 成立. 通过 (2.7) 式和

得到

根据已知得到 (III) 成立.

引理2.3 设 $u\in W_{0}^{s,p}(\Omega )\setminus\left \{ 0 \right \}$, $\varrho \in (0,p_{s}^{\ast }-\theta p )$, 得到

(I) 若 $0< \left [ u \right ]_{s,p}<r(\varrho )$, 则 $I(u)>0$;

(II) 若 $I(u)\le 0$, 则 $\left [ u \right ] _{s,p}\ge r(\varrho )$.

证 根据引理 2.1 和公式 (2.2) 可以得到

从而得到

若 $0< [u]_{s,p}< r(\varrho )$, 则

因此 (I) 成立. 同理可证 (II) 成立.

引理2.4 假设 $s\in (0,1)$, $2< p< N/s$, $\theta \in [1,p_{s}^{\ast }/p)$, $t\in [0,T)$, 并且 $u$ 为方程 (1.5) 的弱解, 则

证 通过定义 1.1 可得到

通过简单的计算得到

因此有

引理2.5 假设 $u_0\in W_{0}^{s,p}(\Omega )$, $u(t)$ 为方程 (1.5) 的弱解, 如果 $J(u_0)\le h$, $I(u_0)<0$, 则对于任意的 $0\le t< T$ 都有 $u(t)\in V$.

证 考虑当 $J(u_0)<h$ 时, 若假设结论不成立, 则根据 $I(u_0)<0$ 以及连续性可知存在 $t_0\in (0,T)$ 使得对于任意 $t\in [0,t_0)$ 有 $I(u(t_0))=0$ 且 $I(u(t))<0$. 根据引理 2.3, 可知对于任意的 $t\in [0,t_0)$ 有 $\left [ u \right ] _{s,p}\ge r(\varrho )$. 根据连续性可得到 $\left [ u(t_{0}) \right ] _{s,p}\ge r(\varrho )$, 从而有 $u(t_0)\ne 0$, 因此 $u(t_0)\in N$, 根据 $h$ 的定义可知 $J(u_0)\ge h$, 这与 $J(u_0)\le h$ 矛盾, 从而有对于任意的 $0\le t< T$ 有 $I(u(t))<0$, 即对于任意的 $0\le t< T$ 都有 $u(t)\in V$. 同理可证 $J(u_0)=h$ 的情况.

在本节中, 我们将证明方程 (1.5) 的解在次临界状态 $J(u_{0})<h$ 下的全局存在性以及在有限时间内的爆破性.

定理3.1 设 $u_0\in W_{0}^{s,p} (\Omega )$ 且 $J(u_0)<h$, $I(u_0)>0$. 则方程 (1.5) 存在一个弱解 $u$ 在 $\Omega \times [0,T)$ 上满足 $u(t)\in L^{\infty } (0,T;W_{0}^{s,p}(\Omega ) )$, $u_t\in L^{2}(0,T;H^{s}_{0}(\Omega ))$ 且对于任意的 $t\in [0,+\infty )$, 有 $u(t)\in W$. 满足能量不等式

证 通过定义 1.1, 选取 $\left \{ {\eta _k(x)} \right \} _{k=1}^{\infty }$ 为 $W_{0}^{s,p}(\Omega )$ 中的正交基. 构造方程 (1.5) 的近似解,

满足下面的常微分方程

设

并且

于 $W_{0}^{s,p}(\Omega )$.

将方程 (3.2) 两边同时乘以 $d_{km}^{'}(t)$, 对 $k$ 从 1 到 $m$ 求和, 再关于 $t$ 积分, 可以得到

根据 (3.3) 式, 可以得到 $d_{km}(0)=b_{km}$, $J(u_m(x,0))\to J(u_0)< h$. 根据已知 $J(u_0)<h$, 得到对于足够大的 $m$ 有

同理可得 $I(u_m(x,0))\to I(u_0)>0$. 由 (2.4) 式可知对于足够大的 $m$ 有 $u_m(x,0)\in W$.

以下将证明 $u_m(t)\in W$. 如果 $m\to +\infty$, $0\le t\le T$, 假设 $u_m(t)\notin W$, 那么存在 $t_0\in (0,T)$ 使得 $u_m(t_0)\in \partial W$, 则有 $I(u_m(t_0))=0$, $u(t_0)\ne 0$, 因此, $u_m(t_0)\in N$. 通过 (2.6) 式得到 $J(u_m(t_0))>h$, 显然是矛盾的. 对于任意的 $0\le t\le T$, 存在足够大的 $m$ 使得 $u_m(t_0)\in W$. 通过 (3.4) 式和

得到

当 $m$ 足够大时有

其中 $C$ 为 $L^{\theta p}(\Omega )\hookrightarrow L^{q}(\Omega )$ 的最佳嵌入常数. 因此得到当 $m$ 足够大时, $T=+\infty$, 对于任意的 $0\le t< +\infty$ 有 $u_m(t)\in W$.

接下来, 将证明 $u(t)\in L^\infty (0,T;W_{0}^{s,p}(\Omega ))$, $u_t\in L^2(0,T;H^{s}_{0}(\Omega ))$. 根据 (3.5), (3.6), (3.7) 式可以得到 $u(t)\in L^\infty (0,T;W_{0}^{s,p}(\Omega ))$ 以及 $u_t\in L^2(0,T;H^{s}_{0}(\Omega ))$. 基于上述证明, 可以得到存在 $u$ 和 $\left \{ u_m \right \}$ 使得

通过 (3.8), (3.9), (3.12) 式以及 Aubin-Lions 紧性定理可得到

则有

通过 (3.13) 式, 可以得到

通过同样的方法可以得到

因为 $\left \{ {\eta _k(x)} \right \} _{k=1}^{\infty }$ 为 $W_{0}^{s,p}(\Omega )\subset H^{s}_{0}(\Omega )$ 中的正交基, 故 $\left \{ {\Delta \eta _k(x)} \right \} _{k=1}^{\infty }\in H^{-s}(\Omega )$. 将方程 (3.2) 乘以任意函数 $f(t)$ 并在 $\left [ 0,T \right ]$ 上积分得到

根据 (3.12) 式有

根据函数 $f(t)$ 的任意性以及 (3.2) 式, 固定一个 $k$, 令 $m\to +\infty$, 则有

因为 $\left \{ {\eta _k(x)} \right \} _{k=1}^{\infty }$ 在 $W_{0}^{s,p}(\Omega )$ 中是稠密的, 则可得到

其中对于 $t\in (0,+\infty )$ 有 $v\in W_{0}^{s,p}(\Omega )$, 并且在 $W_{0}^{s,p}(\Omega )$ 中有 $u(x,0)=u_0(x)$.

最后证明能量不等式 (3.1) 成立. 根据引理 2.1, 可以得到

通过弱下半连续性和上述不等式, 可以得到

因此得到能量不等式

定理3.2 设 $u_0\in W_{0}^{s,p} (\Omega )$, $J(u_0)<h$, $I(u_0)<0$, $u=u(t)$ 是方程 (1.5) 的一个弱解, 则 $u(t)$ 在有限时间 $T$ $(T>0)$ 内爆破, 即

证 假设 $u_0\in W_{0}^{s,p} (\Omega )$, $u=u(t)$ 是方程 (1.5) 的一个弱解. 因为 $J(u_0)<h$, $I(u_0)<0$, 结合引理 2.5 得到 $u(t)\in V$, $0<t<T^{\ast }$, 因此得到 $J(u)<h$, $I(u)<0$, $t\in[0,T^{\ast } )$. $T^ {\ast}=\infty$, 定义辅助函数

通过简单的计算得到

根据 (3.1) 和 (2.3) 式可得

根据 (3.16) 式和 $I(u)<0$ 可得到 $G^{"}(t)>0$, 因此有

根据 Hölder 不等式得到

令 $x=a^{\frac{1}{2} }\left \| \nabla u_{\tau } \right \|_{2}\text{, }y=a^{\frac{1}{2} }\left \| \nabla u \right \|_{2}\text{, }z=\left \| u_{\tau } \right \|_{2}\text{, }w=\left \| u \right \|_{2}$, 根据不等式

可得

根据 (3.15) 和 (3.19) 式有

结合 (3.14), (3.17) 和 (3.20) 式有

因此, 得到

令 $t_{0}>0$, 根据 (3.18) 式有

(a) 如果 $J(u_{0})\le 0$, 则根据 (3.17) 式可得

由 Hölder 不等式得到

其中

根据微分不等式原理, 得到

由于 $G>0$, 根据 (3.21) 式, 当 $t$ 足够大时有

(b) 如果 $0< J(u_{0})< h$, 根据 $I(u)<0$, $t\in (0,+\infty )$ 和引理 2.5 可知存在 $\rho ^{\ast }\in(0,1)$ 使得 $I(\rho ^{\ast }u)=0$. 因此通过 (2.3) 和 (2.6) 式可得

通过计算可得到

设 $T^{\ast }>t_{0}$, 定义函数

根据已知可得到 $H(t)\ge H(0)>0\text{, }H^{'}(t)=G^{'}(t)-G^{'}(0)$, $H^{"}(t)=G^{"}(t)>0$. 通过 (3.16) 式可得到

则有

因此得到 $y(t)$ 是 $t\in[t_{0},T^{\ast }]$ 上的一个凹函数. 由于 $y(t_{0})>0$ 并且 $y^{'}(t_{0})<0$, 则存在 $T<T^{\ast }$ 使得 $\lim\limits_{t \to T^{-}} y(t)=0$, 即

在本节中, 我们将证明方程 (1.5) 的解在临界状态 $J(u_{0})=h$ 下的全局存在性以及在有限时间内的爆破性.

定理4.1 设 $u_0\in W_{0}^{s,p} (\Omega )$ 且 $J(u_0)=h$, $I(u_0)\ge0$. 则方程 (1.5) 存在一个弱解 $u$ 在 $\Omega \times [0,T)$ 上满足 $u(t)\in L^{\infty }(0,T;W_{0}^{s,p}(\Omega ) )$, $u_t\in L^{2}(0,T;H^{s}_{0}(\Omega ))$ 且对于任意的 $t\in [0,+\infty )$ 有 $u(t)\in W$.

证 根据 $J(u_{0})=h$, $[u_{0}]_{s,p}\ne 0$. 设 $\rho _{m}=1-\frac{1}{m}$, $m=2\text{,}\ 3\text{,}\cdots\text{, }$ 并考虑下面的近似问题

根据 $I(u_{0})\ge0$, $\rho _{m}\in (0,1)$ 和引理 2.2 可得

和

通过 (2.2) 式, 可得到

因此得到 $\rho _{m}\frac{\mathrm{d} }{\mathrm{d}\rho _{m}}J(\rho _{m}u_{0}) =I(\rho _{m}u_{0})> 0$. 通过引理 2.2 (III) 可知

根据定理 3.1 可知存在 $m$, 使得方程 (4.1) 存在解 $u_{m}(t)$, 其中 $u_{m}(t)\in L^{\infty }(0,+\infty ;W_{0}^{s,p}(\Omega ))$, $u_{mt}\in L^{2}(0,+\infty ;H^{s}_{0}(\Omega ))$, $u_{m}\in W$, $t\in[0,+\infty )$ 并且满足

其中 $\forall \varphi\in W^{s,p}_{0}(\Omega )$, $t>0$, 另外

根据 (2.3) 式可得

通过 (4.2) 和 (4.3) 式, 可以得到

再根据定理 3.1 的证明, 定理 4.1 得证.

定理4.2 设 $u_0\in W_{0}^{s,p} (\Omega )$, $J(u_0)=h$, $I(u_0)<0$, $u=u(t)$ 是方程 (1.5) 的一个弱解, 则 $u(t)$ 在有限时间 $T$ $(T>0)$ 内爆破.

证 由于 $J(u_{0})=h$, $I(u_{0})<0$, 根据引理 2.5 得到对于任意的 $0\le t<T^{\ast }$ 有 $u(t)\in V$, 从而有

因此对于 $\forall t\ge 0$ 有 $(\left \| u_{t} \right \|^{2}_{2}+\alpha\left \| \nabla u_{t} \right \|^{2}_{2})>0$. 通过能量不等式可知存在 $t_{1}>0$, 使得

设 $t_{1}$ 为初始时间, 根据定理 3.2 的证明, 定理 4.2 得证.