Loading [MathJax]/jax/element/mml/optable/BasicLatin.js

数学物理学报, 2025, 45(2): 512-533

具奇异敏感的趋向性犯罪模型中的警察威慑效应

李彬,1, 谢莉,2,*

1宁波工程学院统计与数据科学学院 浙江宁波 315211

2重庆师范大学数学科学学院 重庆 401331

The Effects of Police Deployment in a Chemotaxis System with Singular Sensitivity for Criminal Activities

Li Bin,1, Xie Li,2,*

1School of Statistics and Data Science, Ningbo University of Technology, Zhejiang Ningbo 315211

2School of Mathematical Sciences, Chongqing Normal University, Chongqing 401331

通讯作者: * 谢莉,E-mail: mathxieli@cqnu.edu.cn

收稿日期: 2023-11-28   修回日期: 2024-10-23  

基金资助: 宁波市自然基金(2022J147)
重庆市科委项目(CSTB2023NSCQ-MSX0411)
重庆市教委科研项目(KJZD-M202000502)
重庆市教委科研项目(CXQT21014)

Received: 2023-11-28   Revised: 2024-10-23  

Fund supported: NSF of Ningbo Municipality(2022J147)
Chongqing Science and Technology Commission Project(CSTB2023NSCQ-MSX0411)
Research Project of Chongqing Education Commission(KJZD-M202000502)
Research Project of Chongqing Education Commission(CXQT21014)

作者简介 About authors

李彬,E-mail:blimath@163.com

摘要

该文在有界光滑域 ΩRn(n3)上研究了一个带奇异灵敏度的两组分非局部模型, 该模型是三组分的 Jones-Brantingham-Chayes 趋向性模型的一个简化模型, 后者被用于模拟在警察威慑下犯罪活动的时空动态. 该文在较大趋化敏感系数范围内证明了相应初边值问题拥有全局经典解. 值得指出的是, 相较于无警察威慑效应的 Short et al 趋向性犯罪模型的相关结果, 警察威慑扩大了确保解全局存在的趋化敏感系数范围, 在某种意义下这也表明了警察威慑效应对模型解性质具有正则化效应. 注意, 先前数值结果 (Jones, Brantingham and Chayes. Math Models Methods Appl Sci, 2010) 表明警察威慑有益于镇压犯罪热点的形成, 因此该文的研究结果也是相应数值结果的一个理论支持.

关键词: 趋向性模型; 奇异敏感; 经典解; 全局存在性; 长时间行为

Abstract

As a simplified version of a three-component chemotaxis system introduced by Jones et al to model the spatio-temporal behavior of criminal activities under the effects of police deployment, a two-component non-local system with singular sensitivity is considered over a bounded domain ΩRn with n3. For all reasonably regular initial data, the existence of classical solution of the corresponding initial-boundary value problem is established globally in time. In particular, we enlarge the range of chemotactic sensitivity χ1, compared to known results on the Short et al model which describes the spatio-temporal behavior of criminal activities without the effects of police deployment, which reveals that the police deployment has a regularization effect on solution. This is also a theoretical supplement to the previously numerical result unveiled by Jones, Brantingham and Chayes (Math Models MethodsAppl Sci, 2010) that the police deployment is beneficial for suppressing the formation of criminal hotspots.

Keywords: chemotaxis system; singular sensitivity; classical solution; global existence; long-time behavior

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

本文引用格式

李彬, 谢莉. 具奇异敏感的趋向性犯罪模型中的警察威慑效应[J]. 数学物理学报, 2025, 45(2): 512-533

Li Bin, Xie Li. The Effects of Police Deployment in a Chemotaxis System with Singular Sensitivity for Criminal Activities[J]. Acta Mathematica Scientia, 2025, 45(2): 512-533

1 引言及主要结果

为了模拟犯罪活动的时空动态, Short 等[35],[36],[37]建立了如下趋向性犯罪模型

{ut=Δuχ1(ulnv)κ1uv+h1,vt=Δvv+uv+h2,
(1.1)

其中 χ1>0, κ1>0, u(x,t) 刻画了犯罪密度, v(x,t) 表示抽象的犯罪吸引力, hi 为相应源项, i=1,2. 模型 (1.1) 的相关背景及推导过程, 可参看文献 [2],[4],[13],[39],[38],[42]. 更多犯罪模型和趋向性模型可参看综述文献 [3],[8].

模型 (1.1) 一被建立, 在有界光滑域 ΩRn上的初边值问题便受到了广泛的关注: 其经典解的局部存在性可参看文献 [31]; 如果进一步假设 n=1[32],[45], 或者 n2χ1<2n[10],[34], 或者初始值和 hi 都小[1],[40], i=1,2, 则该解也是全局存在的. 在广义可解框架下, 趋化敏感系数 χ1 的范围能被进一步扩大: 如果 n=2[48], 或者 n=3χ1(0,3)[16], 则重整化解是全局存在的; 对任意的 χ1>0, 二维的广义下解是全局存在的[22]. 另外, 非线性扩散 (i.e., Δum) 或 logistic 源 (i.e., aubuα) 并入 (1.1) 式的第一个方程都能在适当意义下扩大趋化系数的范围: 假设 m>32[33], 或者 m>1χ1<32[52], 则相应的二维初边值问题具有全局弱解; 当 n=2α=2 时, 相应问题有全局广义解[14], 其最终光滑性参见文献 [30], 经典解的全局存在性则至少需要 α>2[14],[44]. 另一方面, 解的无界性探讨则可参见文献 [11]; 解的长时间行为的研究则可参见文献 [1],[16],[32],[34],[40],[48]; 更多相关稳态问题的研究则可参见文献 [5],[7],[18],[26],[27],[28],[37],[43].

最近, 多种警察策略被并入 Short et al 模型 (1), 具体可参见文献 [17],[36],[53],[29],[37]. 特别地, 假设警察随机移动的同时也定向向高犯罪率区域移动, Jones 等[17]建立了如下带警察策略的趋向性犯罪模型

{ut=Δuχ1(ulnv)κ1uvκ2uκ3uw+h1, xΩ, t>0,[1mm]vt=Δvv+uv+h2, xΩ, t>0,[1mm]τwt=Δwχ(wlnv), xΩ, t>0,[1mm]uν=vν=wν=0, xΩ, t>0,[1mm]u(x,0)=u0(x),v(x,0)=v0(x),τw(x,0)=τw0, xΩ,
(1.2)

其中 τ0, χ1>0, χ>0, κi>0, i=1,2,3. 这里, w(x,t) 表示警察密度, ΩRn 是一个带光滑边界 Ω 的有界域, ν 为边界的外法向量. 值得指出的是, 文献 [17] 中的数值结果表明警察策略是可以减少犯罪活动的, 但目前还没有这方面的理论结果.

1.1 主要结果

为了填补这方面的空白, 本文在 τ=0 情况下探讨警察策略对相关问题解的正则效应. 类似于文献 [41], 在这种情况下三组分的系统 (1.2) 能被转化为两组分的系统. 事实上, 把 (1.2) 式中的第三个方程乘以 lnwχlnv, 并通过分部积分可得

Ωw|(lnwχlnv)|2dx=0.

注意, 最值原理表明 w>0; 据此, 进一步可得

w(x,t):=mwvχΩvχ(,t)dx,

其中 mw:=Ωwdx>0 表示总警察人数. 为了简单, 设 κ1=κ2=1, κ:=κ3mw, 此时三组分的系统 (1.2) 可以简化为如下带非局部阻尼的两组分趋化模型

{ut=Δuχ1(ulnv)uvuκuvχΩvχdx+h1, xΩ, t>0,[1mm]vt=Δvv+uv+h2, xΩ, t>0,[1mm]uν=vν=0, xΩ, t>0,[1mm]u(x,0)=u0(x),v(x,0)=v0(x), xΩ.
(1.3)

相较于 (1.1), (1.3) 式中出现了一个非局部阻尼项 κuvχΩvχdx, 该项潜在地刻画了警察威慑效应. 故接下来我们将集中探讨该非局部阻尼项对解性质的正则效应.

假设源函数满足

0hiC1(¯Ω×[0,))L(Ω×(0,)),i=1,2,
(1.4)

初始值 (u0,v0) 满足

{0u0C0(¯Ω),u00,[1mm]v0W1,(¯Ω),inf
(1.5)

基于假设 (1.4)-(1.5) 式, 设 n\geq3

\widetilde{\chi}_1:=\widetilde{\chi}_1(\chi)=\left\{\begin{split}&\frac{2}{n},&\chi\in\left(0,1\right],\\&\frac{2}{n}\cdot\sqrt{\frac{\chi(n-\chi)}{n-1}},&\chi\in\left(1,\frac{n+2}{4}\right],\\&\max\left\{\frac{2}{n}\cdot\sqrt{\frac{\chi(n-\chi)}{n-1}},\,\sqrt{\frac{4\chi-n}{(2\chi-1)n}}\right\},&\chi\in\left(\frac{n+2}{4},\frac n2\right),\\&\sqrt{\frac{4\chi- n}{(2\chi-1)n}}, &\chi\geq\frac{n}{2}.\end{split}\right.
(1.6)

据此, 初边值问题 (1.3) 的经典解的全局存在性可表述如下

定理1.1 n\geq3, \chi_1>0, \kappa>0, \chi>0 且 (1.4)- (1.5) 式成立. 假设 \widetilde{\chi}_1 由 (1.6) 式确定且 \chi_1<\widetilde{\chi}_1, 则初边值问题 (1.3) 拥有唯一的经典解 (u,v). 该解满足 u,v>0, 且对任意的 r>n 都有

\left\{\begin{split}&u \in \mathcal{C}^0\big(\overline{\Omega}\times[\\ [1mm]&v\in \mathcal{C}^0\big([0,\infty); W^{1,r}(\overline{\Omega})\big)\cap \mathcal{C}^{2,1}\big(\overline{\Omega}\times(0,\infty)\big).\end{split}\right.

注1.1 一个直接的计算可知: \widetilde{\chi}_1:=\widetilde{\chi}_1(\chi) 是关于 \chi 单增的且连续的 (\chi 刻画了警察执法效率). 据此, 若 \chi>1, 则 \widetilde{\chi}_1>\frac2n. 注意, 条件 \chi_1<\frac2n 能确保不带警察威慑的系统具有全局经典解[10]. 因此, 相较于文献 [10] 中结果, 定理 1.1 表明: 在某种意义下警察威慑有利于镇压犯罪热点的形成; 且警察威慑效应越大, 犯罪热点越难形成. 此外, 当 \chi\rightarrow+\infty 时, \sqrt{\frac{2}{n}}>\widetilde{\chi}_1\rightarrow\sqrt{\frac2n}; 条件 \chi_1<\sqrt{\frac2n} 能确保具线性信号产出的对数 Keller-Segel 系统

\left\{\begin{split}&u_t= \Delta u-\chi_1\nabla\cdot\left(u\nabla\ln v\right),\\&v_t=\Delta v- v+ u\end{split}\right.

拥有全局经典解[47]. 这也从数学上表明了当 (1.3) 式中的 \chi 充分大时, 非局部项 -\frac{\kappa uv^\chi}{\int_\Omega v^\chi{\rm d}x} 似乎能消除由非线性产出项 +uv 带来的非光滑效应.

在适当增加关于 h_2 的假设下, 如

\begin{align} h_2(x,t)\geq\delta>0,\quad x\in\Omega,\,\, t>0,\label{eq-iva1} \end{align}
(1.7)

我们能得到经典解的一致有界性. 除此之外, 如果进一步假设

\int_0^\infty\int_\Omega h_1(\cdot,t){\rm d}x{\rm d}t<\infty,
(1.8)
\int_0^\infty\int_\Omega|h_2(\cdot,t)-h_{2,\infty}(\cdot)|{\rm d}x{\rm d}t<\infty,\,\,\, h_{2,\infty}\in \mathcal{C}^1(\overline{\Omega}),
(1.9)

则该有界解将在时间充分大时收敛于相应的平衡态. 具体而言, 经典解的有界性和长时间行为能被刻画如下

定理1.2 令定理 1.1 中条件和 (1.7) 式都成立, 且 (u,v) 源至定理 1.1. 对于任意的 r>n, 存在独立于 tC>0 使得

\begin{align} \label{eq-bedd} \|u(\cdot,t)\|_{L^\infty}+\|v^{-1}(\cdot,t)\|_{L^\infty}+\| v(\cdot,t)\|_{W^{1,r}}\leq C,\quad t>0. \end{align}
(1.10)

更进一步假设 (1.8) 和 (1.9) 式成立, 则当 t\rightarrow\infty 时有

\begin{align}\label{eq-deu} \|u(\cdot,t)\|_{L^\infty}+\|v(\cdot,t)-v_\infty(\cdot)\|_{W^{1,r}}\rightarrow0, \end{align}
(1.11)

其中 v_\infty 为椭圆问题

\left\{\begin{split}&0=\Delta v_\infty-v_\infty+h_{2,\infty},&x\in\Omega,\\ [1mm]&\nabla v_\infty\cdot\nu=0,&x\in\partial\Omega\end{split}\right.
(1.12)

的解. 这里, h_{2,\infty}\in \mathcal{C}^1(\overline{\Omega}) 且满足 (1.9) 式.

1.2 技术策略与文章结构

本文的第一个目的是, 调查初边值问题 (1.3) 经典解的全局存在性. 注意, (1.3) 式中第一个方程的非局部项显然破坏了一些已知方法, 如被用于 (1.1) 式[10], 或具线性信号产出的对数 Keller-Segel 模型[49],[47], 或线性产出项 +u 代替了 (1.1) 式中非线性信号产出项 +uv[23],[25] 的方法. 据此, 为了发展一个可解性理论, 我们的出发点是, 利用

\begin{align*} \int_s^t \int_\Omega u v^\chi{\rm d}x \left\{\int_\Omega v^\chi{\rm d}x\right\}^{-1}{\rm d}\sigma \end{align*}

的有界性 (参见 (2.8) 式), 考察 \|v(\cdot,t)\|_{L^\chi} (\chi>1) 的时间演化且建立其一致有界性, 细节参见引理 3.1. 进而联立第 3 节中建立的一致先验估计, 可找到 p_1>\frac n2 使得 \|u(\cdot,t)\|_{L^{p_1}} 是有界的, 其中难点是怎么确定趋化系数 \chi_1 的范围, 细节可参看引理 3.2 和引理 3.3. 基于这些得到的有界性, 通过迭代策略可得经典解的全局可解性. 值得指出的是, 相较于不具警察威慑效应的 Short et al 模型 (1), 我们扩大了趋化系数 \chi_1 的范围. 因此, 我们的结果是先前数值结果 (警察威慑有益于镇压犯罪热点的形成[17]) 的一个理论补充.

本文的第二个目标是, 探讨解的大时间行为. 为此, 我们首先利用解成分 v 的点态下界建立解的有界性, 参见引理 5.1. 其次, 利用此有界性, 我们可得 \|u(\cdot,t)\|_{L^1}\|v(\cdot,t)-v_\infty(\cdot)\|_{L^2} 的衰减估计, 其中 v_\infty 为椭圆问题 (1.12) 的解. 最后, 我们进一步发展文献 [12],[30] 中的方法, 并利用迭代法得到了期望的解的长时间行为.

本文剩下部分内容安排如下: 第 2 节将证明 (1.3) 式经典解的局部存在性; 一些关键的先验估计将在第 3 节建立; 在第 4 节中, 我们构建了初边值问题 (1.3) 经典解的全局存在性; 第 5 节致力于证明定理 1.2 中期望的解的长时间行为.

符号说明: 除了特别说明之外, 通篇文章中 c_iC_i 都表示某个正常数, 这些常数在不同的地方可能是不同的.

2 局部存在性与准备工作

首先, 我们关注于初边值问题 (1.3) 的局部可解性以及延拓准则.

引理2.1\chi,\chi_1,\kappa>0 以及 \Omega\subset\mathbb{R}^n (n\geq3) 是一有界光滑域. 对任意满足 (1.4)- (1.5) 式的 h_i 和初始值, 存在一个时间 T_{\max}\in(0,\infty] 以及对任意 r>n 存在唯一 (u,v) 满足

\left\{\begin{array}{l} u \in \mathcal{C}^{0}\left(\bar{\Omega} \times\left[0, T_{\max }\right)\right) \cap \mathcal{C}^{2,1}\left(\bar{\Omega} \times\left(0, T_{\max }\right)\right), \\ v \in \mathcal{C}^{0}\left(\left[0, T_{\max }\right) ; W^{1, r}(\bar{\Omega})\right) \cap \mathcal{C}^{2,1}\left(\bar{\Omega} \times\left(0, T_{\max }\right)\right) \end{array}\right.

使得 (u,v) 是初边值问题 (1.3) 在 \Omega\times[0,T_{\max}) 上的经典解. 进一步, 若 T_{\max}<\infty, 则有

\begin{matrix}\label{eq-blowup} \limsup\limits_{t\rightarrow T_{\max}}\big(\|u(\cdot,t)\|_{L^\infty}+\|v^{-1}\|_{L^\infty}+\|\nabla v(\cdot,t)\|_{L^r}\big)=\infty. \end{matrix}
(2.1)

T\in(0,1) 且充分小. 对给定的非负函数 \overline{u} \in \mathcal{C}^{\theta, \frac{\theta}2} (\overline{\Omega}\times(0,T)) (\theta\in(0,1)), 由 (1.4)-(1.5) 式可知初边值问题

\left\{\begin{split}&v_{t}=\Delta v -v+ \overline{u}v +h_2, & x \in \Omega,\, t>0,\\&\frac{\partial v}{\partial \nu}=0, & x \in \partial \Omega, \,t>0,\\&v(x, 0)=v_{0}(x), & x \in \Omega\end{split}\right.

拥有唯一非负经典解 v:=v[\overline{u}]\in\mathcal{C}^0\big([0,T); W^{1,p}(\overline{\Omega})\big)\cap\mathcal{C}^{2+\theta_1,1+\frac{\theta_1}{2}}\big(\overline{\Omega}\times(0,T)\big), 这里, p>n\theta_1\in(0,1). 进一步, 利用最值原理和 Neumann 热半群的性质可得

v(x, t)\geq {\rm e}^{t(\Delta-1)} v_{0} \geq {\rm e}^{-t}\inf _{x \in \overline{\Omega}} v_{0}(x)>0,\quad x\in\Omega,\,\, t\in(0,T).
(2.2)

据此, 我们有

v^{-1}(x, t)\leq {\rm e}^{t}\left(\inf _{x \in \overline{\Omega}} v_{0}(x)\right)^{-1},\quad x\in\Omega,\,\, t\in(0,T),

以及

\left(\int_\Omega v^\chi{\rm d}x\right)^{-1}\leq {\rm e}^{\chi t}\left(\inf _{x \in \overline{\Omega}} v_{0}(x)\right)^{-\chi}|\Omega|^{-1},\quad t\in(0,T).

因此, 对于 v:=v[\overline{u}] 再次利用 (1.4)-(1.5) 式和抛物 Schauder 估计[20]可知

\left\{\begin{split}&u_{t}=\Delta u-\chi_1\nabla \cdot\left(u\nabla\ln v\right)-u v-u-\frac{\kappa u v^\chi}{\int_\Omega v^\chi{\rm d}x}+ h_1, &\ x\in \Omega,\ t>0,\\&\frac{\partial u}{\partial \nu}=0, & x \in \partial \Omega, \,t>0,\\&u(x, 0)=u_{0}(x), & x \in \Omega\end{split}\right.
(2.3)

拥有唯一非负经典解 u:=u[v]\in\mathcal{C}^0\big([0,T);L^\infty(\overline{\Omega})\big)\cap\mathcal{C}^{2+\theta_2,1+\frac{\theta_2}{2}}\big(\overline{\Omega}\times(0,T)\big), 这里 \theta_2\in(0,1). 据此, 对于任意的 \overline{u}\in X_R 我们能定义映射 \Psi(\overline{u})=u, 这里

X_R:=\{f:\,\, 0\leq f\in \mathcal{C}^{\theta, \frac{\theta}2}(\overline{\Omega}\times(0,T)),\,\,\,\|f\|_{L^\infty}\leq R\}

以及 R 仅依赖于初始值 (u_0,v_0). 对于充分小的 T, 易证明 u\in X_R 以及 \Psi 是映 X_R 到自身的压缩映射, 这也表明 (u,v) 是初边值问题 (1.3) 在 \overline{\Omega}\times(0,T) 上的经典解. (2.2) 式确保了 v 的正性. 为了验证 u 的正性, 记 z 是如下问题

\left\{\begin{split}&z_{t}=\Delta z-\nabla \cdot\left(a(x,t)z\right)+b(x,t)z, &\ x\in \Omega,\ t>0,\\&\frac{\partial z}{\partial \nu}=0, & x \in \partial \Omega, \,t>0,\\&z(x, 0)=z_0(x)=u_{0}(x), & x \in \Omega\end{split}\right.
(2.4)

的经典解, 这里 a(x,t):=\chi_1\nabla\ln v 以及 b(x,t):=-v-1-\frac{\kappa v^\chi}{\int_\Omega v^{\chi}{\rm d}x}. 接下来, 我们证明 z>0. 事实上, 据边界条件可得 a\cdot\nu|_{\partial\Omega}=0; 基于假设 (4), 可找到 K>0 使得 0\leq z_0\leq K\int_\Omega z_0 {\rm d}x \geq K^{-1}; 利用 (2.2) 式和 v 的正则性, 存在 p_1\geq2,p_2\geq1, q_1>2, q_2>1 满足

\frac{1}{q_1}+\frac{n}{2p_1}<\frac12,\quad \frac{1}{q_2}+\frac{n}{2p_2}<1,

使得对给定的 T>0 都有

\int_0^T\|a(\cdot,t)\|_{L^{p_1}}^{q_1}{\rm d}t\leq K,\quad \int_0^T\|b(\cdot,t)\|_{L^{p_2}}^{q_2}{\rm d}t\leq K.

基于此, 利用文献 [50,命题 1.1] 可得在 \Omega\times(0,T) 上都有 z(x,t)\geq C(T,K,p_1,p_2,q_1,q_2)>0. 另一方面, 由比较原理和 h_1\geq0 易知 z 是初边值问题 (2.3) 的一个下解; 故 u\geq z>0.

类似于映射 \Psi 的压缩性证明, 易知初边值问题 (1.3) 的经典解是唯一的.

更进一步, 利用迭代策略可将 (u,v) 延拓到最大时间区间 [0,T_{\max}) 上, 这里要么 T_{\max}=\infty 要么 T_{\max}<\infty. 如果 T_{\max}<\infty, 则延拓准则 (2.1) 式成立.

为了展示 T_{\max}=\infty, 我们开始于一些基本的先验估计.

引理2.2 假设 (u,v) 源至引理 2.1, 则有

v(x, t)\geq c_1 {\rm e}^{-t}>0,\quad x\in\Omega,\,\,t \in(0,T_{\max}),
(2.5)
\|u(\cdot,t)\|_{L^1}\leq c_2,\quad t \in(0,T_{\max}),
(2.6)
\|v(\cdot,t)\|_{L^1} \leq c_3,\quad t \in(0,T_{\max}),
(2.7)

且对任意的 t \in(0,T_{\max}) 以及 s\in[t]

\begin{matrix} \int_s^t\int_\Omega u v{\rm d}x{\rm d}\sigma+\kappa \int_s^t \int_\Omega u v^\chi{\rm d}x \left\{\int_\Omega v^\chi{\rm d}x\right\}^{-1}{\rm d}\sigma\leq c_2(1+t-s), \label{eq-0vL-2} \end{matrix}
(2.8)

其中

c_1:=\inf _{x \in \overline{\Omega}} v_{0}(x),
c_2:=\max\{\|u_0\|_{L^1}, \|h_1\|_{L^\infty(0,\infty;L^1(\Omega))}\}

以及

c_3:=\max\{\|u_0\|_{L^1}+\|v_0\|_{L^1}, \|h_1\|_{L^\infty(0,\infty;L^1(\Omega))})+\|h_2\|_{L^\infty(0,\infty;L^1(\Omega))}\}.

事实上, 类似于 (2.2) 式, 直接可得 (2.5) 式. 另外, 对 (1.3) 式中的第一个方程积分可得

\begin{matrix}\label{eq-ul1a1} \frac{\rm d}{{\rm d}t}\int_\Omega u{\rm d}x+ \int_\Omega u v {\rm d}x+ \int_\Omega u{\rm d}x +\kappa \int_\Omega u v^\chi{\rm d}x \left\{\int_\Omega v^\chi{\rm d}x\right\}^{-1}=\int_\Omega h_1{\rm d}x, \end{matrix}
(2.9)

结合 (1.4) 式可知对任意的 t\in(0,T_{\max}) 都有

\begin{matrix}\label{eq-ul1a2} \int_\Omega u{\rm d}x+{\rm e}^{-t}\int_0^t{\rm e}^s\int_\Omega uv {\rm d}x{\rm d}s\leq {\rm e}^{-t}(\|u_0\|_{L^1}-\|h_1\|_{L^\infty(0,\infty;L^1(\Omega))})+\|h_1\|_{L^\infty(0,\infty;L^1(\Omega))}. \end{matrix}
(2.10)

这也表明了 (2.6) 式对 (2.9) 式进行时间积分, 则对任意的 t\in(0,T_{\max}) 0\le s\le t 都有

\begin{align*} \int_s^t \int_\Omega u v {\rm d}x{\rm d}\sigma+\kappa\int_s^t \int_\Omega u v^\chi{\rm d}x \left\{\int_\Omega v^\chi{\rm d}x\right\}^{-1}{\rm d}\sigma \leq \|h_1\|_{L^\infty(0,\infty;L^1(\Omega))}(t-s)+\|u(\cdot, s)\|_{L^1}. \end{align*}

因此, (2.8) 式成立.

类似地, 对 (1.3) 式中第二个方程积分, 则有

\begin{align*} \frac{\rm d}{{\rm d}t}\int_\Omega v{\rm d}x+ \int_\Omega v{\rm d}x = \int_\Omega u v{\rm d}x +\int_\Omega h_2{\rm d}x, \end{align*}

结合 (1.4) 式进一步可得

\begin{align*} \int_\Omega v{\rm d}x\leq {\rm e}^{-t}\int_0^t{\rm e}^s\int_\Omega uv {\rm d}x{\rm d}s+{\rm e}^{-t}(\|v_0\|_{L^1}-\|h_2\|_{L^\infty(0,\infty;L^1(\Omega))})+\|h_2\|_{L^\infty(0,\infty;L^1(\Omega))}. \end{align*}

又因 (2.10) 式, 可得

\begin{align*} \int_\Omega v{\rm d}x\leq \max\{\|u_0\|_{L^1}+\|v_0\|_{L^1}, \|h_1\|_{L^\infty(0,\infty;L^1(\Omega))})+\|h_2\|_{L^\infty(0,\infty;L^1(\Omega))}\},\quad t\in(0,T_{\max}), \end{align*}

因此, (2.7) 式成立.

作为准备, 我们将利用文献 [21],[47] 中策略建立 \int_\Omega u^pv^{-r}{\rm d}x (p>1, r>0) 的有界性.

引理2.3\chi_1<1(u,v) 来自引理 2.1. 令

\begin{matrix} r_{\pm}(p):=\frac{p-1}{2}\left(1\pm\sqrt{1-p\chi_1^2}\right).\label{eq-defr} \end{matrix}
(2.11)

对任意的 p\in(1,\frac{1}{\chi_1^2})r\in(r_{-}(p),r_{+}(p)) 存在 c_4>0 使得

\begin{matrix} \int_\Omega u^pv^{-r}(\cdot,t){\rm d}x\leq c_4{\rm e}^{rt},\quad t\in(0,T_{\max}).\label{eq-Eulnu} \end{matrix}
(2.12)

基于 (1.3) 式, 直接的计算表明

\begin{align*} \frac{\rm d}{{\rm d}t}\int_\Omega u^pv^{-r}{\rm d}x=&-r\int_\Omega u^pv^{-r-1} v_{t}{\rm d}x+p\int_\Omega u^{p-1}v^{-r}u_t {\rm d}x\\ =&-r\int_\Omega u^pv^{-r-1}\left\{\Delta v-v+uv+h_2\right\}{\rm d}x\\ &+p\int_\Omega u^{p-1}v^{-r}\left\{\Delta u-\chi_1\nabla \cdot\left(u\nabla\ln v\right)- u v-u-\frac{\kappa u v^\chi}{\int_\Omega v^\chi{\rm d}x}+ h_1\right\}{\rm d}x. \end{align*}

利用分部积分可得

\begin{align*} \frac{\rm d}{{\rm d}t}\int_\Omega u^pv^{-r}{\rm d}x=&-r(r+1+p\chi_1) \int_\Omega u^pv^{-r-2} |\nabla v|^2{\rm d}x-p(p-1)\int_\Omega u^{p-2}v^{-r} |\nabla u|^2{\rm d}x\\ &+(2pr+\chi_1 p(p-1))\int_\Omega u^{p-1}v^{-r-1} \nabla v\cdot\nabla u{\rm d}x\\ &+r\int_\Omega u^pv^{-r}{\rm d}x-r\int_\Omega u^{p+1}v^{-r}{\rm d}x-r\int_\Omega u^pv^{-r-1} h_2{\rm d}x\\ &-p\int_\Omega u^{p}v^{-r+1}{\rm d}x-p\int_\Omega u^{p}v^{-r}{\rm d}x\\ &-p\kappa \left\{\int_\Omega v^\chi{\rm d}x\right\}^{-1}\int_\Omega u^{p}v^{-r+\chi}{\rm d}x+p\int_\Omega u^{p-1}v^{-r} h_1{\rm d}x. \end{align*}

运用 Hölder 不等式和 Young 不等式可知

\begin{align*} p\int_\Omega u^{p-1}v^{-r} h_1{\rm d}x\leq& p\|h_1\|_{L^\infty(\Omega\times(0,\infty))}\int_\Omega u^{p-1}v^{-r}{\rm d}x\\ \leq&\frac p2 \int_\Omega u^{p}v^{-r}{\rm d}x+C_1\int_\Omega v^{-r}{\rm d}x \end{align*}

其中 C_1:=\left(\frac{2(p-1)}{p}\right)^{p-1}\|h_1\|_{L^\infty(\Omega\times(0,\infty))}^p, 联立 (2.5) 式, 进一步可知

\begin{align*} p\int_\Omega u^{p-1}v^{-r} h_1{\rm d}x\leq \frac p2 \int_\Omega u^{p}v^{-r}{\rm d}x+C_2{\rm e}^{rt}, \end{align*}

这里 C_2:=C_1c_1^{-r}|\Omega|c_1 源至 (2.5) 式. 类似地, 我们有

\begin{align*} r\int_\Omega u^pv^{-r}{\rm d}x\leq \frac r2\int_\Omega u^{p+1}v^{-r}{\rm d}x+C_3{\rm e}^{rt}, \end{align*}

其中

C_3:=\frac{r}{p+1}\left(\frac{2p}{p+1}\right)^p.

由于 u,vh_2 的非负性, 记 C_4:=C_2+C_3 则有

\begin{align*} \frac{\rm d}{{\rm d}t}\int_\Omega u^pv^{-r}{\rm d}x\leq&-r(r+1+p\chi_1) \int_\Omega u^pv^{-r-2} |\nabla v|^2{\rm d}x-p(p-1)\int_\Omega u^{p-2}v^{-r} |\nabla u|^2{\rm d}x\\ &+(2pr+\chi_1 p(p-1))\int_\Omega u^{p-1}v^{-r-1} \nabla v\cdot\nabla u{\rm d}x -\frac p2\int_\Omega u^{p}v^{-r}{\rm d}x+C_4{\rm e}^{rt}\\ =&\left\{\frac{(2pr+\chi_1p(p-1))^2}{4p(p-1)}-r(r+1+p\chi_1)\right\} \int_\Omega u^pv^{-r-2} |\nabla v|^2{\rm d}x\\ &-p(p-1)\int_\Omega\left| u^{\frac p2-1}v^{-\frac r2} \nabla u-\frac{2pr+\chi_1p(p-1)}{2p(p-1)}u^{\frac p2}v^{-\frac r2-1}\nabla v\right| {\rm d}x\\ &-\frac p2\int_\Omega u^{p}v^{-r}{\rm d}x+C_4{\rm e}^{rt}, \end{align*}

这进一步表明

\begin{align*} \frac{\rm d}{{\rm d}t}\int_\Omega u^pv^{-r}{\rm d}x\leq &-g(p,r) \int_\Omega u^pv^{-r-2} |\nabla v|^2{\rm d}x -\frac p2\int_\Omega u^{p}v^{-r}{\rm d}x+C_4{\rm e}^{rt}, \end{align*}

其中

g(p,r):=r(r+1+p\chi_1)-\frac{(2pr+\chi_1p(p-1))^2}{4p(p-1)}.

接下来, 我们还需证明 g(p,r) 的非负性. 事实上, 易知

\begin{align*} 4(p-1)g(p,r)=\,&4(p-1)r(r+1+p\chi_1)-4pr^2-4pr\chi_1(p-1)-\chi_1^2p(p-1)^2 \\ =&-4r^2+4(p-1)r-\chi_1^2p(p-1)^2\\ =&-4\left(r-\frac{p-1}{2}\right)^2+(p-1)^2(1-p\chi_1^2). \end{align*}

据此, 对任意的 r\in(r_{-}(p),r_{+}(p)), 都有 4(p-1)g(p,r)>0, 进一步可得

\begin{align*} \frac{\rm d}{{\rm d}t}\int_\Omega u^pv^{-r}{\rm d}x\leq & -\frac p2\int_\Omega u^{p}v^{-r}{\rm d}x+C_4{\rm e}^{rt},\quad t\in(0,T_{\max}). \end{align*}

由常微分方程比较原理易得 (2.12) 式.

3 一些关键的先验估计

为了探讨警察威慑对解性质的正则效应, 关键步骤是利用 (1.3) 式中非局部项建立一系的先验估计. 这里我们的出发点和创新点是追踪泛函 \int_\Omega v^\chi(\cdot,t){\rm d}x (\chi>1) 的时间演化估计.

引理3.1 \chi>1. 对任意的 t\in(0,T_{\max}) 存在独立于 t 的正常数 c_5 使得

\begin{matrix}\label{eq-uL1vLc} \int_\Omega v^\chi(\cdot,t){\rm d}x \leq c_5. \end{matrix}
(3.1)

基于 (1.3) 式中第二个方程, 由分部积分易得

\begin{align*} & \frac{\rm d}{{\rm d}t}\int_\Omega v^\chi{\rm d}x+\chi(\chi-1)\int_\Omega v^{\chi-2}|\nabla v|^2{\rm d}x+\chi \int_\Omega v^\chi{\rm d}x\\ &= \chi\int_\Omega u v^\chi{\rm d}x +\chi\int_\Omega v^{\chi-1}h_2{\rm d}x. \end{align*}

据此, 利用 Young 不等式可得

\begin{matrix}\label{eq-ul1-r1} \frac{\rm d}{{\rm d}t}\int_\Omega v^\chi{\rm d}x+\frac{4(\chi-1)}{\chi}\left\|\nabla v^{\frac{\chi}2}\right\|^2_{L^2}+ \int_\Omega v^\chi{\rm d}x\leq& \chi\int_\Omega u v^\chi{\rm d}x +\int_\Omega h_2^\chi{\rm d}x. \end{matrix}
(3.2)

再运用 Gagliardo-Nirenberg 不等式可知

\begin{split}\int_\Omega v^\chi{\rm d}x=\|v^{\frac{\chi}2}\|^2_{L^2}\le &C_1\|v^{\frac{\chi}2}\|^{2\theta}_{L^{\frac{2}{\chi}}}\|\nabla v^{\frac{\chi}2}\|^{2(1-\theta)}_{L^2}+C_1\|v^{\frac{\chi}2}\|^{2}_{L^{\frac{2}{\chi}}}\\\le&C_1\|v\|^{\chi\theta}_{L^1}\|\nabla v^{\frac{\chi}2}\|^{2(1-\theta)}_{L^2}+C_1\|v\|^{\chi}_{L^1},\end{split}

这里 \theta=\frac{1}{\frac{\chi n}2+1-\frac{n}2}, C_1>0 依赖于 \chi\Omega. 由于 (2.7) 式, 记 C_2:=C_1\max\{c_3^{\chi\theta},c_3^\chi\}, 其中 c_3 源至 (2.7) 式, 我们有

\int_\Omega v^\chi{\rm d}x\le C_2\|\nabla v^{\frac{\chi}2}\|^{2(1-\theta)}_{L^2}+C_2.
(3.3)

对 (3.3) 式应用 Young 不等式可知: 对来自 (2.6) 式的 c_2

C_3:=\frac{C_2c_2\chi}{\kappa}+\theta\left(\frac{\chi(1-\theta)}{4(\chi-1)}\right)^{\frac{1-\theta}{\theta}}\left(\frac{C_2c_2\chi}{\kappa}\right)^\frac{1}{\theta},

我们有

\frac{4(\chi-1)}{\chi}\left\|\nabla v^{\frac{\chi}2}\right\|^{2}_{L^2}+C_3\ge \frac{c_2\chi}{\kappa}\int_\Omega v^\chi{\rm d}x.
(3.4)

联立 (3.2) 式和记 g(t):=\chi\int_\Omega u v^\chi(\cdot,t){\rm d}x\left\{\int_\Omega v^\chi (\cdot,t){\rm d}x\right\}^{-1}, 则有

\begin{align*} &\frac{\rm d}{{\rm d}t}\int_\Omega v^\chi{\rm d}x+\Big(\frac{c_2\chi}{\kappa}+1\Big)\int_\Omega v^\chi{\rm d}x\\ \leq& g(t)\int_\Omega v^\chi{\rm d}x +\int_\Omega h_2^\chi{\rm d}x+C_3. \end{align*}

由上述微分不等式可得

\begin{align*} \int_\Omega v^\chi{\rm d}x\le&\Big(\int_\Omega v_0^\chi{\rm d}x\Big)\exp\left\{-\int_0^t\frac{c_2\chi}{\kappa}+1-g(s){\rm d}s\right\}\\ &+\int_0^t\exp\left\{-\int_s^t\frac{c_2\chi}{\kappa}+1-g(\sigma){\rm d}\sigma\right\}\Big(\int_\Omega h_2^\chi(\cdot,s){\rm d}x+C_3\Big){\rm d}s. \end{align*}

回顾 (2.8) 式, 对于任意的 t\in(0,T_{\max})s\in[t] 我们有

\begin{align*} &\int_0^tg(s){\rm d}s\leq \frac{c_2\chi}{\kappa}(1+t),\\ &\int_s^tg(\sigma){\rm d}\sigma\leq \frac{c_2\chi}{\kappa}(1+t-s), \end{align*}

进而可得

\begin{align*} \int_\Omega v^\chi{\rm d}x\le&\Big(\int_\Omega v_0^\chi{\rm d}x\Big)\exp\left\{\frac{c_2\chi}{\kappa}-t\right\}\\ &+\Big(\|h_2^\chi\|_{L^\infty((0, \infty; L^1(\Omega))}+C_3\Big)\exp\left\{\frac{c_2\chi}{\kappa}\right\}\Big(1-\exp\left\{-t\right\}\Big), \end{align*}

这表明了 (3.1) 式.

找到 p_1>\frac n2 并得到 \|u\|_{L^{p_1}} 的有界性是证明经典解全局存在的核心. 为了利用引理 2.2, 引理 2.3 和引理 3.1 建立所需有界性, 关键点是确定趋化系数 \chi_1 的范围, 其方法不同于文献 [47],[49] 中用于具线性信号产出的对数 Keller-Segel 模型的方法, 也不同于文献 [10] 中用于 (1.1) 式的方法, 还不同于文献 [23],[25] 中用于 (1.1) 式中 +uv+u 代替情形下的方法.

引理3.2 n\geq3, \chi>1, \chi_1>0 以及 \widetilde{\chi}_1 定义于 (5). 若 \chi_1<\widetilde{\chi}_1, 则存在 p\in(\frac n2,\frac{1}{\chi_1^2}) 使得

\begin{matrix}\label{eq-chi1} \sqrt{1- p\chi_1^2}>1-\frac{2\chi\left(\frac 2n p-1\right)}{ p-1}. \end{matrix}
(3.5)

我们将根据 (3.5) 式右手边的正负性分开证明. 如果右手边非正, 由于p>1, 则有

\begin{matrix}\label{eq-chii} 1-\frac{2\chi\left(\frac 2n p-1\right)}{ p-1}\le 0 \iff (4\chi-n)p\ge (2\chi-1)n. \end{matrix}
(3.6)

\chi>\max\{1,\frac n4\}, 则对任意的 \chi_1<\sqrt{\frac{4\chi-n}{(2\chi-1)n}} 都有 \frac{1}{\chi_1^2}>\frac{(2\chi-1)n}{4\chi-n}. 据此, 存在 p\in(\frac n2,\frac{1}{\chi_1^2}) (例如, p=\frac{(2\chi-1)n}{4\chi-n}) 使得 (3.5) 式成立.

另一方面, 我们考虑 (3.6) 式右手边为正的情形. 事实上, 要么当 n=3,4 时, 对任意的 \chi>1 都有 p\in (\frac n2, \frac{(2\chi-1)n}{4\chi-n}); 要么当 n\geq5 时, 对任意的 \chi>\frac n4 都有 p\in (\frac n2, \frac{(2\chi-1)n}{4\chi-n}); 要么当 n\ge5 时, 对任意的 \chi\in(1,\frac n4] 都有 p>\frac n2, 则有

\begin{align*} 1-\frac{2\chi\left(\frac 2n p-1\right)}{ p-1}>0. \end{align*}

在此情形下, (3.5) 式等价于

\begin{matrix}\label{eq-chi2} \chi_1^2<g(p):=\frac{4\chi(2p-n)[(n-2\chi)p+n(\chi-1)]}{n^2p(p-1)^2}. \end{matrix}
(3.7)

直接计算表明

g(n-1)=\frac{4}{n^2}\cdot\frac{\chi(n-\chi)}{n-1}.

如果 1<\chi<\frac n2, 则 \frac n2<n-1<\frac{(2\chi-1)n}{4\chi-n}. 因此, 我们能取 p=n-1 使得 (3.7) 式成立. 这进一步表明只要 \chi_1<\frac{2}{n}\cdot\sqrt{\frac{\chi(n-\chi)}{n-1}} 就有 (3.5) 式成立.

最后, 我们比较上述两个阀值 \frac{2}{n}\cdot\sqrt{\frac{\chi(n-\chi)}{n-1}}\sqrt{\frac{4\chi-n}{(2\chi-1)n}} 的大小, 就可以得到 (1.6) 中 \widetilde{\chi}_1 的定义并结束引理的证明.

对于 \chi>1 情形运用引理 3.2 且对 \chi\in(0,1] 情形运用文献 [10,定理 1] 可得任意 \chi>0 情形下 \|u\|_{L^{p_1}} (p_1>\frac n2) 的有界性.

引理3.3n\geq3, \chi_1<\widetilde{\chi}_1, \widetilde{\chi}_1 定义于 (1.6) 式. 则存在 p_1>\frac n2c_6=c_6(p_1)>0 使得

\begin{matrix}\label{eq-ulp1} \int_\Omega u^{p_1}(\cdot,t){\rm d}x\leq c_6{\rm e}^{c_6t},\quad t\in(0,T_{\max}). \end{matrix}
(3.8)

r_{\pm}(\cdot) 定义于 (2.11) 式. 根据引理 2.3, 如果 \chi_1<\sqrt{\frac 2n}, 则对任意 p\in(\frac n2,\frac{1}{\chi_1^2})r\in(r_{-}( p),r_{+}( p)) 都有

\begin{matrix} \int_\Omega u^{ p}v^{-r}(\cdot,t){\rm d}x\leq c_4{\rm e}^{rt},\quad t\in(0,T_{\max}),\label{eq-ulp10} \end{matrix}
(3.9)

其中 c_4 源于 (2.12) 式. 对 p_1\in(\frac n2, p), 利用 Hölder 不等式可得

\int_\Omega u^{p_1}{\rm d}x\leq \left\{\int_\Omega u^{ p}v^{-r}{\rm d}x\right\}^{\frac{p_1}{ p}} \left\{\int_\Omega v^{\frac{p_1r}{ p-p_1}}{\rm d}x\right\}^{\frac{ p-p_1}{ p}},

联立 (3.9) 式, 确保存在 C_1:=\frac{p_1r}{p}C_2:=c_4^{\frac{p_1}{p}} 使得

\begin{matrix} \|u(\cdot,t)\|_{L^{p_1}}\leq C_2{\rm e}^{C_1t}\|v\|_{L^\frac{p_1r}{ p-p_1}}^{\frac{p_1r}{ p}},\quad t\in(0,T_{\max}).\label{eq-ulp11} \end{matrix}
(3.10)

如果 \chi\leq1, 据 (2.7) 式, 只要 \chi_1<\frac 2n 就可根据文献 [10,定理 1] 的证明找到 p,p_1,r 使得

\begin{align*} \frac{p_1r}{ p-p_1}\leq1. \end{align*}

这进一步表明对任意的 \chi\leq1 都有 (3.8) 式成立.

对于 \chi>1 的情形, 我们已经在 (3.1) 式中建立了 \|v\|_{L^\chi} 的有界性, 现只需验证在相应参数范围内能得到 \|v\|_{L^\frac{p_1r}{ p-p_1}} 的有界性, 即, 能找到 p,p_1,r 使得

\begin{matrix} \frac{p_1r}{ p-p_1}\leq\chi.\label{eq-ulp12} \end{matrix}

事实上, (3.11) 式等价于 p_1\leq \frac{ p\chi}{r+\chi}, 据此我们仅需找到 pr 使得 \frac n2<\frac{ p\chi}{r+\chi}. 令 r 充分靠近 r_{-}( p), 只需验证 \frac n2<\frac{ p\chi}{r_{-}( p)+\chi} 成立. 因此, 根据 (2.11) 式我们只需验证如下不等式

\begin{align*} \chi\left(\frac2n p-1\right)>\frac{ p-1}{2}\left(1-\sqrt{1- p\chi_1^2}\right), \end{align*}

即验证

\begin{matrix} \sqrt{1- p\chi_1^2}>1-\frac{2\chi\left(\frac 2n p-1\right)}{ p-1}.\label{eq-ulp13} \end{matrix}

回顾引理 3.2, 对任意的 \chi>1 我们能找到 p\in(\frac n2,\frac{1}{\chi_1^2}) (\chi_1<\widetilde{\chi}_1) 使得 (3.12) 式成立, 这里 \widetilde{\chi}_1 定义于 (1.6) 式. 这确保了 (3.11) 式, 再联立 (3.10) 和 (3.1) 式, 可知对任意 p_1\in(\frac n2, p)

\begin{align*} \|u(\cdot,t)\|_{L^{p_1}}\leq C_3{\rm e}^{C_1t},\quad t\in(0,T_{\max}), \end{align*}

其中 C_3:=C_2c_5^{\frac{p_1r}{p}}|\Omega|^{\frac{(p-p_1)\chi-p_1r}{p_1r\chi}}, c_5 确定于 (3.1) 式. 我们因此得到了 (3.8) 式.

4 全局经典解

根据引理 3.3, 利用适当的迭代策略我们能得到更高的正则性估计, 进一步可以证明定理 1.1. 据 Neumann 热半群性质可知有界性 (3.8) 式可确保 \|\nabla v\|_{L^{q_1}} (q_1>n) 的有界性.

引理4.1 设引理 3.3 中条件都成立. 对任意的 \chi>0 存在 q_1>nc_7>0 使得

\begin{matrix} \|\nabla v(\cdot,t)\|_{L^{q_1}}\leq c_7{\rm e}^{c_7t},\quad t\in(0,T_{\max}).\label{eq-vwlq} \end{matrix}
(4.1)

q_1=2p_2, p_2\in(\frac n2,p_1), p_1 来自 (3.8) 式. 据 Neumann 热半群性质 (参见文献 [6,引理 2.1], 文献 [46,引理 1.3]) 及

\begin{matrix}\label{eq-wcov} v(\cdot, t) &={\rm e}^{t(\Delta-1)} v_{0}+\int_{0}^{t} {\rm e}^{(t-s)(\Delta-1)} (uv + h_2){\rm d}s, \end{matrix}
(4.2)

我们能找到 C_1>0 使得

\begin{align*} \|\nabla v(\cdot,t)\|_{L^{q_1}}\leq & \|\nabla {\rm e}^{t(\Delta-1)}v_0\|_{L^{q_1}}+\int_0^t\left\|\nabla {\rm e}^{(t-s)(\Delta-1)}\left(uv+h_2\right)\right\|_{L^{q_1}}{\rm d}s\\ \le&C_1\|\nabla v_0\|_{L^{q_1}}+C_1\int_0^t\left(1+(t-s)^{-\frac12-\frac n2(\frac{1}{p_2}-\frac{1}{q_1})}\right){\rm e}^{-(t-s)}\left\|uv +h_2\right\|_{L^{p_2}}{\rm d}s. \end{align*}

运用 Hölder 不等式可知

\left\|uv +h_2\right\|_{L^{p_2}}\leq \|u\|_{L^{p_1}}\|v\|_{L^{\frac{p_1p_2}{p_1-p_2}}}+\|h_2\|_{L^{p_2}},

联合 (1.5) 和 (3.8) 式, 我们能找到 C_2:=C_2(h_2,c_6)>0 满足

\left\|uv +h_2\right\|_{L^{p_2}}\leq C_2{\rm e}^{C_2t}\|v\|_{L^{\frac{p_1p_2}{p_1-p_2}}}+C_2.

\frac {n}{2}(\frac{1}{p_2}-\frac{1}{q_1})<\frac12, 存在 C_3:=C_3(C_1, v_0)>0 使得

\begin{align*} \|\nabla v(\cdot,t)\|_{L^{q_1}}\leq &C_3+C_2{\rm e}^{C_2t}\int_0^t\left(1+(t-s)^{-\frac12-\frac n2(\frac{1}{p_2}-\frac{1}{q_1})}\right){\rm e}^{-(t-s)}\|v\|_{L^{\frac{p_1p_2}{p_1-p_2}}}{\rm d}s. \end{align*}

利用 Gagliardo-Nirenberg 不等式, 存在 C_4>0 服从

\|v\|_{L^{\frac{p_1p_2}{p_1-p_2}}}\leq C_4\left(\|v\|_{L^1}^{1-\varrho}\|\nabla v\|_{L^{q_1}}^\varrho+\|v\|_{L^1}\right),

其中

0<\varrho:=\frac{nq_1(p_1p_2-p_1+p_2)}{p_1p_2(nq_1-n+q_1)}<\frac{2n(p_1p_2-p_1+p_2)}{2np_1p_2} <1.

联立 (2.7) 式, 可找到 C_5>0 使得

\|v\|_{L^{\frac{p_1p_2}{p_1-p_2}}}\leq C_5\left(\|\nabla v\|_{L^{q_1}}^\varrho+1\right),

据此, 存在 C_i>0, i=6,7,8, 使得

\begin{align*} \|\nabla v(\cdot,t)\|_{L^{q_1}}\leq &C_6+C_7{\rm e}^{C_2t}+C_8{\rm e}^{C_2t}\int_0^t\left(1+(t-s)^{-\frac12-\frac n2(\frac{1}{p_2}-\frac{1}{q_1})}\right){\rm e}^{-(t-s)}\|\nabla v\|_{L^{q_1}}^\varrho {\rm d}s. \end{align*}

K(T):=\sup_{t\in(0, T)}\|\nabla v(\cdot, t)\|_{L^{q_1}}, 对任意的 T\in (0, T_{\max}) 存在 C_9>0 满足

K(T)\leq C_6+C_7{\rm e}^{C_2t}+C_9{\rm e}^{C_2t}K^\varrho(T).

因为 0<\varrho<1, 利用 Young 不等式可知: 存在 C_{10}>0 使得

K(T)\le C_{10}(1+{\rm e}^{C_2t}).

这已足够确保 (4.1) 式.

接下来, 我们将建立 \|u\|_{L^\infty} 的有界性.

引理4.2 设引理 4.1 中条件都成立. 则存在 c_8>0 使得

\begin{matrix} \|u(\cdot,t)\|_{L^\infty}\leq c_8{\rm e}^{c_8t},\quad t\in(0,T_{\max}).\label{eq-ulinf} \end{matrix}
(4.3)

类似于引理 4.1 的证明, 利用最值原理以及 uv\kappa\{\int_\Omega v^\chi{\rm d}x\}^{-1}uv^{\chi} 的非负性, 对任意的 q_2\in(n,q_1) (q_1 源至 (4.1) 式) 存在 C_1>0 使得

\begin{align*} \|u (\cdot,t)\|_{L^\infty} &\leq \|{\rm e}^{t(\Delta-1)}u_0\|_{L^\infty}+\chi_1\int_0^t\left\|{\rm e}^{(t-s)(\Delta-1)}\left\{\nabla\cdot\left(u \nabla\ln v \right)+h_1\right\}\right\|_{L^\infty}{\rm d}s\\ &\le \|u_0\|_{L^\infty}+C_1\int_0^t\left(1+(t-s)^{-\frac12-\frac{n}{2q_2}}\right){\rm e}^{-(t-s)}\|u \nabla\ln v \|_{L^{q_2}}{\rm d}s\\ &\quad +C_1\int_0^t\left(1+(t-s)^{-\frac{n}{2q_2}}\right){\rm e}^{-(t-s)}\|h_1\|_{L^{q_2}}{\rm d}s. \end{align*}

运用 Hölder 不等式, 我们有

\left\|u \nabla \ln v \right\|_{L^{q_2}}\leq \|u\|_{L^{\frac{q_1q_2}{q_1-q_2}}}\|\nabla v\|_{L^{q_1}}\|v^{-1}\|_{L^\infty},

联立 (4.1) 和 (2.5) 式, 进一步可证: 存在 C_2>0 使得

\left\|u \nabla \ln v \right\|_{L^{q_2}}\leq C_2{\rm e}^{C_2t}\|u\|_{L^{\frac{q_1q_2}{q_1-q_2}}}.

据此, 我们能从 (1.4) 式推得: 存在 C_3>0 使得

\begin{align*} \|u (\cdot,t)\|_{L^\infty} &\le C_3+C_3{\rm e}^{C_2t}\int_0^t\left(1+(t-s)^{-\frac12-\frac{n}{2q_2}}\right){\rm e}^{-(t-s)}\|u\|_{L^{\frac{q_1q_2}{q_1-q_2}}}{\rm d}s. \end{align*}

利用插值不等式和 (2.6) 式, 则有

\|u\|_{L^{\frac{q_1q_2}{q_1-q_2}}}\leq\|u\|_{L^1}^{1-\sigma}\|u\|_{L^\infty}^\sigma\leq C_4\|u\|_{L^\infty}^\sigma,

其中 \sigma:=1-\frac{q_1-q_2}{q_1q_2}C_4>0. 综上, 存在 C_5>0 使得

\begin{align*} \|u (\cdot,t)\|_{L^\infty} &\le C_3+C_5{\rm e}^{C_2t}\int_0^t\left(1+(t-s)^{-\frac12-\frac{n}{2q_2}}\right){\rm e}^{-(t-s)}\|u\|_{L^\infty}^\sigma {\rm d}s. \end{align*}

G(T):=\sup_{t\in(0, T)}\|u(\cdot, t)\|_{L^\infty}, 对任意的 T\in (0, T_{\max}) 存在 C_6>0 使得

G(T)\leq C_3+C_6{\rm e}^{C_2t}G^\sigma(T),

这里利用了 \int_0^t\left(1+(t-s)^{-\frac12-\frac{n}{2q_2}}\right){\rm e}^{-(t-s)}{\rm d}s<\inftyq_2>n 的事实. 由于 0<\sigma<1, 利用 Young 不等式可知存在 C_7>0 使得

G(T)\leq C_7(1+{\rm e}^{C_2t}),

这进一步表明了期望的 (4.3) 式.

基于已建立的先验估计, 我们可将 (4.1) 式中的有界性扩展到任意的 r>n.

引理4.3 设引理 4.2 中条件都成立. 则对任意的 r>n 存在 c_9>0 使得

\begin{matrix} \|\nabla v(\cdot,t)\|_{L^{r}}\leq c_9{\rm e}^{c_9t},\quad t\in(0,T_{\max}).\label{eq-0vwlq2} \end{matrix}
(4.4)

q_1 取自 (4.1) 式, 类似于引理 4.1 的证明, 对任意的 r\geq q_1 再次利用 (4.2) 式, 可知存在 C_1>0 使得

\begin{align*} \|\nabla v(\cdot,t)\|_{L^{r}}\leq &C_1\|\nabla v_0\|_{L^{r}}+C_1\int_0^t\left(1+(t-s)^{-\frac12-\frac{n}{2q}}\right){\rm e}^{-(t-s)}\left\|uv +h_2\right\|_{L^{\infty}}{\rm d}s. \end{align*}

基于 Hölder 不等式, (1.4) 和 (4.3) 式, 存在 C_2>0 满足

\begin{align*} \left\|uv +h_2\right\|_{L^{\infty}}\leq&\|u\|_{L^\infty}\|v\|_{L^\infty}+\|h_2\|_{L^\infty}\\ \leq&C_2\left({\rm e}^{C_2t}\|v\|_{L^\infty}+1\right), \end{align*}

联立 Gagliardo-Nirenberg 不等式, (2.7) 和 (4.1) 式, 我们能找到 C_3>0 使得

\begin{align*} \left\|uv +h_2\right\|_{L^{\infty}}\leq &C_2\left\{{\rm e}^{C_2t}(\|v\|_{L^1}^{1-\sigma}\|\nabla v\|_{L^{q_1}}^\sigma+\|v\|_{L^1})+1\right\}\\ \leq&C_3\left(1+{\rm e}^{C_3t}\right), \end{align*}

这里由于 q_1>n0<\sigma:=\frac{q_1n}{q_1n+q_1-n}<1. 进一步, 存在 C_4,C_5>0 使得

\begin{align*} \|\nabla v(\cdot,t)\|_{L^{r}}\leq &C_4+C_5(1+{\rm e}^{C_3t})\int_0^t\left(1+(t-s)^{-\frac12-\frac{n}{2r}}\right){\rm e}^{-(t-s)}{\rm d}s. \end{align*}

注意 r\ge q_1>n, 故对任意的 r\ge q_1 不等式 (4.4) 都成立. 至于 r\in(n,q_1), 利用 (4.1) 式和 Hölder 不等式则可知结论仍成立.

基于以上的准备, 我们现将完成定理 1.1 的证明.

定理 1.1 的证明 对任意的 \chi>0, 基于延拓准则 (2.1), 联立估计 (2.5), (4.3) 和 (4.4) 式则可得初边值问题 (1.3) 经典解的全局存在性, 据此我们完成了定理 1.1 的证明.

5 长时间行为

本节将致力于经典解的有界性和长时间行为. 为此, 我们首先建立如下有界性.

引理5.1 设定理 1.1 的条件和 (1.7) 式都成立. 对任意的 r>n 存在独立于 tC>0 使得

\begin{matrix}\label{eq-0bedd} \|u(\cdot,t)\|_{L^\infty}+\|v^{-1}(\cdot,t)\|_{L^\infty}+\| v(\cdot,t)\|_{W^{1,r}}\leq C,\quad t>0. \end{matrix}
(5.1)

注5.1 条件 (1.7) 被用于获取解成分 v 的点态正下界. 利用 \Omega 的凸性和文献 [51,引理 2.3] (或者文献 [15,引理 3.1]), 通过修正文献 [23,推论 3.1] (或者文献 [24,引理 4.1]) 的证明, 条件 (1.7) 可以被放宽到 \inf_{t>0}\int_\Omega h_2(x,t){\rm d}x>0.

首先建立 v 的点态正下界. 事实上, 设 \widehat{c}_1:=\{\inf _{x \in \overline{\Omega}} v_{0}(x),\delta\} (\delta 取自 (1.7) 式), 对 (1.3) 式中第二个方程运用最值原理, 可知 \widehat{c}_1 是如下初边值问题

\begin{matrix}\label{eq-Vlower} v(x,t)\geq \widehat{c}_1, \quad x\in\Omega,\,\, t>0 \end{matrix}
(5.2)

的一个下解. 基于 (5.2) 式, 类似于引理 2.3, 对任意的 p\in(1,\frac{1}{\chi_1^2})r\in(r_{-}(p),r_{+}(p)) 存在 \widehat{c}_2>0 使得

\begin{matrix} \int_\Omega u^pv^{-r}(\cdot,t){\rm d}x\leq \widehat{c}_2,\quad t\in(0,\infty),\label{eq-0Eulnu} \end{matrix}
(5.3)

其中 r_{\pm}(p) 定义于 (2.11) 式. 基于 (5.3) 式, 类似于引理 3.3 的证明, 对 n\ge3\widetilde{\chi}_1 取自于 (1.6) 式使得 \chi_1<\widetilde{\chi}_1, 则存在 p_1\in(\frac n2,\frac{1}{\chi_1^2})\widehat{c}_3>0 使得

\begin{matrix} \int_\Omega u^{p_1}(\cdot,t){\rm d}x\leq \widehat{c}_3,\quad t\in(0,\infty).\label{eq-0ulp1} \end{matrix}
(5.4)

据此, 适当地修正引理 4.1, 引理 4.2 和引理 4.3 的证明可以得到

\begin{matrix} \|u(\cdot,t)\|_{L^\infty}+\|\nabla v(\cdot,t)\|_{L^{r}}\leq \widehat{c}_4,\quad t\in(0,\infty), \label{eq-vwlq2} \end{matrix}
(5.5)

这里 r>n\widehat{c}_4>0. 联立 (5.2), (2.7) 式以及 Gagliardo-Nirenberg 不等式表明了期望的 (5.1) 式.

为了探讨解的长时间行为, 我们首先关注椭圆问题 (6), 其可解性可直接据文献 [19] 可得.

引理5.2 对给定的 h_{2,\infty}\in \mathcal{C}^1(\overline{\Omega}), 椭圆问题 (1.12) 拥有唯一经典解 v_\infty, 且存在 \theta\in(0,1) 使其属于\mathcal{C}^{2+\theta}(\overline{\Omega}).

作为准备, 我们也需要如下源至微分不等式的衰减估计, 其证明参见文献 [9,引理 4.6] (或文献 [23,引理 2.5]).

引理5.3

0\leq y\in \mathcal{C}^1([0,\infty)).

若存在 0\leq g(t)\in \mathcal{C}([0, \infty))\cap L^\infty([0,\infty)) 满足

\begin{align*} \lim_{t\rightarrow \infty}\int_t^{t+1}g (s){\rm d}s=0,\quad t>0 \end{align*}

以及对某个 \lambda>0 服从

y'(t)+\lambda y (t)\leq g (t),\quad t>0,

则当 t\rightarrow \infty 时有

y(t)\rightarrow0.

引理 5.3 的一个直接的推论是, 可得 \|u\|_{L^1} 的衰减估计.

引理5.4 设引理 5.1 中的条件都成立. 进一步假设 (1.8) 式成立, 则当 t\rightarrow \infty 时有

\int_\Omega u(\cdot,t){\rm d}x\rightarrow0,
(5.6)
\int_t^{t+1}\int_\Omega u v(\cdot,s) {\rm d}x{\rm d}s \rightarrow0.
(5.7)

由 (2.9), (5.2) 式以及 uv 非负性, 可得

\begin{matrix}\label{eq-deL1u} \frac{\rm d}{{\rm d}t}\int_\Omega u{\rm d}x+\int_\Omega u{\rm d}x+\int_\Omega uv{\rm d}x \leq\int_\Omega h_1{\rm d}x, \end{matrix}
(5.8)

联立 (1.8) 式和引理 5.3, 可知 (5.6) 式成立.

把 (5.8) 式关于时间变量在 [t,t+1] 上积分可得

\begin{align*} \int_\Omega u(\cdot,t+1){\rm d}x+ \int_t^{t+1}\int_\Omega u v {\rm d}x{\rm d}s\leq\int_\Omega u(\cdot,t){\rm d}x+\int_t^{t+1}\int_\Omega h_1{\rm d}x{\rm d}s,\quad t>0. \end{align*}

再次利用 (1.8) 和 (5.6) 式, 知期望的 (5.7) 式成立.

再次利用引理 5.3, (5.6) 和 (5.7) 式, 可得 v-v_\inftyL^2 衰减估计, 其中 v_\infty 是椭圆问题 (1.12) 的唯一经典解. 为了方便, 记 \widehat{v}:=v-v_\infty. 由定理 1.1 和引理 5.2, 下列初边值问题

\left\{\begin{split}&\widehat{v}_{t}=\Delta\widehat{v}-\widehat{v}+uv+h_2-h_{2,\infty},&x\in\Omega,\,\,t>0,\\&\nabla \widehat{v}\cdot\nu=0,&x\in\partial\Omega,\,\,t>0,\\&\widehat{v}(x,0)=v_0(x)-v_\infty(x),&x\in\Omega\end{split}\right.
(5.9)

拥有唯一经典解 \widehat{v}.

引理5.5 设引理 5.4 中条件都成立, 且 \widehat{v} 是 (5.9) 式的唯一经典解. 假设 (1.9) 式成立, 则当 t\rightarrow\infty 时有

\begin{matrix}\label{eq-ul1e0} \|\widehat{v}(\cdot,t)\|_{L^2} \rightarrow0. \end{matrix}
(5.10)

把 (5.9) 式中第一个方程乘以 \widehat{v}, 由分部积分可知

\begin{align*} \frac12\frac{\rm d}{{\rm d}t}\int_\Omega\widehat{v}^2{\rm d}x+\int_\Omega|\nabla\widehat{v}|^2{\rm d}x+\int_\Omega\widehat{v}^2{\rm d}x=\int_\Omega uv \widehat{v}{\rm d}x+\int_\Omega \widehat{v}(h_2-h_{2,\infty}){\rm d}x, \end{align*}

进一步结合 Hölder 不等式和 uv 的非负性, 可得

\begin{align*} \frac12\frac{\rm d}{{\rm d}t}\int_\Omega\widehat{v}^2{\rm d}x+\int_\Omega|\nabla\widehat{v}|^2{\rm d}x+\int_\Omega\widehat{v}^2{\rm d}x\leq\|\widehat{v} \|_{L^\infty}\left\{\int_\Omega uv{\rm d}x+\int_\Omega |h_2-h_{2,\infty}|{\rm d}x\right\}. \end{align*}

利用 Sobolev 不等式, 可知对任意的 r>n 存在 C_1>0 使得

\|\widehat{v} \|_{L^\infty}\leq C_1\|\widehat{v} \|_{W^{1,r}},

联立 (5.1) 式和引理 5.2, 可找到 C_2>0 使得

\|\widehat{v} \|_{L^\infty}\leq C_2,\quad t>0.

据此, 有

\begin{align*} \frac12\frac{\rm d}{{\rm d}t}\int_\Omega\widehat{v}^2{\rm d}x+\int_\Omega|\nabla\widehat{v}|^2{\rm d}x+\int_\Omega\widehat{v}^2{\rm d}x\leq C_2\left\{\int_\Omega uv{\rm d}x+\int_\Omega |h_2-h_{2,\infty}|{\rm d}x\right\}. \end{align*}

注意到 (1.9) 和 (5.7) 式可表明: 当 t\rightarrow\infty 时, 则有

C_2\int_t^{t+1}\left\{\int_\Omega uv{\rm d}x+\int_\Omega |h_2-h_{2,\infty}|{\rm d}x\right\}{\rm d}s\rightarrow0.

再次运用引理 5.3, 则可得期望的 (5.10) 式.

基于引理 5.1, 引理 5.4 以及引理 5.5, 我们能进一步发展文献 [12],[30] 中的方法并用此研究解的长时间行为. 我们首先建立 \|u\|_{L^\infty} 的衰减估计.

引理5.6 设引理 5.5 中条件成立. 则当 t\rightarrow\infty 时都有

\begin{matrix}\label{eq-uLide} \|u(\cdot,t)\|_{L^\infty} \rightarrow0. \end{matrix}
(5.11)

基于 u 的常数变易公式, uv\kappa\{\int_\Omega v^\chi{\rm d}x\}^{-1}uv^{\chi} 的非负性以及最值原理, 对于任意的 t>0 我们有

\begin{align*} \|u(\cdot,t)\|_{L^\infty} &\leq\|{\rm e}^{(\Delta-1) t}u_0\|_{L^\infty}+\chi_1\int_{0}^t\left\|{\rm e}^{(\Delta-1)(t-s)} \left\{\nabla\cdot\left( u \nabla \ln v\right)+h_1\right\}\right\|_{L^\infty}{\rm d}s, \end{align*}

再利用 Neumann 热半群的性质可知, 存在 \lambda>0, r>n 以及 C_1>0 使得

\begin{align*} \|u(\cdot,t)\|_{L^\infty} &\leq C_1\big(1+t^{-\frac{n}{2r}}\big){\rm e}^{-\lambda t}\|u_0\|_{L^r}+C_1\int_{0}^t\big(1+(t-s)^{-\frac12-\frac{n}{2r}}\big){\rm e}^{-\lambda(t-s)}\left\|u\nabla\ln v\right\|_{L^r}{\rm d}s\\ &\quad +C_1\int_{0}^t\big(1+(t-s)^{-\frac{n}{2r}}\big){\rm e}^{-\lambda(t-s)}\left\|h_1\right\|_{L^r}{\rm d}s\\ &=:V_1+V_2+V_3. \end{align*}

注意, 当 t\rightarrow\infty 时, 易有

\begin{matrix}\label{eq-deuV1} V_1\rightarrow0. \end{matrix}
(5.12)

对于 V_2, 利用 Hölder 不等式和 (5.1) 式, 存在 C_2>0 满足

\begin{align*} \|u\nabla\ln v\|_{L^r}\leq \|v^{-1}\|_{L^\infty}\|u\|_{L^{2r}}\|\nabla v\|_{L^{2r}}\leq C_2\|u\|_{L^{2r}},\quad t>0, \end{align*}

再结合插值不等式和 (5.1) 式, 可知存在 C_3>0 使得

\begin{align*} \|u\nabla\ln v\|_{L^r}\leq C_2\|u\|_{L^1}^{\frac{1}{2r}}\|u\|_{L^\infty}^{\frac{2r-1}{2r}}\leq C_3\|u\|_{L^1}^{\frac{1}{2r}},\quad t>0. \end{align*}

因此, 我们有

\begin{align*} V_2&\leq C_3\int_{0}^\frac t2\big(1+(t-s)^{-\frac12-\frac{n}{2r}}\big){\rm e}^{-\lambda(t-s)}\|u\|_{L^1}^{\frac{1}{2r}}{\rm d}s+C_3\int_{\frac t2}^t\big(1+(t-s)^{-\frac12-\frac{n}{2r}}\big){\rm e}^{-\lambda(t-s)}\|u\|_{L^1}^{\frac{1}{2r}}{\rm d}s\\ &=: V_{21}+ V_{22}. \end{align*}

对于 V_{21}, 由 (2.6) 式知存在 C_4>0 使得

\begin{align*} V_{21}&=C_3\int_{\frac t2}^{t}\big(1+s^{-\frac12-\frac{n}{2r}}\big){\rm e}^{-\lambda s}\|u(\cdot,t-s)\|_{L^1}^{\frac{1}{2r}}{\rm d}s\\ &\leq C_4\int_{\frac t2}^{t}\big(1+s^{-\frac12-\frac{n}{2r}}\big){\rm e}^{-\lambda s}{\rm d}s, \end{align*}

再结合如下不等式

\int_0^\infty \big(1+s^{-\frac12-\frac{n}{2r}}\big){\rm e}^{-\lambda s}{\rm d}s\leq C_5

知存在 C_5>0 使得当 t\rightarrow\infty 时, 则有

\begin{matrix}\label{eq-deuV21} V_{21}\rightarrow0. \end{matrix}
(5.13)

现在我们估计 V_{22}. 事实上, 不难发现 C_6>0 使得

\begin{align*} V_{22}\leq&C_3\sup\limits_{s>\frac t2}\|u(\cdot,s)\|_{L^1}^{\frac{1}{2r}}\int_{\frac t2}^t\big(1+(t-s)^{-\frac12-\frac{n}{2r}}\big){\rm e}^{-\lambda(t-s)}{\rm d}s\\ \le& C_3\sup\limits_{s>\frac t2}\|u(\cdot,s)\|_{L^1}^{\frac{1}{2r}}\int_{0}^\infty\big(1+s^{-\frac12-\frac{n}{2r}}\big){\rm e}^{-\lambda s}{\rm d}s\\ \leq&C_6\sup\limits_{s>\frac t2}\|u(\cdot,s)\|_{L^1}^{\frac{1}{2r}},\quad t>0. \end{align*}

这联立 (5.6) 式表明, 当 t\rightarrow\infty 时, 则有 V_{22}\rightarrow0, 据此, 进一步结合 (5.13) 式可得

\begin{matrix}\label{eq-deuV2} V_{2}\rightarrow0. \end{matrix}
(5.14)

对于 V_3, 我们类似可得

\begin{align*} V_3&\leq C_1\int_{0}^\frac t2\big(1+(t-s)^{-\frac{n}{2r}}\big){\rm e}^{-\lambda(t-s)}\|h_1\|_{L^r}{\rm d}s+C_1\int_{\frac t2}^t\big(1+(t-s)^{-\frac{n}{2r}}\big){\rm e}^{-\lambda(t-s)}\|h_1\|_{L^r}{\rm d}s\\ &=: V_{31}+ V_{32}. \end{align*}

对于 V_{31}, 由 (1.4) 式知存在 C_7>0 使得

\begin{align*} V_{31}= &C_1\int_{\frac t2}^t\big(1+s^{-\frac{n}{2r}}\big){\rm e}^{-\lambda s}\|h_1(\cdot,t-s)\|_{L^r}{\rm d}s\\ \leq&C_7\int_{\frac t2}^t\big(1+s^{-\frac{n}{2r}}\big){\rm e}^{-\lambda s}{\rm d}s, \end{align*}

故当 t\rightarrow\infty 时, 则有

\begin{matrix}\label{eq-deuV31} V_{31} \rightarrow0. \end{matrix}
(5.15)

对于 V_{32}, 利用 (1.4) 式和分部积分, 我们能找到 C_9>0 使得

\begin{align*} \|h_1\|_{L^r}\leq \|h_1\|_{L^1}^{\frac1r}\|h_1\|_{L^\infty}^{1-\frac1r}\leq C_9\|h_1\|_{L^1}^{\frac1r}, \end{align*}

据此可得

\begin{align*} V_{32}\leq &C_9\int_{\frac t2}^t\big(1+(t-s)^{-\frac{n}{2r}}\big){\rm e}^{-\lambda(t-s)}\|h_1\|_{L^1}^{\frac1r}{\rm d}s. \end{align*}

利用 Hölder 不等式, 可知对 n\geq3r>n 存在 C_{10}>0 使得

\begin{align*} V_{32}\leq &C_9\left\{\int_{\frac t2}^t\big(1+(t-s)^{-\frac{n}{2(r-1)}}\big){\rm e}^{-\frac{r}{r-1}\lambda(t-s)}{\rm d}s\right\}^\frac{r-1}{r}\left\{\int_{\frac t2}^t\left\|h_1(\cdot,s)\right\|_{L^1}{\rm d}s\right\}^\frac{1}{r}\\ \leq&C_{10}\left\{\int_{\frac t2}^t\left\|h_1(\cdot,s)\right\|_{L^1}{\rm d}s\right\}^\frac{1}{r},\quad t>0, \end{align*}

这里我们也利用了 \frac{n}{2(r-1)}<1 的事实. 再结合 (1.8) 式, 知当 t\rightarrow\infty 时则有 V_{32}\rightarrow0, 联立 (5.15) 式表明 V_{3}\rightarrow0. 最后, 再联立 (5.12) 和 (5.14) 式, 可得 (5.11) 式.

引理 5.6 的一个直接推论是, 可得期望的 \widehat{v} 的衰减估计.

引理5.7 设引理 5.6 中条件以及 (1.9) 式都成立, 则当 t\rightarrow\infty 时, 对任意的 r>n 都有

\begin{matrix}\label{eq-wLide} \|\widehat{v}(\cdot,t)\|_{W^{1,r}} \rightarrow0. \end{matrix}
(5.16)

基于 (5.9) 式, Neumann 热半群性质表明, 对任意的 r>n 存在 C_1>0 使得

\begin{align*} \|\nabla \widehat{v}(\cdot,t)\|_{L^{r}}&\leq \|\nabla {\rm e}^{t(\Delta-1)}\widehat{v}_0\|_{L^{r}}+\int_0^t\left\|\nabla {\rm e}^{(t-s)(\Delta-1)}\left(uv+h_2-h_{2,\infty}\right)\right\|_{L^r}{\rm d}s\\ &\le C_1(1+t^{-\frac12}){\rm e}^{- t}\|\widehat{v}_0\|_{L^{r}} \\ &\quad +C_1\int_0^t\left(1+(t-s)^{-\frac12}\right){\rm e}^{- (t-s)}\left\|uv +h_2-h_{2,\infty}\right\|_{L^r}{\rm d}s\\ &=:\Pi_1+\Pi_2. \end{align*}

显然地, 当 t\rightarrow\infty 时, 有

\begin{matrix}\label{eq-deuP1} \Pi_1\rightarrow0. \end{matrix}
(5.17)

另一方面, Hölder 不等式结合 (5.1) 式, 表明存在 C_2>0 使得

\|uv\|_{L^r}\leq \|u\|_{L^\infty}\|v\|_{L^r}\leq C_2 \|u\|_{L^\infty};

插值不等式和 (1.4) 式确保存在 C_3>0 使得

\begin{align*} \|h_2-h_{2,\infty}\|_{L^r}\leq \|h_2-h_{2,\infty}\|_{L^1}^{\frac1r}\|h_2-h_{2,\infty}\|_{L^\infty}^{1-\frac1r}\leq C_3\|h_2-h_{2,\infty}\|_{L^1}^{\frac1r}. \end{align*}

进一步联立 Hölder 不等式表明, 存在 C_4>0 使得

\begin{align*} \Pi_2\leq &C_4\int_0^{\frac{ t} 2}\left(1+(t-s)^{-\frac12}\right){\rm e}^{- (t-s)}\left(\|u\|_{L^\infty}+\left\|h_2-h_{2,\infty}\right\|_{L^1}^\frac1r\right){\rm d}s\\ &+C_4\int_{\frac t2}^t\left(1+(t-s)^{-\frac12}\right){\rm e}^{-(t-s)}\left(\|u\|_{L^\infty}+\left\|h_2-h_{2,\infty}\right\|_{L^1}^\frac1r\right){\rm d}s\\ \leq&C_4\int_{\frac t2}^t\left(1+s^{-\frac12}\right){\rm e}^{-s}\left(\|u(\cdot, t-s)\|_{L^\infty}+\left\|(h_2-h_{2,\infty})(\cdot, t-s)\right\|_{L^1}^\frac1r\right){\rm d}s\\ &+C_4\sup_{s>\frac{t}2}\|u(\cdot, s)\|_{L^\infty}\int_{\frac t2}^t\left(1+(t-s)^{-\frac12}\right){\rm e}^{-(t-s)}{\rm d}s\\ &+C_4\left\{\int_\frac t2^t\left((1+(t-s)^{-\frac12})^{\frac{r}{(r-1)}}\right){\rm e}^{- \frac{r(t-s)}{r-1}}{\rm d}s\right\}^{\frac{r-1}{r}}\left\{\int_\frac t2^t\|h_2-h_{2,\infty}\|_{L^1}{\rm d}s\right\}^\frac1r. \end{align*}

由于 (1.4) 式, (5.1) 式且 \frac{r}{2(r-1)}<1, 我们能找到 C_5>0 使得

\begin{align*} \Pi_2\leq&C_5\int_{\frac t2}^t\left(1+s^{-\frac12}\right){\rm e}^{-s}{\rm d}s +C_5\sup_{s>\frac{t}2}\|u(\cdot, s)\|_{L^\infty}+C_5\left\{\int_\frac t2^t\|h_2-h_{2,\infty}\|_{L^1}{\rm d}s\right\}^\frac1r, \end{align*}

再利用 (5.11)式, (1.9) 式以及

\int_0^\infty\left(1+s^{-\frac12}\right){\rm e}^{-s}{\rm d}s <\infty,

则可知, 当 t\rightarrow\infty 时, 有 \Pi_2\rightarrow0. 联立 (5.17) 式, 表明 \|\nabla\widehat{v}(\cdot,t)\|_{L^r} \rightarrow0, 最后再利用 (5.10) 式和 Sobolev 不等式, 可得 (5.16) 式.

利用引理 5.6 和引理 5.7, 可立即得到定理 1.2 刻画的长时间行为.

定理 1.2 的证明 事实上, 我们仅需联立引理 5.6 和引理 5.7 和回顾 \widehat{v} 的定义就可得到 (1.11) 式, 进而完成定理 1.2 的证明.

参考文献

Ahn J, Kang K, Lee J.

Global well-posedness of logarithmic Keller-Segel type systems

J Differential Equations, 2021, 287: 185-211

[本文引用: 2]

Bellomo N, Colasuonno F, Knopoff D, Soler J.

From a systems theory of sociology to modeling the onset and evolution of criminality

Math Models Methods Appl Sci, 2015, 10: 421-441

[本文引用: 1]

Bellomo N, Outada N, Soler J, et al.

Chemotaxis and cross-diffusion models in complex environments: Models and analytic problems toward a multiscale vision

Math Models Methods Appl Sci, 2022, 32(4): 713-792

[本文引用: 1]

Berestycki H, Nadal J.

Self-organised critical hot spots of criminal activity

Eur J Appl Math, 2010, 21: 371-399

[本文引用: 1]

Berestycki H, Wei J, Winter M.

Existence of symmetric and asymmetric spikes for a crime hotspot model

SIAM J Math Anal, 2014, 46(1): 691-719

[本文引用: 1]

Cao X.

Global bounded solutions of the higher-dimensional Keller-Segel system under smallness conditions in optimal spaces

Discrete Contin Dyn Syst, 2015, 35(5), 1891-1904

[本文引用: 1]

Cantrell R, Cosner C, Manásevich R.

Global bifurcation of solutions for crime modeling equations

SIAM J Math Anal, 2012, 44(3): 1340-1358

[本文引用: 1]

D'Orsogna M, Perc M.

Statistical physics of crime: A review

Phys Life Rev, 2015, 12: 1-21

DOI:10.1016/j.plrev.2014.11.001      PMID:25468514      [本文引用: 1]

Containing the spread of crime in urban societies remains a major challenge. Empirical evidence suggests that, if left unchecked, crimes may be recurrent and proliferate. On the other hand, eradicating a culture of crime may be difficult, especially under extreme social circumstances that impair the creation of a shared sense of social responsibility. Although our understanding of the mechanisms that drive the emergence and diffusion of crime is still incomplete, recent research highlights applied mathematics and methods of statistical physics as valuable theoretical resources that may help us better understand criminal activity. We review different approaches aimed at modeling and improving our understanding of crime, focusing on the nucleation of crime hotspots using partial differential equations, self-exciting point process and agent-based modeling, adversarial evolutionary games, and the network science behind the formation of gangs and large-scale organized crime. We emphasize that statistical physics of crime can relevantly inform the design of successful crime prevention strategies, as well as improve the accuracy of expectations about how different policing interventions should impact malicious human activity that deviates from social norms. We also outline possible directions for future research, related to the effects of social and coevolving networks and to the hierarchical growth of criminal structures due to self-organization. Copyright © 2014 Elsevier B.V. All rights reserved.

Espejo E, Winkler M.

Global classical solvability and stabilization in a two dimensional chemotaxis-Navier-Stokes system modeling coral fertilization

Nonlinearity, 2018, 31: 1227-1259

[本文引用: 1]

Freitag M.

Global solutions to a higher-dimensional system related to crime modeling

Math Meth Appl Sci, 2018, 41(16): 6326-6335

[本文引用: 7]

Fuest M, Heihoff F.

Unboundedness phenomenon in a model of urban crime

Commun Contemp Math, 2021, 26(7): Article 2350032

[本文引用: 1]

Fujie K, Ito A, Winkler M, Yokota T.

Stabilization in a chemotaxis model for tumor invasion

Discrete Contin Dyn Syst, 2016, 36: 151-169

[本文引用: 2]

Gu Y, Wang Q, Yi G.

Stationary patterns and their selection mechanism of urban crime models with heterogeneous near-repeat victimization effect

Eur J Appl Math, 2017, 28(1): 141-178

[本文引用: 1]

Heihoff F.

Generalized solutions for a system of partial differential equations arising from urban crime modeling with a logistic source term

Z Für Angew Math Phys, 2020, 71(3): Article 80

[本文引用: 2]

Hillen T, Painter K, Winkler M.

Convergence of a cancer invasion model to a logistic chemotaxis model

Math Models Methods Appl Sci, 2013, 23(1): 165-198

[本文引用: 1]

Jiang Y, Yang L.

Global solvability and stabilization in a three-dimensional cross-diffusion system modeling urban crime propagation

Acta Appl Math, 2022, 178: Article 11

[本文引用: 2]

Jones P, Brantingham P, Chayes L.

Statistical models of criminal behavior: the effects of law enforcement actions

Math Models Methods Appl Sci, 2010, 20: 1397-1423

[本文引用: 4]

Kolokolnikov T, Ward M, Wei J.

The stability of steady-state hot-spot patterns for a reaction-diffusion model of urban crime

Discrete Contin Dyn Syst, 2014, 19(5): 1373-1410

[本文引用: 1]

Ladyzhenskaya O A, Ural'tseva N W. Linear and Quasilinear Elliptic Equations. New York-London: Academic Press, 1968

[本文引用: 1]

Ladyzhenskaya O A, Solonnikov V A, Ural'tseva N.

Linear and Quasilinear Equations of Parabolic Type. Translations of Mathematical Monographs

Providence, RI: American Mathematical Society, 1968

[本文引用: 1]

Lankeit J, Winkler M.

A generalized solution concept for the Keller-Segel system with logarithmic sensitivity: global solvability for large nonradial data

NoDEA-Nonlinear Differ Equ Appl, 2017, 24(4): Article 49

[本文引用: 1]

Li B, Xie L.

Generalized solution to a 2D parabolic-parabolic chemotaxis system for urban crime: Global existence and large time behavior

Eur J Appl Math, 2024, 35(3): 409-429

[本文引用: 1]

Li B, Xie L.

Generalized solution and eventual smoothness in a logarithmic Keller-Segel system for criminal activities

Math Models Methods Appl Sci, 2023, 33(6): 1281-1330

[本文引用: 4]

Li B, Xie L.

Smoothness effects of a quadratic damping term of mixed type on a chemotaxis-type system modeling propagation of urban crime

Nonliear Anal Real World Appl, 2023, 73: Article 103912

[本文引用: 1]

Li B, Wang Z, Xie L.

Regularization effect of the mixed-type damping in a higher-dimensional logarithmic Keller-Segel system related to crime modeling

Math Biosci Eng, 2023, 24(3): 4532-4559

[本文引用: 2]

Lloyd D, O'Farrell H.

On localised hotspots of an urban crime model

Phys D, 2013, 253: 23-39

[本文引用: 1]

Lloyd D, Santitissadeekorn N, Short M.

Exploring data assimilation and forecasting issues for an urban crime model

Eur J Appl Math, 2016, 27(3): 451-478

[本文引用: 1]

Mei L, Wei J.

The existence and stability of spike solutions for a chemotax is system modeling crime pattern formation

Math Models Methods Appl Sci, 2020, 30: 1727-1764

[本文引用: 1]

Pitcher A.

Adding police to a mathematical model of burglary

Eur J Appl Math, 2010, 21: 401-419

[本文引用: 1]

Qiu Z, Li B.

Eventual smoothness of generalized solutions to a singular chemotaxis system for urban crime in space dimension 2

Electron Res Arch, 2023, 31(6): 3218-3244

[本文引用: 3]

Rodríguez N, Bertozzi A.

Local existence and uniqueness of solutions to a PDE model for criminal behavior

Math Models Methods Appl Sci, 2010, 20(supp1): 1425-1457

[本文引用: 1]

Rodríguez N, Winkler M.

On the global existence and qualitative behavior of one-dimensional solutions to a model for urban crime

Eur J Appl Math, 2022, 33(5): 919-959

[本文引用: 2]

Rodríguez N, Winkler W.

Relaxation by nonlinear diffusion enhancement in a two-dimensional cross-diffusion model for urban crime propagation

Math Models Methods Appl Sci, 2020, 30: 2105-2137

[本文引用: 1]

Shen J, Li B.

Mathematical analysis of a continuous version of statistical models for criminal behavior

Math Meth Appl Sci, 2020, 43(1): 409-426

[本文引用: 2]

Short M, D'Orsogna M, Pasour V, et al.

A statistical model of criminal behavior

Math Models Methods Appl Sci, 2008, 18: 1249-1267

[本文引用: 1]

Short M, Brantingham P, Bertozzi A, Tita G.

Dissipation and displacement of hotspots in reaction-diffusion models of crime

Proc Natl Acad Sci, 2010, 107: 3961-3965

[本文引用: 2]

Short M, Bertozzi A, Brantingham P.

Nonlinear patterns in urban crime: hotspots, bifurcations, and suppression

SIAM J Appl Dyn Syst, 2010, 9(9): 462-483

[本文引用: 3]

Short M, Mohler G, Brantingham P, Tita G.

Gang rivalry dynamics via coupled point process networks

Discrete Contin Dyn Syst, 2014, 19(5): 1459-1477

[本文引用: 1]

Tania N, Vnderleib B, Heathc J, Edelstein-Keshetc L.

Role of social interactions in dynamic patterns of resource patches and forager aggregation

Proc Natl Acad Sci, 2012, 109: 1128-1133

[本文引用: 1]

Tao Y, Winkler M.

Global smooth solutions in a two-dimensional cross-diffusion system modeling propagation of urban crime

Commun Math Sci, 2021, 19(3): 829-849

[本文引用: 2]

Tao Y, Winkler M.

Small-signal solutions to a nonlocal cross-diffusion model for interaction of scroungers with rapidly diffusing foragers

Math Models Methods Appl Sci, 2023, 33(1): 103-138

[本文引用: 1]

Tse W, Ward M.

Asynchronous instabilities of crime hotspots for a 1-D reaction-diffusion model of urban crime with focused police patrol

SIAM J Appl Dyn Syst, 2018, 17(3): 2018-2075

[本文引用: 1]

Tse W, Ward M.

Hotspot formation and dynamics for a continuum model of urban crime

Eur J Appl Math, 2016, 27: 583-624

[本文引用: 1]

Wang D, Feng Y.

Global well-posedness and uniform boundedness of a higher dimensional crime model with a logistic source term

Math Meth Appl Sci, 2022, 45(8): 4727-4740

[本文引用: 1]

Wang Q, Wang D, Feng Y.

Global well-posedness and uniform boundedness of urban crime models: One-dimensional case

J Differential Equations, 2020, 269(7): 6216-6235

[本文引用: 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.

Global solutions in a fully parabolic chemotaxis system with singular sensitivity

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

[本文引用: 4]

Winkler M.

Global solvability and stabilization in a two-dimensional cross-diffusion system modeling urban crime propagation

Ann Inst Henri Poincaré Anal Non Linéaire, 2019, 36: 1747-1790

[本文引用: 2]

Winkler M.

The two-dimensional Keller-Segel system with singular sensitivity and signal absorption: global large-data solutions and their relaxation properties

Math Models Methods Appl Sci, 2016, 26: 987-1024

[本文引用: 2]

Winkler M.

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

Ann Inst Henri Poincaré Anal Non Linéaire, 2024, 41: 95-127

[本文引用: 1]

Xiang T.

Finite time blow-up in the higher dimensional parabolic-elliptic-ODE minimal chemotaxis-haptotaxis system

J Differential Equations, 2022, 336: 44-72

[本文引用: 1]

Yang L, Yang X.

Global existence in a two-dimensional nonlinear diffusion model for urban crime propagation

Nonliear Anal, 2022, 224: 113086

[本文引用: 1]

Zipkin J, Short M, Bertozzi A.

Cops on the dots in a mathematical model of urban crime and police response

Discrete Contin Dyn Syst, 2014, 19(5): 1479-1506

[本文引用: 1]

/