Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations

The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing numerical dissipations for high speed compressible viscous ows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an arti cial compression method (ACM) of Harten (1978) but utilize it in an entirely di erent context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch o the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992 ) used by Gerritsen and Olsson (1996) in an adaptive mesh re nement method, to determine regions where re nement should be done. The other is the modi cation of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand alone options for numerical algorithm other than the Yee et al. scheme.

[1]  Ingrid Daubechies,et al.  Ten Lectures on Wavelets , 1992 .

[2]  Stéphane Mallat,et al.  Characterization of Signals from Multiscale Edges , 2011, IEEE Trans. Pattern Anal. Mach. Intell..

[3]  Randall J. LeVeque,et al.  A study of numerical methods for hyperbolic conservation laws with stiff source terms , 1990 .

[4]  Ami Harten,et al.  Self adjusting grid methods for one-dimensional hyperbolic conservation laws☆ , 1983 .

[5]  S. Mallat A wavelet tour of signal processing , 1998 .

[6]  H. C. Yee,et al.  A class of high resolution explicit and implicit shock-capturing methods , 1989 .

[7]  Stéphane Mallat,et al.  Singularity detection and processing with wavelets , 1992, IEEE Trans. Inf. Theory.

[8]  M. Farge Wavelet Transforms and their Applications to Turbulence , 1992 .

[9]  A. Harten,et al.  The artificial compression method for computation of shocks and contact discontinuities. 3: Self-adjusting hybrid schemes , 1977 .

[10]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[11]  S. Mallat Multiresolution approximations and wavelet orthonormal bases of L^2(R) , 1989 .

[12]  Chi-Wang Shu Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[13]  Neil D. Sandham,et al.  Progress in the Development of a Class of Efficient Low Dissipative High Order Shock-Capturing Methods , 2000 .

[14]  Gilbert Strang,et al.  Wavelets and Dilation Equations: A Brief Introduction , 1989, SIAM Rev..

[15]  H. C. Yee,et al.  Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence , 2002 .

[16]  H. C. Yee,et al.  Entropy Splitting and Numerical Dissipation , 2000 .

[17]  G. Strang Wavelet transforms versus Fourier transforms , 1993, math/9304214.

[18]  Neil D. Sandham,et al.  Low-Dissipative High-Order Shock-Capturing Methods Using Characteristic-Based Filters , 1999 .

[19]  Stéphane Mallat,et al.  Multifrequency channel decompositions of images and wavelet models , 1989, IEEE Trans. Acoust. Speech Signal Process..

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

[21]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[22]  H. C. Yee Upwind and Symmetric Shock-Capturing Schemes , 1987 .