[论文解读] A cure for instabilities due to advection-dominance in POD solution to advection-diffusion-reaction equations

Mejdi Azaïez, Tomás Chacón Rebollo+1 more|arXiv (Cornell University)|Jul 12, 2019

Model Reduction and Neural Networks参考文献 58被引用 18

一句话总结

本文提出一种稳定后处理（SPP）策略，以解决在强对流主导区域中，对流-扩散-反应方程的 Proper Orthogonal Decomposition 降阶模型（POD-ROM）中存在的数值不稳定性问题。通过在离线快照生成和在线解计算过程中应用 SPP，该方法显著减少了虚假振荡，提升了精度，即使在低模态 POD 基下，也能实现接近离线水平的性能。

ABSTRACT

In this paper, we propose to improve the stabilized POD-ROM introduced by S. Rubino in [37] to deal with the numerical simulation of advection-dominated advection-diffusion-reaction equations. In particular, we introduce a stabilizing post-processing strategy that will be very useful when considering very low diffusion coefficients, i.e. in the strongly advection-dominated regime. This strategy is applied both for the offline phase, to produce the snapshots, and the reduced order method to simulate the new solutions. The new process of a posteriori stabilization is detailed in a general framework and applied to advection-diffusion-reaction problems. Numerical studies are performed to discuss the accuracy and performance of the new method in handling strongly advection-dominated cases.

研究动机与目标

• 解决标准 POD-Galerkin ROM 在对流主导流动中因未解析小尺度模态耗散不足而导致的已知不稳定性问题。
• 克服现有稳定 POD-ROM（如 SD-POD-ROM）在扩散系数极低（即强对流主导区域）时的局限性。
• 开发一种通用的后处理稳定框架，适用于离线快照生成和在线 ROM 解算过程。
• 提升在长时间模拟以及具有陡峭层或间断问题中的降阶模型的精度与鲁棒性。
• 证明后处理可有效抑制噪声模态的影响，且无需额外在线计算成本，仅在标准 ROM 求解基础上增加后处理步骤。

提出的方法

• 提出一种稳定后处理（SPP）技术，应用于离线阶段的快照和在线 ROM 解算过程。
• 将 SPP 作为滤波操作应用，以抑制降阶解中的高频振荡，模拟数值稳定性。
• 将 SPP 整合进标准 POD-ROM 框架中，以稳定有限元格式（SD-POD-ROM）作为基线方法。
• 在离线阶段采用滤波对流稳定方法生成更少噪声的快照，从而降低未解析模态在降阶基中的影响。
• 在在线阶段应用相同的 SPP 技术以稳定降阶解，尤其在陡峭层或前缘附近表现更优。
• 在通用框架下构建该方法，适用于对流-扩散-反应问题，可推广至其他 PDE。

实验结果

研究问题

• RQ1 后处理策略是否能有效稳定在标准 ROM 和稳定 ROM 均失效的强对流主导区域中的 POD-ROM 解？
• RQ2 在离线阶段应用后处理是否能提升快照质量，并减少噪声模态在降阶基中的影响？
• RQ3 在离线和在线阶段均结合 SPP 对降阶解的精度和振荡水平有何影响？
• RQ4 所提方法是否能在保持低计算成本的同时，实现与全阶模型相当的精度？
• RQ5 该方法在长时间模拟及具有移动陡峭前缘的问题中表现如何？

主要发现

• 当 ν = 10⁻⁸ 且 r = 90 时，采用在线 SPP 的 SD-ROM 的 L² 误差为 0.0589，接近离线解的精度。
• SPP 方法显著降低了振荡：在 r = 90 时，采用后处理的 SD-ROM 误差（0.0589）低于标准 G-ROM（0.1224），甚至低于 G-ROM 搭配后处理的版本（0.0884）。
• 即使仅使用 r = 30 个模态（捕捉 99.71% 的动能），SPP 增强的 SD-ROM 也表现出比标准 ROM 更优的振荡抑制与稳定性。
• 该方法能有效定位并抑制陡峭层附近的振荡，同时比标准或 G-ROM 方法更好地保持前缘振幅。
• SPP 策略在不同 POD 模态数量下均表现有效，随着 r 增大，振荡逐渐减小，表明其具有鲁棒性与可扩展性。
• 数值结果表明，SPP 可作为现有稳定 POD-ROM 的高效、低成本增强手段，尤其在对流主导场景中效果显著。