数学物理学报, 2026, 46(2): 535-551

趋化系统解的全局存在性与有界性——献给陈化教授 70 寿辰

郭雅平,, 李佳琳,, 吕文斌,*

山西大学数学与统计学院 太原 030006

Global Existence for a Class of Chemotaxis Systems with Signal-Dependent Motility and Generalized Logistic Source

Guo Yaping,, Li Jialin,, Lyu Wenbin,*

School of Mathematics and Statistics, Shanxi University, Taiyuan 030006

通讯作者: *吕文斌, Email:lvwenbin@sxu.edu.cn

收稿日期: 2025-12-11   修回日期: 2026-01-5  

基金资助: 山西省基础研究计划(202503021211056)
山西省科技创新人才团队(202204051002015)

Received: 2025-12-11   Revised: 2026-01-5  

Fund supported: NSF of Shanxi Province(202503021211056)
special fund for Science and Technology Innovation Teams of Shanxi Province(202204051002015)

作者简介 About authors

郭雅平,Email:gyp2016@sxu.edu.cn

李佳琳,Email:245265420@qq.com

摘要

该文主要考虑了在光滑有界区域 $\Omega\subset\mathbb{R}^n\,(n\geqslant1)$ 上具有齐次 Neumann 边界条件的一类运动依赖于信号的趋化模型$$\begin{equation*} \begin{cases} u_t=\Delta(\gamma (v)u)+\rho u-\mu u^\alpha,&x\in\Omega,\,t>0,\\ v_t=\Delta v-v+u^\beta,&x\in\Omega,\,t>0 \end{cases} \end{equation*}$$的解的全局存在性与有界性. 当 $\rho$, $\mu>0$, $\alpha> 1$, $\beta>0$ 满足适当的条件, 且运动函数 $\gamma\in C^3((0,+\infty))$ 满足一定条件时, 该文证明了上述方程对所有足够光滑的初始值都存在一个全局古典解. 这改进了 [Lv W B, Wang Q Y. Proc Roy Soc Edinburgh, 2021, 151(2): 821-841], [Tao X Y, Fang Z B. Z Angew Math Phys, 2022, 73(3): Art 123] 中得到的结果.

关键词: 全局存在性; 趋化性; 有界性; 广义 Logistic 项

Abstract

This paper is concerned with the global existence for a class of Keller-Segel model $$\begin{equation*} \begin{cases} u_t=\Delta(\gamma (v)u)+\rho u-\mu u^\alpha,&x\in\Omega,\,t>0,\\ v_t=\Delta v-v+u^\beta,&x\in\Omega,\,t>0, \end{cases} \end{equation*}$$ under homogeneous Neumann boundary conditions in a smoothly bounded domain $\Omega\subset\mathbb{R}^n\,(n\geqslant1)$. It is proved that for $\rho\in\mathbb{R},\,\mu>0$, $\alpha> 1$, $\beta>0$ satisfying certain additional relations, and under suitable assumptions on the motility function $\gamma$, the system admits a global classical solution for all sufficiently smooth initial data. This result improves recent ones established in [Lv W B, Wang Q Y. Proc Roy Soc Edinburgh, 2021, 151(2): 821-841], [Tao X Y, Fang Z B. Z Angew Math Phys, 2022, 73(3): Art 123].

Keywords: global existence; chemotaxis; boundedness; general logistic source

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

本文引用格式

郭雅平, 李佳琳, 吕文斌. 趋化系统解的全局存在性与有界性——献给陈化教授 70 寿辰[J]. 数学物理学报, 2026, 46(2): 535-551

Guo Yaping, Li Jialin, Lyu Wenbin. Global Existence for a Class of Chemotaxis Systems with Signal-Dependent Motility and Generalized Logistic Source[J]. Acta Mathematica Scientia, 2026, 46(2): 535-551

1 引言和主要结论

为了描述在实验中观察到的细菌运动中形成的条纹斑图, 文献 [15] 提出并分析了如下的密度抑制型模型

$\begin{equation}\label{1.1} \begin{cases} u_t=\Delta(\gamma(v)u)+f(u),&x\in\Omega, t>0,\\ v_t=\Delta v-v+g(u),& x\in \Omega, t>0, \end{cases} \end{equation}$

其中 $u$ 和 $v$ 分别表示大肠杆菌细胞密度和酰基高丝氨酸内酯 (AHL) 的浓度;

函数 $f(u)$ 描述了细菌的生长和死亡, $g(u)$ 表示 AHL 的生成量;$\gamma'(v)<0$ 刻画了化学信号浓度对细胞运动的抑制作用.

事实上, 系统 (1.1) 是一种特殊的 Keller-Segel 模型

$\begin{equation*} \begin{cases} u_t=\nabla\cdot(D(u,v)\nabla u)-\nabla\cdot(\chi(u,v)\nabla v)+f(u),&x\in\Omega, t>0,\\ v_t=\Delta v-v+u,& x\in \Omega, t>0. \end{cases} \end{equation*}$

此模型是由 Keller 和 Segel 在 1970 年建立并做出了开创性工作[11]. 因其在数学、生物学及药理学领域的重要意义, 模型随后被众多学者从多方面展开了深入研究. 众所周知, 对于此模型来说, 一个有趣且具有挑战性的问题是研究解的全局存在性与爆破.

目前已有研究证实,在二维及更高维情形下解会产生爆破[37,39,40].然而,一个合适的 Logistic 项可有效抑制 Keller-Segel 模型解的爆破[3,12,13,23,28,30,32,34,36,42,44,46,47,49].

下文将对系统 (1.1) 的现有研究成果作简要综述. 除非另有说明, 这里总是假定系统满足齐次 Neumann 边界条件.

若 $f(u)\equiv0$ 且 $g(u)=u$. 文献 [27] 首次给出了该系统全局解的存在性结果: 当运动函数 $\gamma(v)$ 存在正的上下界时, 在二维空间中古典解全局有界, 三维及更高维空间中至少存在一个全局弱解. 然而 $\gamma(v)$ 未必存在先验的正下界/上界. 例如, 当 $\gamma(v)=\frac{c_0}{v^k}$ 时, 文献 [48] 证明了对任意 $k>0$ 只要 $c_0>0$ 很小, 该系统在任意维数的空间上都存在全局有界解. 若 $\gamma(v)=\frac{1}{c+v^k}$ (其中 $c\geqslant 0$), 当 $n=1$ 且 $k>0$, $n=2$ 且 $k\in(0,2)$ 或 $n=3$ 且 $k\in(0,\frac{3}{4})$ 时, 文献 [2] 证明了具有大初值的系统 (1.1) 的弱解的全局存在性. 文献 [33] 证明了在 $\gamma(v)$ 满足某些条件的基础上, 当 $2\leqslant n\leqslant4$ 时, 系统存在全局有界的古典解. 当 $\gamma(v)={\rm e}^{-\chi v}$ 时, 文献 [10] 发现了临界质量现象: 若 $n=2$, 则存在一个临界值 $m=4\pi/\chi>0$, 使得当初始细胞质量 $\|u_0\|_{L^1(\Omega)}>m$ 时, 系统 (1.1) 的解可能会发生爆破; 而当 $\|u_0\|_{L^1(\Omega)}<m$ 时, 则系统 (1.1) 的古典解全局存在且关于时间一致有界. 该结果在[4,5] 中得到了进一步的完善, 表明爆破发生在无穷远处.

若 $f(u)=\rho u-\mu u^\alpha$ 且 $g(u)=u$. 文献 [9] 证明了: 当运动函数 $\gamma(v)$ 在 $[0,\infty)$ 上满足 $\gamma\in C^3([0,+\infty))$, $\gamma(v)>0$, $\gamma'(v)<0$, $\lim\limits_{v \to +\infty}\gamma(v)=0$ 且 $\lim\limits_{v \to +\infty}\frac{\gamma'(v)}{\gamma(v)}$ 存在, 以及 $\mu>\frac{K_0}{16}$ (其中 $K_0=\max\limits_{0\leqslant v \leqslant +\infty}\frac{|\gamma'(v)|^2}{\gamma(v)}$) 时, 系统 (1.1) 在二维空间中存在唯一的全局古典解, 而且该解全局渐近收敛于平衡态 $(1,1)$. 随后, 文献 [14] 和 [31] 在 $\mu>0$ 足够大的条件下, 将类似结果推广到更高维 ($n\geqslant 3$) 的空间中. 对于广义 Logistic 源 (即, 把 $\mu u(1-u)$ 替换为 $\mu(u-u^\theta)$), 若 $\theta>\max\{2,\frac{n+2}{2}\}$ 且 $\gamma(v)$ 在 $[0,\infty)$ 上满足 $\gamma\in C^3([0,+\infty))$ 和 $\gamma(v)>0$, 则系统 (1.1) 在任意维数的空间中均存在全局解[19,20].

若 $f(u)=\rho u-\mu u^\alpha$ 且 $g(u)=u^{\beta}$. 文献 [24] 证明了: 如果运动函数满足 $\gamma\in C^3([0,+\infty))$ 是正函数, 且在 $[0,\infty)$ 上 $\gamma'(v)\leqslant0$, 同时参数满足 $\rho$, $\mu>0$, $\alpha>1$ 及 $0<\beta<\frac{2\alpha}{n+2}$, 则系统存在全局有界解. 此外, 该文献还分析了时间趋于无穷大时解的渐近行为.

此外, 相关研究可参考文献[16-18,21,22,41,43,45] 及其参考文献.

本文考虑以下具有广义 Logistic 项的信号依赖于运动的模型

$\begin{equation}\label{1.2} \begin{cases} u_t=\Delta(\gamma(v)u)+\rho u-\mu u^\alpha, &x\in\Omega,\,t>0,\\ v_t=\Delta v-v+u^\beta, &x\in\Omega,\,t>0,\\ \displaystyle\frac{\partial u}{\partial \nu}=\frac{\partial v}{\partial \nu}=0, &x\in\partial\Omega,t>0,\\ u(x,0)=u_0(x),\,v(x,0)=v_0(x), &x\in\Omega. \end{cases} \end{equation}$

对于系统 (1.2), 假设其满足以下条件

• $\Omega\subset \mathbb{R}^n\,(n\geqslant1)$ 是一个具有光滑边界 $\partial\Omega$ 的有界区域, $\nu$ 是 $\partial\Omega$ 的单位外法向量.

• $u$ 和 $v$ 分别表示细菌密度与趋化因子浓度. 初始条件满足

$\begin{equation}\label{1.3} \begin{gathered} u_0\in C^0(\overline{\Omega}),\,v_0\in W^{1,\infty}(\Omega),\\ u_0\geqslant0\,\text{且}\,v_0\geqslant C_0>0, \end{gathered} \end{equation}$

其中 $C_0$ 为正常数.

• $\rho>0$, $\mu>0$ 且 $\alpha>1$, $\beta>0$ 均为常数.

• 假设运动函数 $\gamma\in C^3((0,+\infty))$ 在 $(0,+\infty)$ 上恒正且满足

$\begin{equation}\label{1.4} \gamma(v)>\gamma_1,|\gamma'(v)|\leqslant \gamma_2,\gamma"(v)\leqslant 0. \end{equation}$

我们可以列举出很多满足条件 (1.4) 的函数. 例如

${\gamma(v)=a+\frac{v}{\sqrt{1+v^2}},\,v>0,}$

其中 $a>0$.

基于上述模型设定, 接下来陈述本文的主要结果.

定理 1.1 设 $\Omega\subset \mathbb{R}^n$ 为具有光滑边界的有界区域, $\rho>0$, $\mu>0$ {且正函数 $\gamma\in C^3((0,+\infty))$.} 若满足下列条件之一

$\begin{align*} &{(\mathbf{H1}) \alpha > 1,\ 0<\beta<\frac{2}{n+2}\alpha,} \\ &(\mathbf{H2}) 1<\alpha<\frac{n+2}{n},\ \frac{2}{n+2}\alpha\leqslant\beta<\frac{2}{n}\text{ 且满足条件 (1.4)}, \\ &(\mathbf{H3}) \alpha=\frac{n+2}{n},\ \beta=\frac{2}{n}\text{ 且满足条件 (1.4)}, \\ &(\mathbf{H4}) \alpha>\frac{n+2}{n},\ \frac{2}{n+2}\alpha\leqslant\beta<\alpha-1,\ \text{ 且满足条件 (1.4)}, \\ &(\mathbf{H5}) \alpha>\frac{n+2}{n},\ \beta=\alpha-1,\ \mu>\frac{p_0-1}{p_0}\gamma_2M_{\frac{p_0+\alpha-1}{\alpha-1}}^{\frac{\alpha-1}{p_0+\alpha-1}}\text{ 且满足条件 (1.4)}, \end{align*}$

其中, $p_0>\frac{n\beta}{2}$ 为固定常数, 且 $M_{\frac{p_0+\alpha-1}{\alpha-1}}$ 为由引理 2.2 定义的常数. 则对于任意满足条件 (1.3) 的初值, 系统 (1.2) 存在唯一的非负古典解 $(u,v)$ 满足

$\hspace{-2.3cm} u\in C^0(\overline{\Omega}\times[0,+\infty))\cap C^{2,1}(\overline{\Omega}\times(0,+\infty)),$
$v\in\bigcap_{p>\max\{2,n\}}C^0([0,+\infty);W^{1,p}(\Omega))\cap C^{2,1}(\overline{\Omega}\times(0,+\infty)). $

此外, 若 $\gamma(0)>0$, 则该解关于时间一致有界, 即存在常数 $C > 0$, 使得

$\sup_{t\in[0,+\infty)}\left(\|u(\cdot,t)\|_{L^{\infty}(\Omega)} + \|v(\cdot,t)\|_{W^{1,\infty}(\Omega)} \right)\leqslant C.$

注 1.1 当 $n=2$ 时, 上述五个条件在坐标系中的范围如下

注 1.2 在定理 1.1 中对运动函数 $\gamma$ 的假设条件比文献 [20] 中的假设更弱. 文献[20] 考虑了 $\beta=1$ 的特殊情形. 实际上, 即使无文献 [20] 中的如下条件

$\text{在}\,(0,+\infty) \text{上}\,\gamma'<0\, \text{且}\, \text{在}\, [\varepsilon,+\infty)\, \text{上} \frac{\gamma'}{\gamma}\,\text{有界}, $

其中 $\varepsilon>0$, 我们依然可以得到此类系统的全局解存在.

注 1.3 定理 1.1 的结果表明, 在特定参数范围内, 系统的全局古典解的存在性可得到严格证明. 值得注意的是, 这个参数范围具有随空间维数增大而逐渐受限、范围收窄的特征. 事实上, 我们的分析表明: 当条件 (H1)--(H5) 中有一个成立时, 系统存在全局经典解. 相较于文献 [24] 中的条件

$n\geqslant1, \,\ \rho,\; \mu>0,\ 0<\beta<\frac{2}{n+2}\alpha,$

即仅为 (H1) 的条件, 本文给出的条件更为宽松.这一发现不仅拓展了对系统动力学行为的理解,还为阐释空间维数、参数取值与全局古典解之间的关联提供了更为全面的视角.

下面概述本文方法的核心思路.根据 (1.2) 的结构信息以及扩散与交叉扩散项的特定关系,首先推导 $v$ 的正下界与上界, 这是本文的主要创新点.为实现这一目标, 我们首先通过比较方法得到 $v$ 的下界;同时, 借助经典的半群方法得出 $v$ 的上界. 本文不依赖于文献 [5] 中的条件

$\lim_{s\rightarrow+\infty}\gamma(s)<1.$

换言之, 无需上述条件, 仍可得到此类系统的全局解存在.然后, 在获得 $v$ 的正下界与上界后, 便能处理 $\Delta(\gamma(v)u)$ 项, 并建立 $u$ 的关于时间的一致有界性估计.

2 预备知识

在本节中, 我们给出线性热方程解的若干基本估计.

引理 2.1 设 $T\in(0,\infty]$, $q\geqslant 1$. 对每一个

$r_1 \in \begin{cases} \left[1,\frac{nq}{n - 2q}\right), & q \leqslant \frac{n}{2}, \\ [1,+\infty], & q > \frac{n}{2}, \end{cases}$
$r_2 \in \begin{cases} \left[1,\frac{nq}{n + 2 - 2q}\right), & q \leqslant \frac{n + 2}{2}, \\ [1,+\infty], & q > \frac{n + 2}{2}, \end{cases}$
$r_3 \in \begin{cases} \left[1,\frac{nq}{n - q}\right), & q \leqslant n, \\ [1,+\infty], & q > n \end{cases}$

$r_4 \in \begin{cases} \left[1,\frac{nq}{n + 2 - q}\right), & q \leqslant n+2, \\ [1,+\infty], & q > n+2, \end{cases}$

存在常数 $C_i = C (r_i) > 0 \, (i = 1, 2,3,4)$, 使得当 $z_0\in W^{1,\infty}(\Omega)$ 且

$z \in C^0(\overline{\Omega} \times [0,T)) \cap C^{2,1}(\overline{\Omega} \times (0,T)), \quad w \in C^0(\overline{\Omega} \times [0,T))$

满足如下初边值问题

$ \begin{cases} z_{t}=\Delta z-z+w, &x \in \Omega, t\in (0,T),\\ \displaystyle\frac{\partial z}{\partial \nu} = 0, & x\in\partial\Omega, t\in(0,T),\\ z(x,0)=z_{0}(x),&x \in \Omega.\\ \end{cases}$

时, 则对于所有 $ t\in(0,T)$ 有

$\begin{aligned}\label{2.1} & {\left\| { z(\cdot,t)} \right\|_{{L^{r_1}(\Omega)}}} \leqslant C_1\left({\left\| { z_0} \right\|_{{W^{1,\infty}}(\Omega)}}+\sup_{s\in(0,T)}\|w(\cdot,s)\|_{L^q(\Omega)}\right), \end{aligned}$
$\|z(\cdot, t)\|_{L^{r_{2}}(\Omega)} \leqslant C_{2}\left(\left\|z_{0}\right\|_{W^{1, \infty}(\Omega)}+\sup _{s \in(0, T)} \int_{0}^{t} \mathrm{e}^{-(\lambda+1)(t-s)}\|w(\cdot, s)\|_{L^{q}}^{q} \mathrm{~d} s\right)$
$\|\nabla z(\cdot, t)\|_{L^{r_{3}}(\Omega)} \leqslant C_{3}\left(\left\|z_{0}\right\|_{W^{1, \infty}(\Omega)}+\sup _{s \in(0, T)}\|w(\cdot, s)\|_{L^{q}(\Omega)}\right),$

$\begin{equation}\label{2.4} \begin{aligned} {\left\| {\nabla z(\cdot,t)} \right\|_{{L^{r_4}(\Omega)}}} \leqslant C_4\left({\left\| {z_0} \right\|_{{W^{1,\infty}}(\Omega)}}+ \sup_{s\in(0,T)}\int_0^{t}{\rm e}^{-(\lambda+1)(t-s)} \left\| {w\left(\cdot, s \right)} \right\|_{{L^q}}^q \mathrm{d}s\right). \end{aligned} \end{equation}$

不等式 (2.2) 的证明可参见文献 [21,引理 2.2],不等式 (2.3) 的证明可参见文献 [6,引理 4.1],不等式 (2.4) 的证明可参见文献 [20,引理 2.5].因此, 下面仅证明不等式 (2.1). 利用半群表示法, 对所有 $t\in (0,T)$, 得

$z = {\rm e}^{t(\Delta - 1)}z_0 + \int_{0}^{t}{\rm e}^{(t - s)(\Delta - 1)}w(s) \mathrm{d}s.$

利用带 Neumann 边界条件的 Laplace 半群的光滑性[35]和 Hölder 不等式得,存在常数 $\lambda,\, C_1 > 0$, 对于所有的 $t\in (0,T)$, 有

$\begin{align*} &\|z(\cdot,t)\|_{L^{r_1}(\Omega)}\\ \leqslant& C_1\|z_0\|_{L^{\infty}(\Omega)} + \int_{0}^{t}\|{\rm e}^{(t - s)\Delta}w(s)\|_{L^{r_1}(\Omega)}\mathrm{d}s \\ \leqslant& C_1\|z_0\|_{L^{\infty}(\Omega)} + C_1\int_{0}^{t}{\rm e}^{-\lambda(t - s)}\left(1 + (t - s)^{-\frac{n}{2}(\frac{1}{q} - \frac{1}{r_1})_+}\right)\|w(s)\|_{L^{q}(\Omega)}\mathrm{d}s \\ \leqslant& C_1\|z_0\|_{L^{\infty}(\Omega)} + C_1\sup_{s\in(0,T)}\|w(\cdot,s)\|_{L^{q}(\Omega)}\int_{0}^{t}{\rm e}^{-\lambda(t - s)}\left(1 + (t - s)^{-\frac{n}{2}(\frac{1}{q} - \frac{1}{r_1})_+}\right)\mathrm{d}s. \end{align*}$

当 $q\leqslant \frac{n}{2}$ 时, 有

$\begin{equation}\label{2.5} \begin{aligned} -\frac{n}{2}\left(\frac{1}{q} - \frac{1}{r_1}\right)_+ > -\frac{n}{2}\left(\frac{1}{q} - \frac{n - 2q}{nq}\right) = - 1. \end{aligned} \end{equation}$

当$q>\frac{n}{2}$ 时, 有

$\begin{equation}\label{2.6} \begin{aligned} -\frac{n}{2}\left(\frac{1}{q} - \frac{1}{r_1}\right)_+ \geqslant -\frac{n}{2q} > - 1. \end{aligned} \end{equation}$

结合 (2.5) 和 (2.6) 式得, 存在常数 $C_2 > 0$ 使得对于所有 $t \in (0,T)$ 有

$\int_{0}^{t}{\rm e}^{-\lambda(t - s)}\left(1 + (t - s)^{-\frac{n}{2}(\frac{1}{q} - \frac{1}{r_1})_+}\right) \mathrm{d}s \leqslant C_2.$

因此, 不等式 (2.1) 成立. 该引理得证.

关于线性热方程的极大 Sobolev 正则性, 我们有如下结论 (具体可参考文献 [7,引理 2.1] 与文献 [8]).

引理 2.2 设 $r \in (1, +\infty)$. 则存在常数 $M_r > 0$ 和 $C>0$, 使得对任意 $T \in (0, +\infty]$, $\tau_0 \in [0, T)$,$ z(\cdot,\tau_0 ) \in W^{2,r}(\Omega)$ 满足外法向导数条件$\displaystyle\frac{\partial z(\cdot,\tau_0)}{\partial \nu} = 0$ 且函数

$z \in C^{0}(\overline{\Omega} \times [0,T)) \cap C^{2,1}(\overline{\Omega} \times (0,T)),\quad w \in C^{0}(\overline{\Omega} \times [0,T))$

满足

$ \begin{cases} z_{t}=\Delta z-z+w, &x \in \Omega, t\in (0,T),\\ \displaystyle \frac{\partial z}{\partial \nu}=0, &x \in \partial\Omega, t\in (0,T),\\ z(x,0)=z_{0}(x),&x \in \Omega,\\ \end{cases}$

$\int_{\tau_{0}}^{T} \int_{\Omega} {\rm e}^{r\tau} (|z|^{r} +{ |\Delta z|^{r}}+|z_{t}|^{r}) \mathrm{d}\tau \leqslant M_{r} \int_{\tau_{0}}^{T} {\rm e}^{r\tau} \|w\|_{L^{r}(\Omega)}^{r} \mathrm{d}\tau +M_{r} \|z(\cdot,\tau_{0})\|_{W^{2,r}(\Omega)}^{r},$

$\int_{\tau_{0}}^{T} \int_{\Omega}(|z|^{r} + { |\Delta z|^{r}}+|z_{t}|^{r})\mathrm{d}\tau \leqslant C\int_{\tau_{0}}^{T}\|w\|_{L^{r}(\Omega)}^{r} \mathrm{d}\tau+M_{r} \|z(\cdot,\tau_{0})\|_{W^{2,r}(\Omega)}^{r},$

其中 $M_r$ 和 $C$ 都不依赖于 $T$.

下面给出一个 Grönwall 不等式[29], 主要用于获得解关于时间的一致估计.

引理 2.3 设 $T_{\max}>0$, $\tau\in(0,T_{\max})$. 假设 $a,b,y$ 是定义在 $(0,T_{\max})$ 上的三个正的局部可积函数且$y'$ 在 $(0, T_{\max})$ 上是局部可积的, 若对于所有 $t\in(0,T_{\max})$ 满足

$y'(t) \leqslant a(t)y(t)+b(t)$

且对所有 $t\in[0, T_{\max}-\tau)$,

$\int_{t}^{t+\tau} a\leqslant a_{1}, \ \ \int_{t}^{t+\tau} b \leqslant a_{2}, \ \ \int_{t}^{t+\tau} y \leqslant a_{3},$

其中 $a_{i}(i=1,2,3)$ 是正数. 那么对所有 $t\in[\tau, T_{\max})$,

$y(t) \leqslant\left(\frac{a_{3}}{\tau}+a_{2}\right) {\rm e}^{a_{1}}.$

3 全局存在性与一致有界性

本节, 我们将证明系统 (1.2) 解的全局存在性与有界性.

3.1 局部存在性及延拓准则

引理 3.1 设 $\Omega\subset \mathbb{R}^n$ 是一个具有光滑边界的有界区域,运动函数 $\gamma$ 满足条件 (1.5), 以及$\rho>0,\,\mu>0$ 和 $\alpha>1$, $\beta>0$ 是常数.如果初值满足条件 (1.3),则存在 $T_{max}\in(0,\infty]$ 和一对非负函数 $(u,v)$

$\hspace{-2.3cm} u\in C^0(\overline{\Omega}\times[0,T_{\max}))\cap C^{2,1}(\overline{\Omega}\times(0,T_{\max})),$
$v\in\bigcap_{p>\max\{2,n\}}C^0([0,T_{\max});W^{1,p}(\Omega))\cap C^{2,1}(\overline{\Omega}\times(0,T_{\max})),$

构成系统 (1.2) 在 $\Omega\times(0,T_{\max})$ 上的古典解.此外, 有

$ T_{\max}=\infty \quad \text{或} \quad \limsup\limits_{t\nearrow T_{\max}}(\|u(\cdot,t)\|_{L^\infty(\Omega)}+\|v(\cdot,t)\|_{W^{1,\infty}(\Omega)})=\infty.$

基于文献 [25] 中运用 Schauder 不动点定理的标准方法或者文献 [1]中的抽象理论都可以证明上述解的局部存在性, 这里省略这个引理的证明.

接下来, 我们给出系统 (1.2) 解的一些基本性质.

3.2 $u$ 的有界性

首先, 我们回顾解 $u$ 的有界性.

引理 3.2 如果 $\rho>0,\,\mu>0,\ \alpha>1$ 和$\beta>0$ 是常数, 则存在一个常数 $C>0$, 使得对于所有 $t \in (0,T_{\max})$, 有

$\begin{equation}\label{3.1} \int_\Omega u\leqslant C \end{equation}$

且对任意 $0<\tau<\min\{T_{\max},1\}$, $t \in (0,T_{\max}-\tau)$, 有

$\begin{equation}\label{3.2} \int_t^{t+\tau}\int_\Omega u^\alpha \mathrm{d}s\leqslant C. \end{equation}$

此外, 对于每个 $a>0$, 存在 $C>0$, 使得对任意 $t \in (0,T_{\max})$, 有

$\begin{equation}\label{3.3} \int_0^t{\rm e}^{-a(t-s)}\int_\Omega u^\alpha \mathrm{d}s\leqslant C. \end{equation}$

对 (1.2) 式中的第一个方程在 $\Omega$ 上积分, 任意 $t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.4} \frac{\mathrm{d}}{\mathrm{d}t}\int_{\Omega}u=\int_{\Omega}(\rho u-\mu u^\alpha). \end{equation}$

利用 Young 不等式可得, 存在常数 $C_1>0$, 使得对任意 $ t \in (0,T_{\max})$ 有

$\begin{equation*} \begin{aligned} \frac{\mathrm{d}}{{\mathrm{d}t}}\int_\Omega u + \int_\Omega u &= \left( {\rho + 1} \right)\int_\Omega u - \mu \int_\Omega {{u^\alpha }} \\ &\leqslant \mu \int_\Omega {{u^\alpha }} + {C_1} - \mu \int_\Omega {{u^\alpha }} = C_1. \end{aligned} \end{equation*}$

结合比较原理, 由此可推出所需结论 (3.1).

对 (3.4) 式在 $(t,t+\tau)$ 上积分, 得到 (3.2) 式.

在 (3.4) 式的两端分别乘以 ${\rm e}^{as}$ 并在 $(0, t)$ 上对 $s$ 积分得,

存在常数 $C_2>0$, 使得对于任意 $ t \in (0,T_{\max})$ 有

$\begin{equation*} \begin{aligned} \mu\int_{0}^{t}{\rm e}^{as}\int_{\Omega}u^{\alpha}\mathrm{d}s & = -\int_{0}^{t}{\rm e}^{as}\frac{\mathrm{d}}{\mathrm{d}s}\int_{\Omega}u \mathrm{d}s + \rho\int_{0}^{t}{\rm e}^{as}\int_{\Omega}u \mathrm{d}s\\ &\leqslant \int_{\Omega}u_0 + (a + \rho)\int_{0}^{t}{\rm e}^{as}\int_{\Omega}u \\ &\leqslant C_2{\rm e}^{at}, \end{aligned} \end{equation*}$

即 (3.3) 式成立. 该引理得证.

3.3 $v$ 的下界

接下来, 利用比较原理得 $v$ 的下界估计, 参见文献 [38].

引理 3.3 若 $T_{\max}<+\infty$, 则存在一个常数 $C>0$ 使得对所有 $t \in (0,T_{\max})$ 有

$\inf_{x\in\Omega}v(x,t)\geqslant C.$

通过比较原理并结合 $u$ 的正性, 由系统 (1.2) 的第二个方程得, 对于所有 $(x,t) \in \Omega \times (0, T_{\max})$ 有

$v(x,t) \geqslant {\rm e}^{-t} \inf_{y\in\Omega} v_0(y).$

这意味着对任意 $ (x,t) \in \Omega \times (0, T_{\max})$ 有

$v(x,t) \geqslant {\rm e}^{-T_{\max}} \inf_{y\in\Omega} v_0(y).$

该引理得证.

注 3.1 此引理仅针对解的全局存在性证明起作用, 对于解的一致有界性由于加了额外的条件, 这个引理完全不需要.

3.4 $v$ 的上界

由半群的估计, 立即可得 $v$ 的上界估计.

引理 3.4 如果 $\frac{\alpha}{\beta}>\frac{n+2}{2}$, 则存在一个常数 $C>0$, 使得对于所有 $ t \in (0,T_{\max})$

$\begin{equation}\label{3.5} \begin{aligned} \|v(\cdot,t)\|_{L^\infty}\leqslant C. \end{aligned} \end{equation}$

由引理 2.1 中的 (2.2) 式和引理 3.2 中的 (3.3) 式, 可得到 (3.5) 式.

3.5 改进的两个全局存在性判断准则

引理 3.5 假设 $\mu > 0$, $\rho > 0$. 若对于

$p_1>\frac{n\beta}{2},$

存在 $K > 0$, 使得对任意 $ t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.6} \begin{aligned} \|u(\cdot,t)\|_{L^{p_1}(\Omega)} \leqslant K, \end{aligned} \end{equation}$

或对于

$ p_2>\frac{(n+2)\beta}{2}, $

存在 $K > 0$ 使得对任意 $t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.7} \begin{aligned} \int_{0}^{t }{\rm e}^{-(\lambda+1)(t-s)}\|u(\cdot,s)\|_{L^{p_2}(\Omega)}^{p_2} \mathrm{d}s \leqslant K, \end{aligned} \end{equation}$

那么存在常数 $C > 0$, 使得对任意 $t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.8} \begin{aligned} \|u(\cdot,t)\|_{L^\infty(\Omega)}+\|v(\cdot,t)\|_{W^{1,\infty}(\Omega)} \leqslant C. \end{aligned} \end{equation}$

步骤 1 我们将证明: 存在 $p_0 > n$ 和 $C_1> 0$, 使得对任意 $t \in (0,T_{\max})$, 有

$\|\nabla v(\cdot,t)\|_{L^{p_0}(\Omega)} \leqslant C_1.$

(3.6) 式对应的情形. 由引理 2.1 中 (2.3) 式知, 对任意 $\sigma \in [1,(\frac{np_1}{n\beta-p_1})_+)$, 存在常数 $C_2 > 0$, 使得对任意 $t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.9} \begin{aligned} \|\nabla v(\cdot,t)\|_{L^{\sigma}(\Omega)} \leqslant C_2. \end{aligned} \end{equation}$

鉴于 $(\frac{np_1}{n\beta-p_1})_+ >n$, 由此可以从 (3.9) 式中得到所要证明的结论.

(3.7) 式对应的情形. 根据引理 2.1 中 (2.4) 式得, 对任意 $\kappa\in [1,(\frac{np_2}{(n+2)\beta-p_2})_+)$, 存在常数 $C_3 > 0$, 使得对任意 $t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.10} \begin{aligned} \|\nabla v(\cdot,t)\|_{L^{\kappa}(\Omega)} \leqslant C_{3}. \end{aligned} \end{equation}$

鉴于$(\frac{np_2}{(n+2)\beta-p_2})_+>n$, 由此可从 (3.10) 式中得到所要证明的结论.

步骤 2 我们断言, 对任意 $ p\geqslant2$, 存在常数 $C_4 > 0$, 使得对任意 $t \in (0,T_{\max})$ 有

$\int_\Omega u^p\leqslant C_{4}.$

实际上, 当 $\displaystyle\frac{\alpha}{\beta}>\frac{n+2}{2}$ 时, 结合引理 3.4 和引理 3.5; 当 $\displaystyle\frac{\alpha}{\beta}\leqslant\frac{n+2}{2}$ 时, 结合 (1.4) 式以及引理 3.3, 得存在常数 $\gamma_1, \gamma_2>0$, 使得对任意 $t \in (0,T_{\max})$ 有

$\gamma(v)>\gamma_1,\ |\gamma'(v)|\leqslant \gamma_2.$

利用 (1.2) 式的第一个方程、边界条件、分部积分公式和 Hölder 不等式, 可得对任意 $t \in (0,T_{\max})$

$\begin{equation}\label{3.11} \begin{aligned} \frac{\rm d}{{\rm d}t}\int_{\Omega}u^p &= p\int_{\Omega}u^{p - 1}u_t = p\int_{\Omega}u^{p - 1}\left(\Delta(\gamma(v)u) + \rho u - \mu u^{\alpha}\right) \\ &= -p(p - 1)\int_{\Omega}\gamma(v)u^{p - 2}|\nabla u|^2 - p(p - 1)\int_{\Omega}\gamma'(v)u^{p - 1}\nabla u\cdot\nabla v \\ &\quad + \rho p\int_{\Omega}u^p - \mu p\int_{\Omega}u^{p + \alpha - 1} \\ &\leqslant -\frac{p(p - 1)}{2}\int_{\Omega}\gamma(v)u^{p - 2}|\nabla u|^2 + \frac{p(p - 1)}{2}\int_{\Omega}\frac{|\gamma'(v)|^2}{\gamma(v)}u^{p}|\nabla v|^2 \\ &\quad + \rho p\int_{\Omega}u^p - \mu p\int_{\Omega}u^{p + \alpha - 1} \\ &\leqslant -\frac{p(p - 1)\gamma_1}{2}\int_{\Omega}u^{p - 2}|\nabla u|^2 + \frac{p(p - 1)\gamma_2^2}{2\gamma_1}\left(\int_{\Omega}|\nabla v|^{p_0}\right)^{\frac{2}{p_0}} \\ &\quad \times \left(\int_{\Omega} u^{p\frac{p_0}{p_0 - 2}}\right)^{\frac{p_0 - 2}{p_0}} + \rho p\int_{\Omega}u^p - \mu p\int_{\Omega}u^{p + \alpha - 1} \\ &\leqslant -\frac{p(p - 1)\gamma_1}{2}\int_{\Omega}u^{p - 2}|\nabla u|^2 + \frac{p(p - 1)\gamma_2^2C_1^2}{2\gamma_1}\left\|u^{\frac{p}{2}}\right\|_{L^{\frac{2p_0}{p_0 - 2}}}^2 \\ &\quad + \rho p\int_{\Omega}u^p - \mu p\int_{\Omega}u^{p + \alpha - 1}. \end{aligned} \end{equation}$

用 Gagliardo-Nirenberg 不等式和 Young 不等式得, 对任意 $t \in (0,T_{\max})$, 存在常数 $C_{5},C_{6}>0$ 满足

$\begin{equation}\label{3.12} \begin{aligned} \left\|u^{\frac{p}{2}}\right\|_{L^{\frac{2p_0}{p_0 - 2}}}^2 &\leqslant C_5\left(\left\|\nabla\left(u^{\frac{p}{2}}\right)\right\|_{L^2}^{\frac{2n}{p_0}}\left\|u^{\frac{p}{2}}\right\|_{L^2}^{\frac{2(p_0 - n)}{p_0}} + \left\|u^{\frac{p}{2}}\right\|_{L^{\frac{2}{p}}}^2\right) \\ &\leqslant\frac{\gamma_1^2}{2\gamma^2_2C_1^2}\int_{\Omega}u^{p - 2}|\nabla u|^2 + C_6\int_{\Omega}u^p + C_6. \end{aligned} \end{equation}$

结合 (3.11) 和 (3.12) 式知, 存在常数 $C_{7}>0$, 使得对任意 $t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.13} \begin{aligned} \frac{\rm d}{{\rm d}t}\int_{\Omega}u^p + p\int_{\Omega}u^p \leqslant &\left(p+\rho p+C_7\right)\int_{\Omega}u^p- \mu p\int_{\Omega}u^{p + \alpha - 1}+ C_{7}. \end{aligned} \end{equation}$

然后, 利用 Young 不等式得上述断言.

步骤 3 应用文献 [26] 中的 Moser 迭代, 可以推断出存在常数 $C_{8}>0$, 使得对任意 $t \in (0,T_{\max})$ 有

$\begin{equation*} \|u(\cdot,t)\|_{L^{\infty}(\Omega)}\leqslant C_{8}, \end{equation*}$

相关细节可参见文献 [26]. 将引理 2.1 应用于系统 (1.2), 不难推出存在常数 $C_{9} > 0$, 满足对任意 $t \in (0,T_{\max})$ 有

$\|u(\cdot,t)\|_{L^\infty(\Omega)}+\|v(\cdot,t)\|_{W^{1,\infty}(\Omega)} \leqslant C_{9}.$

至此, 该引理得证.

3.6 $u$ 的 $L^p$ 估计

首先, 给出一个关于 $u$ 的 $L^p$ 估计的不等式, 该不等式对后续引理的证明至关重要.

引理3.6 设 $p>1$, $\rho>0,\,\mu>0,\ \alpha>1$, $\beta>0$ 且运动函数满足条件 (1.4), 那么对任意的 $\varepsilon>0$, 存在常数 $C>0$, 使得对任意 $\tau_0\in (0,T_{\max})$ 且 $t\in (\tau_0,T_{\max})$ 有

$\begin{equation*} \begin{aligned} \int_{\Omega}u^{p} &\leqslant {\rm e}^{-\frac{p+\alpha-1}{\alpha-1}(t - \tau_0)}\int_{\Omega}u^{p}(\tau_0) + \int_{\tau_0}^{t}{\rm e}^{-\frac{p+\alpha-1}{\alpha-1}(t - \tau)} \Bigg[ -\mu p\int_{\Omega}u^{p + \alpha - 1} \\ &\quad + \frac{p(p-1)\gamma_2}{p+\alpha-1}\varepsilon\int_{\Omega}u^{p+\alpha-1} + \frac{(\alpha-1)(p-1)\gamma_2}{(p+\alpha-1)\varepsilon^{\frac{p}{\alpha-1}}} M_{\frac{p+\alpha-1}{\alpha-1}}\int_{\Omega}u^{\frac{\beta(p+\alpha-1)}{\alpha-1}}(\tau) \\ &\quad + \left(\rho p + \frac{p+\alpha-1}{\alpha-1}\right)\int_{\Omega}u^{p} \Bigg] \mathrm{d}\tau \\ &\quad + \frac{(\alpha-1)(p-1)\gamma_2}{(p+\alpha-1)\varepsilon^{\frac{p}{\alpha-1}}} M_{\frac{p+\alpha-1}{\alpha-1}} \|v(\cdot,\tau_0)\|_{W^{2,\frac{p+\alpha-1}{\alpha-1}}(\Omega)}^{\frac{p+\alpha-1}{\alpha-1}}, \end{aligned} \end{equation*}$

其中 $M_{\frac{p+\alpha-1}{\alpha-1}}$ 是引理 2.2 中的常数.

由 (1.4) 式得, 对任意 $t \in (0,T_{\max})$ 有

$|\gamma'(v)| \leqslant \gamma_2.$

利用 (1.2) 式的第一个方程、边界条件和分部积分得, 对任意 $t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.14} \begin{aligned} \frac{\rm d}{{\rm d}t}\int_{\Omega}u^{p} =& p\int_{\Omega}u^{p - 1}u_t = p\int_{\Omega}u^{p - 1}\left(\Delta(\gamma(v)u) + \rho u - \mu u^{\alpha}\right) \\ =& -p(p - 1)\int_{\Omega}\gamma(v)u^{p - 2}|\nabla u|^{2} - p(p - 1)\int_{\Omega}\gamma'(v)u^{p - 1}\nabla u\cdot\nabla v \\ &+ \rho p\int_{\Omega}u^{p} - \mu p\int_{\Omega}u^{p + \alpha - 1}. \end{aligned} \end{equation}$

接下来, 我们估计 (3.14) 式的右侧. 结合 (1.4) 式得, 对任意 $t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.15} \begin{aligned} &\int_{\Omega}\gamma'(v)u^{p - 1}\nabla u\cdot\nabla v= \frac{1}{p}\int_{\Omega}\gamma'(v)\nabla u^p\cdot\nabla v \\ =&-\frac{1}{p}\int_{\Omega}u^p\nabla\cdot(\gamma'(v)\nabla v) \\ =&-\frac{1}{p}\int_{\Omega}\gamma'(v)u^p\Delta v - \frac{1}{p}\int_{\Omega}\gamma"(v)u^p|\nabla v|^2 \\ \geqslant&-\frac{\gamma_2}{p}\int_{\Omega}u^p|\Delta v|. \end{aligned} \end{equation}$

利用 Young 不等式得, 对任意的 $\varepsilon>0$, 有对任意 $t \in (0,T_{\max})$

$\begin{align*} \int_{\Omega}u^p|\Delta v|\leqslant&\frac{p}{p+\alpha-1}\varepsilon\int_{\Omega}u^{p+\alpha-1} +\frac{\alpha-1}{(p+\alpha-1)\varepsilon^{\frac{p}{\alpha-1}}}\int_{\Omega}|\Delta v|^{\frac{p+\alpha-1}{\alpha-1}}. \end{align*}$

利用 (3.14) 和 (3.15) 式可得, 对任意 $t \in (0,T_{\max})$,

$\begin{equation}\label{3.16} \begin{aligned} \frac{\rm d}{{\rm d}t}\int_{\Omega}u^{p} &\leqslant (p-1)\gamma_2\int_{\Omega}u^{p}|\Delta v| + \rho p\int_{\Omega}u^{p} - \mu p\int_{\Omega}u^{p + \alpha - 1} \\ &\leqslant\frac{p(p-1)\gamma_2}{p+\alpha-1}\varepsilon\int_{\Omega}u^{p+\alpha-1} +\frac{(\alpha-1)(p-1)\gamma_2}{(p+\alpha-1)\varepsilon^{\frac{p}{\alpha-1}}}\int_{\Omega}|\Delta v|^{\frac{p+\alpha-1}{\alpha-1}}\\ &\quad+\rho p\int_{\Omega}u^{p} - \mu p\int_{\Omega}u^{p + \alpha - 1}. \end{aligned} \end{equation}$

在 (3.16) 式两边加上 $\frac{p+\alpha-1}{\alpha-1}\int_\Omega u^p$ 可得, 对任意 $t \in (0,T_{\max})$,

$\begin{equation*} \begin{aligned} \frac{\rm d}{{\rm d}t}\int_{\Omega}u^{p}+\frac{p+\alpha-1}{\alpha-1}\int_{\Omega}u^{p}\leqslant &\frac{p(p-1)\gamma_2}{p+\alpha-1}\varepsilon\int_{\Omega}u^{p+\alpha-1}+\frac{(\alpha-1)(p-1)\gamma_2}{(p+\alpha-1)\varepsilon^{\frac{p}{\alpha-1}}}\int_{\Omega}|\Delta v|^{\frac{p+\alpha-1}{\alpha-1}}\\ &+\left(\rho p + \frac{p+\alpha-1}{\alpha-1}\right)\int_{\Omega}u^{p}-\mu p\int_{\Omega}u^{p + \alpha - 1}. \end{aligned} \end{equation*}$

这意味着对于所有 $t \in (\tau_0,T_{\max})$ 有

$\begin{equation}\label{3.17} \begin{aligned} \int_{\Omega}u^{p} &\leqslant {\rm e}^{-\frac{p+\alpha-1}{\alpha-1}(t - \tau_0)}\int_{\Omega}u^{p}(\tau_0) + \int_{\tau_0}^{t}{\rm e}^{-\frac{p+\alpha-1}{\alpha-1}(t - \tau)} \Bigg\{ -\mu p\int_{\Omega}u^{p + \alpha - 1} \\ &\quad + \frac{p(p-1)\gamma_2}{p+\alpha-1}\varepsilon\int_{\Omega}u^{p+\alpha-1} + \frac{(\alpha-1)(p-1)\gamma_2}{(p+\alpha-1)\varepsilon^{\frac{p}{\alpha-1}}}\int_{\Omega}|\Delta v|^{\frac{p+\alpha-1}{\alpha-1}} \\ &\quad + \left(\rho p + \frac{p+\alpha-1}{\alpha-1}\right)\int_{\Omega}u^{p} \Bigg\} \mathrm{d}\tau. \end{aligned} \end{equation}$

利用引理 2.2 可得, 对任意 $t \in (\tau_0,T_{\max})$ 有

$\begin{equation*} \begin{aligned} &\int_{\tau_0}^{t}{\rm e}^{-\frac{p+\alpha-1}{\alpha-1}(t - \tau)} \int_{\Omega}|\Delta v(\tau)|^{\frac{p+\alpha-1}{\alpha-1}}\mathrm{d}\tau \\ \leqslant & M_{\frac{p+\alpha-1}{\alpha-1}}\|v(\cdot,\tau_0)\|_{W^{2,\frac{p+\alpha-1}{\alpha-1}}(\Omega)}^{\frac{p+\alpha-1}{\alpha-1}} + M_{\frac{p+\alpha-1}{\alpha-1}}\int_{\tau_0}^{t}{\rm e}^{-\frac{p+\alpha-1}{\alpha-1}(t - \tau_0)} \int_{\Omega}u^{\frac{\beta(p+\alpha-1)}{\alpha-1}}\mathrm{d}\tau. \end{aligned} \end{equation*}$

将上述结果带入 (3.17) 式, 即得所求结论.

当 $\beta<\alpha-1$ 时, 现在给出 $u$ 的 $L^p$ 估计.

引理 3.7 设 $\rho>0,\,\mu>0,\ \alpha>1$, $0<\beta<\alpha-1$ 且运动函数满足条件 (1.4). 则对于任意的 $p > 1$, 存在常数 $C = C (p) > 0$, 使得对任意 $t \in (0,T_{\max})$ 有

$\int_{\Omega}u^{p} \leqslant C.$

根据引理 3.6, 取 $\displaystyle\varepsilon = \frac{\mu (p+\alpha-1)}{3(p-1)\gamma_2}$ 可以得到, 对任意 $\tau_0 \in (0,T_{\max})$ 和 $ t \in (\tau_0,T_{\max})$ 有

$\begin{equation}\label{3.18} \begin{aligned} \int_{\Omega}u^{p} &\leqslant {\rm e}^{-\frac{p+\alpha-1}{\alpha-1}(t - \tau_0)}\int_{\Omega}u^{p}(\tau_0) + \int_{\tau_0}^{t}{\rm e}^{-\frac{p+\alpha-1}{\alpha-1}(t - \tau)} \Bigg[ -\frac{2\mu p}{3}\int_{\Omega}u^{p + \alpha - 1} \\ &\quad + (\alpha-1)\left(\frac{\mu}{3}\right)^{-\frac{p}{\alpha-1}} \left[\frac{(p-1)\gamma_2}{p+\alpha-1}\right]^{\frac{p+\alpha-1}{\alpha-1}} M_{\frac{p+\alpha-1}{\alpha-1}}\int_{\Omega}u^{\frac{\beta(p+\alpha-1)}{\alpha-1}}(\tau) \\ &\quad + \left(\rho p + \frac{p+\alpha-1}{\alpha-1}\right)\int_{\Omega}u^{p} \Bigg]\mathrm{d}\tau \\ &\quad + (\alpha-1)\left(\frac{\mu}{3}\right)^{-\frac{p}{\alpha-1}} \left[\frac{(p-1)\gamma_2}{p+\alpha-1}\right]^{\frac{p+\alpha-1}{\alpha-1}} M_{\frac{p+\alpha-1}{\alpha-1}} \|v(\cdot,\tau_0)\|_{W^{2,\frac{p+\alpha-1}{\alpha-1}}(\Omega)}^{\frac{p+\alpha-1}{\alpha-1}}. \end{aligned} \end{equation}$

由 Young 不等式得, 存在常数 $C_1,C_2>0$, 使得对任意 $\tau \in (0,T_{\max})$ 有

$\left(\rho p + \frac{p+\alpha-1}{\alpha-1}\right)\int_{\Omega}u^{p} \leqslant \frac{\mu p}{3}\int_{\Omega}u^{p + \alpha - 1} + C_1,$

$\begin{equation*} \begin{aligned} (\alpha-1)\left(\frac{\mu}{3}\right)^{-\frac{p}{\alpha-1}} \left[\frac{(p-1)\gamma_2}{(p+\alpha-1)}\right]^{\frac{p+\alpha-1}{\alpha-1}} M_{\frac{p+\alpha-1}{\alpha-1}}\int_{\Omega}u^{\frac{\beta(p+\alpha-1)}{\alpha-1}}(\tau) \leqslant \frac{\mu p}{3}\int_{\Omega}u^{p + \alpha - 1} + C_2. \end{aligned} \end{equation*}$

该引理得证.

当 $\beta=\alpha-1$ 时, $u$ 的 $L^p$ 估计是在 $\mu > 0$ 足够大的条件下推导出来的.

引理 3.8 设 $\rho>0,\,\alpha>1$, $\beta=\alpha-1$, $p_0>\frac{n\beta}{2}$ 且运动函数满足条件 (1.4). 若

$\mu>\frac{p_0-1}{p_0}\gamma_2M_{\frac{p_0+\alpha-1}{\alpha-1}}^{\frac{\alpha-1}{p_0+\alpha-1}}$

则存在常数 $C = C(p_0) > 0$, 使得对任意 $\tau \in (0,T_{\max})$ 有

$\int_\Omega u^{p_0}{\rm d}x\leqslant C.$

$g(\varepsilon)=\left[p_0\varepsilon+\frac{(\alpha-1)}{\varepsilon^{\frac{p_0}{\alpha-1}}} M_{\frac{p_0+\alpha-1}{\alpha-1}}\right]\frac{(p_0-1)\gamma_2}{p_0+\alpha-1},$

则有

$g'(\varepsilon)=\left[1-\varepsilon^{-\frac{p_0+\alpha-1}{\alpha-1}} M_{\frac{p_0+\alpha-1}{\alpha-1}}\right]\frac{p_0(p_0-1)\gamma_2}{p_0+\alpha-1}.$

故当 $\varepsilon=M_{\frac{p_0+\alpha-1}{\alpha-1}}^{\frac{\alpha-1}{p_0+\alpha-1}}$ 时, $g(\varepsilon)$ 取最小值

$(p_0-1)\gamma_2M_{\frac{p_0+\alpha-1}{\alpha-1}}^{\frac{\alpha-1}{p_0+\alpha-1}}.$

由于

$\mu >\frac{p_0-1}{p_0}\gamma_2M_{\frac{p_0+\alpha-1}{\alpha-1}}^{\frac{\alpha-1}{p_0+\alpha-1}}$

及函数 $g(\varepsilon)$ 的连续性可得, 存在一个常数 $\varepsilon_0 > 0$, 使得

$\begin{equation}\label{3.19} \mu p_0>g(\varepsilon_0)\geqslant(p_0-1)\gamma_2M_{\frac{p_0+\alpha-1}{\alpha-1}}^{\frac{\alpha-1}{p_0+\alpha-1}}. \end{equation}$

由引理 3.6, 取 $p = p_0$ 和 $\varepsilon = \varepsilon_0$, 可推出对任意 $\tau_0 \in (0,T_{\max})$ 和 $t \in (\tau_0,T_{\max})$ 有

$\begin{equation}\label{3.20} \begin{aligned} \int_{\Omega}u^{p_0} & \leqslant {\rm e}^{-\frac{p_0+\alpha-1}{\alpha-1}(t-\tau_0)}\int_{\Omega}u^{p_0}(\tau_0) \\ &\quad+ \int_{\tau_0}^{t}{\rm e}^{-\frac{p_0+\alpha-1}{\alpha-1}(t-\tau)}\left\{-\left[\mu p_0-g(\varepsilon_0)\right]\int_{\Omega}u^{p_0+\alpha- 1}(\tau)\right.\\ &\quad+\left(\rho p_0+ \frac{p_0+\alpha-1}{\alpha-1}\right)\int_{\Omega}u^{p_0}(\tau)\}\mathrm{d}\tau\\ &\quad+\frac{(\alpha-1)(p_0-1)\gamma_2}{(p_0+\alpha-1)\varepsilon_0^{\frac{p_0}{\alpha-1}}} M_{\frac{p_0+\alpha-1}{\alpha-1}}\|v(\cdot,\tau_0)\|_{W^{2,\frac{p_0+\alpha-1}{\alpha-1}}(\Omega)}^{\frac{p_0+\alpha-1}{\alpha-1}}. \end{aligned} \end{equation}$

由 (3.19) 式和 Young 不等式得, 存在一个常数 $C_1>0$, 使得

$\begin{equation*} \begin{aligned} \left(\rho p_0+ \frac{p_0+\alpha-1}{\alpha-1}\right)\int_{\Omega}u^{p_0}(\tau) \leqslant \left[\mu p_0-g(\varepsilon_0)\right]\int_{\Omega}u^{p_0+\alpha- 1}(\tau)+ C_1, \end{aligned} \end{equation*}$

再结合 (3.20) 式, 即得所求结论.

当 $\displaystyle\alpha=\frac{n+2}{n}$ 和 $\displaystyle\beta=\frac{2}{n}$ 时,利用引理 2.3 可导出 $u$ 的 $ L^p$ 估计.

引理3.9 设 $\rho>0,\,\mu>0$, $\displaystyle\alpha=\frac{n+2}{n}$, $\displaystyle\beta=\frac{2}{n}$ 且运动函数满足 (1.4) 式. 那么存在常数 $C>0$, 使得对任意 $t \in (0,T_{\max})$ 有

$\int_\Omega u^\alpha \mathrm{d}x\leqslant C.$

利用系统 (1.2) 的第一个方程、边界条件和分部积分得, 对任意 $t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.21} \begin{aligned} \frac{\rm d}{{\rm d}t}\int_{\Omega}u^{\alpha} =&\alpha\int_{\Omega}u^{\alpha-1}u_t=\alpha\int_{\Omega}u^{\alpha-1}\left(\Delta(\gamma(v)u)+\rho u-\mu u^{\alpha}\right) \\ =&-\alpha(\alpha-1)\int_{\Omega}\gamma(v)u^{\alpha-2}|\nabla u|^{2} -\alpha(\alpha-1)\int_{\Omega}\gamma'(v)u^{\alpha-1}\nabla u\cdot\nabla v \\ &+\rho\alpha\int_{\Omega}u^{\alpha}-\mu\alpha\int_{\Omega}u^{2\alpha-1}. \end{aligned} \end{equation}$

接下来, 估计 (3.21) 式的右侧. 结合 (1.4) 式得, 对任意 $t \in (0,T_{\max})$ 有

$\begin{equation*} \begin{aligned} &\int_{\Omega}\gamma'(v)u^{\alpha-1}\nabla u\cdot\nabla v=\frac{1}{\alpha}\int_{\Omega}\gamma'(v)\nabla u^\alpha\cdot\nabla v \\ =&-\frac{1}{\alpha}\int_{\Omega}u^\alpha\nabla\cdot(\gamma'(v)\nabla v) \\ =&-\frac{1}{\alpha}\int_{\Omega}\gamma'(v)u^\alpha\Delta v-\frac{1}{\alpha}\int_{\Omega}\gamma"(v)u^\alpha|\nabla v|^2 \\ \geqslant&-\frac{\gamma_2}{\alpha}\int_{\Omega}u^\alpha|\Delta v|. \end{aligned} \end{equation*}$

该式与 (3.21) 式和 Young 不等式相结合可得, 对任意 $t \in (0,T_{\max})$ 有

$\begin{equation}\label{3.22} \begin{aligned} \frac{\rm d}{{\rm d}t}\int_{\Omega}u^{\alpha} \leqslant-\alpha(\alpha-1)\gamma_1\int_{\Omega}u^{\alpha-2}|\nabla u|^{2} +(\alpha-1)\gamma_2\int_{\Omega}u^\alpha|\Delta v|+C. \end{aligned} \end{equation}$

再次利用 Young 不等式, 对任意 $t \in (0,T_{\max})$ 有

$\begin{equation*} \begin{aligned} \int_{\Omega}u^\alpha|\Delta v|\leqslant&\left(\int_{\Omega}|\Delta v|^{\frac{\alpha}{\beta}}\right)^{\frac{\beta}{\alpha}} \left(\int_{\Omega}u^{\frac{\alpha^2}{\alpha-\beta}}\right)^{\frac{\alpha-\beta}{\alpha}} = \left(\int_{\Omega}|\Delta v|^{\frac{\alpha}{\beta}}\right)^{\frac{\beta}{\alpha}}\|u^{\frac{\alpha}{2}}\|_{L^{2\alpha}(\Omega)}^2\\ \leqslant& C\left(\int_{\Omega}|\Delta v|^{\frac{\alpha}{\beta}}\right)^{\frac{\beta}{\alpha}} \left(\|\nabla u^{\frac{\alpha}{2}}\|_{L^2(\Omega)}^{\frac{2}{\alpha}}\|u^{\frac{\alpha}{2}}\|_{L^2(\Omega)}^{\frac{2(\alpha-1)}{\alpha}}+\|u^{\frac{\alpha}{2}}\|_{L^2(\Omega)}^2\right)\\ \leqslant&\frac{\alpha\gamma_1}{\gamma_2}\int_{\Omega}u^{\alpha-2}|\nabla u|^{2} +C\int_{\Omega}|\Delta v|^{\frac{\alpha}{\beta}}\int_{\Omega}u^\alpha+C\int_{\Omega}u^{\alpha}. \end{aligned} \end{equation*}$

将上述不等式代入 (3.22) 式得, 对任意 $t \in (0,T_{\max})$ 有

$\frac{\rm d}{{\rm d}t}\int_{\Omega}u^{\alpha}\leqslant C\int_{\Omega}|\Delta v|^{\frac{\alpha}{\beta}}\int_{\Omega}u^\alpha+C\int_{\Omega}u^{\alpha}+C.$

由引理 3.2 和引理 2.2 得, 存在常数 $C>0$, 满足对任意 $t \in (0,T_{\max}-\tau)$ 有

$\int_t^{t+\tau}\int_\Omega u^\alpha \mathrm{d}s\leqslant C $

$\int_t^{t+\tau}\int_{\Omega}|\Delta v|^{\frac{\alpha}{\beta}} \mathrm{d}s\leqslant C.$

利用引理 2.3, 即得所求结论.

3.7 全局存在性: 定理 1.1 的证明

基于上述分析, 证明系统 (1.2) 解的全局存在性和有界性.

[定理 1.1 的证明] 这里将分五种不同情形进行讨论.

情况 1 如果 (H1) 成立, 则

$\alpha>\frac{(n+2)\beta}{2}.$

利用引理 3.2 的 (3.3) 和引理 3.5, 得到 (3.8) 式.

情况 2 如果 (H2) 成立, 则

$1>\frac{n\beta}{2}.$

利用引理 3.2 的 (3.1) 和引理 3.5, 得到 (3.8) 式.

情况 3 如果 (H3) 成立, 则

$\alpha=\frac{(n+2)}{n}>1=\frac{n\beta}{2}.$

利用引理 3.9 和引理 3.5, 得到 (3.8) 式.

情况 4 如果 (H4) 成立, 则

$\beta<\alpha-1.$

利用结合引理 3.7 和引理 3.5, 得到 (3.8) 式.

情况 5 如果 (H5) 成立, 则

$\beta=\alpha-1, \quad \mu>\frac{p_0-1}{p_0}\gamma_2M_{\frac{p_0+\alpha-1}{\alpha-1}}^{\frac{\alpha-1}{p_0+\alpha-1}}.$

利用引理 3.8 和引理 3.5, 得到 (3.8) 式.

由 (3.8) 式得, 系统 (1.2) 的解在任意 $(0,T_{\max})$ 内均满足一致有界. 利用引理 3.1, 即可得到解的全局存在性.

参考文献

Amann H. Nonhomogeneous linear and quasilinear elliptic and parabolic boundary value problems.//Schmeisser H J, Triebel H. Function Spaces, Differential Operators and Nonlinear Analysis. Wiesbaden: Vieweg+Teubner Verlag, 1993, 133: 9--126

[本文引用: 1]

Desvillettes L, Kim Y J, Trescases A, Yoon C W.

A logarithmic chemotaxis model featuring global existence and aggregation

Nonlinear Anal Real World Appl, 2019, 50: 562-582

[本文引用: 1]

Ding M Y, Wang W, Zhou S L, Zheng S N.

Asymptotic stability in a fully parabolic quasilinear chemotaxis model with general logistic source and signal production

J Differential Equations, 2020, 268(11): 6729-6777

[本文引用: 1]

Fujie K, Jiang J.

Boundedness of classical solutions to a degenerate Keller--Segel type model with signal-dependent motilities

Acta Appl Math, 2021, 176(1): Art 3

[本文引用: 1]

Fujie K, Jiang J.

Comparison methods for a Keller-Segel-type model of pattern formations with density-suppressed motilities

Calc Var Partial Differential Equations, 2021, 60(3): Art 92

[本文引用: 2]

Horstmann D, Winkler M.

Boundedness vs. blow-up in a chemotaxis system

J Differential Equations, 2005, 215(1): 52-107

[本文引用: 1]

Ishida S, Seki K, Yokota T.

Boundedness in quasilinear Keller-Segel systems of parabolic-parabolic type on non-convex bounded domains

J Differential Equations, 2014, 256(8): 2993-3010

[本文引用: 1]

Ishida S, Yokota T.

Boundedness in a quasilinear fully parabolic Keller-Segel system via maximal Sobolev regularity.

Discrete Contin Dyn Syst Ser S, 2020, 13(2): 212-232

[本文引用: 1]

Jin H Y, Kim Y J, Wang Z A.

Boundedness, stabilization, and pattern formation driven by density-suppressed motility

SIAM J Appl Math, 2018, 78(3): 1632-1657

[本文引用: 1]

Jin H Y, Wang Z A.

Critical mass on the Keller-Segel system with signal-dependent motility

Proc Amer Math Soc, 2020, 148(11): 4855-4873

[本文引用: 1]

Keller E F, Segel L A.

Initiation of slime mold aggregation viewed as an instability

J Theoret Biol, 1970, 26(3): 399-415

[本文引用: 1]

Lankeit J.

Eventual smoothness and asymptotics in a three-dimensional chemotaxis system with logistic source

J Differential Equations, 2015, 258(4): 1158-1191

[本文引用: 1]

Liu D M, Tao Y S.

Boundedness in a chemotaxis system with nonlinear signal production

Appl Math J Chinese Univ, 2016, 31(4): 379-388

[本文引用: 1]

Liu Z R, Xu J.

Large time behavior of solutions for density-suppressed motility system in higher dimensions

J Math Anal Appl, 2019, 475(2): 1596-1613

[本文引用: 1]

Liu C, Fu X, Liu L, et al.

Sequential establishment of stripe patterns in an expanding cell population

Science, 2011, 334: 238-241

[本文引用: 1]

Lv W B.

Global existence for a class of chemotaxis-consumption systems with signal-dependent motility and generalized logistic source

Nonlinear Anal Real World Appl, 2020, 56: Art 103160

[本文引用: 1]

Lv W B, Wang Q.

A chemotaxis system with signal-dependent motility, indirect signal production and generalized logistic source: Global existence and asymptotic stabilization

J Math Anal Appl, 2020, 488(2): Art 124108

Lv W B, Wang Q Y.

Global existence for a class of chemotaxis systems with signal-dependent motility, indirect signal production and generalized logistic source

Z Angew Math Phys, 2020, 71(2): Art 53

[本文引用: 1]

Lv W B, Wang Q Y.

Global existence for a class of Keller-Segel models with signal-dependent motility and general logistic term

Evol Equ Control Theory, 2021, 10(1): 25-36

[本文引用: 1]

Lv W B, Wang Q Y.

An $n$-dimensional chemotaxis system with signal-dependent motility and generalized logistic source: global existence and asymptotic stabilization

Proc Roy Soc Edinburgh, 2021, 151(2): 821-841

[本文引用: 5]

Lyu W B

Weak solutions to a class of signal-dependent motility Keller-Segel systems with superlinear damping

J Differential Equations, 2022, 312: 209-236

[本文引用: 2]

Lyu W B, Wang Z A.

Global classical solutions for a class of reaction-diffusion system with density-suppressed motility

Electron Res Arch, 2022, 30(3): 995-1015

[本文引用: 1]

Nakaguchi E, Osaki K.

Global existence of solutions to a parabolic-parabolic system for chemotaxis with weak degradation

Nonlinear Anal, 2011, 74(1): 286-297

[本文引用: 1]

Tao X Y, Fang Z B.

Global boundedness and stability in a density-suppressed motility model with generalized logistic source and nonlinear signal production

Z Angew Math Phys, 2022, 73(3): Art 123

[本文引用: 2]

Tao Y S, Winkler M.

A chemotaxis-haptotaxis model: the roles of nonlinear diffusion and logistic source

SIAM J Math Anal, 2011, 43(2): 685-704

[本文引用: 1]

Tao Y S, Winkler M.

Boundedness in a quasilinear parabolic-parabolic Keller-Segel system with subcritical sensitivity

J Differential Equations, 2012, 252(1): 692-715

[本文引用: 2]

Tao Y S, Winkler M.

Effects of signal-dependent motilities in a Keller-Segel-type reaction-diffusion system

Math Models Methods Appl Sci, 2017, 27(9): 1645-1683

[本文引用: 1]

Tello J I, Winkler M.

A chemotaxis system with logistic source

Comm Partial Differential Equations, 2007, 32(4-6): 849-877

[本文引用: 1]

Temam R. Infinite-Dimensional Dynamical Systems in Mechanics and Physics. New York: Springer-Verlag, 1988

[本文引用: 1]

Viglialoro G.

Very weak global solutions to a parabolic-parabolic chemotaxis-system with logistic source

J Math Anal Appl, 2016, 439(1): 197-212

[本文引用: 1]

Wang J P, Wang M X.

Boundedness in the higher-dimensional Keller-Segel model with signal-dependent motility and logistic growth

J Math Phys, 2019, 60(1): Art 011507

[本文引用: 1]

Wang W.

Global boundedness of solutions resulting from both the self-diffusion and the logistic-type source

Z Angew Math Phys, 2019, 70(4): Art 99

[本文引用: 1]

Wang Z A, Zheng J S.

Global boundedness of the fully parabolic Keller-Segel system with signal-dependent motilities

Acta Appl Math, 2021, 171(1): Art 25

[本文引用: 1]

Winkler M.

Chemotaxis with logistic source: very weak global solutions and their boundedness properties

J Math Anal Appl, 2008, 348(2): 708-729

[本文引用: 1]

Winkler M.

Aggregation vs. global diffusive behavior in the higher-dimensional Keller-Segel model

J Differential Equations, 2010, 248(12): 2889-2905

[本文引用: 1]

Winkler M.

Boundedness in the higher-dimensional parabolic-parabolic chemotaxis system with logistic source

Comm Partial Differential Equations, 2010, 35(8): 1516-1537

[本文引用: 1]

Winkler M.

Blow-up in a higher-dimensional chemotaxis system despite logistic growth restriction

J Math Anal Appl, 2011, 384(2): 261-272

[本文引用: 1]

Winkler M.

Global solutions in a fully parabolic chemotaxis system with singular sensitivity

Math Methods Appl Sci, 2011, 34(2): 176-190

[本文引用: 1]

Winkler M.

Finite-time blow-up in the higher-dimensional parabolic-parabolic Keller-Segel system

J Math Pures Appl, 2013, 100(5): 748-767

[本文引用: 1]

Winkler M.

Finite-time blow-up in low-dimensional Keller-Segel systems with logistic-type superlinear degradation

Z Angew Math Phys, 2018, 69(2): Art 69

[本文引用: 1]

Winkler M.

Can simultaneous density-determined enhancement of diffusion and cross-diffusion foster boundedness in Keller-Segel type systems involving signal-dependent motilities?

Nonlinearity, 2020, 33(12): 6590-6623

[本文引用: 1]

Winkler M.

The role of superlinear damping in the construction of solutions to drift-diffusion problems with initial data in $L^1$

Adv Nonlinear Anal, 2020, 9(1): 526-566

[本文引用: 1]

Winkler M.

Global generalized solvability in a strongly degenerate taxis-type parabolic system modeling migration-consumption interaction

Z Angew Math Phys, 2023, 74(1): Art 32

[本文引用: 1]

Winkler M.$L^1$ solutions to parabolic Keller-Segel systems involving arbitrary superlinear degradation. Ann Sc Norm Super Pisa Cl Sci (5), 2023, 24(1): 141-172

[本文引用: 1]

Winkler M.

A quantitative strong parabolic maximum principle and application to a taxis-type migration-consumption model involving signal-dependent degenerate diffusion.

Ann Inst H Poincaré C Anal Non Linéaire, 2024, 41(1): 95-127

[本文引用: 1]

Xiang T.

Sub-logistic source can prevent blow-up in the 2D minimal Keller-Segel chemotaxis system

J Math Phys, 2018, 59(8): Art 081502

[本文引用: 1]

Yan J L, Fuest M.

When do Keller-Segel systems with heterogeneous logistic sources admit generalized solutions?

Discrete Contin Dyn Syst, 2021, 26(8): 4093-4109

[本文引用: 1]

Yoon C, Kim Y J.

Global existence and aggregation in a Keller-Segel model with Fokker-Planck diffusion

Acta Appl Math, 2017, 149: 101-123

[本文引用: 1]

Zhuang M D, Wang W, Zheng S N.

Boundedness in a fully parabolic chemotaxis system with logistic-type source and nonlinear production

Nonlinear Anal Real World Appl, 2019, 47: 473-483

[本文引用: 1]

/