期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 387, 期 -, 页码 245-279出版社
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2018.10.050
关键词
Discontinuous Galerkin schemes; A posteriorisubcell limitation; Arbitrary high-order; DG subcell FV formulation; Flux reconstruction (CPR); Hyperbolic conservation laws
资金
- [ANR-17-CE23-0019 Fast4HHO]
In this paper, we present a new limiter for discontinuous Galerkin (DG) schemes, based on subcell resolution through reconstructed flux correction, for hyperbolic conservation laws. This limitation technique is constructed by means of a subcell Finite Volume (FV) formulation. The algorithm is thus simple, and is able to retain the very high accurate subcell resolution of DG schemes. The main idea of this correction procedure is to preserve as much as possible the high accuracy and the very precise subcell resolution of DG schemes. Consequently, an a posteriori correction will only be applied locally at the subcell scale where it is needed, but still ensuring the scheme conservation. To do so, we first reformulate DG schemes as a subcell FV scheme provided the use of the correct numerical flux. This very simple development provides us with the so-called DG reconstructed flux. This theoretical part will serve as a basis for our limiter framework. Practically, at each time step, we compute a DG candidate solution and check if this solution is admissible (for instance positive, non-oscillating, ...). If it is the case, we go further in time. Otherwise, we return to the previous time step and correct locally, at the subcell scale, the numerical solution. This is why it is refereed to as a posteriori limitation. To this end, each cell is subdivided into subcells. Then, if the solution is locally detected as bad, we substitute the DG reconstructed flux on the subcell boundaries by a robust first-order or second-order TVD numerical flux. And for subcell detected as admissible, we keep the high-order reconstructed flux which allows us to retain the very high accurate resolution and conservation of the DG scheme. Furthermore, only the solution inside troubled subcells and its first neighbors will have to be recomputed, elsewhere the solution remains unchanged. Numerical results on various type problems and test cases will be presented, both in 1D and 2D on Cartesian grids, to assess the very good performance of the design limiting algorithm. (C) 2018 The Author(s). Published by Elsevier Inc.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据