[论文解读] A WENO-type slope-limiter for a family of piecewise polynomial methods

Darren Engwirda, Maxwell Kelley|arXiv (Cornell University)|Jun 27, 2016

Numerical methods for differential equations参考文献 21被引用 12

一句话总结

本文提出了一种适用于分段抛物线法（PPM）和分段四次法（PQM）的WENO型斜率限制器，通过非线性组合自然重建与单调限制重建，实现高阶精度的同时抑制激波附近的振荡。该方法使PPM和PQM格式在光滑问题上分别实现三阶和五阶收敛，相较于传统单调限制器显著降低了数值耗散。

ABSTRACT

A new, high-order slope-limiting procedure for the Piecewise Parabolic Method (PPM) and the Piecewise Quartic Method (PQM) is described. Following a Weighted Essentially Non-Oscillatory (WENO)-type paradigm, the proposed slope-limiter seeks to reconstruct smooth, non-oscillatory piecewise polynomial profiles as a non-linear combination of the natural and monotone-limited PPM and PQM interpolants. Compared to existing monotone slope-limiting techniques, this new strategy is designed to improve accuracy at smooth extrema, while controlling spurious oscillations in the neighbourhood of sharp features. Using the new slope-limited PPM and PQM interpolants, a high-order accurate Arbitrary-Lagrangian-Eulerian framework for advection-dominated flows is constructed, and its effectiveness is examined using a series of one- and two-dimensional benchmark cases. It is shown that the new WENO-type slope-limiting techniques offer a significant improvement in accuracy compared to existing strategies, allowing the PPM- and PQM- based schemes to achieve fully third- and fifth-order accurate convergence, respectively, for sufficiently smooth problems.

研究动机与目标

• 降低大气和海洋模拟中高阶有限体积格式的数值耗散。
• 解决单调斜率限制器在光滑极值点处导致的精度退化问题，该问题会人为压平局部极值。
• 开发一种高阶、非振荡的重构策略，能够在保持单调性的同时保留局部极值。
• 将类似WENO的非线性加权方法拓展至分段多项式方法（PPM和PQM），以在ALE和半拉格朗日框架中提升精度。
• 通过最小化非等熵混合的虚假贡献，实现PPM和PQM的完全三阶和五阶收敛。

提出的方法

• 提出一种WENO型非线性重构方法，通过非线性权重组合自然PPM/PQM插值与单调限制插值。
• 利用目标拐点位置控制重构多项式的形状，确保单调性并避免振荡。
• 通过多项式系数上的二次约束，推导出修正的边坡斜率与边界值估计，以实现一致且单调的轮廓。
• 采用两步限制过程：首先通过系数约束将拐点定位在单元边界；其次通过调整相邻值与斜率来修正不一致的边坡斜率。
• 采用非线性加权策略，在光滑区域优先选择自然PPM/PQM重构，在不连续区域则倾向使用单调限制版本。
• 将新型斜率限制器集成至任意拉格朗日-欧拉（ALE）框架中，用于对流主导流动，实现时间和空间上的高阶精度。

实验结果

研究问题

• RQ1 WENO型非线性重构策略是否能在PPM和PQM中同时提升光滑极值点处的精度，同时在激波附近保持非振荡特性？
• RQ2 所提出的斜率限制器是否能使PPM和PQM在光滑问题中实现其理论上的三阶和五阶收敛速率？
• RQ3 在长时间对流模拟中，新限制器与传统单调限制器相比，在数值耗散和非等熵混合方面表现如何？
• RQ4 类似WENO的非线性加权方法是否能有效适配至分段多项式方法，而不仅限于点值重构？
• RQ5 为确保重构多项式轮廓的单调性与一致性，对边坡斜率与边界值需进行何种修正？

主要发现

• 所提出的WENO型斜率限制器使基于PPM的格式在足够光滑的问题上实现三阶收敛，消除了传统单调限制器中常见的精度退化现象。
• 在相同条件下，基于PQM的格式实现了五阶收敛，表明新限制器能保持其完整的理论精度阶次。
• 与传统单调限制器相比，该方法显著降低了数值耗散，尤其在光滑极值点处，通过避免对局部极值的人为压平实现。
• 重构轮廓保持单调且非振荡，拐点被约束在单元边界，确保了稳定性和一致性。
• 修正后的边坡斜率与边界值校正消除了斜率估计中的不一致性，保证了最终插值函数的精确单调性。
• 一维和二维基准测试验证了新限制器在对流主导流动中更高的精度与更强的鲁棒性。