Ensemble algorithm for parametrized flow problems with energy stable open boundary conditions

技术摘要：具有能量稳定开放边界条件的参数化流问题的集成算法

问题陈述
本文解决了求解受变化的初始条件、体力和粘度系数 ($\nu_j$) 影响的非定常不可压缩纳维-斯托克斯方程集成（Ensemble）计算的挑战。在敏感性分析、降阶建模和集成预报等应用中，集成成员的数量 ($J$) 可能很大。标准的半隐式方法要求在每个时间步内，由于系统矩阵对 $\nu_j$ 和非线性项的依赖性，必须为每个成员 $j$ 组装并求解一个独立的线性系统。这导致了极高的计算成本。虽然全显式方法允许单次系统组装，但它们会施加严格的时间步长限制，尤其是在高雷诺数流动的情况下。以往的集成方法针对恒定粘度或特定的粘度变化（例如 [8]）通常需要对粘度系数进行限制性假设（例如有界的相对差异）以确保稳定性。

方法论
作者提出了一种新型的集成计算方法，该方法将线性系统矩阵与集成成员索引 $j$ 解耦，从而允许在每个时间步求解单个具有多个右端项的线性系统。其核心方法论包括：

1. 分裂策略： 将非线性对流项分解为一个集成平均分量和一个波动分量。平均分量采用半隐式处理，而波动分量采用显式处理。
2. 粘度处理： 与以往要求对粘度系数变化进行特定约束（以确保稳定性）的工作（如 [8]）不同，本文引入了一种不同的扩散项处理方式。通过定义最大粘度 $\nu_\infty = \max_j \nu_j$ 并对扩散项进行分裂，作者推导出了一个使线性系统矩阵与 $j$ 无关的方案，且无需对粘度变化的幅度施加限制。
3. 开放边界条件 (OBC)： 该方法被扩展到具有开放边界 ($\Gamma_N$) 的区域。边界被划分为流出区和回流区。作者利用一种涉及特征长度尺度 $L$ 和类 Heaviside 函数的方法来处理流出的对流特性，以实现能量稳定的边界条件。
4. 离散化： 本文提出了阶一阶（算法 1–3）和二阶（算法 4–6）时间步进方案。空间离散使用保形有限元空间（具体为 Taylor-Hood $P_2/P_1$ 对）结合偏斜对称三线性形式以确保能量稳定性。
5. 稳定性分析： 理论证明确立了在一阶方案下满足特定 CFL 条件（涉及速度波动）时的能量稳定性。对于二阶方案，虽然没有提供不带粘度限制的严格稳定性证明（由于 $\nabla \partial_{tt} u$ 项的复杂性），但数值实验表明其在与一阶方案相同的条件下具有稳定性。

主要贡献

• 无限制粘度稳定性： 主要的理论贡献是提出了一种针对可变粘度系数的集成方案，该方案在无需满足先前文献中存在的限制性假设 $|\nu_j - \nu|/\nu \leq \sqrt{\mu}$ 的情况下保证了稳定性。
• 开放边界扩展： 作者将集成方法扩展到具有开放边界的问题，推导了弱形式，并在特定的时间步长约束下证明了一阶方案的能量稳定性。
• 算法效率： 所提算法导致每个时间步仅有一个系数矩阵，可通过多右端项求解器高效求解，与求解 $J$ 个独立系统相比，显著降低了计算成本。

结果
使用 FreeFem++ 软件包进行了数值实验：

1. 收敛性研究： 在具有精确解的区域进行的收敛性测试证实，所提二阶方案在速度和压力上均达到了预期的二阶收敛率。其精度与独立的顺序模拟相当。
2. 绕圆柱流动： 算法在涉及三种不同粘度值的绕圆柱基准问题上进行了测试。得到的计算结果（阻力、升力和压力降）与使用 Dirichlet 和“无操作”（do-nothing）流出条件的文献参考值相符。开放边界条件成功地允许涡旋自然退出区域，避免了与预设抛物线流出剖面相关的人工反射问题。
3. 带收缩的通道流： 模拟了一个具有两个出口且粘度变化的复杂流动场景。尽管集成模拟是在经过自适应非线性滤波器稳定后的较粗网格上进行的，但集成结果在定性上与独立的直接数值模拟 (DNS) 运行结果相匹配。

意义与主张
本文声称提供了一个稳健且高效的参数化流模拟框架。通过消除先前集成方法中对粘度变化所需的限制性假设，所提算法扩大了集成模拟在物理参数高度变化问题中的适用范围。将集成方法扩展到具有可证明能量稳定性的开放边界条件，是向实现无界流真实模拟迈出的重要一步。作者总结道，虽然集成方案在定性和定量上均表现良好（相对于独立模拟），但仍需进一步研究以充分量化计算优势，特别是针对高雷诺数流动在自适应加密网格上的 CPU 时间对比。