Welcome to Journal of Beijing University of Chemical Technology, Today is
Email Alert  RSS
Management and Mathematics

Well-posedness of the Cauchy problem for compressible Navier-Stokes-Cahn-Hilliard equations in 3D

  • DongJie HOU ,
  • AoMing ZHAO ,
  • YaZhou CHEN , *
Expand
  • College of Mathematics and Physics, Beijing University of Chemical Technology, Beijing 100029, China

Received date: 2022-07-29

  Online published: 2023-10-18

Copyright

All rights reserved, without authorization

Abstract

The well-posedness of the Cauchy problem for the 3D compressible Navier-Stokes-Cahn-Hilliard equation, which describes immiscible two-phase flow with a diffuse interface has been studied. The energy method and the Schauder fixed point theorem were used to prove the existence and uniqueness of the global strong solution for small initial perturbations near the phase separation.

Cite this article

DongJie HOU , AoMing ZHAO , YaZhou CHEN . Well-posedness of the Cauchy problem for compressible Navier-Stokes-Cahn-Hilliard equations in 3D[J]. Journal of Beijing University of Chemical Technology, 2023 , 50(5) : 118 -125 . DOI: 10.13543/j.bhxbzr.2023.05.013

引言

气液两相流或多相流在化学工程、航空航天、微生物等工业领域有着重要的应用价值,吸引了大量学者对其进行研究。与单相流相比,两相流的研究由于存在扩散界面和相互作用而变得更加复杂。1958年,Cahn和Hilliard[1]引入界面自由能,提出著名的Cahn-Hilliard方程来描述两种不混溶流体之间的扩散界面。后来Lowengrub等[2]在Cahn-Hilliard方程中加入了流体运动与扩散界面的相互作用,提出了Navier-Stokes-Cahn-Hilliard(NSCH)方程,该方程可更好地刻画两种可压缩非混相流体流动的物理特性。关于NSCH方程的研究已经有很多。Abels等[3]使用文献[4]引入的框架,在不限制初值大小的情况下,证明了NSCH方程在有界域上初边值问题弱解的存在性。Chen等[5]证明了一维可压缩NSCH方程周期边值和混合边值问题解的适定性,随后Chen等[6]又研究了三维周期边值问题强解的全局存在性及大时间行为。王暐翼等[7]证明了带有van der Waals状态一维NSCH方程周期边值解的适定性。另一种常用的非混相两相流模型是Navier-Stokes-Allen-Cahn系统[8]。以上两个模型的本质区别在于Navier-Stokes-Cahn-Hilliard模型中关于相场的方程是一个四阶方程,它相对于ρφ是守恒的,而Navier-Stokes-Allen-Cahn中的相场方程是一个非守恒的二阶方程。
在前人工作的基础上,本文主要研究三维可压缩NSCH方程组的Cauchy问题的适定性。对于此类问题,我们克服了强非线性项带来的困难,在初值小扰动的条件下通过能量方法证明了全局强解的存在唯一性。

1 模型的构造及主要定理

考虑如下描述三维可压缩非混相两相流体流动的NSCH偏微分方程组。
$\left\{\begin{array}{l}\rho_t+\operatorname{div}(\rho \boldsymbol{u})=0 \\(\rho \boldsymbol{u})_t+\operatorname{div}(\rho \boldsymbol{u} \otimes \boldsymbol{u})=\operatorname{div} \boldsymbol{T} \\(\rho \phi)_t+\operatorname{div}(\rho \phi \boldsymbol{u})=\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \mu \\\mu=\frac{\partial f}{\partial \phi}-\frac{1}{\rho} \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\end{array}\right.$
式中,$\rho=\rho_1+\rho_2$是混合流体总密度;u为流体的速度,且满足$\rho \boldsymbol{u}=\rho_1 \boldsymbol{u}_1+\rho_2 \boldsymbol{u}_2, \boldsymbol{u}_i、\rho_i(i=1, 2)$分别为第i种组分的速度和密度;$\phi=\phi_1-\phi_2$为组分间浓度差,$\phi_i=\rho_i / \rho$μ为化学势;T为Cauchy应力张量;f为界面自由能密度。
$\left\{\begin{array}{l}\boldsymbol{T}=\boldsymbol{S}-p \boldsymbol{I}-\left(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \otimes \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi-\frac{|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi|^2}{2} \boldsymbol{I}\right) \\\boldsymbol{S}=2 \nu(\phi) \boldsymbol{D}(\boldsymbol{u})+\lambda(\phi)(\operatorname{div} \boldsymbol{u}) \boldsymbol{I}\end{array}\right.$
式中, $\nu(\phi)>0, \lambda(\phi) \geqslant 0$为黏性系数,ν(·), λ(·)∈C3(R)且$\lambda+\frac{2}{3} \nu \geqslant 0$S为牛顿黏性张力,D(u)= $\frac{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}+\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{\mathrm{T}} \boldsymbol{u}}{2}$为应变张量;I为单位矩阵;压力p是关于密度ρ的函数,且满足p′(ρ)>0。
f的计算公式为
$f=\frac{\phi^4}{4}-\frac{\phi^2}{2}$
我们要研究的模型满足如下初始条件
$\left\{ {\begin{array}{*{20}{l}}{(\rho , \mathit{\boldsymbol{u}}, \phi )(\mathit{\boldsymbol{x}}, 0) = \left( {{\rho _0}, {\mathit{\boldsymbol{u}}_0}, {\phi _0}} \right)(\mathit{\boldsymbol{x}})}\\{\mathop {\lim }\limits_{\mathit{\boldsymbol{|x}}| \to \pm \infty } \left( {{\rho _0}, {\mathit{\boldsymbol{u}}_0}, \left| {{\phi _0}} \right|} \right) = (\mathit{\boldsymbol{\overline \rho }} , 0, 1)}\end{array}} \right.$
式中,x=(x1, x2, x3), $\bar{\rho}$为给定的一个正实数, $\mathop {\lim }\limits_{|x| \to \pm \infty } \left| {{\phi _0}} \right| = 1$表示两相流的初始时刻浓度差φ0在无穷远处为1或-1。以下是本文的主要结论。
定理1  如果初值(ρ0, u0, φ0)满足
$\begin{aligned}& \quad\left(\rho_0-\bar{\rho}, \boldsymbol{u}_0\right) \in H^3, \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi_0 \in H^2, \phi_0^2-1 \in L^2, \\& \inf _{\boldsymbol{x} \in \mathbb{R}^3} \rho_0(\boldsymbol{x})>0\end{aligned}$
则存在δ>0,当
$\begin{aligned}& \left\|\rho_0-\bar{\rho}\right\|_{H^3}+\left\|\boldsymbol{u}_0\right\|_{H^3}+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi_0\right\|_{H^2}+ \\& \left\|\phi_0^2-1\right\|_{L^2} \leqslant \delta\end{aligned}$
那么方程(1)~(4)存在唯一解(ρ, u, φ)满足
$\begin{aligned}&(\rho-\bar{\rho}, \boldsymbol{u}) \in C\left([0, \infty) ; H^3\right), \phi^2-1 \in C([0, \infty) ;\left.L^2\right) \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \in C\left([0, \infty) ; H^2\right), \rho \in L^2\left([0, \infty) ; H^3\right) \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \in L^2\left([0, \infty) ; H^4\right), \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u} \in L^2\left([0, \infty) ; H^3\right)\end{aligned}$
$\begin{aligned}& \|(\rho-\bar{\rho}, \boldsymbol{u})(t)\|_{H^3}^2+\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi(t)\|_{H^2}^2+ \\& \left\|\phi^2(t)-1\right\|_{L^2}^2+\int_0^t\left(\|\rho-\bar{\rho}\|_{H^3}^2+\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}\|_{H^3}^2+\right. \\& \left.\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|_{H^4}^2\right) \mathrm{d} \tau \leqslant C\left(\left\|\left(\rho_0-\bar{\rho}, \boldsymbol{u}_0\right)\right\|_{H^3}^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi_0\right\|_{H^2}^2+\right. \\& \left.\left\|\phi_0^2-1\right\|_{L^2}^2\right)\end{aligned}$
式中,用Hl, l>0表示Hl(R3), L2表示L2(R3)。

2 主要定理的证明

证明定理1全局解存在性的思路如下:在方程组局部解存在的基础上,运用能量估计方法得到解的一致估计,从而将其延拓到全局解。

2.1 局部解的存在性

$\sigma=\rho-\bar{\rho}, \forall M>0, 0<T<\infty$,构造解空间如下。
$\begin{aligned}& X_M([0, T])=\{(\sigma, \boldsymbol{u}, \phi) \mid(\sigma, \boldsymbol{u}) \in C([0, \\& \left.T] ; H^3\right), \phi^2-1 \in C\left([0, T] ; L^2\right), \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \in C([0, T] ; \\& \left.H^2\right), \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \sigma \in L^2\left([0, T] ; H^2\right), \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u} \in L^2\left([0, T] ; H^3\right), \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \in L^2\left([0, T] ; H^4\right), \inf _{x \in \mathbb{R}^3, t \in[0, T]} \rho(\boldsymbol{x}, t)>0, \\& \sup _{t \in[0, T]}\left(\|(\sigma, \boldsymbol{u})(t)\|_{H^3}+\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|_{H^2}+\| \phi^2-\right. \\& 1 \|) \leqslant M\}\end{aligned}$
对方程组(1)采用线性化方法结合Schauder不动点定理,可以得到如下局部强解的存在唯一性命题,这里略去证明,具体方法可参考文献[9]。
命题1  $\forall M>0$,如果初始值(ρ0, u0, φ0)满足$\left\|\left(\rho_0-\bar{\rho}, \boldsymbol{u}_0\right)\right\|_{H^3}+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi_0\right\|_{H^2}+\left\|\phi_0^2-1\right\| \leqslant M$, $\mathop {\inf }\limits_{x \in {\mathbb{R}^3}} {\rho _0}(\boldsymbol{x}) > 0$, 那么$\exists T^*>0$,方程(1)在X2M([0, T*])内存在唯一局部解(ρ, u, φ)。

2.2 全局解的存在性

设(σ, u, φ)∈XM([0, T])为方程组(10)的局部解,将方程组(1)变形如下。
$\left\{ {\begin{array}{*{20}{l}}{{\sigma _t} + \bar \rho {\mathop{\rm div}\nolimits} {\boldsymbol{u}} = {g_1}}\\{{{\boldsymbol{u}}_t} - \frac{2}{{\bar \rho }}{\mathop{\rm div}\nolimits} [\nu (\phi ){\boldsymbol{D}}({\boldsymbol{u}})] - \frac{1}{{\bar \rho }}\mathit{\boldsymbol{ \boldsymbol{\nabla}}} [(\nu (\phi ) + }\\{\quad \quad \lambda (\phi )){\mathop{\rm div}\nolimits} {\boldsymbol{u}}] + \frac{{{p^\prime }(\bar \rho )}}{{\bar \rho }}\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \sigma + \frac{1}{{\bar \rho }}\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi = {g_2}}\\{\rho {\phi _t} + \rho {\boldsymbol{u}} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi = \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \mu }\\{\rho \mu = \rho \left( {{\phi ^3} - \phi } \right) - \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi }\end{array}} \right.$
式中,
$\left\{\begin{aligned}g_1= & -\operatorname{div}(\sigma \boldsymbol{u}) \\g_2= & -(\boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}) \boldsymbol{u}+h_1(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \sigma-2 h_2(\sigma) \cdot \\& \operatorname{div}[ \nu(\phi) \boldsymbol{D}(\boldsymbol{u})]-h_2(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}[(\nu(\phi)+ \\& \lambda(\phi)) \operatorname{div} \boldsymbol{u}-\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi]\end{aligned}\right.$
其中,${h_1}(\sigma ) = \frac{{{p^\prime }(\bar \rho )}}{{\bar \rho }} - \frac{{{p^\prime }(\rho )}}{\rho }, {h_2}(\sigma ) = \frac{1}{{\bar \rho }} - \frac{1}{\rho }$
由解空间(9),结合Sobolev嵌入定理可得, 存在M0>0足够小,使得$\forall$0<MM0,从而有
$0<\frac{\bar{\rho}}{2} \leqslant \rho(x, t) \leqslant 2 \bar{\rho}$
命题2  假设(σ0, u0)∈H3, φ02-1∈L2, $\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi_0$H2,则存在一个只依赖于初值和T的常数C,使得
$\begin{aligned}& \|(\sigma, \boldsymbol{u})(t)\|_{H^3}^2+\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi(t)\|_{H^2}^2+ \\& \left\|\phi^2(t)-1\right\|_{L^2}^2+\int_0^t\left(\|\sigma\|_{H^3}^2+\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}\|_{H^3}^2+\right. \\& \left.\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|_{H^4}^2\right) \mathrm{d} \tau \leqslant C\left(\left\|\left(\sigma_0, \boldsymbol{u}_0\right)\right\|_{H^3}^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi_0\right\|_{H^2}^2+\right. \\& \left.\left\|\phi_0^2-1\right\|_{L^2}^2\right)\end{aligned}$
命题2的证明由以下4个引理得到。
引理1  若(σ, u, φ)∈XM([0, T])是方程(10)的局部解,则
$\begin{array}{r}\left\|\left(\sigma, \boldsymbol{u}, \phi^2-1, \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\right)(t)\right\|^2+\int_0^t \|(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mu, \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}, \\\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi, \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi)\left\|^2 \mathrm{~d} \tau \leqslant\right\|\left(\sigma, \boldsymbol{u}, \phi^2-1, \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\right)(0) \|^2\end{array}$
$\|\phi\|_{L^{\infty}} \leqslant C$
证明:设
$G(\rho)=\rho \int_{\bar{\rho}}^\rho \frac{p(\eta)-p(\bar{\rho})}{\eta^2} \mathrm{~d} \eta, \rho>0$
由式(15)和方程组(1)中第一式得到
$\begin{gathered}\rho G^{\prime}(\rho)=G(\rho)+(p(\rho)-p(\bar{\rho})), \rho G^{\prime \prime}(\rho)=p^{\prime}(\rho) \\G(\rho)_t+\operatorname{div}(G(\rho) \boldsymbol{u})+(p(\rho)-p(\bar{\rho})) \operatorname{div} \boldsymbol{u}=0\end{gathered}$
方程组(1)中第二式乘以u再积分,结合式(16)可得
$\begin{aligned}& \quad \frac{\mathrm{d}}{\mathrm{d} t} \int\left(\frac{1}{2} \rho \boldsymbol{u}^2+G(\rho)\right) \mathrm{d} x+\frac{1}{2} \int \nu(\phi) \cdot\\& \left|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}+\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{\mathrm{T}} \boldsymbol{u}\right|^2 \mathrm{~d} \boldsymbol{x}+\int \lambda(\phi)|\operatorname{div} \boldsymbol{u}|^2 \mathrm{~d} \boldsymbol{x}+\int \boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \cdot \\& \mathrm{d} \boldsymbol{x}=0\end{aligned}$
方程组(10)中第三式乘以μ, 并运用(10)中第四式,可得
$\begin{aligned}& \frac{1}{2} \frac{\mathrm{d}}{\mathrm{d} t} \int|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi|^2 \mathrm{~d} \boldsymbol{x}+\frac{1}{4} \frac{\mathrm{d}}{\mathrm{d} t} \int \rho\left(\phi^2-1\right)^2 \mathrm{~d} \boldsymbol{x}+ \\& \int|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mu|^2 \mathrm{~d} \boldsymbol{x}=-\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}(\boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi) \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathrm{d} \boldsymbol{x}\end{aligned}$
将式(17)、(18)相加,得到
$\begin{aligned}& \frac{\mathrm{d}}{\mathrm{d} t} \int\left(\frac{1}{2} \rho \boldsymbol{u}^2+G(\rho)+\frac{1}{2}|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi|^2+\frac{\rho}{4} \cdot\right. \\& \left.\left(\phi^2-1\right)^2\right) \mathrm{d} \boldsymbol{x}+\frac{\nu_0}{2}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}\|^2+\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mu\|^2 \leqslant 0\end{aligned}$
由式(9)、(11)、(15)可得
$c_{\bar{\rho}}(\rho-\bar{\rho})^2 \leqslant G(\rho) \leqslant C_{\bar{\rho}}(\rho-\bar{\rho})^2$
其中,$c_{\bar{\rho}}$$C_{\bar{\rho}}$是与$\bar{\rho}$有关的常数。
将式(19)在[0, T]上积分,结合式(20)得
$\begin{aligned}& \left\|\left(\sigma, \boldsymbol{u}, \phi^2-1, \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\right)\right\|^2+\int_0^t\|(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mu, \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u})\|^2 \mathrm{~d} \tau \leqslant \\& \left\|\left(\sigma_0, \boldsymbol{u}_0, \phi_0^2-1, \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi_0\right)\right\|^2\end{aligned}$
$\begin{aligned}& \left\|\phi^2-1\right\|_{L^{\infty}} \leqslant C\left\|\phi^2-1\right\|_{H^2}=C\left(\left\|\phi^2-1\right\|_{L^2}+\right. \\& \left.\|\phi \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|_{L^2}+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^2\left(\phi^2-1\right)\right\|_{L^2}\right) \leqslant M\end{aligned}$
故可得式(14)。
在方程组(10)中第四式两边乘以-Δφ, 关于x积分得
$\begin{array}{r}\|\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\|^2+\int \rho\left(3 \phi^2-1\right)|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi|^2 \mathrm{~d} \boldsymbol{x}= \\-\int \rho \mu \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \mathrm{d} \boldsymbol{x}-\int\left(\phi^2-1\right) \phi \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \sigma \mathrm{d} \boldsymbol{x}\end{array}$
由式(9)、(11)和Hölder不等式得到
$\begin{aligned}& \int \rho\left(3 \phi^2-1\right)|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi|^2 \mathrm{~d} \boldsymbol{x} \geqslant \frac{\bar{\rho} m_0}{2}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|^2- \\& \int \rho \mu \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \mathrm{d} \boldsymbol{x} \leqslant 2 \bar{\rho}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mu\|\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\| \leqslant \frac{\varepsilon}{4}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|^2+ \\& \frac{4 \bar{\rho}^2}{\varepsilon}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mu\|^2 \int\left(\phi^2-1\right) \phi \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \sigma \mathrm{d} \boldsymbol{x} \leqslant\left\|\left(\phi^2-1\right)\right\|_{L^6} \cdot \\& \|\phi \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \sigma\|_{L^3} \leqslant M\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|^2\end{aligned}$
M足够小,可得
$\|\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\|^2+\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|^2 \leqslant\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mu\|^2$
结合式(23)、(21)得到式(13), 引理1得证。
引理2  若(σ, u, φ)∈XM([0, T])是方程(10)的局部解,则
$\begin{aligned}& \quad \frac{\mathrm{d}}{\mathrm{d} t}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \phi\right\|^2+\frac{1}{4 \bar{\rho}^2}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|^2 \leqslant M \cdot\\& \left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^2\right), k=0, \\& 1, 2\end{aligned}$
证明:将方程组(10)中第四式代入(10)中第三式,可得
$\begin{aligned}& \phi_t+\boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi+\frac{1}{\rho} \mathit{\boldsymbol{ \boldsymbol{\Delta}}}\left(\frac{\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi}{\rho}\right)-\frac{3 \phi^2-1}{\rho} \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi=\frac{6 \phi}{\rho} . \\& (\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi)^2\end{aligned}$
对式(25)求$\mathit{\boldsymbol{ \boldsymbol{\nabla}}}$k,然后乘以$-\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \phi$再积分得
$\begin{aligned}& \quad \frac{1}{2} \frac{\mathrm{d}}{\mathrm{d} t} \int\left|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \phi\right|^2 \mathrm{~d} \boldsymbol{x}+\int\left|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}\left(\frac{\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi}{\rho}\right)\right|^2 \mathrm{~d} \boldsymbol{x}- \\& \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k(\boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi) \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \phi \mathrm{d} \boldsymbol{x}=\sum\limits_{1 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \cdot \\& \mathit{\boldsymbol{ \boldsymbol{\Delta}}}\left(\frac{\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \mathrm{d} \boldsymbol{x}+\sum\limits_{1 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\Delta}}}\left[\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\right] . \\& \frac{\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi}{\rho} \mathrm{d} \boldsymbol{x}-\sum\limits_{1 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{3 \phi^2-1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \mathrm{d} \boldsymbol{x}- \\& \int \frac{3 \phi^2-1}{\rho}\left|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\right|^2 \mathrm{~d} \boldsymbol{x}-\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k\left(\frac{6 \phi}{\rho}(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi)^2\right) \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \phi \mathrm{d} \boldsymbol{x}\end{aligned}$
由于
$\begin{array}{*{20}{l}}{\int {{{\left| {\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \left( {\frac{{{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^k}\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi }}{\rho }} \right)} \right|}^2}} {\rm{d}}\mathit{\boldsymbol{x}} = \int {\frac{1}{{{\rho ^2}}}} {{\left| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^{k + 1}}\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi } \right|}^2}{\rm{d}}\mathit{\boldsymbol{x}} + }\\{\int {\frac{{|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \rho {|^2}}}{{{\rho ^4}}}} {{\left| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^k}\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi } \right|}^2}{\rm{d}}\mathit{\boldsymbol{x}} - \int {\frac{{3{\phi ^2} - 1}}{\rho }} {{\left| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^k}\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi } \right|}^2}{\rm{d}}\mathit{\boldsymbol{x}} \le }\\{ - \int {\frac{{{\phi ^2} - 1}}{\rho }} {{\left| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^k}\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi } \right|}^2}{\rm{d}}\mathit{\boldsymbol{x}}}\end{array}$
故可得
$\begin{aligned}& \quad \frac{1}{2} \frac{\mathrm{d}}{\mathrm{d} t}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \phi\right\|^2+\frac{1}{4 \bar{\rho}^2}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\right\|^2 \leqslant I_1^i, i= \\& 1, 2, \cdots, 6\end{aligned}$
根据Hölder不等式和Gagliardo-Nirenberg不等式对I1i(i=1, 2,…,6)进行估计。
$\begin{aligned}& I_1^1=\int \frac{\phi^2-1}{\rho}\left|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\right|^2 \mathrm{~d} \boldsymbol{x} \leqslant\left\|\phi^2-1\right\|_{L^{\infty}} \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\right\|^2 \leqslant M\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2\end{aligned}$
$\begin{aligned}& I_1^2=\sum\limits_{0 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \phi \mathrm{d} \boldsymbol{x} \leqslant \sum\limits_{0 \leqslant l \leqslant k} \| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{u} \cdot \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi\|\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi \|\end{aligned}$
$l \leqslant\left[\frac{k+1}{2}\right]$时,有
$\begin{aligned}& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi\right\| \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{u}\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi\right\|_{L^6} \leqslant \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \boldsymbol{u}\right\|^{1-\frac{l}{k+1}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^{\frac{l}{k+1}}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|^{\frac{l}{k+1}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^{1-\frac{l}{k+1}} \leqslant \\& M\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^{\frac{l}{k+1}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^{1-\frac{l}{k+1}} \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|\right)\end{aligned}$
$\frac{{l - 1}}{3} = \left( {\frac{\alpha }{3} - \frac{1}{2}} \right)\left( {1 - \frac{l}{{k + 1}}} \right) + \left( {\frac{{k + 2}}{3} - \frac{1}{2}} \right){\rm{\cdot}}\frac{l}{{k + 1}}$,可确定$\alpha=\frac{1}{2}-\frac{l}{2(k+1-l)} \in\left[0, \frac{1}{2}\right]$
$\left[\frac{k+1}{2}\right]+1 \leqslant l \leqslant k$时,有
$\begin{aligned}& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi\right\| \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{u}\right\|_{L^6}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi\right\|_{L^3} \leqslant \\& \|\boldsymbol{u}\|^{1-\frac{l+1}{k+2}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^{\frac{l+1}{k+2}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \phi\right\|^{\frac{l+1}{k+2}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^{1-\frac{l+1}{k+2}} \leqslant \\& M\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^{\frac{l+1}{k+2}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^{1-\frac{l+1}{k+2}} \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|\right)\end{aligned}$
$\frac{{k - l}}{3} = \left( {\frac{\alpha }{3} - \frac{1}{2}} \right)\frac{{l + 1}}{{k + 2}} + \left( {\frac{{k + 2}}{3} - \frac{1}{2}} \right)(1 - \left. {\frac{{l + 1}}{{k + 2}}} \right)$
可确定$\alpha=\frac{k+2}{2(l+1)} \in\left[\frac{1}{2}, 1\right]$
因此,可得到
$I_1^2 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2\right)$
$\begin{aligned}\;\;\;I_1^3 & =\sum\limits_{1 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \mathit{\boldsymbol{ \boldsymbol{\Delta}}}\left(\frac{\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \mathrm{d} \boldsymbol{x} \leqslant \\\sum\limits_{1 \leqslant l} & \left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2}\left(\frac{\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi}{\rho}\right)\right\|\right)\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|_{L^2}\end{aligned}$
其中
$\begin{aligned}& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2}\left(\frac{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi}{\rho}\right)\right\| \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l(\sigma)\right\|_{L^6} \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2}\left(\frac{\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi}{\rho}\right)\right\|_{L^3} \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l(\sigma)\right\|_{L^6}\left(\|\left(\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2} \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\right) \cdot\right. \\& \left.\frac{1}{\rho}+\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2}\left(\frac{1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \|_{L^3}\right) \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l(\sigma)\right\|_{L^6} \cdot \\& \left(\left\|\left(\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+4} \phi\right)\right\|_{L^3}\left\|\frac{1}{\rho}\right\|_{L^{\infty}}+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2}\left(\frac{1}{\rho}\right)\right\|_{L^3} \cdot\right. \\& \left.\|\mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\|_{L^{\infty}}\right) \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|\right)\end{aligned}$
$I_1^3 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2\right)$
$\begin{aligned}& I_1^4=-\sum\limits_{1 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}\left[\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\right] . \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}\left(\frac{\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi}{\rho}\right) \mathrm{d} \boldsymbol{x} \leqslant \sum\limits_{1 \leqslant l \leqslant k}\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{1}{\rho}\right)\right\|_{L^3} \cdot\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+3} \phi\right\|_{L^6}+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{l+1}\left(\frac{1}{\rho}\right)\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2} \phi\right\|_{L^6}\right) \times \\& \left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}\left(\frac{1}{\rho}\right)\right\|_{L^{\infty}}\right) \leqslant \\& \sum\limits_{1 \leqslant l \leqslant k}\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+3} \phi\right\|_{L^6}+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{l+1} \sigma\right\|_{L^3} \cdot\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2} \phi\right\|_{L^6}\right)\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \sigma\right\|\right) \leqslant \\& M \sum\limits_{1 \leqslant l \leqslant k}\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|\right)\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|\right) \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|^2\right)\end{aligned}$
$\begin{aligned}& I_1^5=\sum\limits_{1 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{3 \phi^2-1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \mathrm{d} \boldsymbol{x} \leqslant \\& \sum\limits_{1 \leqslant l \leqslant k}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{3 \phi^2-1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2} \phi\right\|\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\| \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{3 \phi^2-1}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2} \phi\right\| \leqslant\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma\right\|_{L^3}+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \phi\right\|_{L^3}\right) \cdot\\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2} \phi\right\|_{L^6} \leqslant\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \sigma\right\|^{1-\frac{l}{k+2}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^{\frac{l}{k+2}}+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \phi\right\|^{1-\frac{l}{k+2}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^{\frac{l}{k+2}}\right) \quad \times \quad\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|^{\frac{l}{k+2}} \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^{1-\frac{l}{k+2}} \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|\right)\end{aligned}$
其中α满足
$\begin{aligned}& \quad \frac{l}{3}-\frac{1}{3}=\left(\frac{\alpha}{3}-\frac{1}{2}\right)\left(1-\frac{l}{k+2}\right)+\left(\frac{k+1}{3}-\right. \\& \left.\frac{1}{2}\right) \frac{l}{k+2}\end{aligned}$
$I_1^5 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2\right)$
$\begin{aligned}& I_1^6=\sum\limits_{0 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{6 \phi}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l}(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi) 2 \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \phi \mathrm{d} \boldsymbol{x} \leqslant \\& \sum\limits_{1 \leqslant l \leqslant k}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{6 \phi}{\rho}\right) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l}(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi)^2\right\|\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\| \| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l\left(\frac{6 \phi}{\rho}\right) . \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l}(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi)^2 \| \leqslant\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma\right\|_{L^3}+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \phi\right\|_{L^3}\right) . \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi\right\|_{L^6}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi\|_{L^{\infty}} \leqslant\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma\right\|_{L^3}+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \phi\right\|_{L^3}\right)\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi\right\|_{L^6} \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|\right)\end{aligned}$
$I_1^6 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2\right)$
I1i(i=1, 2, …, 6)的估计式代入式(27),可得式(24),引理2得证。
引理3  若(σ, u, φ)∈XM([0, T])是方程(10)的局部解,则
$\begin{aligned}& \frac{\mathrm{d}}{\mathrm{d} t}\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}\right\|+\frac{p^{\prime}(\bar{\rho})}{\bar{\rho}^2}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma\right\|^2\right)+\frac{\nu_0}{\bar{\rho}} \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^2 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2\right), k=0, \\& 1, 2, \cdots, 3\end{aligned}$
证明:对式(10)求$\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k$得到
$\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma_t+\bar{\rho} \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}=\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k g_1$
$\begin{aligned}& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}_t-\frac{2}{\bar{\rho}} \operatorname{div}\left(\nu(\phi) \boldsymbol{D}\left(\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}\right)\right)-\frac{1}{\bar{\rho}} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}\{[\boldsymbol{\nu}(\phi)+ \\& \left.\lambda(\phi)] \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}\right\}+\frac{p^{\prime}(\bar{\rho})}{\bar{\rho}} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \sigma+\frac{1}{\bar{\rho}} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi)= \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k g_2+\frac{2}{\bar{\rho}} \operatorname{div}\left\{\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k[\nu(\phi) \boldsymbol{D}(\boldsymbol{u})]-\nu(\phi) \boldsymbol{D}\left(\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}\right)\right\}+ \\& \frac{1}{\bar{\rho}} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}\left\{\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k[\nu(\phi)+\lambda(\phi)] \operatorname{div} \boldsymbol{u}-[\nu(\phi)+\lambda(\phi)] \cdot\right. \\& \left.\operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}\right\}\end{aligned}$
对式(29)的第二式乘以$\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}$,并利用(29)的一式可得
$\begin{aligned}& \frac{1}{2} \frac{\mathrm{d}}{\mathrm{d} t}\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}\right\|^2+\frac{p^{\prime}(\bar{\rho})}{\bar{\rho}^2}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma\right\|^2\right)+\frac{1}{\bar{\rho}} \cdot \\& \int \nu(\phi)\left|\boldsymbol{D}\left(\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}\right)\right|^2 \mathrm{~d} \boldsymbol{x}+\frac{1}{\bar{\rho}} \int[\nu(\phi)+\lambda(\phi)] \cdot \\& \left|\operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u}\right|^2 \mathrm{~d} \boldsymbol{x}=I_2^1+I_2^2+I_2^3+I_2^4+I_2^5+I_2^6+I_2^7\end{aligned}$
其中,
$\begin{aligned}& I_2^1=\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \operatorname{div}(\sigma \boldsymbol{u}) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma \mathrm{d} \boldsymbol{x} \leqslant \sum\limits_{0 \leqslant l \leqslant k+1} \| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma \cdot \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\|\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma \|\end{aligned}$
$l \leqslant\left[\frac{k}{2}\right]$时,有
$\begin{aligned}& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\| \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{\sigma}\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\|_{L^6} \leqslant \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \sigma\right\|^{1-\frac{l-1}{k}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma\right\|^{\frac{l-1}{k}}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}\|^{\frac{l-1}{k}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^{1-\frac{l-1}{k}} \leqslant \\& M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma\right\|+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|\right)\end{aligned}$
其中α$\frac{{l - 1}}{3} = \left( {\frac{\alpha }{3} - \frac{1}{2}} \right)\left( {1 - \frac{{l - 1}}{k}} \right) + \left( {\frac{k}{3} - \frac{1}{2}} \right) \cdot \frac{{l - 1}}{k}$确定。
$l \geqslant\left[\frac{k}{2}\right]+1$时,有
$\begin{aligned}& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\| \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma\right\|_{L^6}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi\right\|_{L^3} \leqslant \\& \|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \sigma\|^{1-\frac{l-1}{k}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma\right\|^{\frac{l+1}{k}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \boldsymbol{u}\right\|^{\frac{l+1}{k}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^{1-\frac{l+1}{k}} \leqslant \\& M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma\right\|+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|\right)\end{aligned}$
其中α$\frac{{k - l}}{3} = \left( {\frac{\alpha }{3} - \frac{1}{2}} \right)\left( {\frac{{l + 1}}{k}} \right) + \left( {\frac{{k + 1}}{3} - \frac{1}{2}} \right) \cdot \left( {1 - \frac{{l + 1}}{k}} \right)$确定。
$I_2^1 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \sigma\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^2\right)$
$\begin{array}{*{20}{l}} {I_2^2 = \int {{\nabla ^k}} \left[ { - (\boldsymbol{u},\nabla )\boldsymbol{u} + {h_1}(\sigma )\nabla \sigma } \right]{\nabla ^k}\boldsymbol{u}{\text{d}}\boldsymbol{x} = } \\ { - \int {{\nabla ^{k - 1}}} \left[ { - (\boldsymbol{u},\nabla )\boldsymbol{u} + {h_1}(\sigma )\nabla \sigma } \right]{\nabla ^{k + 1}}\boldsymbol{u}{\text{d}}\boldsymbol{x} \leqslant } \\ {\left( {\parallel \boldsymbol{u}{\parallel _{{L^3}}}{{\left\| {{\nabla ^k}\boldsymbol{u}} \right\|}_{{L^6}}} + {{\left\| {{\nabla ^{k - 1}}\boldsymbol{u}} \right\|}_{{L^6}}}\parallel \nabla \boldsymbol{u}{\parallel _{{L^3}}}} \right) \cdot } \\ {\left\| {{\nabla ^{k + 1}}{\mathbf{u}}} \right\| + \left( {{{\left\| {{h_1}(\sigma )} \right\|}_{{L^3}}}{{\left\| {{\nabla ^k}\sigma } \right\|}_{{L^6}}} + \parallel \nabla \sigma {\parallel _{{L^6}}} \cdot } \right.} \\ {\left. {{{\left\| {{\nabla ^{k - 1}}{h_1}(\sigma )} \right\|}_{{L^3}}}} \right)\left\| {{\nabla ^{k + 1}}\boldsymbol{u}} \right\| \leqslant M\left( {\left\| {{\nabla ^{k + 1}}\boldsymbol{u}} \right\| + } \right.} \\ {\left. {\left\| {{\nabla ^k}\sigma } \right\|} \right)\left\| {{\nabla ^{k + 1}}\boldsymbol{u}} \right\| \leqslant M\left( {{{\left\| {{\nabla ^k}\sigma } \right\|}^2} + {{\left\| {{\nabla ^{k + 1}}\boldsymbol{u}} \right\|}^2}} \right)} \end{array}$
$\begin{aligned}& I_2^3=\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-1}\left\{2 h_2(\sigma) \operatorname{div}[\nu(\phi) \boldsymbol{D}(\boldsymbol{u})]\right\} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u} \mathrm{d} \boldsymbol{x}+ \\& \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-1}\left\{h_2(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}[(\nu(\phi)+\lambda(\phi)) \operatorname{div} \boldsymbol{u}]\right\} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u} \mathrm{d} \boldsymbol{x}= \\& \sum\limits_{0 \leqslant l \leqslant k-1} C_{k-1}^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l h_2(\sigma) \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l-1}(\nu(\phi) \boldsymbol{D}(\boldsymbol{u})) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u} \cdot \\& \mathrm{d} \boldsymbol{x}+\sum\limits_{0 \leqslant l \leqslant k-1} C_{k-1}^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l h_2(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l}((\nu(\phi)+\lambda(\phi)) \cdot \\& \operatorname{div} \boldsymbol{u}) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u} \mathrm{d} \boldsymbol{x}\end{aligned}$
$\begin{aligned}& \sum\limits_{0 \leqslant l \leqslant k-1} C_{k-1}^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l h_2(\sigma) \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l-1}(\nu(\phi) \boldsymbol{D}(\boldsymbol{u})) \cdot \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u} \mathrm{d} \boldsymbol{x} \leqslant \sum\limits_{0 \leqslant l \leqslant k-1}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l h_2(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l}(\nu(\phi) \boldsymbol{D}(\boldsymbol{u}))\right\| \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\| \leqslant \sum\limits_{0 \leqslant l \leqslant k-1}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma\right\|_{L^3}\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\|_{L^6} \cdot\right. \\& \left.\|\nu(\phi)\|_{L^{\infty}}+\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}\|_{L^{\infty}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \phi\right\|_{L^6}\right)\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\| \leqslant \\& \sum\limits_{0 \leqslant l \leqslant k-1}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\|_{L^6}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|+ \\& \sum\limits_{0 \leqslant l \leqslant k-1}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \sigma\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \phi\right\|_{L^6}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\| \leqslant \\& M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^2\right)\end{aligned}$
用同样的方法可以对
$\begin{aligned}& \sum\limits_{0 \leqslant l \leqslant k-1} C_{k-1}^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l h_2(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l}((\nu(\phi)+\lambda(\phi)) \cdot\\& \operatorname{div} \boldsymbol{u}) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u} \mathrm{d} \boldsymbol{x} \\&\end{aligned}$
进行估计,然后可得
$I_2^3 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^2\right)$
$I_2^4=\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u} \mathrm{d} \boldsymbol{x}$
k=0时,$I_2^4 = \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi {\rm{\cdot}}\mathit{\boldsymbol{u}}{\rm{d}}\mathit{\boldsymbol{x}} \le \parallel \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi {\parallel _{{L^3}}}\parallel \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \parallel {\rm{\cdot}}\parallel {\boldsymbol{u}}{\parallel _{{L^6}}} \le \parallel \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi {\parallel ^{\frac{1}{2}}}{\left\| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^2}\phi } \right\|^{\frac{1}{2}}}\left\| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^2}\phi } \right\|\parallel \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mathit{\boldsymbol{u}}\parallel \le M\left( {{{\left\| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^2}\phi } \right\|}^2} + \parallel \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mathit{\boldsymbol{u}}{\parallel ^2}} \right)$
k=1时,$I_2^4 = \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}} (\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi ) \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mathit{\boldsymbol{u}}{\rm{d}}\mathit{\boldsymbol{x}} = - \int {(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi )} \cdot {\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^2}\mathit{\boldsymbol{u}}{\rm{d}}\mathit{\boldsymbol{x}} \le \parallel \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi {\parallel _{{L^3}}}\parallel \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi {\parallel _{{L^6}}}\left\| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^2}\mathit{\boldsymbol{u}}} \right\| \le \parallel \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi {\parallel ^{\frac{1}{2}}}{\left\| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^2}\phi } \right\|^{\frac{1}{2}}}\left\| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^3}\phi } \right\|\left\| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^2}\mathit{\boldsymbol{u}}} \right\| \le M \cdot \left( {{{\left\| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^2}\phi } \right\|}^2} + \parallel \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \mathit{\boldsymbol{u}}{\parallel ^2}} \right)$
k≥2时,$I_2^4 = \int {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^k}} (\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi ) \cdot {\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^k}\mathit{\boldsymbol{u}}{\rm{d}}\mathit{\boldsymbol{x}} = - \int {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^{k - 1}}} \cdot (\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi ) \cdot {\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^{k + 1}}\mathit{\boldsymbol{u}}{\rm{d}}\mathit{\boldsymbol{x}} = - \sum\limits_{0 \le l \le k - 1} {C_{k - 1}^l} \int {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^{l + 1}}} \phi \cdot {\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^{k - l + 1}}\phi \cdot {\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^{k + 1}}\mathit{\boldsymbol{u}}{\rm{d}}\mathit{\boldsymbol{x}} \le \sum\limits_{0 \le l \le k - 1} {\left\| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^{l + 1}}\phi \cdot {\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^{k - l + 1}}\phi } \right\|} \cdot \left\| {{\mathit{\boldsymbol{ \boldsymbol{\nabla}}} ^{k + 1}}\mathit{\boldsymbol{u}}} \right\|$
$\begin{aligned}& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{l+1} \phi \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi\right\| \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{l+1} \phi\right\|_{L^3} \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \phi\right\|_{L^6} \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \phi\right\|^{1-\frac{l}{k}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^{\frac{l}{k}} . \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^2 \phi\right\|^{\frac{l}{k}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^{1-\frac{l}{k}} \leqslant M\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|\end{aligned}$
其中α$\frac{l}{3} = \left( {\frac{\alpha }{3} - \frac{1}{2}} \right)\left( {1 - \frac{l}{k}} \right) + \left( {\frac{{k + 2}}{3} - \frac{1}{2}} \right){\rm{\cdot}}\frac{l}{k}$确定。
$I_2^4 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^2\right)$
$\begin{array}{*{20}{l}} {I_2^5 = \sum\limits_{0 \leqslant l \leqslant k} {C_k^l} \int {{\nabla ^l}} \nu (\phi )\boldsymbol{D}\left( {{\nabla ^{k - l}}\boldsymbol{u}} \right):\nabla {\nabla ^k}\boldsymbol{u}{\text{d}}\boldsymbol{x} \leqslant } \\ {\sum\limits_{0 \leqslant l \leqslant k} {{{\left\| {{\nabla ^l}\phi } \right\|}_{{L^3}}}} \left\| {{\nabla ^{k + 1}}\boldsymbol{u}} \right\|{{\left\| {{\nabla ^{k - l + 1}}\boldsymbol{u}} \right\|}_{{L^6}}}} \end{array}$
其中
$\begin{aligned}& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \phi\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\|_{L^6} \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \phi\right\|^{1-\frac{l-1}{k}} \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^{\frac{l-1}{k}}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}\|^{\frac{l-1}{k}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^{1-\frac{l-1}{k}} \leqslant M \cdot \\& \left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|\right)\end{aligned}$
其中α$\frac{{l - 1}}{3} = \left( {\frac{\alpha }{3} - \frac{1}{2}} \right)\left( {1 - \frac{{l - 1}}{k}} \right) + \left( {\frac{{k + 2}}{3} - } \right.\left. {\frac{1}{2}} \right)\frac{{l - 1}}{k}$确定。
$I_2^5 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \phi\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^2\right)$
用同样的方法可得
$\begin{aligned}& I_2^6=\frac{1}{\bar{\rho}} \sum\limits_{1 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l(\nu(\phi)+\lambda(\phi)) \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \cdot \\& \boldsymbol{u} \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u} \mathrm{d} \boldsymbol{x} \leqslant M\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^2\end{aligned}$
$\begin{aligned}& I_2^7=\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k\left(h_2(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\right) \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u} \mathrm{d} \boldsymbol{x} \leqslant M\left(\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \cdot\right. \\& \left.\phi\left\|^2+\right\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u} \|^2\right)\end{aligned}$
$I_2^i(i=1, 2, \cdots, 7)$代入式(30)可得式(28),引理3得证。
引理4  若(σ, u, φ)∈XM([0, T])是方程(10)的局部解,则
$\begin{aligned}& \quad \frac{\mathrm{d}}{\mathrm{d} t} \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x}+\frac{p^{\prime}(\bar{\rho})}{2 \bar{\rho}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2 \leqslant \\& M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|^2\right)+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^2, k= \\& 0, 1, 2\end{aligned}$
证明:对方程组(29)的第二式乘以$\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma$,并利用(29)的第一式可得
$\begin{array}{*{20}{l}} {\frac{{\text{d}}}{{{\text{d}}t}}\int {{\nabla ^k}} \boldsymbol{u} \cdot {\nabla ^{k + 1}}\sigma {\text{d}}\boldsymbol{x} + \frac{{{p^\prime }(\bar \rho )}}{{\bar \rho }}{{\left\| {{\nabla ^{k + 1}}\sigma } \right\|}^2} - \bar \rho \cdot } \\ {\int {{{\left( {\operatorname{div} {\nabla ^k}\boldsymbol{u}} \right)}^2}} {\text{d}}\boldsymbol{x} \leqslant I_3^1 + I_3^2 + I_3^3 + I_3^4 + I_3^5 + I_3^6 + I_3^7} \end{array}$
式中,
$\begin{aligned}& I_3^1=\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \operatorname{div}(\boldsymbol{u} \sigma) \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k \boldsymbol{u} \mathrm{d} \boldsymbol{x} \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\| \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1}(\boldsymbol{u} \sigma)\right\| \leqslant\left(\|\boldsymbol{\sigma}\|_{L^{\infty}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\| \cdot\right. \\& \left.\|\boldsymbol{u}\|_{L^{\infty}}\right)\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\| \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{u}\right\|^2\right)\end{aligned}$
$\begin{aligned}& I_3^2=\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k\left[-(\boldsymbol{u}, \mathit{\boldsymbol{ \boldsymbol{\nabla}}}) \boldsymbol{u}+h_1(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{\sigma}\right] \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x}= \\& \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k[-(\boldsymbol{u}, \mathit{\boldsymbol{ \boldsymbol{\nabla}}}) \boldsymbol{u}] \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x}+\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k\left[h_1(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \sigma\right] \cdot \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x}\end{aligned}$
其中
$\begin{aligned}& \quad \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k[-(\boldsymbol{u}, \mathit{\boldsymbol{ \boldsymbol{\nabla}}}) \boldsymbol{\mu}] \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant \sum\limits_{0 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{u} \cdot \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant M \sum\limits_{0 \leqslant l \leqslant k}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\| \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \boldsymbol{\sigma}\right\|\end{aligned}$
$\begin{aligned}& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\| \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \boldsymbol{u}\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\|_{L^6} \leqslant\\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \boldsymbol{u}\right\|^{1-\frac{l}{k+1}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^{\frac{l}{k+1}}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}\|^{\frac{l}{k+1}} \| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \cdot\\& \boldsymbol{u}\left\|^{1-\frac{l}{k+1}} \leqslant M\right\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u} \|\end{aligned}$
其中α$\frac{{l - 1}}{3} = \left( {\frac{\alpha }{3} - \frac{1}{2}} \right)\left( {1 - \frac{l}{{k + 1}}} \right) + \left( {\frac{{k + 2}}{3} - } \right.\left. {\frac{1}{2}} \right)\frac{l}{{k + 1}}$确定。
由同样的方法可得
$\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k\left[h_1(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \sigma\right] \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant M\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2$$I_3^2 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^2\right)$
$\begin{aligned}& I_3^3=-\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k\left\{2 h_2(\sigma) \operatorname{div}[\nu(\phi) \boldsymbol{D}(\boldsymbol{u})]\right\} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \cdot \\& \sigma \mathrm{d} \boldsymbol{x}-\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k\left\{h_2(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}[(\nu(\phi)+\lambda(\phi)) \operatorname{div} \boldsymbol{u}]\right\} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \cdot \\& \mathrm{d} \boldsymbol{x}=-\sum\limits_{0 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l h_2(\sigma) \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1}[\nu(\phi) \boldsymbol{D}(\boldsymbol{u})] \cdot \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x}-\sum\limits_{0 \leqslant l \leqslant k} C_{k-1}^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l h_2(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l}[(\nu(\phi)+ \\& \lambda(\phi)) \operatorname{div} \boldsymbol{u}] \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|^2\right)\end{aligned}$
$\begin{aligned}& I_3^4=-\frac{1}{\bar{\rho}} \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k(\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi) \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant \\& -\sum\limits_{0 \leqslant l \leqslant k} C_k^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{l+1} \phi \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant \sum\limits_{0 \leqslant l \leqslant k} \| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{l+1} \phi \cdot \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2} \phi\|\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \|\end{aligned}$
其中
$\begin{aligned}& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{l+1} \phi \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2} \phi\right\| \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{l+1} \phi\right\|_{L^3} \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2} \phi\right\|_{L^6} \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \phi\right\|^{1-\frac{l+1}{k+2}} \mid\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|^{\frac{l+1}{k+2}} \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^2 \phi\right\|^{\frac{l+1}{k+2}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|^{1-\frac{l+1}{k+2}} \leqslant M\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|\end{aligned}$
其中α$\frac{l}{3} = \left( {\frac{\alpha }{3} - \frac{1}{2}} \right)\left( {1 - \frac{{l + 1}}{{k + 2}}} \right) + \left( {\frac{{k + 3}}{3} - } \right.\left. {\frac{1}{2}} \right)\frac{{l + 1}}{{k + 2}}$确定。
$I_3^4 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2\right)$
$\begin{aligned}& I_3^5=\frac{2}{\bar{\rho}} \int \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k(\nu(\phi) \boldsymbol{D}(\boldsymbol{u})) \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant \\& \int \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k(\nu(\phi) \boldsymbol{D}(\boldsymbol{u})) \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x}=\sum\limits_{0 \leqslant l \leqslant k} C_{k+1}^l \int\left(\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \cdot\right. \\& \left.\nu(\phi) \operatorname{div} \boldsymbol{D}\left(\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \boldsymbol{u}\right)+\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{l+1} \nu(\phi) \boldsymbol{D}\left(\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \boldsymbol{u}\right)\right) \cdot \\& \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant \sum\limits_{1 \leqslant l \leqslant k}\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \phi\right\|_{L^3} \quad\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+2} \boldsymbol{u}\right\|_{L^6}+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{l+1} \phi\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\|_{L^6}\right)\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\| \leqslant M\left(\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \cdot\right. \\& \left.\boldsymbol{u}\left\|^2+\right\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \|^2\right)\end{aligned}$
$\begin{aligned}& I_3^6=\frac{1}{\bar{\rho}} \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1}\{[\nu(\phi)+\lambda(\phi)] \operatorname{div} \boldsymbol{u}\} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant \\& \sum\limits_{0 \leqslant l \leqslant k} C_{k+1}^l \int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l[\nu(\phi)+\lambda(\phi)] \cdot \operatorname{div} \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l} \boldsymbol{u} \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant \\& \sum\limits_{0 \leqslant l \leqslant k}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \phi\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\|_{L^6}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|\end{aligned}$
其中
$\begin{aligned}& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^l \phi\right\|_{L^3}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k-l+1} \boldsymbol{u}\right\|_{L^6} \leqslant\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^\alpha \phi\right\|^{1-\frac{l}{k}} \cdot \\& \left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|^{\frac{l}{k}}\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}} \boldsymbol{u}\|^{\frac{l}{k}}\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^{1-\frac{l}{k}} \leqslant M\left(\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \cdot\right. \\& \left.\phi\|+\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u} \|\right)\end{aligned}$
其中α$\frac{{l - 1}}{3} = \left( {\frac{\alpha }{3} - \frac{1}{2}} \right)\left( {1 - \frac{l}{k}} \right) + \left( {\frac{{k + 3}}{3} - } \right.\left. {\frac{1}{2}} \right)\frac{l}{k}$确定。
$\begin{aligned}& I_3^6 \leqslant M\left(\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \phi\right\|^2+\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma\right\|^2+\right. \\& \left.\left\|\mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+2} \boldsymbol{u}\right\|^2\right)\end{aligned}$
$\begin{aligned}& I_3^7=\int \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^k\left[h_2(\sigma) \mathit{\boldsymbol{ \boldsymbol{\nabla}}} \phi \mathit{\boldsymbol{ \boldsymbol{\Delta}}} \phi\right] \cdot \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \mathrm{d} \boldsymbol{x} \leqslant M\left(\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+3} \cdot\right. \\& \left.\phi\left\|^2+\right\| \mathit{\boldsymbol{ \boldsymbol{\nabla}}}^{k+1} \sigma \|^2\right)\end{aligned}$
$I_3^i(i=1, 2, \cdots, 7)$代入式(32)可得式(31),引理4得证。结合引理1~4,命题2得证,进一步定理1得证。

3 结束语

本文研究了三维Navier-Stokes-Cahn-Hilliard方程组Cauchy问题的适定性。在初值小扰动的条件下,我们解决了相场φ在估计过程中带来的一些困难,证明了该方程全局解的存在唯一性。该结果表明,在小扰动情况下,相分离状态保持不变。本文结果可为非混相两相流的模拟计算提供理论基础。
1
CAHN J W , HILLIARD J E . Free energy of a nonuniform system. I. Interfacial free energy[J]. The Journal of Chemical Physics, 1958, 28, 258- 267.

DOI

2
LOWENGRUB J , TRUSKINOVSKY L . Quasi-incompres-sible Cahn-Hilliard fluids and topological transitions[J]. Proceedings of the Royal Society of London A, 1998, 454, 2617- 2654.

DOI

3
ABELS H , FEIREISL E . On a diffuse interface model for a two-phase flow of compressible viscous fluids[J]. Indiana University Mathematics Journal, 2008, 57 (2): 659- 698.

DOI

4
OZ·AN·SKI W S, POOLEY B C. Leray's fundamental work on the Navier-Stokes equations: a modern review of "Sur le mouvement d'un liquide visqueux emplissant l'espace"[EB/OL]. (2017-08-31). arXiv: 1708.09787v1 [math. AP].

5
CHEN Y Z , HE Q L , MEI M , et al. Asymptotic stability of solutions for 1-D compressible Navier-Stokes-Cahn-Hilliard system[J]. Journal of Mathematical Analysis and Applications, 2018, 467, 185- 206.

DOI

6
CHEN Y Z, HONG H, SHI X D. Asymptotic stability of phase separation states for compressible immiscible two-phase flow with periodic boundary condition in 3D[J/OL]. (2021-09-02). arXiv: 2105.13552v3 [math. AP].

7
王暐翼, 童天骄, 陈亚洲. 一维Navier-Stokes-Cahn-Hilliard方程组解的适定性分析[J]. 北京化工大学学报(自然科学版), 2019, 46 (6): 101- 107.

WANG W Y , TONG T J , CHEN Y Z . Well-posedness of solutions for Navier-Stokes-Cahn-Hilliard system in one dimension[J]. Journal of Beijing University of Chemical Technology (Natural Science), 2019, 46 (6): 101- 107.

8
HEIDA M , MÁLEK J , RAJAGOPAL K R . On the development and generalizations of Allen-Cahn and Stefan equations within a thermodynamic framework[J]. Zeitschrift für Angewandte Mathematik und Physik, 2012, 63, 759- 776.

DOI

9
VALLI A , ZAJACZKOWSKI W M . Navier-Stokes equations for compressible fluids: global existence and qualitative properties of the solutions in the general case[J]. Communications in Mathematical Physics, 1986, 103, 259- 296.

Outlines

/