On the non-stationary Navier-Stokes flows and reiterated homogenization

技术摘要：关于非定常 Navier-Stokes 型流动与重迭代均化

问题陈述
本文研究了定义在固定有界区域 $\Omega \subset \mathbb{R}^N$（具体为 $N=2$）上、时间区间 $]0, T[$ 内的非定常（时间依赖）Navier-Stokes 型方程的确定性重迭代均化问题。控制方程描述了具有周期性快速变化粘度系数的流体流动。这些系数的特定标度涉及两个不同的微观尺度 $\varepsilon$ 和 $\varepsilon^2$，导致如下定义：
$$a^\varepsilon_{ij}(x) = a_{ij}\left(\frac{x}{\varepsilon}, \frac{x}{\varepsilon^2}\right)$$
其中 $a_{ij}$ 是定义在两个周期性单元 $Y \times Z$ 的乘积上的连续、有界、对称且一致椭圆函数。本研究关注速度 $u^\varepsilon$ 和压力 $p^\varepsilon$ 的柯西 - 狄利克雷边值问题，其特征方程为：
$$\frac{\partial u^\varepsilon}{\partial t} + P^\varepsilon u^\varepsilon + \sum_{j=1}^N u^\varepsilon_j \frac{\partial u^\varepsilon}{\partial x_j} + \nabla p^\varepsilon = f, \quad \text{div } u^\varepsilon = 0$$
并配有零初始条件和边界条件。主要目标是确定当参数 $\varepsilon \to 0$ 时解 $(u^\varepsilon, p^\varepsilon)$ 的宏观行为。

方法论
分析依赖于重迭代双尺度收敛方法，这是标准双尺度收敛技术的扩展，旨在处理多重振荡尺度。该方法按以下步骤进行：

1. 函数框架：问题在适当的索伯列夫空间中表述，特别是 $W(0, T)$，该空间连续嵌入到 $C(0, T; H)$ 中，并紧嵌入到 $L^2(0, T; H)$ 中。利用空间 $V$（无散度向量场）及其对偶空间来处理不可压缩约束。
2. 先验估计：利用与算子 $P^\varepsilon$ 相关的双线性形式的强制性以及代表对流项的三线性形式的性质，作者在 $W(0, T)$ 中为解序列 $(u^\varepsilon)$ 以及在 $L^2(0, T; L^2(\Omega))$ 中为 $(p^\varepsilon)$ 建立了均匀有界性。
3. 紧性与收敛性：利用重迭代双尺度收敛的紧性结果（定理 1 和定理 2），作者提取了一个收敛于极限三元组 $(u_0, u_1, u_2)$ 的子序列。其中，$u_0$ 代表宏观极限，而 $u_1$ 和 $u_2$ 分别代表与尺度 $y = x/\varepsilon$ 和 $z = x/\varepsilon^2$ 相关的一阶和二阶微观修正项。
4. 均化模型的推导：通过使用特定的振荡测试函数对变分形式进行检验并取极限，作者推导出一组耦合方程。这涉及在周期性单元 $Y$ 和 $Z$ 上求解局部单元问题，以定义修正函数 $\chi_{ik}$ 和 $\eta_{ik}$。
5. 修正项分析：在关于区域和数据的额外正则性假设下，本文推导了修正项结果，证明了原始解的梯度可以通过宏观解的梯度与缩放后的微观修正项梯度之和来近似。

主要贡献与结果

• 宏观均化模型：本文推导了一个唯一的 Navier-Stokes 型宏观极限问题（定理 5）。宏观速度 $u_0$ 和压力 $p_0$ 的极限方程为：
  $$\frac{\partial u_0}{\partial t} + Q u_0 + \sum_{j=1}^N u^0_j \frac{\partial u_0}{\partial x_j} + \nabla p_0 = f$$
  其中 $Q$ 是一个具有常数系数 $q_{ijkh}$ 的均化微分算子。这些系数通过涉及原始系数 $a_{ij}$ 以及在 $Y$ 和 $Z$ 上单元问题解的积分显式定义。
• 唯一性：本文证明了所推导的均化问题解的唯一性（命题 8），确保了宏观模型的适定性。
• 收敛定理：证明了 $u^\varepsilon$ 在 $W(0, T)$ 中弱收敛于 $u_0$，且梯度以双尺度意义重迭代收敛于宏观梯度与微观梯度之和（定理 4）。
• 修正项结果：建立了梯度的强收敛结果（定理 6）。具体而言，在适当的正则性假设下，实际梯度 $\nabla u^\varepsilon$ 与近似值 $\nabla u_0 + \varepsilon \nabla_y u_1 + \varepsilon^2 \nabla_z u_2$ 之间的误差在 $L^2$ 范数下随 $\varepsilon \to 0$ 而消失。

意义与范围
作者将这项工作定位为对具有复杂多尺度粘度结构的 Navier-Stokes 型流动数学分析的贡献。本文声称是首个在确定性设定下，针对特定标度 $a_{ij}(x/\varepsilon, x/\varepsilon^2)$ 解决非定常（非稳态）非线性流体流动的重迭代均化问题。

该研究的动机源于对具有空间变化粘度的多相流进行建模的需求。虽然作者指出该方法可适应更一般（非周期性）的粘度行为，但当前工作严格局限于周期性设定。

本文在时间依赖系数方面保持了适度的范围。作者明确指出，研究假设粘度系数 $a_{ij}$ 仅依赖于空间变量，而不依赖于时间。结论承认，将这些结果扩展到时间依赖粘度是一项困难的任务，因为文献中缺乏此类情形的既定唯一性和估计；同时指出，当前的均化过程只要能够建立这些估计和唯一性结果，仍然有效。

这项工作建立在作者此前关于稳态流动和单尺度均化的研究，以及其他研究人员（引用为 [5]、[6]、[8]）发展的基础重迭代双尺度收敛理论之上。这些结果为用更简单的有效宏观模型近似复杂、快速振荡的流体动力学提供了严格的数学基础。