[论文解读] CIP/multi-moment finite volume method with arbitrary order of accuracy

Feng Xiao, Satoshi Satoshi|arXiv (Cornell University)|Jul 30, 2012

Electromagnetic Simulation and Numerical Methods被引用 7

一句话总结

该论文提出了一种高阶CIP/多矩有限体积方法，通过在单元边界引入多个导数矩，实现了任意阶精度。该方法使用体积平均值和点态导数作为矩，对每个单元重构高阶多项式，并通过通量形式有限体积法和局部黎曼问题更新，随着矩阶数的提高，在激波和间断处表现出更优的分辨率。

ABSTRACT

This paper presents a general formulation of the CIP/multi-moment finite volume method (CIP/MM FVM) for arbitrary order of accuracy. Reconstruction up to arbitrary order can be built on single cell by adding extra derivative moments at the cell boundary. The volume integrated average (VIA) is updated via a flux-form finite volume formulation, whereas the point-based derivative moments are computed as local derivative Riemann problems by either direct interpolation or approximate Riemann solvers.

研究动机与目标

• 为双曲守恒律开发一种具有任意阶精度的保守有限体积方法。
• 通过在单元边界引入多个导数矩，克服传统单矩格式（如MUSCL、WENO）的局限性。
• 通过使用体积积分平均值（VIA）和高阶导数矩，在紧凑模板上实现高阶空间重构。
• 通过通量形式有限体积法和局部黎曼问题更新矩，确保计算效率和守恒性。
• 通过提高导数矩的阶数，展示该方法在解析激波和间断方面更高的精度。

提出的方法

• 利用单元内的体积积分平均值（VIA）和单元边界处的最高K阶导数矩，在每个单元内重构高阶多项式。
• 通过2(K+1)+1个约束条件，包括单元中心处的一阶导数，强制多项式及其导数在单元界面处连续。
• 通过标准通量形式有限体积格式，对VIA进行更新，积分单元边界上的通量差。
• 通过守恒律导出的控制方程演化导数矩，其中通量的空间导数通过多矩重构计算。
• 利用通量雅可比矩阵对最高阶导数矩求解线性化导数黎曼问题。
• 应用简单的斜率限制器，在保持高阶精度的同时抑制大梯度区域的数值振荡。

实验结果

研究问题

• RQ1 仅使用单元内矩而不依赖全局模板，有限体积方法能否实现任意阶精度？
• RQ2 如何在保持守恒性的同时，利用单元边界处的多个导数矩构建高阶多项式重构？
• RQ3 提高导数矩的阶数对激波和间断分辨率有何影响？
• RQ4 在非结构化或复杂网格上，通过导数矩的局部黎曼问题求解，该方法能否保持稳定性和精度？
• RQ5 与现有高阶格式（如WENO或ADER）相比，所提出的多矩方法在效率和精度方面表现如何？

主要发现

• 通过在单元边界引入最高K阶导数矩，该方法实现了任意阶精度，可在紧凑模板上实现高阶重构。
• 数值结果表明，提高导数矩的阶数可显著改善对激波和接触间断的分辨率，结合斜率限制后振荡极小。
• 该格式在1D激波管问题中表现出高精度，其锐度和位置精度与五阶WENO格式相当或更优。
• 对最高阶导数矩采用线性化导数黎曼问题，确保了稳定性和一致性，无需在单元间允许导数不连续。
• 该方法计算高效且基于单单元，避免了全局重构或每阶段多时间步的复杂性。
• 该公式可通过分量重构和基于通量雅可比矩阵的导数黎曼求解器，自然推广至守恒律系统（如欧拉方程）。