一类具有低正则初值的趋化–流体耦合模型解的整体存在性 Global Existence of Solutions for a Chemotaxis-Fluid Coupling Model with Low Regular Initial Data

doi:10.12677/PM.2023.1311341

Pure Mathematics
Vol. 13 No. 11 ( 2023 ), Article ID: 76236 , 14 pages
10.12677/PM.2023.1311341

一类具有低正则初值的趋化–流体耦合模型解的整体存在性

陈娟，韩永杰

●How to Cite this Article

西华大学理学院，四川成都

收稿日期：2023年10月15日；录用日期：2023年11月16日；发布日期：2023年11月24日

摘要

考虑在三维有界区域上带有logistic源的具有信号消耗机制的趋化–流体耦合方程组：

${\begin{cases} n_{t} + u \cdot \nabla n = Δ n - \nabla \cdot (n \nabla c) + r n - μ n^{α}, x \in Ω, t > 0, \\ c_{t} + u \cdot \nabla c = Δ c - n c, x \in Ω, t > 0, \\ u_{t} + (u \cdot \nabla) u = Δ u + \nabla P + n \nabla Φ, x \in Ω, t > 0, \\ \nabla \cdot u = 0, x \in Ω, t > 0 . \end{cases}$

的初边值问题，其中 $Ω \subset ℝ^{3}$ 是一个具有光滑边界的有界区域；n和c满足齐次Neumann边界条件， $u$ 满足Dirichlet边界条件； $Φ \in W^{2, \infty} (Ω)$ ； $r > 0$ ， $μ > 0$ ， $α > 1$ 是给定的参数。此前的结果表明：当初值满足 $n_{0} \in C^{0} (\bar{Ω})$ 时，该模型在三维有界凸区域上存在整体弱解。本文进一步研究了当初值条件正则性更低时，该模型弱解的整体存在性。具体而言，初值 $n_{0}$ 满足 $n_{0} \in L^{1} (Ω)$ ，该模型在三维有界凸区域上存在整体弱解。

关键词

趋化–流体耦合模型，Logistic源，弱解

Global Existence of Solutions for a Chemotaxis-Fluid Coupling Model with Low Regular Initial Data

Juan Chen, Yongjie Han

School of Science, Xihua University, Chengdu Sichuan

Received: Oct. 15^th, 2023; accepted: Nov. 16^th, 2023; published: Nov. 24^th, 2023

ABSTRACT

This paper mainly studies the initial-boundary value problem of chemotaxis-fluid coupling equations

with signal consumption mechanism on a three-dimensional bounded domain with logistic source where $Ω \subset ℝ^{3}$ is bounded domain with smooth boundary; n and c satisfy the homogeneous Neumann boundary condition and $u$ satisfy the Dirichlet boundary condition; where $Φ \in W^{2, \infty} (Ω)$ ; where $r > 0$ , $μ > 0$ , and $α > 1$ are given parameters. Previous results show that the initial value satisfies $n_{0} \in C^{0} (\bar{Ω})$ , the model has a globally weak solution on the three-dimensional bounded convex region. This paper further examines when the initial value condition regularity is lower, the globally existence of weak solutions in this model. Specifically, the initial value $n_{0}$ satisfies $n_{0} \in L^{1} (Ω)$ , the model has a global weak solution on a three-dimensional bounded convex region.

Keywords:Chemotaxis-Fluid Coupling Model, Logistic Source, Weak Solution

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

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

1. 引言

趋化性是指细胞(菌)或生物体对化学物质浓度做出反应的定向运动，是细胞(菌)迁移的重要手段之一。20世纪70年代，Keller和Segel提出了一类抛物–抛物型偏微分方程组用来研究基网柄菌聚合现象 [1] ，其具体形式为

${\begin{cases} n_{t} = Δ n - \nabla \cdot (n S (n) \nabla c), \\ c_{t} = Δ c - c + n, \end{cases}$ (1.1)

其中未知函数 $n (x, t), c (x, t)$ 分别表示细胞密度和化学信号浓度，S表示趋化灵敏度函数。当S恒等于常值函数1时，所得模型被称为Keller-Segel极小模型 [2] [3] 。该模型自提出以来引起了许多数学家的兴趣，主要集中于解的爆破，整体存在性及渐近性等问题 [4] [5] [6] [7] 。在过去的40年里，很多文献致力于研究该模型的齐次Neumann初边值问题的有界性和爆破 [8] [9] [10] 。由于模型同时存在着扩散和迁移两种机制，其动力学行为比较复杂，某些形式会在有限时刻出现爆破，即细胞密度在某一刻趋于无穷大，这与实验观测的生物事实不符。因此，从生物学的角度出发，考虑模型在一定情况下的整体可解性具有现实意义。

当枯草杆菌悬浮在固着的水滴中时，该菌群会呈现复杂的时空行为。具体的说，枯草杆菌细胞会聚集成羽状，同时周围还会自发的产生大规模流体运动并形成对流斑图。为了用数学的方式描述这种现象。Tuval等人提出了如下趋化–流体耦合模型：

${\begin{cases} n_{t} + u \cdot \nabla n = Δ n - \nabla \cdot (n S (x, n, c) \nabla c), x \in Ω, t > 0, \\ c_{t} + u \cdot \nabla c = Δ c - n c, x \in Ω, t > 0, \\ u_{t} + κ (u \cdot \nabla) u = Δ u + \nabla P + n \nabla Φ, x \in Ω, t > 0, \\ \nabla \cdot u = 0, x \in Ω, t > 0, \end{cases}$ (1.2)

其中 $u$ 和P分别表示不可压缩流体的速度以及相应的压力， $Φ$ 是给定的重力势函数， $κ$ 表示非线性对流的强度。建立此模型基于两方面的假设：第一，菌体的运动要受随机扩散(“ $n_{t} = Δ n$ ”)，流体运动(“ $- u \cdot \nabla n$ ”)，氧气量递增梯度(“ $- \nabla (n \nabla c)$ ”)三者的共同作用；第二，单细胞菌体会消耗氧气(“ $- n c$ ”)并借助浮力的反作用对流体的运动(“ $- u \cdot \nabla c$ ”)产生影响。

在自然界中，描述生物种群模型的微分方程通常涉及超线性退化项。例如：经典的人口动态模型最初是由Pierre-Francois Verhulst在1838年提出，表明了一个生物种群在一定时间内的自限性增长，可通

过logistic方程 $\frac{d P}{d t} = r P - μ P^{2}$ 来描述，这里r为生长速率， $\frac{r}{μ}$ 是物种的可容能力，这也与同一物种的死

亡率有关。如果我们考虑种群在空间中的扩散和迁移。例如，在最简单的演化方程中，具有以下形式

$n_{t} = Δ n + r n - μ {| n |}^{α - 1}, α > 1, r \in R, μ > 0.$

由于这种类型的衰减机制通常提供额外的耗散，从而相应地增强了发散的特征。这种超线性阻尼可以用来扩展著名的热方程解理论，从而为具有远低于可积性的初值构造解。于是，在趋化–流体耦合模型中，细菌自身的增殖与死亡也应该被考虑。Mimura等人研究了基于细菌趋化扩散和增长的模型，他们在模型中增加了logistic源项 $r n - μ n^{2} (r \geq 0, μ > 0)$ ，通过考虑细胞的死亡来抵消由细胞数量的无限增长引起的细胞聚集。

我们知道logistic源对方程组解的爆破具有抑制作用，于是考虑具有logistic源的趋化–流体耦合模型：

${\begin{cases} n_{t} + u \cdot \nabla n = Δ n - \nabla \cdot (n S (x, n, c) \nabla c) + r n - μ n^{α}, x \in Ω, t > 0, \\ c_{t} + u \cdot \nabla c = Δ c - n c, x \in Ω, t > 0, \\ u_{t} + κ (u \cdot \nabla) u = Δ u + \nabla P + n \nabla Φ, x \in Ω, t > 0, \\ \nabla \cdot u = 0, x \in Ω, t > 0, \end{cases}$ (1.3)

其中 $Φ \in W^{2, \infty} (Ω)$ ，该模型有很多重要的结果 [11] [12] ，其中：Wang考虑了 $1 < α < 2$ 的情况，得到了模型

弱解的整体存在性和 $\frac{6}{5} \leq α < 2$ 时弱解的最终光滑性 [13] 。

根据上述趋化–流体耦合方程组的发展状况，在Wang的研究基础上，本论文研究在三维有界凸区域上的如下初边值问题：

(1.4)

在文献 [13] 中，由定理1.1可知，当 $r \geq 0, μ > 0, α > 1$ 时，对于初值 $n_{0}$ 满足 $n_{0} \in C^{0} (\bar{Ω})$ ，该模型在三维有界凸区域中存在整体弱解。本文针对初值 $n_{0} \in L^{1} (Ω)$ 的情况，对问题解的存在性进行了研究。

考虑其初值满足：

${\begin{cases} n_{0} \in L^{1} (Ω), n_{0} \geq 0, \\ c_{0} \in W^{1, \infty} (Ω), c_{0} \geq 0, \\ u_{0} \in D (A^{β}), \end{cases}$ (1.5)

其中 $β \in (\frac{3}{4}, 1)$ ， $A : = - P Δ$ 是Stokes算子，其定义域为 $D (A) : = W^{2, 2} (Ω; ℝ^{3}) \cap W_{0}^{1, 2} (Ω; ℝ^{3}) \cap L_{σ}^{2} (Ω)$ ，其中

$L_{σ}^{2} (Ω) : = {φ \in L^{2} (Ω; ℝ^{3}) | \nabla \cdot φ = 0}$ 。基于以上条件的假设，主要的结论如下：

定理1.1. 假设 $Ω \subset ℝ^{3}$ 是一个具有光滑边界的有界凸区域，当 $α > 1, r > 0, μ > 0$ ，且初值 $(n_{0}, c_{0}, u_{0})$ 满足条件(1.5)，则存在函数

${\begin{cases} n \in L_{l o c}^{1} ([0, \infty), W^{1, 1} (Ω)), \\ c \in L_{l o c}^{1} ([0, \infty), W^{1, 1} (Ω)), \\ u \in L_{l o c}^{1} ([0, \infty), W_{0, σ}^{1, 1} (Ω); ℝ^{3}), \end{cases}$ (1.6)

使得 $(n, c, u)$ 在命题2.1 意义下是初边值问题(1.4)的整体弱解。

2. 预备知识

首先给出弱解的定义。在下面的表述中，对于给定向量 $v \in ℝ^{3}$ 以及 $w \in ℝ^{3}$ ，令 $v \otimes w$ 表示矩阵 ${(v \otimes w)}_{i j} : = v_{i} w_{j}, i, j \in {1, 2, 3}$ 。

命题2.1. 若函数组 $(n, c, u)$ 满足下列条件，则称 $(n, c, u)$ 为初边值问题(1.4)的整体弱解：

(i)

其中 $n \geq 0, c \geq 0$ 在 $Ω \times (0, \infty)$ 上几乎处处成立，且 $n c \in L_{l o c}^{1} (\bar{Ω} \times [0, \infty))$ ， $n^{α} \in L_{l o c}^{1} (\bar{Ω} \times [0, \infty))$ ，

$u \otimes u \in L_{l o c}^{1} (\bar{Ω} \times [0, \infty); ℝ^{3 \times 3})$ ；且 $n \nabla c$ ， $n u$ ， $c u \in L_{l o c}^{1} (\bar{Ω} \times [0, \infty); ℝ^{3})$ ；

(ii) 对于任意 $ξ \in C_{0}^{\infty} (Ω \times [0, \infty))$ ，有

$\begin{array}{l} - \int_{0}^{\infty} \int_{Ω} n ξ_{t} - \int_{Ω} n_{0} ξ (\cdot, 0) - \int_{0}^{\infty} n u \cdot \nabla ξ \\ = - \int_{0}^{\infty} \int_{Ω} \nabla n \cdot \nabla ξ + \int_{0}^{\infty} \int_{Ω} n \nabla c \cdot \nabla ξ + r \int_{0}^{\infty} \int_{Ω} n ξ - μ \int_{0}^{\infty} \int_{Ω} n^{α} ξ, \end{array}$ (2.2)

和

$\begin{array}{l} - \int_{0}^{\infty} \int_{Ω} c ξ_{t} - \int_{Ω} c_{0} ξ (\cdot, 0) - \int_{0}^{\infty} \int_{Ω} c u \cdot \nabla ξ \\ = - \int_{0}^{\infty} \int_{Ω} \nabla c \cdot \nabla ξ - \int_{0}^{\infty} \int_{Ω} n c ξ, \end{array}$ (2.3)

对于任意的 $ψ \in C_{0, σ}^{\infty} (Ω \times [0, \infty))$ ，有

$\begin{array}{l} - \int_{0}^{\infty} \int_{Ω} u \cdot ψ_{t} - \int_{Ω} u_{0} \cdot ψ (\cdot, 0) - \int_{0}^{\infty} \int_{Ω} u \otimes u \cdot \nabla ψ \\ = - \int_{0}^{\infty} \int_{Ω} \nabla u \cdot \nabla ψ + \int_{0}^{\infty} \int_{Ω} n \nabla ψ \cdot \nabla Φ . \end{array}$ (2.4)

接下来，为了证明初边值问题(1.4)解的整体存在，我们首先引入一个函数序列 ${n_{0 ε}}_{ε \in (0, 1)}$ 满足

${n_{0 ε}}_{ε \in (0, 1)} \subset C^{0} (\bar{Ω}), n_{0 ε} \geq 0.$ (2.5)

在 $L^{1} (Ω)$ 中，当 $ε \to 0$ 时，有

$n_{0 ε} \to n_{0},$ (2.6)

以及对于任意的 $ε \in (0, 1)$ ，有

$\int_{Ω} n_{0 ε} \leq 2 \int_{Ω} n_{0} .$ (2.7)

其次，考虑下列正则化问题：

${\begin{cases} n_{ε t} + u_{ε} \cdot \nabla n_{ε} = Δ n_{ε} - \nabla \cdot (n_{ε} {F^{'}}_{ε} (n_{ε}) \nabla c_{ε}) + r n_{ε} - μ n_{ε}^{α}, x \in Ω, t > 0, \\ c_{ε t} + u_{ε} \cdot \nabla c_{ε} = Δ c_{ε} - c_{ε} \cdot F_{ε} (n_{ε}), x \in Ω, t > 0, \\ u_{ε t} + (Y_{ε} u_{ε} \cdot \nabla) u_{ε} = Δ u_{ε} + \nabla P_{ε} + n_{ε} \nabla Φ, x \in Ω, t > 0, \\ \nabla \cdot u_{ε} = 0, x \in Ω, t > 0, \\ \partial_{ν} n_{ε} = \partial_{ν} c_{ε} = 0, u_{ε} = 0, x \in \partial Ω, t > 0, \\ n_{ε} (x, 0) = n_{0 ε} (x), c_{ε} (x, 0) = c_{0} (x), u_{ε} (x, 0) = u_{0} (x), x \in Ω, \end{cases}$ (2.8)

其中 $ε \in (0, 1)$ ， $Y_{ε} = {(1 + ε A)}^{- 1}$ ， $F_{ε} (n_{ε}) : = \frac{1}{ε} ln (1 + ε n_{ε})$ 。

此正则化问题的经典解整体存在。

引理2.2. 假设 $Ω \in ℝ^{3}$ 是一个具有光滑边界的有界凸区域， $Φ \in W^{2, \infty} (Ω)$ ，且 $r \geq 0$ ， $μ > 0$ ， $α > 1$ ， $n_{0 ε}$ ， $c_{0}$ 和 $u_{0}$ 满足(1.5)，(2.5)~(2.7)，则对于任意的 $ε \in (0, 1)$ ，存在 $T_{\max, ε} \in (0, \infty]$ 使得初边值问题(2.8)在 $Ω \times (0, T_{\max, ε})$ 上有经典解 $(n_{ε}, c_{ε}, u_{ε}, P_{ε})$ ，且满足

${\begin{cases} n_{ε} \in C^{0} (\bar{Ω} \times [0, T_{\max, ε})) \cap C^{2, 1} (\bar{Ω} \times (0, T_{\max, ε})), \\ c_{ε} \in C^{0} (\bar{Ω} \times [0, T_{\max, ε})) \cap C^{2, 1} (\bar{Ω} \times (0, T_{\max, ε})) \cap \underset{p > 1}{\cap} C^{0} ([0, T_{\max, ε}); W^{1, p} (Ω)), \\ u_{ε} \in C^{0} (\bar{Ω} \times [0, T_{\max, ε}); ℝ^{3}) \cap C^{2, 1} (\bar{Ω} \times (0, T_{\max, ε}); ℝ^{3}), \\ P_{ε} \in C^{1, 0} (\bar{Ω} \times (0, T_{\max, ε})), \end{cases}$ (2.9)

对于任意 $q > 3$ ， $β \in (\frac{3}{4}, 1)$ 成立，且当 $T_{\max, ε} > 0$ 时，

$\lim_{t ↗ T_{\max, ε}} \sup ({‖ n_{ε} (\cdot, t) ‖}_{L^{\infty} (Ω)} + {‖ c_{ε} (\cdot, t) ‖}_{W^{1, q} (Ω)} + {‖ A^{β} u_{ε} (\cdot, t) ‖}_{L^{2} (Ω)}) = \infty .$ (2.10)

证明局部存在性的类似证明过程可以参考文献 [12] 的引理2.1和2.2。

3. 一些基本估计

引理3.1. 设 $(n_{ε}, c_{ε}, u_{ε})$ 是问题(2.8)的经典解，对于所有 $t \in (0, T_{\max, ε})$ 和任意 $ε \in (0, 1)$ ，满足

$\int_{Ω} n_{ε} (t) \leq \max {2 \int_{Ω} n_{0} (x), {(\frac{r}{μ})}^{\frac{1}{α - 1}} | Ω |} : = m^{*} .$ (3.1)

证明在方程组(2.8)的第一个方程两边同时积分，利用Hölder不等式和 $\nabla \cdot u = 0$ 以及边界条件可得

${(\int_{Ω} n_{ε})}_{t} = r \int_{Ω} n_{ε} - μ \int_{Ω} n_{ε}^{α} \leq r \int_{Ω} n_{ε} - μ {| Ω |}^{1 - α} {(\int_{Ω} n_{ε})}^{α},$ (3.2)

解这个常微分方程，即可证得(3.1)式。

引理3.2. 在引理3.1的假设下，存在 $C (τ) > 0$ 使得

$\int_{t}^{t + τ} \int_{Ω} n_{ε}^{α} \leq C (τ),$ (3.3)

对于所有 $t \in (0, T_{\max, ε} - τ)$ 和任意 $ε \in (0, 1)$ 成立，其中 $τ : = \min {1, \frac{1}{2} T_{\max, ε}}$ 。

证明对(3.1)式两端在时间上进行积分即可得证。

引理3.3. 设 $(n_{ε}, c_{ε}, u_{ε})$ 是模型(2.8)的经典解，对于所有 $t \in (0, T_{\max, ε})$ 和任意 $ε \in (0, 1)$ ，有

$n_{ε} \geq 0, c_{ε} \geq 0.$ (3.4)

证明对方程组(2.8)的第一个方程和第二个方程分别利用抛物比较原理，即可得 $n_{ε}$ 和 $c_{ε}$ 的非负性。

引理3.4. 设 $(n_{ε}, c_{ε}, u_{ε})$ 是模型(2.8)的经典解，且初值满足(1.5)，则对于所有 $t \in (0, T_{\max, ε})$ 和任意 $ε \in (0, 1)$ ，有

${‖ c_{ε} (\cdot, t) ‖}_{L^{\infty} (Ω)} \leq {‖ c_{0} (x) ‖}_{L^{\infty}} : = c_{\infty} .$ (3.5)

证明由抛物比较原理即可证得结论成立。

引理3.5. 对于所有 $t \in (0, T_{\max, ε})$ 和任意 $ε \in (0, 1)$ ，有

$\int_{0}^{T_{\max, ε}} \int_{Ω} {| \nabla c_{ε} |}^{2} \leq m_{*},$ (3.6)

其中 $m_{*} = \frac{1}{2} \int_{Ω} c_{0}^{2}$ 。

证明在方程组(2.8)的第二个方程两端同时乘以 $c_{ε}$ ，再在 $Ω$ 和时间上积分，对于所有 $t \in (0, T_{\max, ε})$ ，有

$\frac{1}{2} \int_{Ω} c_{ε}^{2} (t) + \int_{0}^{t} \int_{Ω} {| \nabla c_{ε} |}^{2} = \frac{1}{2} \int_{Ω} c_{0}^{2} - \int_{0}^{t} \int_{Ω} c_{ε}^{2} \cdot \frac{1}{ε} \ln (1 + ε n_{ε}) = m_{*},$

最后取 $t ↗ T_{\max, ε}$ 可证得(3.6)式。

下面，我们希望得到能量泛函

$\int_{Ω} n_{ε} \ln n_{ε} + \frac{1}{2} \int_{Ω} \frac{{| \nabla c_{ε} |}^{2}}{c_{ε}} + K \int_{Ω} {| u_{ε} |}^{2}$ (3.7)

的估计。

引理3.6.存在 $C > 0$ ，使得对于所有 $t \in (0, T_{\max, ε})$ 和任意 $ε \in (0, 1)$ ，有

$\frac{d}{d t} \int_{Ω} n_{ε} \ln n_{ε} + \frac{μ}{2} \int_{Ω} n_{ε}^{α} \ln n_{ε} + \int_{Ω} \frac{| \nabla n_{ε} |}{n_{ε}} \leq \int_{Ω} \frac{\nabla n_{ε} \cdot \nabla c_{ε}}{1 + ε n_{ε}} + C .$ (3.8)

证明对于 $\int_{Ω} n_{ε} \ln n_{ε}$ ，由方程组(2.8)的第一个方程，利用分部积分有

$\begin{matrix} \frac{d}{d t} \int_{Ω} n_{ε} \ln n_{ε} = \int_{Ω} n_{ε t} + \int_{Ω} n_{ε t} \ln n_{ε} \\ = r \int_{Ω} n_{ε} - μ \int_{Ω} n_{ε}^{α} - \int_{Ω} \frac{{| \nabla n_{ε} |}^{2}}{n_{ε}} + \int_{Ω} \frac{\nabla c_{ε} \cdot \nabla n_{ε}}{1 + ε n_{ε}} \\ + r \int_{Ω} n_{ε} \ln n_{ε} - μ \int_{Ω} n_{ε}^{α} \ln n_{ε}, \end{matrix}$

当 $α > 1$ 时，函数 $s \mapsto r s - μ s^{α}, s \in [0, \infty)$ 和 $s \mapsto (r s - \frac{μ}{2} s^{α}) \ln s$ 是有界的，所以存在 $C_{1} > 0$ ，使得对于

所有 $t \in (0, T_{\max, ε})$ ，有

$\frac{d}{d t} \int_{Ω} n_{ε} \ln n_{ε} + \int_{Ω} \frac{{| \nabla n_{ε} |}^{2}}{n_{ε}} + \frac{μ}{2} \int_{Ω} n_{ε}^{α} \ln n_{ε} \leq \int_{Ω} \frac{\nabla c_{ε} \cdot \nabla n_{ε}}{1 + ε n_{ε}} + 2 C_{1},$

则(3.8)式得证。

引理3.7. 存在 $K, k, c > 0$ ，使得对于所有 $t \in (0, T_{\max, ε})$ 和任意 $ε \in (0, 1)$ ，有

$\frac{d}{d t} \int_{Ω} \frac{{| \nabla c_{ε} |}^{2}}{c_{ε}} + k \int_{Ω} c_{ε} {| D^{2} \ln c_{ε} |}^{2} + k \int_{Ω} \frac{{| \nabla c_{ε} |}^{4}}{c_{ε}^{3}} \leq K \int_{Ω} {| \nabla u_{ε} |}^{2} - 2 \int_{Ω} \frac{\nabla n_{ε} \cdot \nabla c_{ε}}{1 + ε n_{ε}} + C .$ (3.9)

证明我们从计算 $\frac{d}{d t} \int_{Ω} \frac{{| \nabla c_{ε} |}^{2}}{c_{ε}}$ 开始，对于所有 $t \in (0, T_{\max, ε})$ 和任意 $ε \in (0, 1)$ ，有

$\begin{matrix} \frac{d}{d t} \int_{Ω} \frac{{| \nabla c_{ε} |}^{2}}{c_{ε}} = \int_{Ω} \frac{c_{ε} 2 \nabla c_{ε} \cdot \nabla c_{ε t} - c_{ε t} {| \nabla c_{ε} |}^{2}}{c_{ε}^{2}} \\ = - 2 \int_{Ω} \frac{{| Δ c_{ε} |}^{2}}{c_{ε}} + 2 \int_{Ω} Δ c_{ε} \cdot \frac{1}{ε} \ln (1 + ε n_{ε}) + 2 \int_{Ω} \frac{Δ c_{ε}}{c_{ε}} u_{ε} \cdot \nabla c_{ε} \\ + \int_{Ω} \frac{{| \nabla c_{ε} |}^{2} Δ c_{ε}}{c_{ε}^{2}} - \int_{Ω} \frac{{| \nabla c_{ε} |}^{2} \cdot \frac{1}{ε} \ln (1 + ε n_{ε})}{c_{ε}} - \int_{Ω} \frac{{| \nabla c_{ε} |}^{2}}{c_{ε}^{2}} u_{ε} \cdot \nabla c_{ε}, \end{matrix}$ (3.10)

由参考文献 [12] 的引理2.7，利用Young不等式和分部积分公式可得，存在 $k_{1} > 0$ ， $k_{2} > 0$ ， $k_{3} > 0$ 使得对于所有 $t \in (0, T_{\max, ε})$ 和任意 $ε \in (0, 1)$ ，有

$\frac{d}{d t} \int_{Ω} \frac{{| \nabla c_{ε} |}^{2}}{c_{ε}} + k_{1} \int_{Ω} c_{ε} {| D^{2} \ln c_{ε} |}^{2} + \frac{k_{1}}{2} \int_{Ω} \frac{{| \nabla c_{ε} |}^{4}}{c_{ε}^{3}} \leq {k_{3} \int_{Ω} \nabla u_{ε} |}^{2} - 2 \int_{Ω} \frac{\nabla n_{ε} \cdot \nabla c_{ε}}{1 + ε n_{ε}} + k_{2} \int_{Ω} c_{ε},$

则(3.9)得证。

引理3.8. 对于所有 $t \in (0, T_{\max, ε})$ 和任意 $ε \in (0, 1)$ ，有

$\frac{1}{2} \frac{d}{d t} \int_{Ω} {| u_{ε} |}^{2} + \int_{Ω} {| \nabla u_{ε} |}^{2} = \int_{Ω} n_{ε} u_{ε} \cdot \nabla Φ .$ (3.11)

证明在方程组(2.8)的第三个方程两端同时乘以 $u_{ε}$ ，再在 $Ω$ 上分部积分，利用分部积分公式和 $\int_{Ω} u_{ε} \cdot \nabla P = 0$ 以及 $\int_{Ω} (Y_{ε} u_{ε} \cdot \nabla) u_{ε} \cdot u_{ε} = 0$ 可证得结论成立。

引理3.9. 存在 $k, K, C > 0$ 使得对于所有 $t \in (0, T_{\max, ε})$ 和任意 $ε \in (0, 1)$ ，有

$\begin{array}{l} \frac{d}{d t} (\int_{Ω} n_{ε} \ln n_{ε} + \frac{1}{2} \int_{Ω} \frac{{| \nabla c_{ε} |}^{2}}{c_{ε}} + K \int_{Ω} {| u_{ε} |}^{2}) + \frac{1}{2} \int_{Ω} \frac{{| \nabla n_{ε} |}^{2}}{n_{ε}} \\ + \frac{k}{2} \int_{Ω} c_{ε} {| D^{2} \ln c_{ε} |}^{2} + \frac{μ}{2} \int_{Ω} n_{ε}^{α} \ln n_{ε} + \frac{k}{2} \int_{Ω} \frac{{| \nabla c_{ε} |}^{4}}{c_{ε}^{3}} + \frac{K}{2} \int_{Ω} {| \nabla u_{ε} |}^{2} \leq C . \end{array}$ (3.12)

证明由引理3.6、3.7和引理3.8可得，存在 $C_{1} > 0$ 使得对于任意的 $ε \in (0, 1)$ ，有

$\begin{array}{l} \frac{d}{d t} (\int_{Ω} n_{ε} \ln n_{ε} + \frac{1}{2} \int_{Ω} \frac{{| \nabla c_{ε} |}^{2}}{c_{ε}} + K \int_{Ω} {| u_{ε} |}^{2}) + \int_{Ω} \frac{{| \nabla n_{ε} |}^{2}}{n_{ε}} + \frac{k}{2} \int_{Ω} c_{ε} {| D^{2} \ln c_{ε} |}^{2} \\ + \frac{μ}{2} \int_{Ω} n_{ε}^{α} \ln n_{ε} + \frac{k}{2} \int_{Ω} \frac{{| \nabla c_{ε} |}^{4}}{c_{ε}^{3}} + K \int_{Ω} {| \nabla u_{ε} |}^{2} \leq 2 K \int_{Ω} n_{ε} u_{ε} \cdot \nabla Φ + C_{1} \end{array}$ (3.13)

对于 $t \in (0, T_{\max, ε})$ 成立。利用 $\nabla Φ$ 的有界性，Hölder不等式以及 $W^{1, 2} (Ω)$ 嵌入到 $L^{6} (Ω)$ ，有

$K \int_{Ω} n_{ε} u_{ε} \cdot \nabla Φ \leq C {(\int_{Ω} n_{ε}^{\frac{6}{5}})}^{\frac{5}{6}} \cdot {(\int_{Ω} u_{ε}^{\frac{6}{5}})}^{\frac{5}{6}} \leq C {‖ n_{ε} ‖}_{L^{5} (Ω)} {‖ u_{ε} ‖}_{L^{6} (Ω)} \leq C_{2} {‖ \nabla n_{ε} ‖}_{L^{\frac{6}{5}} (Ω)} {‖ \nabla u_{ε} ‖}_{L^{2} ( Ω )}$

对于所有的 $t \in (0, T_{\max, ε})$ 成立。对于上式的 ${‖ n_{ε} ‖}_{L^{\frac{6}{5}} (Ω)}$ ，利用Gagliardo-Nirenberg不等式、Young不等式

和(3.1)式，存在 $C_{5} > 0$ ，使得

$K \int_{Ω} n_{ε} u_{ε} \cdot \nabla Φ \leq \frac{K}{2} {‖ \nabla u_{ε} ‖}_{L^{2} (Ω)}^{2} + \frac{1}{2} {‖ \nabla n_{ε}^{\frac{1}{2}} ‖}_{L^{2} (Ω)}^{2} + C_{5} .$

对于所有的 $t \in (0, T_{\max, ε})$ 成立。将这个不等式代入(3.13)，结论得证。

引理3.10. 存在 $C > 0$ 使得对于所有 $t \in (0, T_{\max, ε})$ 和任意的 $ε \in (0, 1)$ ，有

$\int_{Ω} n_{ε} \ln n_{ε} + \int_{Ω} \frac{| \nabla c_{ε} |}{c_{ε}} + K \int_{Ω} {| u_{ε} |}^{2} \leq C .$ (3.14)

对于所有 $t \in [0, T_{\max, ε} - τ)$ ，有

$\int_{t}^{t + τ} \int_{Ω} \frac{{| \nabla n_{ε} |}^{2}}{n_{ε}} + \int_{t}^{t + τ} \int_{Ω} n_{ε}^{α} \ln n_{ε} + \int_{t}^{t + τ} \int_{Ω} \frac{{| \nabla c_{ε} |}^{4}}{c_{ε}} + \int_{t}^{t + τ} \int_{Ω} {| \nabla u_{ε} |}^{2} \leq C,$ (3.15)

$\int_{t}^{t + τ} \int_{Ω} n_{ε}^{\frac{5}{3}} + \int_{t}^{t + τ} \int_{Ω} {| \nabla n_{ε} |}^{\frac{5}{4}} + \int_{t}^{t + τ} \int_{Ω} {| u_{ε} |}^{\frac{10}{3}} + \int_{t}^{t + τ} \int_{Ω} {| \nabla c_{ε} |}^{4} + \int_{Ω} {| \nabla c_{ε} |}^{2} \leq C,$ (3.16)

其中 $τ : = \min {1, \frac{1}{2} T_{\max, ε}}$ 。

证明令

$H_{ε} (t) = \int_{Ω} n_{ε} \ln n_{ε} + \frac{1}{2} \int_{Ω} \frac{{| \nabla c_{ε} |}^{2}}{c_{ε}} + K \int_{Ω} {| u_{ε} |}^{2} .$

首先利用初等不等式 $s \ln s \leq \frac{3}{2} s^{\frac{5}{3}}$ 、Gagliardo-Nirenberg不等式、Young不等式和poincare不等式，结

合(3.1)式可得到，存在 $C_{1}, C_{2} > 0$ ，对于所有 $t \in (0, T_{\max, ε})$ ，有

${H^{'}}_{ε} (t) + C_{1} H_{ε} (t) \leq C_{2},$

于是(3.14)式得证。由(3.12)式可直接证得(3.15)式。对于(3.16)式，同样利用Gagliardo-Nirenberg不等

式，Young不等式，可得 $\int_{t}^{t + τ} \int_{Ω} n_{ε}^{\frac{5}{3}} \leq C_{3}$ ， $\int_{t}^{t + τ} \int_{Ω} {| \nabla n_{ε} |}^{\frac{5}{4}} \leq C_{4}$ ， $\int_{t}^{t + τ} \int_{Ω} {| u_{ε} |}^{\frac{10}{3}} \leq C_{5}$ 对于所有的 $t \in (0, T_{\max, ε} - τ)$ 成

立，其中 $C_{3}, C_{4}, C_{5} > 0$ 。进一步结合(3.15)式的有界性和(3.14)式可推得(3.16)式的最后一项 $\int_{Ω} {| \nabla c_{ε} |}^{2} < C_{6}$ ，其中 $C_{6} > 0$ ，则(3.16)式得证。

引理3.11. 对于任意 $ε \in (0, 1)$ ，问题(2.8)的解是整体存在的，即 $T_{\max, ε} = \infty$ 。

证明假设存在 $ε \in (0, 1)$ ，使得 $T_{\max, ε}$ 是有界的。由引理3.10，存在 $C_{1} > 0, C_{2} > 0$ ，使得 $\int_{t}^{t + τ} \int_{Ω} {| \nabla c_{ε} |}^{4} \leq C_{1}$ ，

$\int_{Ω} {| u_{ε} |}^{2} \leq C_{2}$ 对于所有的 $t \in (0, T_{\max, ε})$ 成立。在方程组(2.8)的第一个方程两端乘以 $n_{ε}^{3}$ ，再在 $Ω$ 上积分，利用分部积分公式可得：

$\frac{1}{4} \frac{d}{d t} \int_{Ω} n_{ε}^{4} + 3 \int_{Ω} n_{ε}^{2} {| \nabla n_{ε} |}^{2} = 3 \int_{Ω} \frac{n_{ε} \cdot n_{ε}^{3}}{1 + ε n_{ε}} \nabla c_{ε} \cdot \nabla n_{ε} + r \int_{Ω} n_{ε}^{4} - μ \int_{Ω} n_{ε}^{3 + α} .$

利用Young不等式和 $\frac{n_{ε}}{1 + ε n_{ε}} \leq n_{ε}$ 可知，对于所有 $t \in (0, T_{\max, ε})$ ，有

$\int_{Ω} n_{ε}^{4} (\cdot, t) \leq C_{3},$ (3.17)

其中 $C_{3} > 0$ 。对于 ${‖ u_{ε} ‖}_{L^{\infty} (Ω)}$ ，因为 $D (1 + ε A)$ 嵌入到 $L^{\infty} (Ω)$ ，对于所有 $t \in (0, T_{\max, ε})$ ，存在 $C_{4}, C_{5} > 0$ ，使得

${‖ Y_{ε} u_{ε} (\cdot, t) ‖}_{L^{\infty} (Ω)} = {‖ {(1 + ε A)}^{- 1} u_{ε} (\cdot, t) ‖}_{L^{\infty} (Ω)} \leq C_{4} {‖ u_{ε} (\cdot, t) ‖}_{L^{2} (Ω)} \leq C_{5} .$ (3.18)

令Stokes方程组 $u_{ε t} + A u_{ε} = P [- (Y_{ε} u_{ε} \cdot \nabla) u_{ε} + n_{ε} \nabla Φ] : = H_{ε} (x, t)$ ，在方程两端同时乘以 $A u_{ε}$ ，再在 $Ω$

上积分，由 ${‖ P φ ‖}_{L^{2} (Ω)} \leq {‖ φ ‖}_{L^{2} (Ω)}$ 对于所有 $φ \in L^{2} (Ω)$ 成立， $\int_{Ω} {| A^{\frac{1}{2}} φ |}^{2} = \int_{Ω} {| \nabla φ |}^{2}$ 对于所有的 $φ \in D (A)$ 成立可得

$\int_{Ω} {| \nabla u_{ε} |}^{2} \leq C_{6},$ (3.19)

于是对于所有 $t \in (0, T_{\max, ε})$ ，有

${‖ H_{ε} ‖}_{L^{2} (Ω)} \leq C_{7} .$ (3.20)

接下来，利用Stokes半群的估计可知 [14] ，存在 $λ > 0$ ， $C_{8} > 0$ ，对于所有 $t \in (0, T_{\max, ε})$ ，有

${‖ A^{β} u_{ε} ‖}_{L^{2} (Ω)} = {‖ A^{β} e^{- t A} u_{0} + \int_{Ω} A^{β} e^{- (t - s) A} H_{ε} (\cdot, s) d s ‖}_{L^{2} (Ω)} \leq C_{8} .$ (3.21)

由于 $D (A^{β})$ 嵌入到 $L^{\infty} (Ω)$ 可知，存在 $C_{9} > 0$ ，使得对于所有 $t \in (0, T_{\max, ε})$ ，有

${‖ u_{ε} ‖}_{L^{\infty} (Ω)} \leq C_{9} .$

同样，利用Neumann热半群的估计可知 [15] ，存在常数 $C_{10} > 0$ ， $C_{11} > 0$ ，使得对于所有 $t \in (0, T_{\max, ε})$ ，有

${‖ \nabla c_{ε} ‖}_{L^{4} (Ω)} \leq C_{10}$

和

${‖ n_{ε} ‖}_{L^{\infty} (Ω)} \leq C_{11}$

成立。证明细节可参考文献 [16] 中的引理4.1。结合(3.21)式，这与引理2.1的局部存在性准则矛盾，于是， $T_{\max, ε} = \infty$ 得证。

引理3.12. 对于所有 $T > 0$ 和任意 $ε \in (0, 1)$ ，存在 $C > 0$ ，有

${‖ n_{ε t} ‖}_{L^{\frac{10}{9}} ((0, T), {(W^{1, 10} (Ω))}^{*})} \leq C .$

证明取任意 $φ \in W^{1, 10} (Ω)$ ，对方程组(2.8)的第一个方程两端同时乘以 $φ$ ，对于任意 $t > 0$ ，有

$\begin{matrix} | \int_{Ω} n_{ε t} \cdot φ | \leq | \int_{Ω} \nabla n_{ε} \cdot \nabla φ | + | \int_{Ω} \frac{n_{ε}}{1 + ε n_{ε}} \nabla c_{ε} \cdot \nabla φ | + r | \int_{Ω} n_{ε} φ | + μ | \int_{Ω} n_{ε}^{α} φ | + | \int_{Ω} n_{ε} u_{ε} \cdot φ | \\ \leq {‖ \nabla n_{ε} ‖}_{L^{9} (Ω)}^{\frac{10}{9}} \cdot {‖ \nabla φ ‖}_{L^{10} (Ω)} + {‖ \frac{n_{ε}}{1 + ε n_{ε}} ‖}_{L^{\frac{10}{9}} (Ω)} \cdot {‖ \nabla φ ‖}_{L^{10} (Ω)} + r {‖ φ ‖}_{L^{10} (Ω)} \cdot {‖ n_{ε} ‖}_{L^{\frac{10}{9}} (Ω)} \\ + μ {‖ φ ‖}_{L^{\infty} (Ω)} \int_{Ω} n_{ε}^{α} + {‖ n_{ε} u_{ε} ‖}_{L^{\frac{10}{9}} (Ω)} \cdot {‖ φ ‖}_{L^{10} (Ω)} . \end{matrix}$

存在常数 $C_{1} > 0$ ，对于所有 $T > 0$ ，有

$\begin{array}{l} {\int_{0}^{T} ‖ n_{ε t} ‖}_{{(W^{1, 10} (Ω))}^{*}}^{\frac{10}{9}} \leq C_{1} \int_{0}^{T} \int_{Ω} {| \nabla n_{ε} |}^{\frac{10}{9}} + C_{1} \int_{0}^{T} \int_{Ω} {| \frac{n_{ε}}{1 + ε n_{ε}} |}^{\frac{10}{9}} \\ + C_{1} r \int_{0}^{T} \int_{Ω} {| n_{ε} |}^{\frac{10}{9}} + C_{1} μ \int_{0}^{T} \int_{Ω} n_{ε}^{α} + C_{1} \int_{0}^{T} \int_{Ω} {| n_{ε} u_{ε} |}^{\frac{10}{9}} . \end{array}$ (3.22)

利用Young不等式，结合引理3.2、3.10和(3.22)式可知，存在常数 $C > 0$ ，使得

$\int_{0}^{T} {‖ n_{ε t} ‖}_{{(W^{1, 10} (Ω))}^{*}}^{\frac{10}{9}} \leq C$

成立，则引理得证。

引理3.13. 对于所有 $T > 0$ 和任意 $ε \in (0, 1)$ ，存在 $C > 0$ ，有

${‖ c_{ε t} ‖}_{L^{2} ((0, T), {(W^{1, 10} (Ω))}^{*})} \leq C .$

证明取任意 $φ \in W_{0}^{1, 2} (Ω)$ ，对方程组(2.8)的第二个方程两端同时乘以 $φ$ ，利用Holder不等式和分部积分公式可知，对于任意的 $T > 0$ ，有

$\begin{array}{l} \int_{0}^{T} \int_{Ω} c_{ε t} φ \leq {(\int_{0}^{T} \int_{Ω} {| \nabla c_{ε} |}^{2})}^{\frac{1}{2}} {‖ \nabla φ ‖}_{L^{2} ((0, T), L^{2} (Ω))} + c_{\infty} {(\int_{0}^{T} \int_{Ω} {| n_{ε} |}^{\frac{5}{3}})}^{\frac{3}{5}} \\ \times {(\int_{0}^{T} \int_{Ω} {| φ |}^{\frac{5}{2}})}^{\frac{2}{5}} + c_{\infty} {(\int_{0}^{T} \int_{Ω} {| u_{ε} |}^{2})}^{\frac{1}{2}} {‖ \nabla φ ‖}_{L^{2} ((0, T), L^{2} ( Ω ))} \end{array}$

成立。因此，结合引理3.10，结论得证。

引理3.14. 对于所有 $T > 0$ 和任意 $ε \in (0, 1)$ ，存在 $C > 0$ ，有

${‖ u_{ε t} ‖}_{L^{2} ((0, T), {(W^{1, 3} (Ω))}^{*})} \leq C .$

证明取任意 $ψ \in L^{2} ((0, T), {(W^{1, 3} (Ω))}^{*})$ 且 ${‖ ψ ‖}_{L^{2} ((0, T), W^{1, 3} (Ω))} = 1$ ，则对于任意的 $T > 0$ ，有

$\begin{matrix} \int_{0}^{T} \int_{Ω} u_{ε t} \cdot ψ \leq | - \int_{0}^{T} \int_{Ω} \nabla u_{ε} \cdot \nabla ψ + \int_{0}^{T} \int_{Ω} Y_{ε} u_{ε} \otimes u_{ε} \nabla ψ + \int_{0}^{T} \int_{Ω} n_{ε} \nabla Φ \cdot ψ | \\ \leq {‖ \nabla u_{ε} ‖}_{L^{2} ((0, T), L^{2} (Ω))} {‖ \nabla ψ ‖}_{L^{2} ((0, T), L^{2} (Ω))} \\ + {‖ Y_{ε} u_{ε} ‖}_{L^{2} ((0, T), L^{6} (Ω))} {‖ u_{ε} ‖}_{L^{\infty} ((0, T), L^{2} (Ω))} {‖ \nabla ψ ‖}_{L^{2} ((0, T), L^{3} (Ω))} \\ + {‖ \nabla Φ ‖}_{L^{\infty} (Ω)} ({‖ n_{ε} ‖}_{L^{2} ((0, T), L^{6 / 5} (Ω))} + {‖ ψ ‖}_{L^{2} ((0, T), L^{6} (Ω))}) . \end{matrix}$

由 $W^{1, 2} (Ω)$ 嵌入到 $L^{6} (Ω)$ 和 $Y_{ε}$ 的性质可知，存在 $k_{1} > 0$ ，使得对于所有 $t > 0$ ，有

${‖ Y_{ε} u_{ε} ‖}_{L^{6} (Ω)} \leq k_{1} {‖ \nabla Y_{ε} u_{ε} ‖}_{L^{2} (Ω)} = k_{1} {‖ A^{1 / 2} Y_{ε} u_{ε} ‖}_{L^{2} (Ω)} = k_{1} {‖ \nabla u_{ε} ‖}_{L^{2} (Ω)} .$

进一步，利用Gagliardo-Nirenberg不等式和Young不等式推导可得

${‖ n_{ε} ‖}_{L^{6 / 5} (Ω)}^{2} = {‖ n_{ε}^{1 / 2} ‖}_{L^{12 / 5} (Ω)}^{4} \leq k_{2} {‖ \nabla n_{ε}^{1 / 2} ‖}_{L^{2} (Ω)}^{2} + k_{3} .$

又由引理3.10中包含的 ${‖ \nabla u_{ε} ‖}_{L^{2} ((0, T), L^{2} (Ω))}, {‖ u_{ε} ‖}_{L^{\infty} ((0, T), L^{2} (Ω))}, {‖ \nabla n_{ε}^{1 / 2} ‖}_{L^{2} ((0, T), L^{2} (Ω))}$ 的有界性，引理3.14得证。

命题3.15. 设 $r > 0, μ > 0$ ， $(n_{ε}, c_{ε}, u_{ε}, P_{ε})$ 是正则化问题(2.8)的经典解，存在函数 $(n, c, u)$ 满足(2.1)及一个子列 ${(ε_{j})}_{j \in ℕ}$ ，使得当 $j \to \infty$ ， $ε_{j} ↘ 0$ 并且当 $ε = ε_{j} ↘ 0$ 时有

$n_{ε} \to n$ 在 $L_{l o c}^{\frac{5}{3}} (\bar{Ω} \times [0, \infty))$ 并在 $Ω \times (0, \infty)$ 上几乎处处收敛， (3.23)

$\nabla n_{ε} ⇀ \nabla n$ 在 $L_{l o c}^{\frac{5}{4}} (\bar{Ω} \times [0, \infty))$ ， (3.24)

$n_{ε}^{α} \to n^{α}$ 在 $L_{l o c}^{1} (\bar{Ω} \times [0, \infty))$ ，(3.25)

$c_{ε} \to c$ 在 $Ω \times (0, \infty)$ 上几乎处处收敛， (3.26)

$c_{ε} \overset{⋆}{⇀} c$ 在 $L^{\infty} ((0, \infty); L^{p} (Ω))$ 对于所有的 $p \in (1, \infty]$ ，(3.27)

$\nabla c_{ε} ⇀ \nabla c$ 在 $L_{l o c}^{4} ([0, \infty); L^{4} (Ω))$ ， (3.28)

$u_{ε} \to u$ 在 $L^{\infty} ((0, \infty); L^{p} (Ω))$ $p \in [1, 6)$ 并且在 $Ω \times (0, \infty)$ 上几乎处处收敛，(3.29)

$u_{ε} ⇀ u$ 在 $L_{l o c}^{\frac{10}{3}} (\bar{Ω} \times [0, \infty))$ ， (3.30)

$\nabla u_{ε} ⇀ \nabla u$ 在 $L_{l o c}^{2} ((0, \infty); L^{2} (Ω))$ ， (3.31)

$Y_{ε} u_{ε} \otimes u_{ε} \to u \otimes u$ 在 $L_{l o c}^{1} (\bar{Ω} \times [0, \infty))$ 。 (3.32)

证明由(3.16)式可直接证得(3.24)式。由引理3.12和(3.16)式，利用Aubin-Lions引理可得 ${n_{ε}}_{ε \in (0, 1)}$ 在 $L_{l o c}^{\frac{5}{3}} (\bar{Ω} \times [0, \infty))$ 上是相对列紧的，因此存在子列 ${(ε_{j})}_{j \in ℕ}$ 使得(3.23)式成立。又由引理3.10可以得到，对于任意 $T > 0$ ， $\int_{0}^{T} \int_{Ω} n_{ε}^{α} \ln n_{ε}$ 的是有界的，则 ${n_{ε}^{α}}_{ε \in (0, 1)}$ 是一致可积的，于是通过(3.23)和Vitali收敛定理可知，存在子列使得(3.25)式成立。对于 $c_{ε}$ 的收敛性，由引理3.4和(3.16)式可直接证得(3.27)式和(3.28)式。由引理3.13和(3.16)式，利用Aubin-Lions引理可得 ${c_{ε}}_{ε \in (0, 1)}$ 在 $L^{\infty} ([0, T), W^{1, 2} (Ω))$ 上相对列紧的，因此存在子列 ${(ε_{j})}_{j \in ℕ}$ 使得式成立。最后，关于 $u_{ε}$ 的收敛性，由引理3.14和(3.16)式，利用Aubin-Lions引理可得 ${u_{ε}}_{ε \in (0, 1)}$ 在 $L^{2} ((0, \infty); L^{p} (Ω))$ 上是相对列紧的，因此存在子列 ${(ε_{j})}_{j \in ℕ}$ 使得(3.29)式成立。又由引理3.10

可直接证得(3.30)式和(3.31)式。对于(3.32)式，由Yosida逼近的性质结合(3.29)式可证得。

4. 主要结果的证明

定理1.1的证明现在的目的是证明 $(n, c, u)$ 是问题(1.6)的形如命题2.1所定义的整体弱解。显然 $n, c$ 有和 $n_{ε}, c_{ε}$ 一样的非负性。又由(3.31)可知 $\nabla \cdot u = 0$ 于 $Ω \times (0, \infty)$ 几乎处处成立。为了得到(2.1)~(2.4)的结论，我们有以下准备：

由(3.23)式、(3.28)式，当 $ε = ε_{j} ↘ 0$ 时有

$\frac{n_{ε} \nabla c_{ε}}{1 + ε n_{ε}} = \frac{n_{ε}}{1 + ε n_{ε}} \cdot \nabla c_{ε} ⇀ n \nabla c$ 于 $L_{l o c}^{1} (\bar{Ω} \times [0, \infty)) .$ (4.1)

实际上，对于 $T > 0$ ，有

$\begin{matrix} {‖ F_{ε} (n_{ε}) - n ‖}_{L^{\frac{5}{3}} (Ω \times (0, T))} \leq {‖ F_{ε} (n_{ε}) - F_{ε} (n) ‖}_{L^{\frac{5}{3}} (Ω \times (0, T))} + {‖ F_{ε} (n) - n ‖}_{L^{\frac{5}{3}} (Ω \times (0, T))} \\ \leq {‖ {F^{'}}_{ε} ‖}_{L^{\infty} ((0, \infty))} {‖ n_{ε} - n ‖}_{L^{\frac{5}{3}} (Ω \times (0, T))} + {‖ F_{ε} (n) - n ‖}_{L^{\frac{5}{3}} (Ω \times (0, T))} . \end{matrix}$ (4.2)

由(3.23)可知，当 $ε = ε_{j} ↘ 0$ 时， ${‖ n_{ε} - n ‖}_{L^{\frac{5}{3}} (Ω \times (0, T))} \to 0$ ，因此，对于 $t > 0$ 有

${‖ F_{ε} (n (\cdot, t)) - n (\cdot, t) ‖}_{L^{\frac{5}{3}} (Ω)}^{\frac{5}{3}} \leq 2^{\frac{5}{3}} {‖ n (\cdot, t) ‖}_{L^{\frac{5}{3}} (Ω)}^{\frac{5}{3}}$ 几乎处处成立。由控制收敛定理知，当 $ε = ε_{j} ↘ 0$ 时，也可得到

$\int_{0}^{T} {‖ F_{ε} (n (\cdot, t)) - n (\cdot, t) ‖}_{L^{\frac{5}{3}} (Ω)}^{\frac{5}{3}} d t \to 0.$ (4.3)

由(3.23)和(3.30)可知， $\frac{3}{5} + \frac{3}{10} = \frac{9}{10} < 1$ ，因此，当 $ε = ε_{j} ↘ 0$ 时，有

$n_{ε} u_{ε} \to n u$ 于 $L_{l o c}^{1} (\bar{Ω} \times [0, \infty)) .$ (4.4)

由(3.26)和(3.5)可知，利用Lebesgue收敛定理，当 $ε = ε_{j} ↘ 0$ 时，有 $c_{ε} \to c$ 于 $L_{l o c}^{\frac{5}{2}} (\bar{Ω} \times [0, \infty))$ 。结合

(4.3)式可得，当 $ε = ε_{j} ↘ 0$ 时，有

$F (n_{ε}) c_{ε} \to n c$ 于 $L_{l o c}^{1} (\bar{Ω} \times [0, \infty)) .$ (4.5)

接下来，由(3.26)和(3.5)式进一步可知，当 $ε = ε_{j} ↘ 0$ 时，有 $c_{ε} \to c$ 于 $L_{l o c}^{\frac{10}{7}} (\bar{Ω} \times [0, \infty))$ 。结合(3.30)式

可得，当 $ε = ε_{j} ↘ 0$ 时，有

$c_{ε} u_{ε} \to c u_{ε}$ 于 $L_{l o c}^{1} (\bar{Ω} \times [0, \infty)) .$ (4.6)

现在，由(4.4)，(4.1)，(4.5)，(4.6)，(3.32)，(3.25)式可以得到命题2.1中(i)的可积性，结合(3.23)，(3.31)，(2.5)和(2.6)式可以保证对方程组(2.8)取极限成立。

事实上，对于 $ε \in (0, 1)$ ，我们对方程组(2.8)中的第一个方程乘以任意的 $ξ \in C_{0}^{\infty} (Ω \times [0, \infty))$ ，然后分部积分，得到恒等式

$\begin{array}{l} - \int_{0}^{\infty} \int_{Ω} n_{ε} ξ_{t} - \int_{Ω} n_{0 ε} ξ (\cdot, 0) = - \int_{0}^{\infty} \int_{Ω} \nabla n_{ε} \cdot \nabla ξ + \int_{0}^{\infty} \int_{Ω} n_{ε} {F^{'}}_{ε} (n_{ε}) \nabla c_{ε} \cdot \nabla ξ \\ + \int_{0}^{\infty} \int_{Ω} n_{ε} u_{ε} \cdot \nabla ϕ + r \int_{0}^{\infty} \int_{Ω} n_{ε} ξ - μ \int_{0}^{\infty} \int_{Ω} n_{ε}^{α} ξ . \end{array}$

对于方程取 $ε = ε_{j} ↘ 0$ ，应用(3.23)，(2.6)，(3.24)，(4.1)，(4.4)式和(3.25)式在上述方程的积分式中，可得到(2.2)式。

同理，对于 $ε \in (0, 1)$ ，我们对方程组(2.8)中的第二个方程乘以任意的 $ξ \in C_{0}^{\infty} (Ω \times [0, \infty))$ ，然后分部积分，得到恒等式

$- \int_{0}^{\infty} \int_{Ω} c_{ε} ξ_{t} - \int_{Ω} c_{0} ξ (\cdot, 0) = - \int_{0}^{\infty} \int_{Ω} \nabla c_{ε} \cdot \nabla ξ - \int_{0}^{\infty} \int_{Ω} F (n_{ε}) c_{ε} ξ + \int_{0}^{\infty} \int_{Ω} c_{ε} u_{ε} \cdot \nabla ξ .$

对于方程取 $ε = ε_{j} ↘ 0$ ，应用(3.26)，(3.28)，(4.5)式和(4.6)式在上述方程的积分式中，可得到(2.3)式。

最后，对于 $ε \in (0, 1)$ ，我们对方程组(2.8)中的第三个方程乘以任意的 $ψ \in C_{0, σ}^{\infty} (Ω \times [0, \infty))$ ，然后分部积分，得到恒等式

$- \int_{0}^{\infty} \int_{Ω} u_{ε} \cdot ψ_{t} - \int_{Ω} u_{0} \cdot ψ (\cdot, 0) = - \int_{0}^{\infty} \int_{Ω} \nabla u_{ε} \cdot \nabla ψ + \int_{0}^{\infty} \int_{Ω} Y_{ε} u_{ε} \otimes u_{ε} \cdot \nabla ψ + \int_{0}^{\infty} \int_{Ω} n_{ε} \nabla Φ \cdot ψ .$

对于方程利用Holder不等式和 $\nabla Φ$ 的有界性，取 $ε = ε_{j} ↘ 0$ ，应用(3.29)，(3.31)，(3.32)和(3.23)式在上述方程的积分式中，可得到(2.4)式。于是，定理1.1得证。

5. 结论

本文主要是集中在对模型初边值问题整体可解性的分析上，至于整体解在时间趋于无穷大时的渐近行为和在某一时刻的光滑性没有做出分析，未来可以对解的这些性质进行研究讨论。对于带有logistic源的趋化–流体耦合模型的初边值问题，本文以及现有的文献在分析其解的整体存在性时总是考虑模型中的一种生物动力学机制的主要作用，未来可以通过对模型中两种或多种生物动力学机制的作用不分主次的同时进行考虑而得到结果。除了logistic源，细胞的非线性扩散也是阻碍其聚集的一种重要机制，在充分利用logistic源的基础上，可以更加细致的研究体积填充与logistic源之间的相互作用对模型初边值问题解的整体有界性的影响。

文章引用

陈娟,韩永杰. 一类具有低正则初值的趋化–流体耦合模型解的整体存在性
Global Existence of Solutions for a Chemotaxis-Fluid Coupling Model with Low Regular Initial Data[J]. 理论数学, 2023, 13(11): 3272-3285. https://doi.org/10.12677/PM.2023.1311341

参考文献

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

2. 王玉兰. 趋化-流体耦合模型研究进展[J]. 西华大学学报(自然科学版), 2016, 35(4): 30-38.

3. Bellomo, N., Bellouquid, A., Tao, Y.S., et al. (2005) Toward a Mathematical Theory of Keller-Segel Models of Pattern Formation in Biological Tissues. Mathematical Models in Application Sciences, 25, 1663-1763. https://doi.org/10.1142/S021820251550044X

4. Biler, P. (1998) Local and Global Solvability of Some Parabolic Systems Modelling Chemotaxis. Advances in Mathematical Sciences and Application, 8, 715-743.

5. Herrero, M.A. and Velazquez, J.J.L. (1997) A Blow-Up Mechanism for Chemotaxis Model. Annali Della Scuola Normale Superiore Di Pisa Classe Di Scienze, 24, 633-683.

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

7. Lankeit, J. (2005) Eventual Smoothness and Asymptotics in a Three-Dimensional Chemotaxis System with Logistic Source. Journal of Differential Equations, 258, 1158-1191. https://doi.org/10.1016/j.jde.2014.10.016

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

9. 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.

10. 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

11. Vorotnikov, D. (2014) Weak Solutions for a Bioconvection Model Related to Bacillus subtillis. Communications in Mathematical Sciences, 12, 545-563. https://doi.org/10.4310/CMS.2014.v12.n3.a8

12. Lankeit, J. (2016) Long-Term Behavior in a Chemotaxis-Fluid System with Logistic Source. Mathematical Models and Methods in Applied Sciences, 26, 2071-2109. https://doi.org/10.1142/S021820251640008X

13. Wang, Y.L. (2020) Global Solvability and Eventual Smoothness in a Chemotaxis-Fluid System with Weak Logistic-Type Degradation. Mathematical Models and Methods in Applied Sciences, 6, 1217-1252. https://doi.org/10.1142/S0218202520400102

14. Giga, Y. (1986) Solution for Semilinear Parabolic Equations in Lp and Regularity of Weak Solutions of the Navier-Stokes System. Journal of Differential Equations, 61, 186-212. https://doi.org/10.1016/0022-0396(86)90096-3

15. Henry, D. (1981) Geometric Theory of Semilinear Parabolic Equations. Lecture Notes in Mathematics, Vol. 840. https://doi.org/10.1007/BFb0089647

16. Zhang, Q. and Li, Y. (2015) Global Weak Solutions for the Three-Dimensional Chemotaxis Navier-Stokes System with Nonlinear Diffusion. Journal of Differential Equations, 259, 3730-3754. https://doi.org/10.1016/j.jde.2015.05.012

期刊菜单