Brinkman's term for the steady-state Navier-Stokes equations with prescribed flux rate or pressure drop in perforated pipes

技术摘要：带混合边界条件的穿孔管道中稳态 Navier-Stokes 方程的均质化研究

问题陈述
本研究探讨了在含有大量微小固体颗粒的扭曲管道内，粘性不可压缩流体的渐近行为。该区域 $\Omega_\varepsilon$ 是一个有界容许区域 $\Omega \subset \mathbb{R}^3$ 的穿孔子集，通过移除 $N(\varepsilon)$ 个直径为 $\varepsilon^3$、间距为 $\varepsilon$ 的颗粒得到。这种配置对应于均质化理论中的临界状态（$\alpha = 3$），即颗粒尺寸随颗粒间距的立方进行缩放。

流体流动受带有混合边界条件的稳态 Navier-Stokes 方程支配。具体而言：

1. 无滑移条件施加在物理管壁（$\Gamma_W$）以及颗粒表面（$\partial K_n^\varepsilon$）上，其中颗粒可能以给定的常数速度 $\mu_n^\varepsilon$ 运动。
2. 伯努利压力条件应用于入口（$\Gamma_I$）和出口（$\Gamma_O$）。这些条件涉及速度的法向导数以及伯努利压力（$\Phi = p + \frac{1}{2}|u|^2$）。
3. 研究分析了两种不同的情景：
   • 给定净通量： 出口处的总通量 $F$ 是固定的，而压降是一个由解决定的未知常数。
   • 给定压降： 入口与出口之间的压降是固定的，而由此产生的通量是一个未知量。

主要的数学挑战在于，如何在不施加对外部力、通量率或压降之微小性假设的情况下，推导出当 $\varepsilon \to 0$ 时的有效（均质化）方程。由于在这些混合边界值问题中，标准的能量估计往往无法提供一致的界限，这使得问题变得更加复杂。

方法论
作者在 Tartar 的一般框架内采用了能量法。证明过程分为三个主要阶段：

1. 一致界限（$\varepsilon$-无关估计）：
   建立解 $(u_\varepsilon, \Phi_\varepsilon)$ 的一致界限是最关键且最具创新性的步骤。由于混合边界条件和未知的压降，作者不能像标准情况那样仅仅依靠用解本身作为测试函数来进行能量估计。

   • 相反，他们利用了基于稳态 Euler 方程解的伯努利定律，采用了一种矛盾论证法。
   • 他们假设解的范数在 $\varepsilon \to 0$ 时趋于无穷大，对序列进行归一化，并证明其极限满足 Euler 方程。
   • 通过分析极限伯努利压力在管壁上的边界行为（利用局部正则性结果和特定的入口/出口几何结构），他们证明了极限压力在壁面上必须是常数。这与通量或压降条件产生了矛盾，从而证明了序列的一致有界性。
   • 对于给定通量问题，提出了一个额外的几何假设，即颗粒不会在入口和出口圆柱附近聚集，以便于进行边界压力分析。

2. 构建修正子与限制算子：
   为了处理穿孔区域，作者构造了位于颗粒外部及周围环形区域的 Stokes 修正子 $(\omega_{n,1}^\varepsilon, \omega_{n,2}^\varepsilon)$。这些修正子被用于定义一个限制算子 $R_\varepsilon$。

   • 不同于以往在亚临界状态下仅需标量电容修正子的工作，临界状态需要这些矢量值的 Stokes 修正子，以解释颗粒群非零的相对电容。
   • 限制算子 $R_\varepsilon(\phi)$ 通过修改测试函数在颗粒附近的性质来将来自均质化区域 $\Omega$ 的测试函数映射到穿孔区域 $\Omega_\varepsilon$，使其满足零边界条件，同时保持散度为零的特性。

3. 取极限过程：
   利用一致界限和限制算子，作者在 Navier-Stokes 方程的弱形式中取极限。

   • 非线性对流项通过 $L^4$ 中的强收敛性实现收敛（由一致界限和紧嵌入保证）。
   • 流体与颗粒之间的相互作用产生了一个极限方程中的附加项。该项源于颗粒运动产生的拖曳力的极限行为，其特征由 Stokes 阻力矩阵 $G_n$ 表征。

主要贡献与结果
本文建立了描述系统均质极限的两个主要定理（定理 2.4 和 2.5）：

• Brinkman 项的出现： 极限区域 $\Omega$ 中的均质方程是增加了 Brinkman 项后的 Navier-Stokes 方程。有效方程的形式为：
  $$-\Delta u + (\nabla \times u) \times u + \nabla \Phi + Gu - J = f$$
  其中：

  • $Gu$代表与流体速度$u$成正比的**摩擦项**，它源于颗粒的集体拖曳作用（其中$G$ 是缩放后的 Stokes 阻力矩阵之和的极限）。
  • $J$ 是一个源项，源于运动颗粒自身的极限通量。
  • 入口和出口的边界条件仍为涉及伯努利压力的混合类型。

• 非零数据下的临界状态扩展： 研究结果将以往关于亚临界状态的均质化发现（特别是文献 [59]）扩展到了临界状态（$\alpha=3$）。至关重要的是，这些结果在不对外部力、通量率或压降做任何微小性假设的情况下依然成立。这是对许多依赖小数据假设来确保唯一性或存在有界性的 Navier-Stokes 均质化结果的重要推进。

• 处理运动颗粒： 该框架可以容纳具有非零常数速度 $\mu_n^\varepsilon$ 的颗粒，前提是其动能是有界一致的（假设 2.11）。与固定颗粒模型相比，这允许对颗粒-流体相互作用进行更真实的建模。

意义
本文声称其最重要的贡献是在不受限数据下，在临界状态下严格推导了稳态 Navier-Stokes 方程的 Brinkman 定律。作者强调，通过涉及伯努利定律的矛盾论证来建立一致界限是一个极其精细的过程，特别是必须考虑到穿孔的复杂结构，以确保该论证在 $\varepsilon$ 层级上依然成立。

通过成功构造一个不会修改在颗粒处消失的测试函数的真正限制算子，并处理了在具有非零颗粒速度的多分散设置下压力的收敛问题，这项工作为建模穿孔管道中的流体流动提供了坚实的数学基础。在这种情况下，颗粒密度足够高以诱发 Brinkman 型拖曳，但又不足以导致达西定律（Darcy's law）。其结果可直接应用于血液动力学、石油工程和空气动力学等领域，在这些领域中，流体流经固体障碍物阵列是一个基本现象。