IGA-LBM: Isogeometric lattice Boltzmann method

技术摘要：IGA-LBM：等几何格子玻尔兹曼方法

问题陈述
格子玻尔兹曼方法（LBM）是一种强大的中观尺度计算流体力学框架，因其内在的并行性和简单的边界处理而备受重视。然而，传统的 LBM 依赖于均匀的笛卡尔晶格，通过“阶梯状”表示来近似曲线边界。正如论文中所阐述的，这种几何离散化引入了显著误差：锯齿状的边界刻面会扭曲壁面法线和有效物体形状，从而传播伪造力并降低边界层预测精度。虽然存在各种缓解策略（例如，基于插值的处理、切割单元法、浸没边界技术），但它们往往会损害数值稳定性、引入插值伪影或产生高昂的计算成本。核心挑战在于开发一种既能保持 LBM 的简洁与高效，又能消除在弯曲定义域上几何离散化误差的框架。

方法论
作者提出了等几何-配置格子玻尔兹曼方法（IGA-LBM），该方法将等几何分析（IGA）与离散速度 BGK 系统相结合。该方法在两个基本方面背离了经典的“流向-碰撞”（stream-collide）模板：

1. 精确几何表示： 方法不再使用笛卡尔网格，而是使用 B 样条或非均匀有理 B 样条（NURBS）对物理域进行参数化。这使得能够直接从 CAD 定义中解析地精确表示曲线边界（如圆、翼型）。分布函数使用与几何形状相同的 NURB 基函数进行近似（等参数化构造）。
2. 强形式配置法： 该方法不再沿晶格链路传输粒子群，而是在贴体参数化网格上求解连续时间的 Boltzmann-BGK 方程。
   • 空间离散化： 分布函数在 Greville 点处进行配置。变换后的平流项通过高阶微分算子进行评估。在均匀网格内部，三次样条配置算子在代数上等同于四阶紧致有限差分（Padé）格式。
   • 边界处理： 在需要解析边界层而采用的壁面聚集网格上，朴素的样条配置法在壁面处会丧失精度。作者使用直接构建在物理节点上的、具有宽度为 5 的一致性边界 Fornberg 有限差分算子，取代了靠近边界处的配置导数。这恢复了直至壁面的均匀四阶精度。
   • 度量项： 逆变晶格速度通过从 NURBS 映射中获得的解析精确度量项计算得出。这确保了度量项不存在网格差分误差。
   • 时间积分： 半离散系统使用显式四阶 Runge–Kutta（RK4）积分器进行推进。
   • 稳定性控制： 由于中心配置算子是非耗散的，因此每步应用一次尺度选择性的隐式六阶低通滤波器，以抑制网格尺度的不稳定性，而不影响已解析的范围。
   • 边界条件： 该方法采用基于矩的一致性边界处理方法（无滑壁面、远场），而不是覆盖粒子群，因为链路级的反弹规则（bounce-back）与配置公式不兼容。

核心贡献与理论结果
论文为该方案提供了完整的理论分析支撑：

• 流体动力学极限： 该方案恢复了具有粘度 $\nu = c_s^2 \tau$ 的不可压缩 Navier–Stokes 方程。值得注意的是，由于该方法直接对连续时间方程进行积分（而非使用精确流向），它避免了经典的 $-\Delta t/2$ 粘度修正伪影。
• 精确自由流保持： 定理 7 证明，平流离散化可以在任意贴体网格上将均匀流保持至机器精度，而无需要求离散几何守恒律（GCL）。这是因为无论度量精度如何，微分算子都会使常数项为零。
• 精度阶数： 内部配置算子被证明是四阶精确的（引理 4）。定理 6 确立了物理节点有限差分闭合在聚集网格上可恢复直至壁面的均匀四阶精度，修正了朴素配置导致的阶数退化。
• 度量精确性： 命题 9 表明，使用解析精确的度量项可以保持映射过程中方向导数离散化的形式阶数。消融研究显示，使用较低阶的差分度量（即使是四阶）如果不能完美匹配，会悄悄降低方案性能，而精确的 NURB 度量消除了这一误差源。
• 稳定性约束： 定理 10 确定了显式时间步长受限于刚性碰撞项，即 $\Delta t \lesssim 2.785 \tau$，这意味着对于固定的流动参数，时间步长缩放为 $\Delta t = \mathcal{O}(\text{Re}^{-1})$。这是非分裂 BGK 系统的必要约束。

数值结果与验证
该方法通过以下基准测试进行了验证：

• 自由流测试： 在弯曲网格（圆柱、椭圆、NACA0012）上，该方法将均匀自由流保持至机器精度（$\sim 10^{-15}$），而传统的阶梯状方法在壁面法线处表现出 $\mathcal{O}(1)$ 的几何误差。
• Taylor–Green 涡流： 2D 和 3D 的收敛性研究证实了平流算子的设计阶数（2、4 和 6）。该方法成功模拟了 3D Taylor–Green 涡流向湍流的转变。
• 驱动顶盖流（Re=100）： 该方法以极高的精度重现了基准速度剖面（Ghia 等人），其精度仅受限于非热模型中的 $\mathcal{O}(\text{Ma}^2)$ 可压缩性底限。
• 绕流圆柱（Re=40 与 Re=100）：
  • 定常流 (Re=40)： 再循环长度和阻力系数与参考数据分别在 $\sim 9\%$ 和 $6\text{--}12\%$ 以内匹配。作者将微小的阻力高估归因于壁面应力评估和有限定义域效应，而非内部精度问题。
  • 非定常流 (Re=100)： Strouhal 数的预测值与参考值偏差在 $\sim 2.5\%$ 以内。尾迹在 Re=1000 以下保持稳定。
• NACA0012 翼型（Re=500, $\alpha=10^\circ$）： 该方法重现了参考表面压力和皮肤摩擦分布（Hafez 等人），误差在百分之几以内。升力和阻力系数均已实现网格收敛。

意义与主张
论文明确界定了 IGA-LBM 的意义及其特定限制：

1. 几何保真度，而非内部优越性： 作者强调，等几何路径的主要优势在于精确的几何和精确的度量，而非更优越的内部算子。在光滑网格上，内部算子等效于标准的四阶紧致有限差分。其收益在于消除了阶梯状误差和度量差分误差，从而允许在弯曲边界上保持形式上的精度阶数。
2. 自由流保持： 虽然精确的自由流保持是一个关键特性，但论文指出这是平流形式在任何光滑网格上的属性，并不需要离散几何守恒律（GCL）。这是一个鲁棒性特征，而非等几何特有的唯一属性。
3. 权衡： 该方法本质上不是保守型的（它是一种强形式点值方案），并且需要一个可调的滤波器参数以保证稳定性。显式时间步长限制随 $\text{Re}^{-1}$ 缩放，这对于高雷诺数模拟可能在计算上较为昂贵。
4. 适用范围： 该方法被呈现为一种基于精确解析几何的高阶、保持自由流的 LBM，用于处理弯曲边界流动。它并不声称在无需进一步发展（例如，开发保守型弱形式离散化）的情况下，就能解决高雷诺数非定常分离流的局限性。

总之，本文提出了一种经过严格分析的高阶 LBM，成功地将等几何分析的优势——特别是精确几何和度量评估——转移到了动力学求解器中，解决了长期以来基于格子的流体模拟中存在的几何离散化误差问题。