数学物理学报, 2019, 39(6): 1381-1404 doi:

论文

具渗流扩散和间接信号产生的Prey-Taxis模型解的整体存在性与稳定性

张利民, 徐海燕, 金春花,

Global Existence and Stability to a Prey-Taxis Model with Porous Medium Diffusion and Indirect Signal Production

Zhang Limin, Xu Haiyan, Jin Chunhua,

通讯作者: 金春花, E-mail: jinchhua@126.com

收稿日期: 2018-09-27  

基金资助: 国家自然科学基金.  11871230
国家自然科学基金.  11571380
广东省自然科学基金杰出青年基金.  2015A030306029

Received: 2018-09-27  

Fund supported: the NSFC.  11871230
the NSFC.  11571380
the Distinguished Young Peopleundefineds Fund of Guangdong Natural Science Fund.  2015A030306029

摘要

该文考虑如下prey-taxis模型

在三维有界区域上的零流边值问题.该文证明了对任意的m1>1,m2>1,对任意大的初值,模型存在一个全局弱解.并在一致有界性的基础上,研究了解的大时间行为,建立了定常状态的全局渐近稳定性理论.确切地说,该文证明了当λ=0,α ≥ 0时,全局弱解强收敛到(ū0,0,0);当λ>0,α=0时,如果λ < F0ū),全局弱解强收敛到(ū0,0,0),如果λ>F0ū),全局弱解强收敛到$\left( {{{\bar u}_0}, 0, k\left( {1 - \frac{{{F_0}(\bar u)}}{\lambda }} \right)} \right) $.

关键词: 趋食性 ; 渗流扩散 ; 全局弱解 ; 一致有界 ; 稳定性

Abstract

In this paper, we consider the following prey-taxis model with nonlinear diffusion and indirect signal production

in a bounded domain of $ {{\mathbb{R}}^{3}}$ withzero-flux boundary condition. It is shown that for any m1>1, m2>1, there exists a global bounded weak solution for any large initial datum. Based on the uniform boundedness property, we also studied the large time behavior of solutions, and the global asymptotically stability of the constant steady states are established. More precisely, we showed that when λ=0, α ≥ 0, the global weak solution converges to (ū0, 0, 0) in the large time limit; when λ>0, α=0, the global weak solution converges to (ū0, 0, 0) if λ < F0(ū), and the global weak solution converges to $\left( {{{\bar u}_0}, 0, k\left( {1 - \frac{{{F_0}(\bar u)}}{\lambda }} \right)} \right) $ if λ > F0(ū).

Keywords: Prey-Taxis ; Porous Medium Diffusion ; Global Weak Solution ; Uniform Boundedness ; Stabilization

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

本文引用格式

张利民, 徐海燕, 金春花. 具渗流扩散和间接信号产生的Prey-Taxis模型解的整体存在性与稳定性. 数学物理学报[J], 2019, 39(6): 1381-1404 doi:

Zhang Limin, Xu Haiyan, Jin Chunhua. Global Existence and Stability to a Prey-Taxis Model with Porous Medium Diffusion and Indirect Signal Production. Acta Mathematica Scientia[J], 2019, 39(6): 1381-1404 doi:

1 引言

本文考虑以下prey-taxis模型

$\begin{eqnarray}\label{1-1}\left\{ \begin{array}{ll}u_t=\Delta u^{m_1}-\chi\nabla\cdot(u\nabla w), & \mbox{in}\ Q, \\ w_t=\Delta w-\mu w+\alpha vF_0(u), & \mbox{in}\ Q, \\ v_t=\Delta v^{m_2}+\lambda v\left(1-\frac{v}{k}\right)-vF_0(u), & \mbox{in}\ Q, \\ \left.(\nabla u^{m_1}-\chi u\cdot\nabla w)\cdot{\bf n}\right|_{\partial\Omega} =\nabla w\cdot{\bf n}|_{\partial\Omega}=\left.\nabla v^{m_2}\cdot{\bf n}\right|_{\partial\Omega}=0, \\ u(x, 0)=u_0(x), w(x, 0)=w_0(x), v(x, 0)=v_0(x), & x\in\Omega, \end{array}\right. \end{eqnarray}$

其中${{m}_{1}}>1,{{m}_{2}}>1,Q=\Omega \times {{\mathbb{R}}^{+}},\Omega \subset {{\mathbb{R}}^{3}}$是一个有界的光滑区域. $u, w, v$分别代表捕食者的密度,猎物释放的趋化因子浓度和猎物种群密度, $\chi$是趋化灵敏度系数, $\mu$是趋化因子的衰减率, $\lambda$猎物种群的增长率和$k$代表环境的承载容量. $F_0$是如下的Holling Ⅱ型函数

其中, $F_0(u)$为由于捕食者活动而导致的猎物死亡率, $A$为猎物被捕食者捕杀的半饱和常数, $F_m$为最大死亡率.

Prey-taxis模型首先由Karevia和Odell在文献[6]中提出.不同于经典的捕食者-食饵模型,捕食者和猎物除了受随机扩散的影响,捕食者的运动也受猎物密度梯度的影响,即prey-taxis项$-\nabla\cdot(\chi(u, v)u\nabla v)$.此后,许多研究者对具有prey-taxis项的捕食-食饵模型进行了广泛的研究.在2008年,文献[1]研究了如下系统

$ \begin{eqnarray}\label{eq-2} \left\{ \begin{array}{ll} &u_t=d_1\Delta u-\nabla\cdot(\chi(u)\nabla v)-au+bg(v)u, \\ &v_t=d_2\Delta v++k(v)-g(v)u \end{array}\right. \end{eqnarray} $

弱解的存在唯一性,并且通过数值模拟显示了班图形成的有趣现象.随后,文献[10]得到了空间维数为1, 2和3时古典解的整体存在性和唯一性.对于多个物种的存在性结果,请参考文献[6, 17].在2009年, Lee等在文献[7]对于如下模型

$ \begin{eqnarray}\label{eq-3} \left\{ \begin{array}{ll} &u_t=\Delta u-\nabla\cdot(\chi(v)u\nabla v)+\gamma u(h(v, u)-\delta(u)), \\ &v_t=\varepsilon\Delta v+vf(v)-uh(v, u), \end{array}\right. \end{eqnarray}$

研究了一维空间中捕食者和猎物的运动对捕食者-食饵系统空间班图形成的作用,并针对不同的情况建立了班图形成的必要条件.更多关于这类模型的研究,我们也可以参考文献[2, 5, 13-14]等.近年来,一种带有间接信号产生的捕食模型被提出,其中捕食者跟随猎物释放的某些化学物质的梯度移动.例如,猎物在捕获过程中分泌的化学引诱物或受伤猎物血液的气味.在文献[11]中, Tello和Wrzosek提出如下模型

其中,猎物的随机运动被忽略.在对系数施加适当的假设条件下,对任意空间维数作者证明了解的整体存在性和局部渐近稳定性,特别是在1维空间中,证明了非常数稳态解的存在性.

然而,从生物学的角度来看,捕食者和猎物的迁移应该被认为是非线性扩散,而渗流扩散是最常用的模型.本文受文献[11]的启发,我们考虑了一个具有渗流慢扩散的捕食模型,其中,捕食者或猎物的迁移率依赖于它们自身的密度.对任意的$m_1>1, m_2>1$,我们首先在3维空间中给出了弱解的全局存在和有界性结果.在此基础上,进一步研究了不同系数假设条件下的平衡态的稳定性.

我们首先给出了本文的假设.

$\left\{ \begin{array}{ll} & u_0, v_0\in L^\infty(\Omega), \nabla u_0^{m_1}, \nabla v_0^{m_2}\in L^2(\Omega), w_0\in W^{2, \infty}(\Omega), \\ & u_0, w_0, v_0\ge 0. \end{array}\right. \$

下面我将给出本篇文章的主要定理.

定理1.1  假设(H)成立, $m_1>1, m_2>1$.那么对任意的$\chi>0, \mu>0, \alpha\ge 0, \lambda\ge 0$,问题(1.1)存在一个非负弱解$(u, w, v)$满足$u\in{\cal X}_1, w\in{\cal X}_2, v\in {\cal X}_3$,其中

$ \begin{equation} \label{1-2} \sup\limits_{t\in(0, \infty)}\left(\|u(\cdot, t)\|_{L^\infty}+\|w\|_{W^{1, \infty}}+\|v(\cdot, t)\|_{L^\infty}\right)\leq M_1, \end{equation}$

$\begin{equation} \label{1-3} \sup\limits_{t\in (0, +\infty)}\int_\Omega(|\nabla u^{m_1}|^2+|\nabla v^{m_2}|^2){\rm d}x+ \sup\limits_{t\in (0, +\infty)}\|(u^{\frac{m_1+1}2}, v^{\frac{m_2+1}2}) \|_{W_2^{1, 1}(Q_1(t))} \le M_2, \end{equation}$

$ \begin{equation} \label{1-4} \sup\limits_{t\in (0, +\infty)}\|w\|_{W_p^{2, 1}(Q_1(t))} \le M_3, \quad \forall p>1, \end{equation} $

其中$M_i$$(i=1, 2, 3)$仅依赖于$\Omega$, $\chi, \mu, \alpha, \lambda, u_0, w_0, v_0$.

定理1.2  假设$u_0\not\equiv0, v_0\not\equiv0$.$(u, w, v)\in {\mathcal X}_1\times{\cal X}_2\times{\cal X}_3$是问题(1.1)的一个整体解.那么对任意的初值,对任意的$\chi>0, \mu>0, \lambda=0, \alpha\geq0$,我们有

定理1.3  假设$u_0\not\equiv0, v_0\not\equiv0$.$(u, w, v)\in {\mathcal X}_1\times{\cal X}_2\times{\cal X}_3$是问题(1.1)的一个整体解.那么对任意的初值,对任意的$\chi>0, \mu>0, \lambda>0, \alpha=0$,如果$\lambda<F_0(\bar{u})$,我们有

如果$\lambda>F_0(\bar{u})$,对任意的$p\ge 1$,我们有

2 预备知识

为了读者方便,我们首先引入一些本文中使用的记号.

记号:$\|\cdot\|_{L^p}=\|\cdot\|_{L^p(\Omega)}$, $Q_{\tau}(t)=\Omega\times(t, t+\tau)$, $Q_T:=Q_T(0)=\Omega\times (0, T)$, $\bar{\varphi}:=\frac{1}{|\Omega|}\displaystyle\int_\Omega\varphi {\rm d}x$.

接下来,我们给出弱解的定义.

定义2.1  称$(u, w, v)$为问题(1.1)的弱解,如果$(u, v, \omega)\in{\cal X}_{1}\times{\cal X}_{2}\times{\cal X}_{3}$$u, v, \omega\geq0$,使得对任意给的$T>0$,我们有

其中$\varphi\in C^{\infty}(\overline{Q_{T}})$$\left.\frac{\partial\varphi}{\partial{\bf n}}\right|_{\partial\Omega}=0$$\varphi(x, T)=0$.

根据文献[4],我们有下面两个引理.

引理2.1  对常数$T>0$,假设$f\geq0, f(t)\in L^1(T, \infty)$并且

那么

引理2.2  对常数$T\geq0$,假设$f(t), g(t) \ge 0, \mathop {\lim }\limits_{t \to \infty } g(t) = 0$$f(t)\in L^1(T, \infty)$.$F(t)=f(t)-g(t)$,有

那么$ \mathop {\lim }\limits_{t \to \infty } f(t) = 0.$

3 弱解的整体存在性及有界性

我们首先考虑以下的正则化问题

$\begin{eqnarray}\left\{\begin{array}{ll}u_{\varepsilon t}=\Delta((u_\varepsilon+\varepsilon)^{m_1-1}u_\varepsilon)-\chi\nabla\cdot(u_\varepsilon\nabla w_\varepsilon), \ \mbox{in}\ Q, \\w_{\varepsilon t}=\Delta w_\varepsilon-\mu w_\varepsilon+\alpha v_\varepsilon F_0(u_\varepsilon), \ \mbox{in}\ Q, \\v_{\varepsilon t}=\Delta((v_\varepsilon+\varepsilon)^{m_2-1} v_\varepsilon)+\lambda v_\varepsilon\left(1-\frac{v_\varepsilon}{k}\right)-v_\varepsilon F_0(u_\varepsilon), \ \mbox{in}\ Q, \\\left.\frac{\partial u_{\varepsilon}}{\partial {\bf n}}\right|_{\partial\Omega}=\left.\frac{\partial w_{\varepsilon}}{\partial {\bf n}}\right|_{\partial\Omega}=\left.\frac{\partial v_{\varepsilon}}{\partial {\bf n}}\right|_{\partial\Omega}=0, \\ u_{\varepsilon}(x, 0)=u_{\varepsilon 0}(x), w_{\varepsilon}(x, 0)=u_{\varepsilon 0}(x), v_{\varepsilon}(x, 0)=v_{\varepsilon 0}(x), \quad x\in\Omega, \end{array}\right.\end{eqnarray}$

其中$u_{\varepsilon 0}, w_{\varepsilon 0}, v_{\varepsilon 0}$满足

$\left\{\begin{array}{ll}u_{\varepsilon 0}, w_{\varepsilon 0}, v_{\varepsilon 0} \in C(\overline\Omega)\cap C^{2+\alpha}(\Omega), \left.\frac{\partial u_{\varepsilon 0}}{\partial{\bf n}}\right|_{\partial\Omega}=0, \left.\frac{\partial w_{\varepsilon 0}}{\partial{\bf n}}\right|_{\partial\Omega}=0, \left.\frac{\partial v_{\varepsilon 0}}{\partial{\bf n}}\right|_{\partial\Omega}=0, \\\|u_{\varepsilon 0}\|_{L^\infty}+\|\nabla u_{\varepsilon 0}^{m_1}\|_{L^2}+\|w_{\varepsilon 0}\|_{W^{2, \infty}}+\|v_{\varepsilon 0}\|_{L^\infty}+\|\nabla v_{\varepsilon 0}^{m_2}\|_{L^2}\\\leq2(\|u_0\|_{L^\infty}+\|\nabla u_0^{m_1}\|_{L^2}+\|w_0\|_{W^{2, \infty}}+\|v_0\|_{L^\infty}+\|\nabla v_0^{m_2}\|_{L^2}), \\u_{\varepsilon 0}\to u_0, w_{\varepsilon 0}\to w_0, v_{\varepsilon 0}\to v_0, \quad \mbox{一致收敛}.\end{array}\right.$

运用不动点定理,我们可以证明以下局部存在结果,参考文献[9, 12, 16].如下,我们给出了问题(3.1)古典解的局部存在性结果.

引理3.1  假设$u_{\varepsilon 0}$, $w_{\varepsilon 0}$, $v_{\varepsilon 0}$满足(H1).那么存在$T_{\max}\in(0, +\infty]$使得问题(3.1)存在唯一的古典解$(u_\varepsilon, w_\varepsilon, v_\varepsilon)\inC(\overline\Omega \times[0, T_{\max}))\cap C^{2+\alpha, 1+\alpha/2}(\Omega \times(0, T_{\max}))$满足

$\begin{equation}u_\varepsilon\ge 0, \quad w_\varepsilon\ge 0, \quad v_\varepsilon\ge 0, \ \forall (x, t)\in \Omega\times (0, T_{\max}), \end{equation}$

使得要么$T_{max} = \infty$,要么

在本文中,我们证明了正则化问题古典解的整体存在性.我们给出以下命题.

命题3.1  假设$u_{\varepsilon 0}$, $w_{\varepsilon 0}$, $v_{\varepsilon 0}$满足(H1).那么对任意的$\varepsilon>0$,问题(3.1)存在唯一的古典解$(u_\varepsilon, w_\varepsilon, v_\varepsilon)$,满足

$\begin{equation}\sup\limits_{t\in(0, \infty)}\left(\|u_\varepsilon(\cdot, t)\|_{L^\infty}+\|w_\varepsilon\|_{W^{1, \infty}}+\|v_\varepsilon(\cdot, t)\|_{L^\infty}\right)\leq M_1, \end{equation}$

$\begin{eqnarray}&&\sup\limits_{t\in (0, +\infty)}\int_\Omega(|\nabla((u_\varepsilon+\varepsilon)^{m_1-1}v_\varepsilon)|^2+|\nabla((v_\varepsilon+\varepsilon)^{m_2-1}v_\varepsilon)|^2) {\rm d}x\nonumber\\&&+\sup\limits_{t\in (0, +\infty)}\|(u_\varepsilon^{\frac{m_1+1}2}, v_\varepsilon^{\frac{m_2+1}2})\|_{W_2^{1, 1}(Q_1(t))}\le M_2, \end{eqnarray}$

$\begin{equation}\label{add3}\sup\limits_{t\in (0, +\infty)}\|w_\varepsilon\|_{W_p^{2, 1}(Q_1(t))}\le M_3, \quad \forall p>1, \end{equation}$

其中$M_i (i=1, 2, 3)$都不依赖于$\varepsilon$.

接下来,我们给出$(u_\varepsilon, w_\varepsilon, v_\varepsilon)$的一些先验估计,我们令$\tau=\min\left\{1, \frac{T_{\max}}{2}\right\}$.容易发现$\tau\le 1$.在下文中,所有的常数$C$, $C_i$, $\tilde C$, $\tilde C_i$, $M$都不依赖于$\tau$$T_{\max}$.

引理3.2  假设$u_{\varepsilon 0}$, $w_{\varepsilon 0}$, $v_{\varepsilon 0}$满足(H1).令$(u_\varepsilon, w_\varepsilon, v_\varepsilon)$是模型(3.1)在$(0, T_{\max})$上的古典解.那么

$\begin{equation} \|u_\varepsilon(\cdot, t)\|_{L^1}=\|u_{\varepsilon0}\|_{L^1}, \end{equation}$

$\begin{equation} \|w_\varepsilon(\cdot, t)\|_{L^\infty}+\|v_\varepsilon(\cdot, t)\|_{L^\infty}\leq C, \end{equation}$

其中$C$不依赖于$\varepsilon$.

  模型(3.1)的第一个方程直接积分,可得(3.6)式.

$\tilde{V}:=\max\{k, \|v_0\|_{L^\infty}\}$,显然, $\tilde{V}$是模型(3.1)第三个方程的一个上解.根据比较原理,我们有$v_\varepsilon\le \tilde{V}$.类似地,令

那么$\tilde{W}$是模型(3.1)第二个方程的一个上解.那么,容易发现(3.7)式成立.证毕.

引理3.3  假设$u_{\varepsilon 0}$, $w_{\varepsilon 0}$, $v_{\varepsilon 0}$满足(H1).令$(u_\varepsilon, w_\varepsilon, v_\varepsilon)$是模型(3.1)在$(0, T_{\max})$上的古典解.那么

$\begin{equation}\|w_\varepsilon\|_{W^{1, \infty}}\leq C, \end{equation}$

其中$C$不依赖于$\varepsilon$.

  根据Duhamel's原理,我们发现解$w_\varepsilon$可以被如下表示

其中$\{e^{t\Delta}\}_{t\geq0}$是在$\Omega$上的Neumann热半群,更多关于Neumann热半群的知识,请参考文献[15].根据(3.7)式,对所有的$t\in(0, T_{\max})$我们有

结合(3.7)式可得引理3.3.证毕.

引理3.4  假设$u_{\varepsilon 0}$, $w_{\varepsilon 0}$, $v_{\varepsilon 0}$满足(H1),并且$m_1>1, m_2>1$.$(u_\varepsilon, w_\varepsilon, v_\varepsilon)$是模型(3.1)在$(0, T_{\max})$上的古典解.则

$\begin{equation}\sup\limits_{t\in(0, T_{max})}\|u_\varepsilon\|_{L^r}\leq C(r), \end{equation}$

其中$C(r)$依赖于$r$.

  对任意的$r>m_1$,在模型(3.1)的第一个方程两边同时乘以$r u_\varepsilon^{r-1}$,然后在$\Omega$上积分,并且结合(3.8)式,可得

根据Gagliardo-Nirenberg插值不等式,对任意小的$\eta>0$

结合上面三个不等式,同时利用(3.7)式,我们可得

这意味着(3.9)式成立.证毕.

下面,我们证明$u_\varepsilon$$L^\infty$估计.

引理3.5  假设$u_{\varepsilon 0}$, $w_{\varepsilon 0}$, $v_{\varepsilon 0}$满足(H1),并且$m_1>1, m_2>1$.$(u_\varepsilon, w_\varepsilon, v_\varepsilon)$是模型(3.1)在$(0, T_{\max})$上的古典解,则

$\begin{equation}\sup\limits_{t\in(0, T_{max})}\|u_\varepsilon\|_{L^\infty}\leq C.\end{equation}$

  对任意的$r\ge2m_1$,在模型(3.1)的第一个方程两边同时乘以$r u_\varepsilon^{r-1}$,然后在$\Omega$上积分,并且结合(3.8)式,可得

这意味着

$\begin{eqnarray}&&\frac{\rm d}{{\rm d}t}\int_\Omega u_\varepsilon^{r} {\rm d}x+\frac{2r(r-1)}{(m_1+r-1)^2}\int_\Omega |\nabla u_\varepsilon^{\frac{r+m_1-1}2}|^2 {\rm d}x+\int_\Omega u_\varepsilon^{r} {\rm d}x \nonumber\\\label{3-8}&\le &Cr^2\int_\Omega u_\varepsilon^{r+1-m_1}{\rm d}x+\int_\Omega u_\varepsilon^{r} {\rm d}x.\end{eqnarray}$

根据Gagliardo-Nirenberg插值不等式,可得

将上面两个不等式代入(3.11)式,可得

其中$C_i$$(i=3, 4, 5, 6)$都不依赖于$r$.注意到

进一步可得

$\begin{eqnarray}\frac{\rm d}{{\rm d}t}\int_\Omega u_\varepsilon^r{\rm d}x+\int_\Omega u_\varepsilon^r{\rm d}x\leq C_7r^5\|u_\varepsilon\|_{L^\frac{r}{2}}^\frac{r(2(m_1-1)+r)}{6(m_1-1)+r}+C_8r^2\|u_\varepsilon\|_{L^\frac{r}{2}}^r.\end{eqnarray}$

$r_j=2r_{j-1}=2^jr_0, r_0=2m_1, M_j=\max\Big\{1, \|u_0\|_{L^\infty}, \sup\limits_{t\in(0, T_{\max})}\|u_\varepsilon\|_{L^{r_j}}\Big\}$,那么有

其中$C$不依赖于$j$.$j\to\infty$,则(3.10)式得证.证毕.

命题3.1的证明  根据引理3.1和引理3.5,对任意的$\varepsilon>0$, $T_{\max}=+\infty$,模型(3.1)存在一个整体的古典解.这意味着$\tau=1$.

注意到$v_\varepsilon F_0(u_\varepsilon)$有界.那么根据线性抛物型方程的$L^p$理论,对任意的$p>1$,我们可得

$\begin{equation}\sup\limits_{t\in(1, \infty)}\int_{t-1}^t\int_\Omega|\Delta w|^p{\rm d}x{\rm d}t\leq C.\end{equation}$

在模型(3.1)的第一个方程两边同时乘以$u_\varepsilon$,然后在$\Omega$上积分,可得

根据(3.3)和(3.13)式,易得

$\begin{equation}\sup\limits_{t\in(1, \infty)}\int_{t-1}^t\int_\Omega (u_\varepsilon+\varepsilon)^{m_1-1}|\nabla u_\varepsilon|^2 {\rm d}x {\rm d}s\le C.\end{equation}$

在模型(3.1)的第三个方程两边同时乘以$v_\varepsilon$,然后在$\Omega$上积分,可得

根据(3.3)式,易得

$\begin{equation}\sup\limits_{t\in(1, \infty)}\int_{t-1}^t\int_\Omega (v_\varepsilon+\varepsilon)^{m_2-1}|\nabla v_\varepsilon|^2 {\rm d}x {\rm d}s\le C.\end{equation}$

在模型(3.1)的第一个方程两边同时乘以$\frac{\partial( (u_\varepsilon+\varepsilon)^{m_1-1}u_\varepsilon)}{\partial t}$,然后在$\Omega$上积分可得

这意味着

根据(3.13)和(3.14)式,进一步有

$\begin{equation}\sup\limits_{t\in (0, +\infty)}\int_\Omega|\nabla( (u_\varepsilon+\varepsilon)^{m_1-1}u_\varepsilon)|^2 {\rm d}x+\sup\limits_{t\in (1, +\infty)}\int_{t-1}^t\int_\Omega(u_\varepsilon+\varepsilon)^{m_1-1} \left|\frac{\partial u_\varepsilon}{\partial t}\right|^2 {\rm d}x\le C.\end{equation}$

在模型(3.1)的第三个方程两边同时乘以$\frac{\partial( (v_\varepsilon+\varepsilon)^{m_2-1}v_\varepsilon)}{\partial t}$,然后在$\Omega$上积分可得

根据(3.15)式,进一步有

$\begin{equation}\sup\limits_{t\in(0, \infty)}\int_\Omega|\nabla((v_\varepsilon+\varepsilon)^{m_2-1}v_\varepsilon)|^2{\rm d}x+\sup\limits_{t\in(1, \infty)}\int_{t-1}^t\int_\Omega(v_\varepsilon+\varepsilon)^{m_2-1}\left|\frac{\partial v_\varepsilon}{\partial t}\right|^2{\rm d}x{\rm d}s\leq C.\end{equation}$

结合(3.16)式和(3.17)式,可得(3.4)式.

定理1.1的证明  由于$(u_\varepsilon, w_\varepsilon, v_\varepsilon, )$是模型(3.1)的古典解,对任意的$\varphi\in C^\infty(\overline Q_T)$,有

其中$\left.\frac{\partial\varphi}{\partial {\bf n}}\right|_{\partial\Omega}=0$$\varphi(x, T)=0$.

根据Aubin-Lions定理,利用命题3.1,对任意$T>0$,我们有

并且对于$u, w, v$,有(1.4)-(1.6)式成立.注意到$u_\varepsilon$$u$都是一致有界的,那么

因此,在$L^\infty(Q_T)$中,对任意的$p\in (1, +\infty)$

类似地,在$L^\infty(Q_T)$中,对任意的$p\in (1, +\infty)$也有

$\varepsilon\to 0$,那么可得

注意到$\nabla u^{m_1}\in L^2(Q_T)$, $\nabla v^{m_2}\in L^2(Q_T)$,那么也有

这意味着满足$u\in{\cal X}_1$, $w\in{\cal X}_2$, $v\in{\cal X}_3$$(u, w, v)$是模型(1.1)的一个整体弱解.证毕.

4 $\lambda=0, \alpha\geq0$时解的大时间行为

在这部分,我们基于$\lambda=0$, $\alpha\geq 0$研究弱解的大时间行为.首先有以下引理.

引理4.1  令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么有

$\begin{equation} \int_0^\infty\int_\Omega vF_0(u){\rm d}x{\rm d}t\leq C_1, \end{equation}$

$\begin{equation} \int_0^\infty\int_\Omega w^2{\rm d}x{\rm d}t+\int_0^\infty\int_\Omega |\nabla w|^2{\rm d}x{\rm d}t\leq C_2, \end{equation}$

$\begin{equation} \int_0^\infty\int_\Omega u^{r+m_1-2}|\nabla u|^2{\rm d}x{\rm d}t\leq C_3, \quad \forall r>(m_1-2)_+, \end{equation}$

其中$C_1$, $C_2$, $C_3$依赖于$u_0, w_0, v_0, \mu, \alpha, \chi, \Omega$.

  模型(1.1)的第三个方程直接在$\Omega\times (0, t)$上积分,对所有的$t>0$,可得

从中可得(4.1)式成立.此外,模型(1.1)的第二个方程两端同时乘以$w$,然后在$\Omega\times(0, t)$上积分,结合(4.1)式,对所有的$t>0$,可得

从而验证(4.2)式成立.

对任意的$(r+1)u^r$,在模型(1.1)的第一个方程两边同时乘以$(r+1)u^r$,可得

利用(1.4)和(4.2)式,并且上面的不等式从$0$$\infty$上积分,可得(4.3)式.证毕.

引理4.2  令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么有

$\begin{eqnarray} \lim\limits_{t\rightarrow\infty}\int_\Omega(|\nabla w(x, t)|^2+w^2(x, t)){\rm d}x=0.\end{eqnarray}$

  在模型(1.1)中第二个方程两边分别乘以$-\Delta w$, $w$,然后它们在$\Omega$上积分,可得

这意味着

$\begin{eqnarray} \frac{\rm d}{{\rm d}t}\int_\Omega(|\nabla w|^2+w^2){\rm d}x+\int_\Omega |\Delta w|^2{\rm d}x+2(\mu+1)\int_\Omega|\nabla w|^2{\rm d}x+\mu\int_\Omega w^2{\rm d}x\le C.\end{eqnarray}$

那么有

根据引理2.1和(4.2)式,可得(4.4)式.证毕.

引理4.3  令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么有

$\begin{eqnarray} \lim\limits_{t\rightarrow\infty}\|w(\cdot, t)\|_{L^\infty}=0.\end{eqnarray}$

  根据Gagliardo-Nirenberg插值不等式,我们有

利用引理4.2和(1.4)式,我们有(4.6)式.证毕.

引理4.4  令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么对任意的$p>1$

$\begin{eqnarray}\lim\limits_{t\rightarrow\infty}\|u(x, t)-\bar u_0\|_{L^p}{\rm d}x=0.\end{eqnarray}$

  令$b^{m_1}(t)=\frac{1}{|\Omega|}\int_\Omega u^{m_1}{\rm d}x$.利用(1.4)和(1.5)式,可得

结合(4.5)式,可得

这意味着对任意的$t>s>0$,有

$\begin{eqnarray}&&\int_\Omega(|u^{m_1}(x, t)-b^{m_1}(t)|^2-|\nabla w(x, t)|^2-|w(x, t)|^2){\rm d}x\nonumber\\&&-\int_\Omega(|u^{m_1}(x, s)-b^{m_1}(s)|^2-|\nabla w(x, s)|^2-| w(x, s)|^2){\rm d}x\geq-C(t-s).\end{eqnarray}$

进一步,根据Poincaré不等式,并且令(4.3)式中$r=m_1$,我们发现

根据(4.4)式和(4.8)式,结合引理2.2,可得

$\begin{eqnarray} \lim\limits_{t\rightarrow\infty}\int_\Omega|u^{m_1}(x, t)-b^{m_1}(t)|^2{\rm d}x=0.\end{eqnarray}$

根据Hölder不等式,对所有的$t>0$,我们发现

这意味着$b(t)\geq \bar u_0$.像对所有的$\xi\geq0$$\eta\geq0$$|\xi^{m_1}-\eta^{m_1}|\geq\eta^{m_1-1}\cdot|\xi-\eta|$一样,对所有的$t>0$,在$ \Omega$上处处有

因此从(4.9)式可得

$\begin{equation} \lim\limits_{t\rightarrow\infty}\int_\Omega|u(x, t)-b(t)|^2{\rm d}x=0.\end{equation}$

注意到

因此

因此

结合(4.10)式,可得

$\begin{equation}\lim\limits_{t\rightarrow\infty}\|u(x, t)-\bar u_0\|_{L^2}{\rm d}x=0.\end{equation}$

$1\le p\le 2$时,根据这个不等式我们显然有(4.7).注意到对任意的$p\ge 2$

并且利用(4.11)式,我们完成了定理的证明.证毕.

引理4.5  令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么有

$\begin{eqnarray}\label{4-12} \lim\limits_{t\rightarrow\infty}\int_\Omega v{\rm d}x=0.\end{eqnarray}$

  由模型(1.1)的第三个方程,可得

根据上面的方程,对任意的$\eta>0$,存在$t>T$,使得

解这个不等式,当$t\rightarrow\infty$时,可得

由于$\eta>0$是任意的,我们可得(4.12)式.证毕.

定理1.2的证明  根据Hölder不等式,当$t\to\infty$时,我们可得

类似于文献[4,引理3.4]的证明或者本文中引理5.7的证明,根据Morse迭代技术,我们进一步得到

利用引理4.3,引理4.4和引理4.5,我们完成了定理1.2的证明.

5 $\lambda>0, \alpha=0$时解的大时间行为

在这部分,我们基于$\lambda>0, \alpha=0$研究弱解的大时间行为.首先有以下引理.

引理5.1  令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么有

$\begin{equation} \int_0^\infty\int_\Omega w^2{\rm d}x{\rm d}t+\int_0^\infty\int_\Omega |\nabla w|^2{\rm d}x{\rm d}t\leq C_1, \end{equation}$

$\begin{equation} \int_0^\infty\int_\Omega u^{r+m_1-2}|\nabla u|^2{\rm d}x{\rm d}t\leq C_2, \quad \forall r>(m_1-2)_+, \end{equation}$

其中$C_1, C_2$依赖于$u_0, w_0, v_0, \chi, \mu, \lambda, \Omega$.利用上面的不等式,进一步有

$\begin{equation}\lim\limits_{t\rightarrow\infty}\int_\Omega (w^2(x, t)+|\nabla w(x, t)|^2){\rm d}x=0, \end{equation}$

$\begin{equation}\lim\limits_{t\rightarrow\infty}\int_\Omega|u(x, t)-\bar u_0|^2{\rm d}x=0, \end{equation}$

这意味着

$\begin{equation}\lim\limits_{t\rightarrow\infty}\|w(\cdot, t)\|_{L^\infty}=0, \end{equation}$

$\begin{equation} \lim\limits_{t\rightarrow\infty}\|u(\cdot, t)-\bar u_0\|_{L^p}=0, \ \forall \ p>1.\end{equation}$

  在模型(1.1)的第二个方程两边同时乘以$w$,然后在$\Omega\times(0, t)$上积分,可得

从中可验证(5.1)式成立.

对任意的$r>(m_1-2)_+$,在模型(1.1)的第一个方程两边同时乘以$u^r$,可得

利用(1.4)和(5.1)式,并且上面的不等式从0到$\infty$上积分,可得(5.2)式.那么类似于引理4.2和引理4.4的证明,进一步有(5.3)-(5.6)式.证毕.

由模型(1.1)的第三个方程,可得

$\begin{eqnarray}\label{5-8} v_t=\Delta v^{m_2}+v\left((\lambda-F_0(\bar{u}))-\frac{\lambda v}{k}\right)+v(F_0(\bar{u})-F_0(u)).\end{eqnarray}$

接下来,我们将分两种情况证明$v$的大时间行为.

首先,考虑$\lambda<F_0(\bar{u})$的情形.

引理5.2  假设$\lambda<F_0(\bar{u})$,并且令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么有

$\begin{equation}\lim\limits_{t\rightarrow\infty}\|v(\cdot, t)\|_{L^\infty}=0.\end{equation}$

  根据(5.7)式,对任意的$\eta>0$,存在$t>T$,使得

其中$A_1:=F_0(\bar{u})-\lambda>0$.解这个微分不等式,当$t\rightarrow\infty$时可得

由于$\eta>0$是任意的,我们有

这意味着对任意的$p\ge 1$

那么,类似于引理5.7的证明,进一步有(5.8)式.证毕.

接下来,考虑$\lambda>F_0(\bar{u})$的情况.令$A_2:=1-\frac{F_0(\bar{u})}{\lambda}$,那么(5.7)式等价于

$\begin{equation}\label{5-11} v_t=\Delta v^{m_2}+\lambda v\left(A_2-\frac{v}{k}\right)+v(F_0(\bar{u})-F_0(u)).\end{equation}$

引理5.3  假设$\lambda>F_0(\bar{u})$,并且令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么有

$\begin{eqnarray}\label{5-12}\int_0^\infty\int_\Omega|u-\bar{u}|^2{\rm d}x{\rm d}t\leq C, \end{eqnarray}$

其中$C$仅依赖于$u_0, w_0, v_0, \chi, \mu, \lambda, \Omega$.

  令$a(t)=\frac{1}{|\Omega|}\int_\Omega u(x, t){\rm d}x, b^\frac{m_1+r}{2}(t)=\frac{1}{|\Omega|}\int_\Omega u^\frac{m_1+r}{2}{\rm d}x$.注意到

这意味着$a(t)\le b(t)$.

根据Poincaré不等式,并且利用(5.2)式,我们发现

$\begin{eqnarray}\label{5-13} \int_0^\infty\int_\Omega|u^\frac{m_1+r}{2}-b^\frac{m_1+r}{2}|^2{\rm d}x{\rm d}t\leqC\int_0^\infty\int_\Omega|\nabla u^\frac{m_1+r}{2}|^2{\rm d}x{\rm d}t\leq \tilde{C}.\end{eqnarray}$

注意到$\frac{u^\frac{m_1+r}{2}-b^\frac{m_1+r}{2}}{u-b}\geq b^\frac{m_1+r-2}{2}$,那么

结合(5.11)式,可得

由于$a(t)=\frac{1}{|\Omega|}\|u_0\|_{L^1}>0$是一个常数,有

$\begin{eqnarray} \int_0^\infty\int_\Omega|u-b|^2{\rm d}x\leq C.\end{eqnarray}$

注意到

$\begin{eqnarray} |a(t)-b(t)|^2=\left|\frac{1}{|\Omega|}\int_\Omega(u-b){\rm d}x\right|^2\leq\frac{1}{|\Omega|}\int_\Omega|u-b|^2{\rm d}x, \end{eqnarray}$

结合(5.12)式和(5.13)式,可得

证毕.

引理5.4  假设$\lambda>F_0(\bar{u})$,并且令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么对任意的$\beta\ge 1$

$\begin{eqnarray}\label{5-16}\int_0^\infty\int_\Omega v^{\beta+m_2-2}|\nabla v|^2{\rm d}x{\rm d}t+\int_0^\infty\int_\Omega v\left(\frac{v}{k}-A_2\right)(v^\beta-(A_2k)^\beta){\rm d}x{\rm d}t \le C, \end{eqnarray}$

其中$C$仅依赖于$u_0, w_0, v_0, \chi, \mu, \lambda, \Omega$.

  对任意的$\beta\ge 1$, (5.9)式的两边同时乘以$v^\beta-(A_2k)^\beta$,那么有

$\varepsilon$足够小,使得$C\varepsilon<\frac \lambda 2$,那么有

结合(5.10)式,对上面的不等式从0到$\infty$上积分,我们完成了定理的证明.证毕.

引理5.5  假设$\lambda>F_0(\bar{u})$,并且令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么有

$\begin{eqnarray}\lim\limits_{t\rightarrow\infty}\int_\Omega v\left(\frac{v}{k}-A_2\right)(v^\beta-(A_2k)^\beta){\rm d}x=0, \quad \forall \beta\ge m_2.\end{eqnarray}$

  根据(1.4)式,我们发现

在上面的不等式中令$\beta\ge m_2$,并且结合(1.4)式和(1.5)式,我们有

那么$\int_\Omega v\left(\frac{v}{k}-A_2\right)(v^\beta-(A_2k)^\beta){\rm d}x$是lipschitz连续.根据(5.14)式,利用2.1,我们有(5.15)式.证毕.

接下来,我们致力于证明

首先证明下面的引理.

引理5.6  假设$\lambda>F_0(\bar{u})$,令$\tilde a(t)=\frac{1}{|\Omega|}\int_\Omega v(x, t){\rm d}x$, ${\tilde b}^{m_2}(t)=\frac{1}{|\Omega|}\int_\Omega v^{m_2}{\rm d}x$.那么

$\begin{eqnarray} \lim\limits_{t\rightarrow\infty}\int_\Omega|v^{m_2}(x, t)-{\tilde b}^{m_2}(x, t)|^2{\rm d}x=0.\end{eqnarray}$

并且如果存在一个序列$t_j$满足$t_j\nearrow\infty$和一个正常数$\sigma_0$,使得

那么$\sigma_0=A_2k$,这意味着

$\begin{eqnarray} \lim\limits_{t\rightarrow\infty}\int_\Omega|v(x, t)-A_2k|^2{\rm d}x=0.\end{eqnarray}$

  注意到

这意味着$\tilde a(t)\le \tilde b(t)$.接下来,根据一个直接计算,利用(1.4)和(1.5)式,可得

那么有

另一方面,根据Poincaré不等式,并且在(5.14)式中令$\beta=m_2$,我们发现

$\begin{eqnarray}\label{5-20} \int_0^\infty\int_\Omega|v^{m_2}-{\tilde b}^{m_2}|^2{\rm d}x{\rm d}t\leq C\int_0^\infty\int_\Omega|\nabla v^{m_2}|^2{\rm d}x{\rm d}t\leq \tilde{C}, \end{eqnarray}$

根据引理2.1,得到(5.16)式.注意到$\frac{v^{m_2}-{\tilde b}^{m_2}}{v-\tilde b}\geq {\tilde b}^{m_2-1}$$\tilde a(t)\le \tilde b(t)$,那么有

$\begin{eqnarray}\label{5-21} \tilde a(t)^{2m_2-2}\int_\Omega|v-\tilde b|^2{\rm d}x\leq {\tilde b}^{2m_2-2}(t)\int_\Omega|v-\tilde b|^2{\rm d}x\leq\int_\Omega|v^{m_2}-{\tilde b}^{m_2}|^2{\rm d}x.\end{eqnarray}$

如果$\lim_{j\rightarrow\infty}\tilde a(t_j)=\sigma_0>0$,注意到

$\begin{eqnarray}\label{5-22} |\tilde a(t)-\tilde b(t)|^2=\left|\frac{1}{|\Omega|}\int_\Omega(v-\tilde b){\rm d}x\right|^2\leq\frac{1}{|\Omega|}\int_\Omega|v-\tilde b|^2{\rm d}x, \end{eqnarray}$

结合(5.19)式和(5.20)式,可得

$\begin{eqnarray}\label{5-23} \lim\limits_{j\rightarrow\infty}(\tilde b(t_j)-\tilde a(t_j))=0.\end{eqnarray}$

我们注意到

根据(5.16), (5.19)和(5.21)式,可得

$\begin{eqnarray}\label{5-24} \lim\limits_{j\rightarrow\infty}|v(x, t_j)-\sigma_0|^2{\rm d}x=0.\end{eqnarray}$

另一方面,我们发现

这意味着

根据(5.15), (5.16)和(5.19)式,可得

由于$\mathop {\lim }\limits_{j \to \infty } \tilde b({t_j}) = {\sigma _0} > 0$.进一步,我们发现对任意的$m_2>1$,必存在一个正整数$i\ge 0$使得$\beta=2^{i}\ge m_2$,那么

因此,我们有

$\begin{eqnarray} \lim\limits_{j\rightarrow\infty}\int_\Omega(v(x, t_j)-A_2k)^2{\rm d}x=0.\end{eqnarray}$

根据(5.22)和(5.23)式,并且注意到

那么我们有$\sigma_0=A_2k$.证毕.

从上面引理可知,如果$\displaystyle\lim_{j\rightarrow\infty}a(t_j)$存在,那么等于$A_2k$$0$.接下来,我们证明$\displaystyle\lim_{j\rightarrow\infty}a(t_j)=0$是不可能的.利用文献[4]的思路,我们有

引理5.7  假设$\lambda>F_0(\bar{u})$,并且令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.如果$\mathop {\lim }\limits_{t \to \infty } {\left\| {v( \cdot ,t)} \right\|_{{L^1}}} = 0$,那么

  根据(1.4)式,并注意到

那么对任意固定的$r>1$,有

因此,对任意的足够小的常数$\varepsilon_0>0$,存在一个常数$T>0$,使得

$\begin{eqnarray} \|v(\cdot, t)\|_{L^{3m_2}}<\varepsilon_0.\end{eqnarray}$

令一个切断函数$\eta(t)\in C^1[T, +\infty)$,对于$t>T+1$满足$\eta(T)=0$, $0\le\eta(t)\le 1$, $|\eta'(t)|\leq C$$\eta(t)=1$.对任意的$r>3m_2$,在模型(1.1)的第三个方程两边同时乘以$r\eta^rv^{r-1}$,然后在$\Omega$上积分,利用(1.4)式,可得

$\begin{eqnarray}&& \frac{\rm d}{{\rm d}t}\displaystyle\int_\Omega\eta^rv^r{\rm d}x+m_2r(r-1)\int_\Omega\eta^rv^{m_2+r-3}|\nabla v|^2{\rm d}x+\int_\Omega\eta^rv^r{\rm d}x \nonumber\\& \le& (\lambda r+1)\displaystyle\int_\Omega\eta^rv^r{\rm d}x+r\int_\Omega\eta'\eta^{r-1}v^r{\rm d}x \nonumber\\& \le& Cr\displaystyle\int_\Omega\eta^{r-1}v^r{\rm d}x.\end{eqnarray}$

注意到对任意的足够小的$\sigma>0$,

使$\sigma$适当地小,并且将上面的两个不等式代入到(5.25)式.我们有

$\begin{eqnarray} \frac{\rm d}{{\rm d}t}\int_\Omega\eta^rv^r{\rm d}x+\int_\Omega\eta^rv^r{\rm d}x\leq C_3r^\frac{5}{2}\|\eta v\|_{\frac{r}{2}}^{r-3m_2}+C_2r\|\eta v\|_\frac{r}{2}^{r-1}.\end{eqnarray}$

$r_j=2r_{j-1}=2^jr_0, r_0=4m_2, M_j=\sup\limits_{t\in(T, \infty)}\|\eta v\|_{L^{r_j}}$,并且注意到$\eta(T)=0$.那么根据一个直接计算,得到如下结论

根据一个迭代过程,我们发现

接下来证明$S_1=\prod\limits_{k=1}^\infty\frac{{r_0}^{2^k}-3m_2}{{r_0}^{2^k}}>0$, $S_2=\prod\limits_{k=1}^\infty\frac{{r_0}^{2^k}-1}{{r_0}^{2^k}}>0$.注意到$\ln\frac{1}{S_1}=\sum\limits_{k=1}^\infty\ln\left(\frac{{r_0}^{2^k}}{{r_0}^{2^k}-3m_2}\right)$, $\ln\frac{1}{S_2}=\sum\limits_{k=1}^\infty\ln\left(1+\frac{1}{{r_0}^{2^k}-1}\right)$,显然, $\sum\limits_{k=1}^\infty\ln\left(1+\frac{3m_2}{{r_0}^{2^k}-3m_2}\right)$$\sum\limits_{k=1}^\infty\ln\left(1+\frac{1}{{r_0}^{2^k}-1}\right)$收敛,这意味着$S_1, S_2>0$.那么

$j\rightarrow\infty$,我们有

证毕.

引理5.8   假设$\lambda>F_0(\bar{u})$,并且令$(u, w, v)\in{\cal X}_1\times{\cal X}_2\times{\cal X}_3$是模型(1.1)的一个解.那么

$\begin{equation} \lim\limits_{t\rightarrow\infty}\int_\Omega|v(x, t)-A_2k|^2{\rm d}x=0.\end{equation}$

  我们首先证明

$\begin{equation} \lim\limits_{t\rightarrow\infty}\tilde a(t)\neq0. \end{equation}$

另一方面,根据引理5.7,我们有

那么存在$T>0$,使得对任意的$t>T$,在$\Omega$上处处有$v\leq\frac{A_2k}{2}$,那么(5.9)式两端同时乘以$\frac{1}{v}$,然后在$\Omega$上积分,对任意的$t>T$,我们有

这意味着当$t\rightarrow\infty$

根据Jensen's不等式,我们发现当$t\rightarrow\infty$

与前面矛盾.因此,我们有(5.28)式.接下来,我们证明

假设相反,根据(5.28)式和引理5.6,存在两个序列$\{t_j^{(1)}\}_j$$\{t_j^{(2)}\}_j$满足当$j\rightarrow\infty$$t_j^{(1)}$, $t_j^{(2)}\rightarrow\infty$,使得

$ \begin{equation}\label{5-31} \lim\limits_{t\rightarrow\infty}\tilde a(t_j^{(1)})=0, \ \ \lim\limits_{t\rightarrow\infty}\tilde a(t_j^{(2)})=A_2k, \end{equation} $

由于$\tilde a(t)$是有界的.注意到

这意味着$\tilde a(t)$是Lipschitz连续.那么根据(5.29)式和中值定理,那么存在$\{t_j^{(3)}\}_j$满足当$j\rightarrow\infty$$t_j^{(3)}\rightarrow\infty$,使得

这与引理5.6矛盾.因此(5.27)式是引理5.6的一个直接推论.证毕.

定理1.3的证明  如果$p<2$,那么

注意到对任意的$p\ge 2$

并且利用引理5.1,引理5.2和引理5.8,我们完成了定理1.3的证明.

参考文献

Aiseba B , Bendahmane M , Noussair A .

A reaction-diffusion system modeling predator-prey with pre-taxis

Nonlinear Anal Real World Anal, 2008, 9: 2086- 2105

DOI:10.1016/j.nonrwa.2007.06.017      [本文引用: 1]

Gao X , Zhou J , Tian M .

Global boundedness and asymptotic behavior for an attraction-repulsion chemotaxis system with logistic source

Acta Math Sci, 2017, 37A: 113- 121

[本文引用: 1]

Jin C .

Global classical solution and boundedness to a chemotaxis-haptotaxis model with re-establishment mechanism

Bull London Math Soc, 2018, 50: 598- 618

DOI:10.1112/blms.12160     

Jin C .

Large time behavior of solutions to a chemotaxis model with porous medium diffusion

J Math Anal Appl, 2019, 478: 195- 211

DOI:10.1016/j.jmaa.2019.05.027      [本文引用: 3]

Jin H , Wang Z .

Global stability of prey-taxis systems

J Differential Equations, 2017, 262: 1257- 1290

DOI:10.1016/j.jde.2016.10.010      [本文引用: 1]

Kareiva P , Odell G .

Swarms of predators exhibit "prey-taxis" if individual predators use area-restricted search

Amer Nat, 1987, 130: 233- 270

DOI:10.1086/284707      [本文引用: 2]

Lee J M , Hillen T , Lewis M A .

Pattern formation in prey-taxis systems

J Biol Dynam, 2009, 3: 551- 573

DOI:10.1080/17513750802716112      [本文引用: 1]

Stinner C , Surulescu C , Winkler M .

Global weak solutions in a PDE-ODE system modeling multiscale cancer cell invasion

SIAM J Math Anal, 2014, 46: 1969- 2007

DOI:10.1137/13094058X     

Sugiyama Y .

Time global existence and asymptotic behavior of solutions to degenerate quasilinear parabolic systems of chemotaxis

Differential Integral Equations, 2007, 20: 133- 180

[本文引用: 1]

Tao Y .

Global existence of classical solutions to a predator-prey model with nonlinear prey-taxis

Nonlinear Anal Real World Anal, 2010, 11: 2056- 2064

DOI:10.1016/j.nonrwa.2009.05.005      [本文引用: 1]

Tello J , Wrzosek D .

Predator-prey model with diffusion and indirect prey-taxis

Mathematical Models and Methods in Applied Science, 2016, 11: 2129- 2162

URL     [本文引用: 2]

Tao Y , Winkler M .

Global existence and boundedness in a Keller-Segel-Stokes model with arbitrary porous medium diffusion

Discrete Contin Dyn Syst, 2012, 32: 1901- 1914

DOI:10.3934/dcds.2012.32.1901      [本文引用: 1]

Wu S , Shi J , Wu B .

Global existence of solutions and uniform persistence of a diffusive predator-prey model with prey-taxis

J Differential Equations, 2016, 260: 5847- 5874

DOI:10.1016/j.jde.2015.12.024      [本文引用: 1]

Wang X , Wang W , Zhang G .

Global bifurcation of solutions for a predator-prey model with prey-taxis

Math Methods Appl Sci, 2015, 38: 431- 443

DOI:10.1002/mma.3079      [本文引用: 1]

Winkler M. Aggregation vs. global diffusive behavior in the higher-dimensional Keller-Segel model. J Differential Equations, 2010, 248: 2889-2905

[本文引用: 1]

Winkler M .

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

J Math Anal Appl, 2008, 348: 708- 729

DOI:10.1016/j.jmaa.2008.07.071      [本文引用: 1]

Zheng P , Mu C , Hu X .

Persistence property in a two-species chemotaxis system with two signals

J Math Phys, 2017, 58 (11): 111501

DOI:10.1063/1.5010681      [本文引用: 1]

Wu Z Q , Yin J X , Wang C P . Introduction to Elliptic and Parabolic Equations. Beijing: Science Press, 2003

/