Coherent vorticity extraction in resistive drift-wave turbulence: Comparison of orthogonal wavelets versus proper orthogonal decomposition

We assess two techniques for extracting coherent vortices out of turbulent flows: the wavelet based Coherent Vorticity Extraction (CVE) and the Proper Orthogonal Decomposition (POD). The former decomposes the flow field into an orthogonal wavelet representation and subsequent thresholding of the coefficients allows one to split the flow into organized coherent vortices with non-Gaussian statistics and an incoherent random part which is structureless. POD is based on the singular value decomposition and decomposes the flow into basis functions which are optimal with respect to the retained energy for the ensemble average. Both techniques are applied to direct numerical simulation data of two-dimensional drift-wave turbulence governed by Hasegawa–Wakatani equation, considering two limit cases: the quasi-hydrodynamic and the quasi-adiabatic regimes. The results are compared in terms of compression rate, retained energy, retained enstrophy and retained radial flux, together with the enstrophy spectrum and higher order statistics.

RésuméCet article compare deux méthodes permettant dʼextraire les contributions cohérentes dans les écoulements turbulents : la méthode CVE (Coherent Vorticity Extraction), basée sur la représentation en ondelettes, et la méthode POD (Proper Orthogonal Decomposition). La première méthode, CVE, décompose lʼécoulement en base dʼondelettes orthogonales puis, grâce au filtrage des coefficients dʼondelettes, permet de séparer lʼécoulement entre une contribution cohérente ayant un comportement statistique non Gaussien et un écoulement résiduel aléatoire ne présentant pas de structures. La seconde méthode, POD, est basée sur la décomposition en valeurs singulières. Elle décompose lʼécoulement dans la base de fonctions qui retient le mieux lʼénergie dʼun ensemble de réalisations de lʼécoulement. Ces deux méthodes sont ensuite utilisées pour analyser les résultats de simulations numériques directes calculant un écoulement turbulent bidimensionnel dʼondes de dérive, gouverné par lʼéquation dʼHasegawa–Wakatani, ceci pour deux cas limites : le régime quasi-hydrodynamique et le régime quasi-adiabatique. Les résultats sont comparés en fonction des taux de compression, des proportions dʼenstrophie et de flux radial retenues, ainsi que du spectre dʼenstrophie (statistique dʼordre deux) et de statistiques dʼordre plus élévé.