Update-Magnitude State Redistribution (UM-SRD): A Shut-off Extension of Weighted SRD for Cut-Cell Methods

技术摘要：更新幅度状态重分布（UM-SRD）

问题陈述
切割网格方法允许有限体积格式在不使用贴体网格的情况下处理复杂几何形状，但由于任意小的切割网格，它们引入了严重的稳定性挑战。状态重分布（SRD）方法，特别是 Berger & Giuliani (2024) 开发的加权 SRD 框架，通过将小网格合并到邻域中、在合并网格上演化平均值，并将结果重新分布回原始网格，解决了这些稳定性问题。

然而，现有加权 SRD 公式的一个关键局限性在于缺乏“关闭”机制。与通量重分布方法不同（后者在通量不平衡为零时后处理自然消失），加权 SRD 在每个时间步都应用非平凡的重分布算子，无论有限体积更新的幅度如何。这种持续激活导致了两个实际缺陷：

1. 无法保持稳态：即使通量平衡恰好为零，SRD 也会用体积加权平均值替换逐点值，阻止收敛到机器精度，并在静态区域引入 $O(1)$ 误差。
2. 不必要的耗散：在具有微小但非零更新的平滑区域中，连续应用重分布会引入与邻域梯度成正比的数值扩散，这与实际物理动力学无关。

方法论
本文介绍了更新幅度状态重分布（UM-SRD），这是一种加权 SRD 的关闭扩展，旨在根据有限体积更新的局部幅度将 SRD 算子与恒等算子进行混合。

该方法在标准的三阶段时间步进框架内运行：

1. 有限体积更新：使用基础格式计算重分布前的状态 $U^*$。
2. 邻域聚合：在合并的邻域 $M_j$ 上计算加权邻域平均值 $\hat{Q}^*_j$。
3. 混合重分布：UM-SRD 不直接应用 SRD 算子 $S_j$，而是应用凸组合：
   $$R_j(s_j) = (1 - s_j)I + s_j S_j$$
   其中 $I$ 是恒等算子，$s_j \in [0, 1]$ 是混合参数。

混合参数 $s_j$ 由平滑的、局部定义的关闭函数 $h(\eta_j)$ 确定：
$$s_j = \frac{\eta_j^p}{\eta_j^p + \tau^p}$$
这里，$\eta_j$ 是邻域 $M_j$ 内最大绝对更新幅度的归一化指示器（具体为 $\eta_j = \frac{\max \|\Delta U_i\|}{\epsilon + \max \|\Delta U_k\|}$）。

• 零更新机制：如果邻域中的所有更新都处于机器精度水平（$\Delta U_i \approx 0$），则 $\eta_j = 0$，$s_j = 0$，该方法精确退化为恒等变换（$U^{n+1} = U^*$）。
• 激活机制：如果发生显著更新（例如在间断附近或当小网格刚度激活时），$\eta_j \approx 1$，$s_j \approx 1$，该方法表现为完全加权的 SRD。

关键分析贡献
本文对具有单个小切割网格的一维线性平流模型中的 UM-SRD 提供了严格分析：

• 守恒性：通过构造证明，对于任何 $s_j \in [0, 1]$，混合算子都是守恒的。
• 总变差缩小（TVD）：在与基础一阶迎风格式相同的 CFL 条件下，证明了 UM-SRD 在原始切割网格上是 TVD 的。这依赖于恒等算子与 TVD SRD 算子的凸组合。
• 稳态保持：证明了 UM-SRD 精确保持基础格式的稳态。如果有限体积更新为零，则 $s_j=0$，且不发生重分布，从而避免了标准 SRD 固有的虚假漂移。
• 精度：在平滑区域，UM-SRD 被证明是基础格式的扰动。如果混合参数随网格尺寸适当缩放（或通过实践中使用的归一化指示器），则对局部截断误差的额外贡献的阶数高于基础一阶格式。因此，UM-SRD 保留了底层迎风方法的一阶收敛率。

数值结果
在一维和二维模型问题上的数值实验证实了理论预测：

• 收敛性：对于平滑平流，UM-SRD 在 $L_1$ 和 $L_\infty$ 范数下表现出一阶收敛，其误差与基础格式和标准加权 SRD 无法区分（差异 < 0.01%）。
• 关闭行为：在具有紧凑脉冲的测试中，当脉冲不存在时，混合参数 $s_j$ 精确降至零，恢复为基础格式；当脉冲与小网格相互作用时，$s_j$ 升至约 0.99。
• 稳定性：在全网格时间步长下（小网格处的局部 CFL 是稳定性极限的 10 倍），UM-SRD 稳定了基础格式，防止了基础迎风格式中观察到的指数发散。
• 非振荡特性：具有间断数据的测试表明，UM-SRD 保持非振荡，在小网格附近与标准 SRD 的行为相匹配。
• 稳态保持：实验证实，虽然标准 SRD 在静态剖面中引入 $O(1)$ 误差，但 UM-SRD 将初始剖面保持到机器精度。
• 二维扩展：该算法成功应用于具有嵌入边界的二维笛卡尔网格，在不修改混合逻辑的情况下展示了一阶收敛。

意义
本文声称，UM-SRD 为现有的加权 SRD 基础设施提供了一种实用且最小化的修改。通过重用既定的邻域和权重结构，并仅添加局部混合机制，它赋予了 SRD 消失更新属性。这将状态重分布针对任意小切割网格的稳定性优势与通量重分布方法的精度保持特性相结合，具体实现了精确的稳态保持，并减少了平滑区域中不必要的耗散。作者指出，虽然当前的分析仅限于一维标量线性平流，但该算法结构与维度无关，并准备好集成到生产级切割网格代码中。