高维空间中具有间接信号产出生物趋化模型解的全局有界性 Global Boundedness of Classical Solutions to a High-Dimensional Chemotaxis Model with Indirect Signal Production

doi:10.12677/AAM.2019.84088

Advances in Applied Mathematics
Vol. 08 No. 04 ( 2019 ), Article ID: 29985 , 10 pages
10.12677/AAM.2019.84088

Global Boundedness of Classical Solutions to a High-Dimensional Chemotaxis Model with Indirect Signal Production

Dongdong Zhang, Qiao Xin^*, Jiangang Tang

●How to Cite this Article

College of Mathematics and Statistics, Yili University, Yining Xinjiang

Received: Apr. 6^th, 2019; accepted: Apr. 21^st, 2019; published: Apr. 28^th, 2019

ABSTRACT

This paper deals with the chemotaxis system of diffusion and aggregation of Mountain Pine Beetle with indirect attractant production and generalized logistic source  ${\begin{cases} u_{t} = Δ u - \nabla \cdot (u χ (v) \nabla v) + f (u) x \in Ω, t > 0, \\ v_{t} = Δ v - v + u g (u), x \in Ω, t > 0, \\ τ w_{t} + w = u, x \in Ω, t > 0, \end{cases}$ in a smoothly bounded domain with homogeneous Neumann boundary conditions and nonnegative initial values, the chemotactic sensitivity function satisfies $χ (v) \leq χ_{0} / {(1 + β v)}^{δ} (δ > 1),$ the consumption function $g (u) = h_{0} / {(1 + h u)}^{υ} (υ \geq 0)$ and the logistic source satisfies $f (u) \leq a u - b u^{2} (a, b > 0) .$ Moreover, $τ, β, χ_{0}, h_{0}$ are given positive parameter. Firstly, the energy estimation method and Gagliardo-Nirenberg inequality are used to establish the local prior estimate of u and w, and then Moser iteration is used. This problem admits a unique global classical solution that is uniformly in-time bounded.

Keywords:Chemoaxis, Globl Existence, Boundedness, Indirect Signal Production

高维空间中具有间接信号产出生物趋化模型解的全局有界性

张冬冬，辛巧^*，汤建钢

伊犁师范大学数学与统计分院，新疆伊宁

收稿日期：2019年4月6日；录用日期：2019年4月21日；发布日期：2019年4月28日

摘要

考虑一个高维空间中描述具有间接信号产出的山地松甲虫扩散和聚集模型 ${\begin{cases} u_{t} = Δ u - \nabla \cdot (u χ (v) \nabla v) + f (u) x \in Ω, t > 0, \\ v_{t} = Δ v - v + u g (u), x \in Ω, t > 0, \\ τ w_{t} + w = u, x \in Ω, t > 0, \end{cases}$ 在齐次Neumann初值边界和非负初值条件下经典解的整体性态。假设 $Ω \in R^{n} (n \geq 1)$ 是一个光滑有界区域， $χ (v) = χ_{0} / {(1 + β v)}^{2}$ 表示趋化敏感函数且满足 $χ (v) \leq χ_{0} / （ 1 + β v ）^{δ} （ δ > 1 ）$ ，消耗函数为 $g (u) = h_{0} / {(1 + h u)}^{υ} (υ \geq 0)$ ，Logistic源满足 $f (u) \leq a u - b u^{2} (a, b > 0)$ ， $τ, β, χ_{0}, h_{0}$ 是正的参数。先利用能量估计方法及Gagliardo-Nirenberg不等式建立u和w的先验估计，再运用Moser迭代证明了此问题存在经典解且一致有界，不存在坍塌现象。

关键词 :趋化性，全局存在性，有界性，间接信号产出

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

1. 引言

生物趋化性为生物种群在某种特定的生物环境中在化学信号物质的影响下做定向运动，关于生物趋化性的研究对生物学领域有很重要的意义。在现实生活中运用也很广泛，如生物除污、伤口的愈合等等。最早的生物趋化模型是由Patlak和Keller-Segel提出，最简单的K-S趋化模型 [1] 如下

${\begin{cases} u_{t} = Δ u - χ \nabla \cdot (u \nabla v), x \in Ω, t > 0, \\ v_{t} = Δ v - v + u, x \in Ω, t > 0, \end{cases}$ (1.1)

其中 $Ω \in R^{n}$ 表示一个光滑有界区域， $u = u (x, t)$ 表示细胞的种群密度， $v = v (x, t)$ 表示化学物质的浓度， $χ \in R$ 表示趋化敏感系数， $- χ \nabla \cdot (u \nabla v)$ 表示趋化灵敏度，体现趋化模型最本质的一项。数学家最注重研究趋化模型解的性态是否具有坍塌行为，即模型解在有限时间内或无限时间内是否爆破，此类问题的研究已经有了许多成果。当 $n = 1$ 时，Yagi等人证明了此模型解全局存在，不存在坍塌现象 [2] 。当 $n = 2$ 时，存在临界常数 $m_{c} = 4 π / χ$ ，Nagai和Senba证明了如果初值满足 ${‖ u_{0} ‖}_{L^{1} (Ω)} < m_{c}$ 时，则此模型解存在且一致有界 [3] ，Horstmann和Wang证明了如果初值满足 ${‖ u_{0} ‖}_{L^{1} (Ω)} > m_{c}$ ，则此模型解在有限时间内或无限时间内发生爆破 [4] 。进一步，Nagai证明了满足 ${‖ u_{0} ‖}_{L^{1} (Ω)} = m_{c}$ 时，此模型经典解也在有限时间内发生爆破 [5] 。当 $n \geq 3$ 时，Winkler在文献 [6] 中证明了若对于任意的 $δ > 0$ ，存在 $ε > 0$ 使得 ${‖ u_{0} ‖}_{L^{n / 2 + δ} (Ω)} < ε$ 和 ${‖ \nabla v_{0} ‖}_{L^{n + δ} (Ω)} < ε$ 成立，此模型解存在且一致有界，而当初值满足 ${‖ u_{0} ‖}_{L^{1} (Ω)} > 0$ 时，文献 [7] 利用Lyapunov泛函证明了此模型解会在有限时间内存在坍塌现象。

受文献 [8] [9] 工作的启发，本文研究具有间接信号产出的带有一般的趋化敏感函数的趋化模型，具体模型如下：

${\begin{cases} u_{t} = Δ u - \nabla \cdot (u χ (v) \nabla v) + f (u) x \in Ω, t > 0, \\ v_{t} = Δ v - v + u g (u), x \in Ω, t > 0, \\ τ w_{t} + w = u, x \in Ω, t > 0, \\ \frac{\partial u}{\partial ν} = \frac{\partial v}{\partial ν} = \frac{\partial w}{\partial ν} = 0, x \in \partial Ω, t > 0, \\ u (x, 0) = u_{0} (x), v (x, 0) = v_{0} (x), w (x, 0) = w_{0} (x), x \in Ω, \end{cases}$ (1.2)

其中是一个光滑有界区域， $u = u (x, t)$ 表示飞行的山地松甲虫的种群密度， $w = w (x, t)$ 表示做窝的山地松甲虫的种群密度， $v = v (x, t)$ 表示化学信号的浓度。当 $n = 2$ 且 $χ (v) = 1$ 时，王笑丹等人在文献 [10] 中证明了模型解全局存在且一致有界。当 $n \geq 2$ 且 $χ (v) = 1$ ，logistic源项为 $u - u^{α} (α > n / 2)$ 时，Li和Tao证明了模型解全局存在且一致有界 [11] 。但是当 $χ (v) = 1, g (u) = 1$ ，此模型的解的性态不明确。本文，主要考虑一般的趋化敏感函数情况，假设趋化敏感函数 $χ (v)$ 和函数 $f (u), g (u)$ 满足

$χ (v) \leq χ_{0} / {(1 + β v)}^{δ}, δ > 1, χ_{0} > 0,$ (1.3)

$f (u) \leq a u - b u^{2}, a, b > 0, f (0) = 0,$ (1.4)

$g (u) = h_{0} / {(1 + h u)}^{υ}, h_{0}, h, υ \geq 0.$ (1.5)

与高维K-S模型解存在有限时间坍塌相比，模型(1.2)的解不会在有限时间内爆破。本文对于具有一般的趋化敏感函数的趋化模型(1.2)经典解的全局有界性结论如下：

定理1：假设非负初值 $(u_{0}, v_{0}, w_{0})$ 满足 $u_{0} \in C^{0} (\bar{Ω})$ 、 $v_{0} \in W^{1, q} (\bar{Ω}) (q > n)$ 和 $w_{0} \in C^{0} (\bar{Ω})$ ，设 $χ (v), f (u), g (u)$ 满足条件(1.3)~(1.5)，则生物趋化模型(1.2)有唯一、非负经典解且关于时间一致有界。

2. 预备知识

为了证明定理1，先给出一些相关的引理。

引理2.1：假设非负初值 $(u_{0}, v_{0}, w_{0})$ 满足 $u_{0} \in C^{0} (\bar{Ω})$ 、 $v_{0} \in W^{1, q} (\bar{Ω}) (q > n)$ 和 $w_{0} \in C^{0} (\bar{Ω})$ 。则存在 $T_{\max} \in [0, \infty)$ 和非负函数 $(u, v, w)$ 满足

$\begin{array}{l} u \in C^{0} (\bar{Ω} \times [0, T_{\max})) \cap C^{2, 1} (\bar{Ω} \times [0, T_{\max})), \\ v \in C^{0} (\bar{Ω} \times [0, T_{\max})) \cap L_{l o c}^{\infty} ([0, T_{\max}), W^{1, q} (Ω)) \cap C^{2, 1} (\bar{Ω} \times [0, T_{\max})), \\ w \in C^{0} (\bar{Ω} \times [0, T_{\max})) \cap C^{2, 1} (\bar{Ω} \times [0, T_{\max})), \end{array}$

使得 $(u, v, w)$ 是模型(1.2)在 $\bar{Ω} \times [0, T_{\max})$ 上唯一的经典解。而且 $T_{\max} < + \infty$ ，则

${‖ u (x, t) ‖}_{L^{\infty} (Ω)} + {‖ v (x, t) ‖}_{w^{1, q} (Ω)} + {‖ w (x, t) ‖}_{L^{\infty} (Ω)} \to \infty, t \to T_{\max} .$

证明：模型(1.2)经典解的局部存在性可利用不动点定理得到，详见文献 [12] [13] [14] 。

引理2.2：假设条件(1.3)~(1.5)成立，则存在常数 $M_{0} > 0$ ，使得模型(1.2)的解 $(u, v, w)$ 满足

$\int_{Ω} u (\cdot, t) d x \leq M_{0}, t \in (0, T_{\max}),$ (2.1)

$\int_{Ω} w (\cdot, t) d x \leq M_{1}, t \in (0, T_{\max}),$ (2.3)

证明：对于模型(1.2)的第一个方程两边在 $Ω$ 上积分，容易得出(2.1)式。对于模型(2.2)的第三个方程两边在 $Ω$ 上积分，可得

$\frac{d}{d t} \int_{Ω} τ w d x + \int_{Ω} w d x = \int_{Ω} u d x,$

那么，存在常数 $M_{1} > 0$ ，解得

$\int_{Ω} w (\cdot, t) d x \leq M_{1}, t \in (0, T_{\max}) .$

接下来，引进一个有用的不等式，详见证明过程参考文献 [3] [14] 。

引理2.3：假设则对所有的 $l^{*} \in (2, \infty)$ 满足 $l^{*} \leq α^{*}$ ， $h \in (0, 2)$ 和 $α \in [h, l^{*}]$ ，则存在常数 $C_{G N} > 0$ 使得

${‖ ψ ‖}_{L^{α} (Ω)} \leq C_{G N} ({‖ \nabla ψ ‖}_{L^{2} (Ω)}^{λ^{*}} {‖ ψ ‖}_{L^{h} (Ω)}^{1 - λ^{*}} + {‖ ψ ‖}_{L^{h} (Ω)}),$

其中 $ψ \in W^{1, 2} (Ω)$ 且

$λ^{*} = \frac{\frac{n}{h} - \frac{n}{α}}{1 - \frac{n}{2} + \frac{n}{h}} .$

3. 解的全局有界性

为了证明定理1，需要先建立 ${‖ u (\cdot, t) ‖}_{L^{k} (Ω)}$ 与 ${‖ w (\cdot, t) ‖}_{L^{k} (Ω)}$ 的有界性，下面给出一些相关的引理及其证明。

引理3.1：假设条件(1.3)~(1.5)成立，则存在常数 $M_{2}, M_{3} > 0$ ， $k > 1$ ，使得对所有的 $t \in (0, T_{\max})$ 有

${‖ u (\cdot, t) ‖}_{L^{k} (Ω)} \leq M_{2}, t \in (0, T_{\max}),$ (3.1)

${‖ w (\cdot, t) ‖}_{L^{k} (Ω)} \leq M_{3}, t \in (0, T_{\max}),$ (3.2)

证明：首先定义函数

$ϕ (s) = \exp ({(1 + β^{*} s)}^{- r}), s > 0.$

设 $k > 1$ ， $r > 1$ 且

$r \leq \min {k - 1 / 8 k, 2 δ - 2},$ (3.3)

和

$r \leq k (\frac{k - 1}{e C_{G N}^{2} k^{2} M^{k (1 - λ)}} - a),$ (3.4)

以及

$β^{*} > χ_{0} \sqrt{2 p (p - 1) / r (r + 1)},$ (3.5)

通过直接计算，容易得到

$\begin{array}{l} = - (k - 1) \int_{Ω} u^{k - 2} ϕ (v) {| \nabla u |}^{2} d x - \int_{Ω} u^{k - 1} ϕ^{'} (v) \nabla u \cdot \nabla v d x \\ + (k - 1) \int_{Ω} u^{k - 1} ϕ (v) χ (v) \nabla u \cdot \nabla v d x + \int_{Ω} u^{p} ϕ^{'} (v) χ (v) {| \nabla v |}^{2} d x \\ + \frac{1}{k} \int_{Ω} u^{k} ϕ (v) f (u) d x - \int_{Ω} u^{k - 1} ϕ^{'} (v) \nabla u \cdot \nabla v d x - \frac{1}{k} \int_{Ω} u^{k} ϕ^{″} (v) {| \nabla v |}^{2} d x \\ - \frac{1}{k} \int_{Ω} u^{k} v ϕ^{'} (v) d x + \frac{1}{k} \int_{Ω} u^{k} ϕ^{'} (v) g (u) d x \end{array}$

因为 $f (s) \leq a s - b s^{2}$ ， $g (s) \geq 0$ 和 $χ (s) \geq 0$ 以及 $ϕ^{'} (s) \leq 0, s \geq 0$ 。那么可以得出

$\begin{array}{l} \frac{1}{k} \frac{d}{d t} \int_{Ω} u^{k} ϕ (v) d x + (k - 1) \int_{Ω} u^{k - 2} ϕ (v) {| \nabla u |}^{2} d x + \frac{1}{k} \int_{Ω} u^{k} ϕ^{″} (v) {| \nabla v |}^{2} d x \\ \leq - 2 \int_{Ω} u^{k - 1} ϕ^{'} (v) \nabla u \cdot \nabla v d x + (k - 1) \int_{Ω} u^{k - 1} ϕ (v) χ (v) \nabla u \cdot \nabla v d x \\ - \frac{1}{k} \int_{Ω} u^{k} v ϕ^{'} (v) d x + a \int_{Ω} u^{k} ϕ (v) d x, \end{array}$ (3.6)

其中

$- \frac{1}{k} \int_{Ω} u^{k} v ϕ^{'} (v) d x \leq \frac{r}{k} \int_{Ω} u^{k} ϕ (v) d x,$

此外，利用杨不等式估计和(1.3)式得出

$- 2 \int_{Ω} u^{k - 1} ϕ^{'} (v) \nabla u \cdot \nabla v d x \leq \frac{k - 1}{4} \int_{Ω} u^{k - 2} ϕ (v) {| \nabla u |}^{2} d x + \frac{4}{k - 1} \int_{Ω} u^{k} \frac{{ϕ^{'}}^{2} (v)}{ϕ (v)} {| \nabla v |}^{2} d x,$

和

$\begin{array}{l} (k - 1) \int_{Ω} u^{k - 1} ϕ (v) χ (v) \nabla u \cdot \nabla v d x \\ \leq \frac{k - 1}{4} \int_{Ω} u^{k - 2} ϕ (v) {| \nabla u |}^{2} d x + (k - 1) \int_{Ω} u^{k} ϕ (v) χ^{2} (v) {| \nabla v |}^{2} d x \\ \leq \frac{k - 1}{4} \int_{Ω} u^{k - 2} ϕ (v) {| \nabla u |}^{2} d x + (k - 1) χ_{0}^{2} \int_{Ω} u^{k} {(1 + β v)}^{- 2 δ} ϕ (v) {| \nabla v |}^{2} d x, \end{array}$

那么，(3.6)式可以写成

$\begin{array}{l} \frac{1}{k} \frac{d}{d t} \int_{Ω} u^{k} ϕ (v) d x + \frac{k - 1}{2} \int_{Ω} u^{k - 2} ϕ (v) {| \nabla u |}^{2} d x + \frac{1}{k} \int_{Ω} u^{k} ϕ^{″} (v) {| \nabla v |}^{2} d x \\ \leq \frac{4}{k - 1} \int_{Ω} u^{k} \frac{{ϕ^{'}}^{2} (v)}{ϕ (v)} {| \nabla v |}^{2} d x + (k - 1) χ_{0}^{2} \int_{Ω} u^{k} {(1 + β v)}^{- 2 δ} ϕ (v) {| \nabla v |}^{2} d x \\ + (\frac{r}{k} + a) \int_{Ω} u^{p} ϕ (v) d x, \end{array}$ (3.7)

根据题意，接下来处理(3.7)式中不等式右端的前两项，我们计算

$Z_{1} (s) : = \frac{4}{k - 1} \frac{{ϕ^{'}}^{2} (s)}{ϕ (s)} = \frac{4}{k - 1} \cdot {(β^{*})}^{2} r^{2} {(1 + β^{*} s)}^{- 2 r - 2} e^{{(1 + β^{*} s)}^{- r}},$

$Z_{2} (s) : = (k - 1) χ_{0}^{2} {(1 + β s)}^{- 2 δ} ϕ (s) = (k - 1) χ_{0}^{2} {(1 + β s)}^{- 2 δ} e^{{(1 + β^{*} s)}^{- r}}$

和

$Z_{3} (s) : = \frac{1}{k} ϕ^{″} (s) = \frac{1}{k} \cdot {(β^{*})}^{2} r (r + 1) {(1 + β^{*} s)}^{- r - 2} e^{{(1 + β^{*} s)}^{- r}} + \frac{1}{k} \cdot {(β^{*})}^{2} r^{2} {(1 + β^{*} s)}^{- 2 r - 2} e^{{(1 + β^{*} s)}^{- r}},$

那么有

$\begin{matrix} \frac{Z_{1} (s)}{\frac{1}{2} Z_{3} (s)} \leq \frac{\frac{4}{k - 1} \cdot {(β^{*})}^{2} r^{2} {(1 + β^{*} s)}^{- 2 r - 2} e^{{(1 + β^{*} s)}^{- r}}}{\frac{1}{2} \cdot \frac{1}{k} \cdot {(β^{*})}^{2} r (r + 1) {(1 + β^{*} s)}^{- r - 2} e^{{(1 + β^{*} s)}^{- r}}} \\ = \frac{8 k r}{(k - 1) (r - 1)} \cdot {(1 + β^{*} s)}^{- r} \leq \frac{8 k r}{k - 1} \leq 1, \end{matrix}$

根据(3.3)式，可知上式显然成立。

$\begin{matrix} \frac{Z_{2} (s)}{\frac{1}{2} Z_{3} (s)} \leq \frac{(k - 1) χ_{0}^{2} {(1 + β s)}^{- 2 δ} e^{{(1 + β^{*} s)}^{- r}}}{\frac{1}{2} \cdot \frac{1}{k} \cdot {(β^{*})}^{2} r (r + 1) {(1 + β^{*} s)}^{- r - 2} e^{{(1 + β^{*} s)}^{- r}}} \\ = \frac{2 k (k - 1) χ_{0}^{2}}{{(β^{*})}^{2} r (r + 1)} \cdot {(1 + β s)}^{- 2 δ} {(1 + β^{*} s)}^{r + 2}, \end{matrix}$

定义函数 $φ (s) : = {(1 + β s)}^{- 2 δ} {(1 + β^{*} s)}^{r + 2}, s \geq 0$ 。

$\begin{matrix} φ^{'} (s) = {(1 + β s)}^{- 2 δ - 1} \cdot {(1 + β^{*} s)}^{r + 1} \cdot {- 2 k β (1 + β^{*} s) + (r + 2) (1 + β s)} \\ \leq {(1 + β s)}^{- 2 δ - 1} \cdot {(1 + β^{*} s)}^{r + 1} \cdot {- (2 δ β - r - 2) - [2 δ β β^{*} - (r + 2) β] s} \\ \leq 0, s \geq 0, \end{matrix}$

所以 $φ (s) \leq φ (0) = 1, s \geq 0$ 。

根据(3.5)式，有

$\frac{Z_{2} (s)}{\frac{1}{2} Z_{3} (s)} \leq \frac{2 k (k - 1) χ_{0}^{2}}{{(β^{*})}^{2} r (r + 1)} \leq 1, s \geq 0,$

则(3.7)式可改写为

$\frac{1}{k} \frac{d}{d t} \int_{Ω} u^{k} ϕ (v) d x + \frac{k - 1}{2} \int_{Ω} u^{k - 2} ϕ (v) {| \nabla u |}^{2} d x \leq (a + \frac{r}{k}) \int_{Ω} u^{k} ϕ (v) d x ， t \in (0, T_{\max}),$

那么根据引理2.3不等式和 $ϕ (s) \leq e, s \geq 0$ 得

$\begin{matrix} \int_{Ω} u^{k} ϕ (v) d x \leq e \int_{Ω} u^{k} d x = e {‖ u^{\frac{k}{2}} ‖}_{L^{2} (Ω)}^{2} \leq e \cdot C_{G N} {‖ u^{\frac{k}{2}} ‖}_{W^{1 ， 2} (Ω)}^{2 h} \cdot {‖ u^{\frac{k}{2}} ‖}_{L^{\frac{2}{k}} (Ω)}^{2 (1 - h)} \\ \leq e \cdot c_{G N} \cdot {(c (\frac{2}{k}))}^{2 h} {({‖ \nabla u^{\frac{k}{2}} ‖}_{L^{2} (Ω)} + {‖ u^{\frac{k}{2}} ‖}_{L^{\frac{2}{k}} (Ω)})}^{2 h} \cdot {‖ u^{\frac{k}{2}} ‖}_{L^{\frac{2}{k}} (Ω)}^{2 (1 - h)}, \end{matrix}$

其中 $C_{G N} > 0$ 和

$h = \frac{\frac{2 n}{k} - \frac{n}{2}}{\frac{2 n}{k} + 1 - \frac{n}{2}} \in (0, 1),$

根据引理2.2的性质可得

$\int_{Ω} u^{k} ϕ (v) d x \leq c_{1} {({‖ \nabla u^{\frac{k}{2}} ‖}^{2}_{L^{2} (Ω)} + 1)}^{h}, c_{1} > 0,$

因为 $ϕ (s) \geq 1, s \geq 0$ ，可知

$\begin{matrix} \frac{k - 1}{2} \int_{Ω} u^{k - 2} ϕ (v) {| \nabla u |}^{2} d x \geq \frac{k - 1}{2} \int_{Ω} u^{k - 2} {| \nabla u |}^{2} d x = \frac{2 (k - 1)}{k^{2}} \int_{Ω} {| \nabla u^{\frac{k}{2}} |}^{2} d x \\ \geq \frac{2 (k - 1)}{k^{2} c_{1}^{\frac{1}{h}}} {(\int_{Ω} u^{k} ϕ (v) d x)}^{\frac{1}{h}} - \frac{2 (k - 1)}{k^{2}}, \end{matrix}$

那么可以得出

$\begin{array}{l} \frac{1}{k} \frac{d}{d t} \int_{Ω} u^{k} ϕ (v) d x \\ \leq - \frac{2 (k - 1)}{k^{2} c_{1}^{\frac{1}{h}}} {(\int_{Ω} u^{k} ϕ (v) d x)}^{\frac{1}{h}} + \frac{2 (k - 1)}{k^{2}} + (a + \frac{r}{k}) \int_{Ω} u^{k} ϕ (v) d x, t \in (0, T_{\max}), \end{array}$

容易解得这个常微分方程

${‖ u (\cdot, t) ‖}_{L^{k} (Ω)} \leq M_{2}, t \in (0, T_{\max}),$

下面证明(3.2)式，模型(1.2)中第三个方程同乘 $w^{k - 1}$ 再在 $Ω$ 上积分得

$\frac{τ}{p} \frac{d}{d t} \int_{Ω} w^{k} d x + \int_{Ω} w^{k} d x = \int_{Ω} u w^{k - 1} d x,$

利用杨不等式可得

$\int_{Ω} u w^{k - 1} d x \leq \frac{1}{4} \int_{Ω} w^{k} d x + C_{k} \int_{Ω} u^{k} d x, C_{k} > 0,$

那么存在常数 $M_{3} > 0$ ，使得下式成立

$\frac{τ}{k} \frac{d}{d t} \int_{Ω} w^{k} d x + \frac{3}{4} \int_{Ω} w^{k} d x \leq M_{3},$

容易得出

${‖ w (\cdot, t) ‖}_{L^{k} (Ω)} \leq M_{3}, t \in (0, T_{\max}),$

4. 定理1的证明

由引理2.1可知，对于 $τ \in (0, T_{\max})$ ，存在 $c (τ) > 0$ 可得

${‖ u (x, t) ‖}_{L^{\infty} (Ω)} + {‖ v (x, t) ‖}_{L^{\infty} (Ω)} + {‖ w (x, t) ‖}_{L^{\infty} (Ω)} \leq c (τ), t \in (τ, T_{\max}),$

根据热半群理论 [14] ，设 $ρ \in (\frac{1}{2}, 1)$ ， $q > 1$ 使得 $2 ρ - \frac{n}{q} > 1 - \frac{n}{2 p}$ ，容易得出

$v (\cdot, t) = e^{t (Δ - 1)} v_{0} + \int_{0}^{t} e^{(t - s) (Δ - 1)} g (u (\cdot, t)) d s, t \in (0, T_{\max}),$

由算子理论和引理3.1，可以得到

$\begin{matrix} {‖ v (\cdot, t) ‖}_{W^{1, 2 p} (Ω)} \leq c {‖ {(Δ + 1)}^{ρ} v (\cdot, t) ‖}_{L^{q} (Ω)} \\ \leq c t^{- ρ} e^{- υ t} {‖ v_{0} ‖}_{L^{q} (Ω)} + c \int_{0}^{t} {(t - s)}^{- ρ} e^{- υ (t - s)} {‖ u (\cdot, t) ‖}_{L^{q} (Ω)} {‖ g (u (\cdot, t)) ‖}_{L^{q} (Ω)} d s \\ \leq c t^{- ρ} + c \int_{0}^{t} {(t - s)}^{- ρ} e^{- υ (t - s)} d s \leq c t^{- ρ} + c \int_{0}^{\infty} σ^{- ρ} e^{- υ σ} d σ \\ \leq c (t^{- ρ} + 1), t \in (0, T_{\max}), \end{matrix}$

其中c为一个正的常数，那么存在常数，由嵌入定理可得

${‖ v (\cdot, t) ‖}_{L^{\infty} (Ω)} \leq M_{4}, t \in (0, T_{\max}) .$

对于任意的 $m \geq 2$ ，通过模型(1.2)的第一个方程计算，容易得到

$\begin{matrix} \frac{d}{d t} \int_{Ω} u^{m} d x = m \int_{Ω} u^{m - 1} u_{t} d x \\ = m \int_{Ω} u^{m - 1} Δ u d x - m \int_{Ω} u^{m - 1} \nabla \cdot (u χ (v) \nabla v) d x + m \int_{Ω} u^{m - 1} f (u) d x \\ = - m (m - 1) \int_{Ω} u^{m - 2} {| \nabla u |}^{2} d x + m \int_{Ω} u^{m - 1} f (u) d x + m (m - 1) \int_{Ω} u^{m - 1} χ (v) \nabla u \cdot \nabla v d x \\ \leq \frac{- 4 (m - 1)}{m} \int_{Ω} {| \nabla u^{\frac{m}{2}} |}^{2} d x + c m (m - 1) \int_{Ω} u^{m - 1} | \nabla u | d x + a m \int_{Ω} u^{m} d x \\ \leq \frac{- 2 (m - 1)}{m} \int_{Ω} {| \nabla u^{\frac{m}{2}} |}^{2} d x + (\frac{c^{2} m (m - 1)}{2} + a m) \int_{Ω} u^{m} d x, \end{matrix}$

其中 $c > 0$ ，那么有

$\frac{d}{d t} \int_{Ω} u^{m} d x + m (m - 1) \int_{Ω} u^{m} d x \leq \frac{- 2 (m - 1)}{m} \int_{Ω} {| \nabla u^{\frac{m}{2}} |}^{2} d x + c_{1} m (m - 1) \int_{Ω} u^{m} d x, c_{1} > 0,$

其中

$c_{1} m (m - 1) \int_{Ω} u^{m} d x \leq \frac{2 (m - 1)}{m} \int_{Ω} {| \nabla u^{\frac{m}{2}} |}^{2} d x + c_{2} m (m - 1) {(\int_{Ω} u^{\frac{m}{2}} d x)}^{2}, c_{2} > 0,$

那么，可以得到

$\frac{d}{d t} \int_{Ω} u^{m} d x + m (m - 1) \int_{Ω} u^{m} d x \leq c_{2} m (m - 1) {(\int_{Ω} u^{\frac{m}{2}} d x)}^{2}, c_{2} > 0,$

解得

$\int_{Ω} u^{m} d x \leq \int_{Ω} u_{0}^{m} d x + c_{2} \sup_{0 \leq t \leq T_{\max}} {(\int_{Ω} u^{\frac{m}{2}} d x)}^{2}, c_{2} > 0,$

定义函数

$G (m) = \max {{‖ u_{0} ‖}_{L^{\infty} (Ω)}, \sup_{0 \leq t \leq T_{\max}} {(\int_{Ω} u^{m} d x)}^{\frac{1}{m}}} \leq c_{3}^{\frac{1}{m}} G (m / 2),$

设 $m = 2^{i}, i = 1, 2, \dots$ ，那么可知

$G (2^{i}) \leq c_{3}^{2^{- i}} G (2^{i - 1}) \leq c_{3}^{2^{- i} + 2^{- i + 1}} G (2^{i - 2}) \leq \dots \leq c_{3} G (1),$

令 $i \to \infty$ ，可以得到

${‖ u (\cdot, t) ‖}_{L^{\infty} (Ω)} \leq M_{5}, t \in (0, T_{\max}) .$

对于 $w (\cdot, t)$ ，根据 ${‖ u (\cdot, t) ‖}_{L^{\infty} (Ω)} \leq M_{5}, t \in (0, T_{\max})$ ，解得

$\begin{matrix} w (x, t) = \frac{1}{τ} \int_{0}^{t} e^{\frac{1}{τ} (s - t)} u (x, s) d s + w_{0} e^{- \frac{1}{τ} t} \\ \leq \frac{1}{τ} {‖ \int_{0}^{t} e^{\frac{1}{τ} (s - t)} u (x, s) d s ‖}_{L^{\infty} (Ω)} + {‖ w_{0} e^{- \frac{1}{τ} t} ‖}_{L^{\infty} (Ω)} \\ \leq C_{1} \cdot τ (1 - e^{- \frac{1}{τ} t}) + C_{2} e^{- \frac{1}{τ} t} \\ \leq C_{1} \cdot τ + C_{2} : = M_{6}, t \in (0, T_{\max}), \end{matrix}$

即

${‖ w (x, t) ‖}_{L^{\infty} (Ω)} \leq M_{6}, t \in (0, T_{\max}) .$

那么可以得到

${‖ u (x, t) ‖}_{L^{\infty} (Ω)} + {‖ v (x, t) ‖}_{L^{\infty} (Ω)} + {‖ w (x, t) ‖}_{L^{\infty} (Ω)} \leq C, t \in (τ, T_{\max}) .$

5. 结论

综上所述，证明了高维空间中具有间接信号产出生物趋化模型解的全局有界性，但模型(1.2)有数值模拟等问题仍可继续研究。

基金项目

伊犁师范学院研究生科研创新项目《一类生物趋化模型解的存在性和渐近性》 (YLSF2017031)；自治区青年科技创新人才培养项目“偏微分方程理论及其在图像处理中的应用” (2017Q081)。

文章引用

张冬冬,辛巧,汤建钢. 高维空间中具有间接信号产出生物趋化模型解的全局有界性
Global Boundedness of Classical Solutions to a High-Dimensional Chemotaxis Model with Indirect Signal Production[J]. 应用数学进展, 2019, 08(04): 780-789. https://doi.org/10.12677/AAM.2019.84088

参考文献

1. Keller, E.F. and Sege, L.A. (1970) Initiation of Slime Mold Aggregation Viewed as an Instability. Journal of Theoretical Biology, 26, 399-415. https://doi.org/10.1016/0022-5193(70)90092-5

2. Osaki, K. and Yagi, A. (2001) Finite Dimensional Attractors for One-Dimensional Keller-Segel Equations. Funkcialaj Ekvacioj, 44, 441-469.

3. Nagai, T., Senba, T. and Yoshida, K. (1997) Application of the Trudinger-Moser Inequality to a Parabolic System of Chemotaxis. Funkcialaj Ekvacioj, 40, 411-433.

4. Horstmann, D. and Wang, G. (2001) Blow-Up in a Chemotaxis Model without Symmetry Assumptions. European Journal of Applied Mathematics, 12, 159-177. https://doi.org/10.1017/S0956792501004363

5. Nagai, T. (2001) Blow up of Nonradial Solutions to Parabol-ic-Elliptic Systems Modeling Chemotaxis in Two-Dimensional Domains. Journal of Inequalities and Applications, 6, 37-55. https://doi.org/10.1155/S1025583401000042

6. Tao, Y.S. and Winkler, M. (2012) Boundedness in a Quasilinear Parabolic-Parabolic Keller-Segel System with Subcritical Sensitivity. Journal of Differential Equations, 252, 692-715. https://doi.org/10.1155/S1025583401000042

7. Winkler, M. (2013) Finite-Time Blow-Up in the Higher-Dimensional Parabolic-Parabolic Keller-Segel System. Journal de Mathématiques Pures et Appliquées, 100, 748-767. https://doi.org/10.1016/j.matpur.2013.01.020

8. Winkler, M. (2010) Absence of Collapse in a Parabolic Chemotaxis System with Signal-Dependent Sensitivity. Mathematische Narchrichten, 283, 1664-1673.

9. Hu, B. and Tao, Y.S. (2016) To the Exclusion of Blow-Up in a Three-Dimensional Chemotaxis-Growth Model with Indirect Attractant Production. Mathematical Models and Methods in Applied Sciences, 26, 2111-2128. https://doi.org/10.1142/S0218202516400091

10. 王笑丹, 陶有山. 具间接信号产出的二维趋向性模型古典解的整体存在性和有界性[J]. 东华大学学报(自然科学版), 2016, 42(6): 922-930.

11. Li, H.Y. and Tao, Y.S. (2018) Boundedness in a Chemotaxis System with Indirect Signal Production and Generalized Logistic Source. Applied Mathematics Letters, 77, 108-113. https://doi.org/10.1016/j.aml.2017.10.006

12. Li, Y. (2015) Global Bounded Solutions and Their Asymptotic Properties under Small Initial Data Condition in a Two-Dimensional Chemotax Is System for Two Species. Journal of Mathematical Analysis and Applications, 429, 1291-1304. https://doi.org/10.1016/j.jmaa.2015.04.052

13. Cho, Y.S. and Wang, Z.A. (2010) Prevention of Blow-Up by Fast Diffusion in Chemotaxis. Journal of Mathematical Analysis and Applications, 362, 553-564. https://doi.org/10.1016/j.jmaa.2009.08.012

14. Delgado, M., Gayte, I., Morales-rodrigo, C. and Suarez, A. (2010) An Angiogenesis Model with Nonlinear Chemotactic Response and Flux at the Tumor Boundary. Nonlinear Analysis: Theory, Methods & Applications, 72, 330-347. https://doi.org/10.1016/j.na.2009.06.057

NOTES

^*通讯作者。

期刊菜单