We employ a wavelet representation in the inviscid Burgers equation to anlyse and treat the resonance phenomenon that appears in the Galerkin-truncated numerical solution to this problem. We use the CVS (Coherent Vorticity Simulation) method to avoid the formation of the resonances, obtaining an appropriate regularized solution. The previously developed method with a complex wavelet frame is applied and expanded to embrace the use of a real orthogonal wavelet basis, which we show to exhibit proper results only under the condition of adding a safety zone in wavelet space. We also apply the complex wavelet based method to the 2D Euler equation problem, showing it is able to filter the resonances in this case either.


This work is done in collaboration with Rodrigo Pereira (Instituto di Fisica, Universidade Federale do Rio de Janeiro, Brazil), Romain Nguyen van yen (Fachbereich Mathematik und Informatik, Freie Universit ̈at Berlin, Germany) and Kai Schneider (CMI, Universit́e d’Aix-Marseille, France). It has been submitted to Physical Review E.
Our papers can be downloaded from http://wavelets.ens.fr

Penalization approaches are nowadays commonly employed to solve boundary or initial-boundary value problems. They consist in embedding the original, possibly complex spatial domain inside a bigger domain having a simpler geometry, for example a torus, while keeping the boundary conditions approximately enforced thanks to new terms that are added to the equations. One particular example is the volume penalization method which, inspired by the physical intuition that a solid wall is similar to a vanishingly porous medium, uses the Brinkman-Darcy drag force as penalization term. The main advantage of such penalized equations is that they can be discretized independently of the geometry of the original problem, since the latter has been encoded into the penalization terms. Such a simplifcation permits a massive reduction in solver development time, since it avoids the issues associated to the design and management of the grid, allowing for example the use of simple spectral solvers in Cartesian geometries. The gain becomes even more substantial when the geometry is time-dependent, as in the case of moving obstacles, or when fluid-structure interaction is taken into account. We present results of a detailed study of the spectral properties of Laplace and Stokes operators, modifed with a volume penalization term designed to approximate Dirichlet conditions in the limit when the penalization parameter, η, tends to zero. The eigenvalues and eigenfunctions are determined either analytically or numerically as functions of η, both in the continuous case and after applying Fourier or finite difference discretization schemes. For fixed η, we find that only the part of the spectrum corresponding to eigenvalues λ ≤ η−1 approaches Dirichlet boundary conditions, while the remainder of the spectrum is made of uncontrolled, spurious wall modes. The penalization error for the controlled eigenfunctions is estimated as a function of λ and η. Surprisingly, in the Stokes case, we show that the eigenfunctions approximately satisfy, with a precision O(η), Navier slip boundary conditions with slip length equal to √η. Moreover, for a given discretization, we show that there exists a value of η, corresponding to a balance between penalization and discretization errors, below which no further gain in precision is achieved. These results shed light on the behavior of volume penalization schemes when solving the Navier-Stokes equations, outline the limitations of the method, and gives indications on how to choose the penalization parameter in practical cases. Possible extensions how to deal with Neumann boundary conditions will also be presented. Finally, different illustrations will be given for vortex-dipole wall interactions, flapping insect wings, fluid-structure interaction and a dynamical mixer.

This work is joint work with Romain Nguyen van yen and Dmitry Kolomenskiy.
Ref: arXiv:1206.0002