Multiresolution analysis as a criterion for effective dynamic mesh adaptation - A case study for Euler equations in the SAMR framework AMROC

Autor(es)

Deiterding, Ralf ^[1] ; Domingues, Margarete Oliveira ^[2] ; Schneider, Kai ^{[3, 4]}

Afiliação do(s) autor(es)

^[1] Univ Southampton, Aerodynam & Flight Mech Res Grp, Boldrewood Innovat Campus, Southampton SO16 7QF, Hants - England

^[2] Natl Inst Space Res INPE, Associate Lab Comp & Appl Math, Av Astronautas 1758, BR-12227010 Sao Jose Dos Campos, SP - Brazil

^[3] Aix Marseille Univ, Inst Math Marseille I2M, CNRS, 39 Rue Joliot Curie, F-13453 Marseille 13 - France

^[4] Cent Marseille, 39 Rue Joliot Curie, F-13453 Marseille 13 - France

• Tipo de documento: Artigo Científico
• Fonte: COMPUTERS & FLUIDS; v. 205, JUN 15 2020.
• Citações Web of Science: 0

Resumo

Dynamic mesh adaptation methods require suitable refinement indicators. In the absence of a comprehensive error estimation theory, adaptive mesh refinement (AMR) for nonlinear hyperbolic conservation laws, e.g. compressible Euler equations, in practice utilizes mainly heuristic smoothness indicators like combinations of scaled gradient criteria. As an alternative, we describe in detail an easy to implement and computationally inexpensive criterion built on a two-level wavelet transform that applies projection and prediction operators from multiresolution analysis. The core idea is the use of the amplitude of the wavelet coefficients as smoothness indicator, as it can be related to the local regularity of the solution. Implemented within the fully parallelized and structured adaptive mesh refinement (SAMR) software system AMROC (Adaptive Mesh Refinement in Object-oriented C++), the proposed criterion is tested and comprehensively compared to results obtained by applying the scaled gradient approach. A rigorous quantification technique in terms of numerical adaptation error versus used finite volume cells is developed and applied to study typical two- and three-dimensional problems from compressible gas dynamics. It is found that the proposed multiresolution approach is considerably more efficient and also identifies besides discontinuous shock and contact waves - in particular smooth rarefaction waves and their interaction as well as small-scale disturbances much more reliably. Aside from pathological cases consisting solely of planar shock waves, the majority of realistic cases show reductions in the number of used finite volume cells between 20 to 40%, while the numerical error remains basically unaltered. (C) 2020 Elsevier Ltd. All rights reserved. (AU)

• Processo FAPESP: 16/50016-9 - Advanced criteria for dynamic mesh adaption in computacional space weather forecasting
• Beneficiário: Margarete Oliveira Domingues
• Modalidade de apoio: Auxílio à Pesquisa - Regular
• Processo FAPESP: 15/25624-2 - Desenvolvimento de modelagem multiescala para instabilidades locais não-lineares em Astrofísica e Geofísica Espacial
• Beneficiário: Margarete Oliveira Domingues
• Modalidade de apoio: Auxílio à Pesquisa - Regular