Are adaptive wavelet discretizations dissipative?
Adaptive wavelet schemes for solving partial differential equations offer an attractive possibility to introduce locally refined grids, which dynamically track the evolution of the solution in scale and space. Automatic error control of the adaptive discretization, with respect to a uniform grid solution, is an advantageous feature. For a review in the context of computational fluid dynamics we refer to [4]. In this work we focus on dynamical Galerkin schemes, where the projection operator changes over time. When selecting a subset of basis functions, the projection operator is non-differentiable and an integral formulation has to be used. We will analyze the projected equations with respect to existence and uniqueness of the solution and prove that non-smooth projection operators introduce dissipation, a result which is crucial for adaptive discretizations of PDEs. For the Burgers equation we will illustrate numerically that thresholding the wavelet coefficients, and thus changing the projection space, will indeed introduce dissipation of energy [3]. We discuss consequences for the so-called ``pseudo-adaptive" simulations, where time evolution and de-alising are done in Fourier space whilst thresholding is carried out in wavelet space. Applications will deal with the inviscid Burgers equation using a simple small-scale periodic filtering scheme, as well as CVS filtering [2], the latter has been applied for three-dimensional turbulence in [1].
Joint work with Rodrigo Pereira, Romain Nguyen van yen and Marie Farge.
[1] N. Okamoto, K. Yoshimatsu, K. Schneider, M. Farge and Y. Kaneda. Coherent vortex simulation of three- dimensional decaying homogeneous isotropic turbulence. SIAM Multiscale Model. Simul., 9(3), 1144--1161, 2011.
[2] R. M. Pereira, R. Nguyen van yen, M. Farge and K. Schneider. Wavelet methods to eliminate resonances in the Galerkin-truncated Burgers and Euler equations. Phys. Rev. E, 87, 033017, 2013.
[3] R. M. Pereira, R. Nguyen van yen, K. Schneider and M. Farge. Adaptive solution of initial value problems by a dynamical Galerkin scheme. Preprint 12/2015, to be submitted.
[4] K. Schneider and O. Vasilyev. Wavelet methods in computational fluid dynamics. Annu. Rev. Fluid Mech., 42, 473--503, 2010.