无滑移边界与大温度变化下低马赫数极限如何表现？

极限系统保留散度约束而非不可压缩性

论文研究三维光滑有界区域中的完全可压缩Navier-Stokes方程，未知量为密度、速度、内能、温度和压力，并采用理想多方气体状态方程。马赫数 ε 是特征流速与特征声速之比，目标是描述当 ε 趋于零时会发生什么。

形式渐近展开表明，解收敛到一个满足低马赫数系统的极限 (r, ω, π)。其第三个方程为 γCv div ω = div(κ∇(1/r))，其中 γ 是比热比，κ 是热导率。当不存在热传导，或极限温度在空间上为常数时，右端项消失，约束退化为 div ω = 0，即熟悉的不可压缩极限。在存在热传导且温度变化很大时，约束保持非平凡：论文的收敛定理给出 div ω = (κ/2)Δθ，同时还有 (1+θ)^{-1}(θ_t + ω·∇θ) = (κ/2)Δθ。因此低马赫数近似滤除了声学可压缩性，却保留了温度与热传导效应，这正是该一般模型比不可压缩极限更难的原因。

大温度变化破坏斜对称消去结构

对于等熵流，压力只依赖于密度，诸如 ρ^ε = 1 + ερ^ε_1 或 ρ^ε = e^{ερ^ε_1} 的变量替换会产生形如 ε^{-1}(div, ∇) 的斜对称首阶算子。该结构使奇异项在 L² 能量恒等式中相互抵消，从而给出 ε 一致估计。

当 κ = 0 时，同样的斜对称形式 ε^{-1}(div, ∇) 仍然可用，奇异项依然抵消。但当 κ > 0 且温度变化很大时，惩罚算子在原始变量下不再斜对称，因为奇异压力项同时与密度和温度发生非线性耦合。论文引入变量替换 p^ε = 1 + εq^ε 与 Θ^ε = 1 + θ^ε，随后引入修正速度 v^ε := u^ε − (κ/2)∇θ^ε，从而恢复奇异声学部分的斜对称性。该变换并未使 v 与温度完全解耦：v 的方程仍包含三阶项 −β∇Δθ^ε，而在无滑移条件下 v 的切向迹等于 −(κ/2)∇_τθ^ε，因此依赖于切向温度梯度。

无滑移壁面移除了涡量边界关系

在有界区域中，声波与边界相互作用，因此全空间或周期情形中使用的基于Fourier的方法不能直接使用，而且法向求导会产生基本能量估计无法控制的边界项。因此高阶法向导数无法通过直接求导来估计。

论文明确指出，无滑移条件无法提供滑移边界情形中所用的涡量边界关系。在涡量滑移和Navier滑移情形中，高阶估计通过 v 的 div-curl 分解闭合，而 curl 部分的估计关键依赖于 curl v 所满足的边界条件。在无滑移条件下，该机制不可用，因此那些情形中得到的高阶法向估计无法直接迁移到本文设定。论文指出，无滑移边界条件与大温度变化下完全Navier-Stokes系统的低马赫数问题此前仍是未解决的。

无散度Stokes提升闭合了估计

核心技术手段是对修正速度的时间导数 v_t 的边界迹作无散度Stokes提升。将 v_t 作为微分后动量方程的检验函数会产生一个涉及 (n·∇v_t)·v_t 的边界项，该项无法控制。为消除它，作者求解一个边界数据为 g = −(κ/2)∇_τθ_t、ℓ_0 = 0 的Stokes问题，得到一个提升 ℓ，其存在性、唯一性和估计由所引用的Stokes理论给出。

随后他们使用 Z := v_t − ℓ 作为新的检验函数。按构造，Z 在边界上为零，且在区域内满足 div Z = div v_t。由于 ℓ 无散度且法向迹为零，该修正不会改变奇异压力项的消去。将所得能量估计与 ℓ 的界结合，即得到所需的 v_t 的 ε 一致估计。论文还指出，这种构造特殊检验函数的方法也适用于其他问题，并引用了别处基于Stokes和基于Bogovskii的构造。

两个斜对称结构把估计分成两块

证明结合了两个斜对称结构。一个将温度与加权未知量耦合，给出温度时间导数的一致估计；另一个将压力与修正速度耦合，控制速度的时间导数与高阶空间导数。具体而言，由 (q, v, θ) 满足的系统用于同时估计 ∇²div v 与 ∇²Δθ 的高阶范数，以及 v_t 在 L^∞_t(L²) 与 L²_t(H¹) 中的低阶范数。第二个系统用于估计 θ_t 在 L^∞_t(H¹) 与 L²_t(H²) 中的范数。

作者引入指数变换 p = e^{εσ}、Θ = e^{ϑ}，它通过 εq = e^{εσ} − 1 与 θ = e^{ϑ} − 1 等价于线性变换，并令 η := (1/2)(εσ − ϑ)、u^ε := εu。在变量 (η, u^ε, ϑ) 下，主奇异部分变为斜对称，且与小温度变化情形具有相同形式。改写后的系统还有一个优点：ϑ_t 的系数为常数。估计通过量 M₁(t) 与 M₂(t) 组织，它们分别控制 (η_t, u^ε_t, ϑ_t) 以及按 ε 的幂加权的二阶时间导数。

一致界在独立于马赫数的时间区间上成立

一致估计定理假设初始数据在 (ρ^ε_0, u^ε_0, q^ε_0) 的 H² 范数、θ^ε_0 的 H³ 范数以及最高二阶时间导数的 L² 或 H¹ 范数上有界，并对 (div u^ε_t, ∇q^ε_t)(0) 与 (ε²u^ε_tt, ε²q^ε_tt, εθ^ε_tt)(0) 附加 ε 加权条件。它还假设 ρ^ε_0 ≥ M_0^{-1} > 0、Θ^ε_0 ≥ A > 0，其中正常数 M_0 与 A 独立于 ε，并要求 β > 0。

在这些假设下，存在正常数 ε_0、T 与 C（在所述意义下均独立于 ε），使得唯一强解对每个 ε ∈ (0, ε_0] 与 t ∈ (0, T] 满足 M^ε(t) ≤ C，且在时空柱上 Θ^ε(x,t) ≥ C^{-1}。能量泛函 M^ε(t) 汇集了 (ρ^ε, q^ε) 的 L^∞_t(H²) 范数、u^ε 的 L^∞_t(H¹) 范数、θ^ε 的 L^∞_t(H²) 范数、(ρ^ε)^{-1} 的 L^∞ 界、(ρ^ε_t, u^ε_t, q^ε_t) 的 L^∞_t(L²) 范数、θ^ε_t 的 L^∞_t(H¹) 范数，以及 (u^ε, ∇θ^ε) 的 L²_t(H³) 范数与 (u^ε_t, ∇θ^ε_t) 的 L²_t(H¹) 范数。关键在于 T 不随 ε 缩小，这正是后续紧性论证得以进行的原因。

紧性把极限过渡到低马赫数系统

收敛定理取满足一致有界性的解，并令 ε 趋于零。结果是 (ρ^ε, u^ε, θ^ε) ⇀ (r, ω, θ) 在 L^∞(0,T; H²(Ω)) 中弱-* 收敛，且 r(1+θ) = 1。证明记录了 ρ^ε、u^ε、θ^ε 及其时间导数在相应空间中的弱-* 极限，并使用Aubin-Lions-Simon引理将 (ρ^ε, θ^ε) 提升为 C([0,T], H^{2−δ}) 中的强收敛，将 u^ε 提升为 C([0,T], H^{1−δ}) 中的强收敛，对任意 δ ∈ (0,1) 成立。

在改写后的系统中取极限，得到低马赫数Navier-Stokes方程：(1+θ)^{-1}(ω_t + ω·∇ω) + ∇π = μΔω + (μ+λ)∇div ω，以及 (1+θ)^{-1}(θ_t + ω·∇θ) = (κ/2)Δθ 与 div ω = (κ/2)Δθ，边界条件为 ω = 0 与 ∂θ/∂n = 0。对于动量方程，取旋度可消去奇异项，得到 curl[r(ω_t + ω·∇ω) − μΔω − (μ+λ)∇div ω] = 0，由此可恢复出满足 ∇π ∈ L²(0,T; L²(Ω)) 的压力 π。极限系统是针对良准备初始数据陈述的，而局部存在性定理要求相容性条件 u^ε_0 = u^ε_t(0) = u^ε_tt(0) = 0，以及 Θ^ε_0、Θ^ε_t(0) 与 Θ^ε_tt(0) 的法向导数在边界上为零。

决定性差异在边界条件而非马赫数

论文的贡献最好被理解为一项边界条件结果，而非小马赫数结果。在所有这些问​​题中马赫数都被驱向零；变化的是区域是否有边界以及速度在边界上的行为。在全空间或周期设定中，频率局部化处理了大温度变化情形。在具有滑移型条件的区域中，涡量边界关系闭合了高阶估计。在无滑移条件下，这两条路径都不可用，而 v_t 边界迹的无散度Stokes提升正是替代它们的手段。

因此实际解读是有条件的。一致估计与向低马赫数系统的收敛是在三维光滑有界区域中、针对强解、在良准备初始数据、β > 0 的条件下、在长度独立于 ε 但并未断言可任意长的时间区间上建立的。极限系统本身保留了一个非平凡的散度约束，因此任何将该结果与经典不可压缩极限作比较的人，都应检查是否存在热传导以及极限温度在空间上是否为常数；若其中任一条件使右端项消失，约束就退化为 div ω = 0，问题也就不同了。该结果本身也不覆盖此无滑移、大温度变化设定下的非良准备初始数据，论文将其留作由其他方法处理的另一类情形。

本文引用的文献