[论文解读] On the global existence and qualitative behavior of one-dimensional solutions to a model for urban crime

Nancy Rodríguez, Michael Winkler|arXiv (Cornell University)|Mar 15, 2019

Mathematical and Theoretical Epidemiology and Ecology Models被引用 9

一句话总结

本文建立了对一维交叉扩散PDE模型在城市犯罪动力学中的经典解的全局存在性与长时间行为，其中犯罪个体被吸引至高吸引力区域。对于任意 χ > 0，当源项平均值为正时，解始终保持全局有界；在犯罪源项衰减假设下，u 衰减至零，而 v 趋近于稳态分布，且在 χ < √(6√3 + 9)/2 ≈ 2.17 范围内存在显式估计，涵盖物理上相关的 χ = 2 情况。

ABSTRACT

We consider the no-flux initial-boundary value problem for the cross-diffusive evolution system \begin{eqnarray*} \left\{ \begin{array}{ll} u_t = u_{xx} - χ\big(\frac{u}{v} \partial_x u \big)_x - uv +B_1(x,t), \qquad &amp; x\in Ω, \ t&gt;0, \\[1mm] v_t = v_{xx} +uv - v + B_2(x,t), \qquad &amp; x\in Ω, \ t&gt;0, \end{array} ight. \end{eqnarray*} which was introduced by Short et al. in [Short2008] with $χ=2$ to describe the dynamics of urban crime In bounded intervals $Ω\subset\mathbb{R}$ and with prescribed suitably regular nonnegative functions $B_1$ and $B_2$, we first prove the existence of global classical solutions for any choice of $χ&gt;0$ and all reasonably regular nonnegative initial data. We next address the issue of determining the qualitative behavior of solutions under appropriate assumptions on the asymptotic properties of $B_1$ and $B_2$. Indeed, for arbitrary $χ&gt;0$ we obtain boundedness of the solutions given strict positivity of the average of $B_2$ over the domain; moreover, it is seen that imposing a mild decay assumption on $B_1$ implies that $u$ must decay to zero in the long-term limit. Our final result, valid for all $χ\in\left(0,\frac{\sqrt{6\sqrt{3}+9}}{2} ight),$ which contains the relevant value $χ=2$, states that under the above decay assumption on $B_1$, if furthermore $B_2$ appropriately stabilizes to a nontrivial function $B_{2,\infty}$, then $(u,v)$ approaches the limit $(0,v_\infty)$, where $v_\infty$ denotes the solution of \begin{eqnarray*} \left\{ \begin{array}{l} -\partial_{xx}v_\infty + v_\infty = B_{2,\infty}, \qquad x\in Ω, \\[1mm] \partial_x v_{\infty}=0, \qquad x\in\partialΩ. \end{array} ight. \end{eqnarray*}

研究动机与目标

• 为一个一维交叉扩散系统建立经典解的全局存在性，该系统用于建模城市犯罪动力学，且 χ > 0 任意。
• 在源项 B₁ 和 B₂ 的渐近条件下，分析解的长时间行为。
• 确定犯罪密度 u 衰减至零且吸引力 v 趋近于稳态分布的条件。
• 通过数值模拟验证理论结果，探索大 χ 值对短期解增长的影响。

提出的方法

• 分析一个耦合犯罪密度 u 的交叉扩散PDE 与吸引力 v 的反应-扩散PDE 的无通量初边值问题。
• 采用能量估计与最大值原理技术，证明对任意 χ > 0，经典解的全局存在性与有界性。
• 利用辅助ODI引理的衰减估计与平均化方法，分析在 B₁ 和 B₂ 的衰减与正性假设下的长时间行为。
• 应用 Arzelà-Ascoli 定理，建立 v 对稳态解 v∞ 的预紧性与一致收敛性。
• 推导出关于 χ 的显式估计（上限为 √(6√3 + 9)/2 ≈ 2.17），使得在 B₁ 与 B₂ 衰减并稳定时，(u,v) → (0, v∞) 成立。
• 使用 MATLAB 的 pdepe 进行数值模拟，探索大 χ 情况下 u 的瞬态增长，显示其最终趋于稳态。

实验结果

研究问题

• RQ1 在何种条件下，系统对任意 χ > 0 存在全局经典解？
• RQ2 当源项 B₁ 随时间衰减时，犯罪密度 u 的长时间行为如何？
• RQ3 当 B₂ 稳定至非平凡极限时，吸引力场 v 是否收敛至稳态分布？
• RQ4 当 χ 超过分析中推导出的理论界限时，特别是短期演化中，解的动力学行为如何？
• RQ5 数值模拟能否揭示大 χ 情况下 u 的瞬态增长现象，即使长期行为稳定？

主要发现

• 对所有 χ > 0 及足够光滑的非负初值与源项 B₁、B₂，全局经典解存在。
• 若 B₂ 在区域上的平均值严格为正，则解保持有界。
• 若 B₁ 衰减足够快，则当 t → ∞ 时，犯罪密度 u 衰减至零。
• 当 χ < √(6√3 + 9)/2 ≈ 2.17（包含 χ = 2 的情况）时，有 (u,v) → (0, v∞) 成立，其中 v∞ 是一个稳态椭圆问题的解。
• 数值模拟显示，对于大 χ，u 的 L∞ 范数在短时间尺度内急剧上升，随后最终衰减至稳态。
• 当 χ 较大时，稳态解变为非均匀且呈驼峰形（在原点处峰值），偏离于均匀稳态 u ≡ 1/2, v ≡ 2。