Numerical symmetry-preserving techniques for low-dissipation shock-capturing schemes

被引:57
作者
Fleischmann, Nico [1 ]
Adami, Stefan [1 ]
Adams, Nikolaus A. [1 ]
机构
[1] Tech Univ Munich, Chair Aerodynam & Fluid Mech, Dept Mech Engn, D-85748 Garching, Germany
基金
欧洲研究理事会;
关键词
Symmetry breaking; Floating-point arithmetic; Low-dissipation schemes; High-resolution schemes; WENO; Shock-capturing; ESSENTIALLY NONOSCILLATORY SCHEMES; TARGETED ENO SCHEMES; HIGH-ORDER; WENO SCHEMES; RESOLUTION; EULER; IMPLEMENTATION; SOLVER; FLOWS; CODE;
D O I
10.1016/j.compfluid.2019.04.004
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Modern applications of computational fluid dynamics involve complex interactions across scales such as shock interactions with turbulent structures and multiphase interfaces. Such phenomena, which occur at very small physical viscosity, require high-resolution and low-dissipation compressible flow solvers. Many recent publications have focused on the design of high-order accurate numerical schemes and provide e.g. weighted essentially non-oscillatory (WENO) stencils up to 17th order for this purpose. As shown in detail by different authors, such schemes tremendously decrease adverse effects of numerical dissipation. However, such schemes are prone to numerically induced symmetry breaking which renders validation for the targeted problem range problematic. In this paper, we show that symmetry-breaking relates to vanishing numerical viscosity and is driven systematically by algorithmic floating-point effects which are no longer hidden by numerical dissipation. We propose a systematic procedure to deal with such errors by numerical and algorithmic formulations which respect floating-point arithmetic. We show that by these procedures inherent symmetries are preserved for a broad range of test cases with high-order shock-capturing schemes in particular in the high-resolution limit for both 2D and 3D. (C) 2019 The Authors. Published by Elsevier Ltd.
引用
收藏
页码:94 / 107
页数:14
相关论文
共 45 条
[21]   Computation and flow visualization in high-speed aerodynamics [J].
Hadjadj, A ;
Kudryavtsev, A .
JOURNAL OF TURBULENCE, 2005, 6 (16)
[22]   ON UPSTREAM DIFFERENCING AND GODUNOV-TYPE SCHEMES FOR HYPERBOLIC CONSERVATION-LAWS [J].
HARTEN, A ;
LAX, PD ;
VAN LEER, B .
SIAM REVIEW, 1983, 25 (01) :35-61
[23]   Uniformly high order accurate essentially non-oscillatory schemes .3. (Reprinted from Journal of Computational Physics, vol 71, pg 231, 1987) [J].
Harten, A ;
Engquist, B ;
Osher, S ;
Chakravarthy, SR .
JOURNAL OF COMPUTATIONAL PHYSICS, 1997, 131 (01) :3-47
[24]   Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points [J].
Henrick, AK ;
Aslam, TD ;
Powers, JM .
JOURNAL OF COMPUTATIONAL PHYSICS, 2005, 207 (02) :542-567
[25]   An adaptive central-upwind weighted essentially non-oscillatory scheme [J].
Hu, X. Y. ;
Wang, Q. ;
Adams, N. A. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2010, 229 (23) :8952-8965
[26]   Positivity-preserving method for high-order conservative schemes solving compressible Euler equations [J].
Hu, Xiangyu Y. ;
Adams, Nikolaus A. ;
Shu, Chi-Wang .
JOURNAL OF COMPUTATIONAL PHYSICS, 2013, 242 :169-180
[27]   Efficient implementation of weighted ENO schemes [J].
Jiang, GS ;
Shu, CW .
JOURNAL OF COMPUTATIONAL PHYSICS, 1996, 126 (01) :202-228
[28]   Implementation of WENO schemes in compressible multicomponent flow problems [J].
Johnsen, Eric ;
Colonius, Tim .
JOURNAL OF COMPUTATIONAL PHYSICS, 2006, 219 (02) :715-732
[29]   An in-cell reconstruction finite volume method for flows of compressible immiscible fluids [J].
Kannan, Karthik ;
Kedelty, Dominic ;
Herrmann, Marcus .
JOURNAL OF COMPUTATIONAL PHYSICS, 2018, 373 :784-810
[30]   Comparison of several difference schemes on 1D and 2D test problems for the Euler equations [J].
Liska, R ;
Wendroff, B .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2003, 25 (03) :995-1017