Point Value Multiscale Algorithms for 2D Compressible Flows

The numerical simulation of physical problems modeled by systems of conservation laws is difficult due to the presence of discontinuities in the solution. High-order shock capturing schemes combine sharp numerical profiles at discontinuities with a highly accurate approximation in smooth regions, but usually their computational cost is quite large. Following the idea of A. Harten [Comm. Pure Appl. Math., 48 (1995), pp. 1305--1342] and Bihari and Harten [SIAM J. Sci. Comput., 18 (1997), pp. 315--354], we present in this paper a method to reduce the execution time of such simulations. It is based on a point value multiresolution transform that is used to detect regions with singularities. In these regions, an expensive high-resolution shock capturing scheme is applied to compute the numerical flux at cell interfaces. In smooth regions a cheap polynomial interpolation is used to deduce the value of the numerical divergence from values previously obtained on lower resolution scales. This method is applied to solve the two-dimensional compressible Euler equations for two classical configurations. The results are analyzed in terms of quality and efficiency.

[1]  J. Ibáñez,et al.  Numerical Relativistic Hydrodynamics , 1993 .

[2]  Sidi Mahmoud Kaber,et al.  Finite volume schemes on triangles coupled with multiresolution analysis , 1999 .

[3]  Antonio Marquina,et al.  Morphology and Dynamics of Relativistic Jets , 1997 .

[4]  A. Harten.  Multiresolution representation of data: a general framework , 1996 .

[5]  A. Harti.  Discrete multi-resolution analysis and generalized wavelets , 1993 .

[6]  Wim Sweldens,et al.  The lifting scheme: a construction of second generation wavelets , 1998 .

[7]  Gabi Ben-Dor,et al.  Shock wave reflection phenomena , 1992 .

[8]  S.S.M. Wong,et al.  Relativistic Hydrodynamics and Essentially Non-oscillatory Shock Capturing Schemes , 1995 .

[9]  J. Quirk.  A Contribution to the Great Riemann Solver Debate , 1994 .

[10]  Wolfgang Dahmen,et al.  Multiresolution schemes for conservation laws , 2001, Numerische Mathematik.

[11]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[12]  A. Harten,et al.  Multiresolution Based on Weighted Averages of the Hat Function I: Linear Reconstruction Techniques , 1998 .

[13]  Björn Sjögreen,et al.  Numerical experiments with the multiresolution scheme for the compressible Euler equations , 1995 .

[14]  Antonio Marquina,et al.  Computing Strong Shocks in Ultrarelativistic Flows: A Robust Alternative , 1999 .

[15]  Wolfgang Dahmen,et al.  Multiresolution schemes for conservation laws , 2001 .

[16]  Barna L. Bihari,et al.  Multiresolution Schemes for the Numerical Solution of 2-D Conservation Laws I , 1997, SIAM J. Sci. Comput..

[17]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[18]  Antonio Marquina,et al.  Capturing Shock Reflections , 1996 .

[19]  Ronald Fedkiw,et al.  An Isobaric Fix for the Overheating Problem in Multimaterial Compressible Flows , 1999 .

[20]  Chaowei Hu,et al.  No . 98-32 Weighted Essentially Non-Oscillatory Schemes on Triangular Meshes , 1998 .

[21]  Antonio Marquina,et al.  Local Piecewise Hyperbolic Reconstruction of Numerical Fluxes for Nonlinear Scalar Conservation Laws , 1994, SIAM J. Sci. Comput..