Hostname: page-component-cd9895bd7-7cvxr Total loading time: 0 Render date: 2024-12-25T02:51:23.941Z Has data issue: false hasContentIssue false

Wavelet tools to study intermittency: application to vortex bursting

Published online by Cambridge University Press:  25 September 2009

JORI RUPPERT-FELSOT
Affiliation:
LMD–IPSL–CNRS, Ecole Normale Supérieure, 24 rue Lhomond, 75231 Paris Cedex 05, France PMMH, ESPCI, 10 rue Vauquelin, 75231 Paris Cedex 05, France
MARIE FARGE*
Affiliation:
LMD–IPSL–CNRS, Ecole Normale Supérieure, 24 rue Lhomond, 75231 Paris Cedex 05, France
PHILIPPE PETITJEANS
Affiliation:
PMMH, ESPCI, 10 rue Vauquelin, 75231 Paris Cedex 05, France
*
Email address for correspondence: [email protected]

Abstract

This paper proposes statistical tools adapted to study highly unsteady and inhomogeneous flows, such as vortex bursting. For this, we use the wavelet representation in which each coefficient keeps track of both location and scale, in contrast to Fourier representation which requires keeping the phase of all coefficients to preserve the spatial structure of the flow. Based on the continuous wavelet transform, we propose several diagnostics, such as the local spectrum and the local intermittency measure. We also use the orthogonal wavelet transform to split each flow realization into coherent and incoherent contributions, which are then analysed independently and from which we define the coherency measure. We apply these wavelet tools to analyse the bursting of a three-dimensional stretched vortex immersed in a steady laminar channel flow. The time evolution of the velocity field is measured by particle image velocimetry during several successive bursts.

Type
Papers
Copyright
Copyright © Cambridge University Press 2009

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

REFERENCES

Antonia, R. A. 1981 Conditional sampling in turbulence measurement. Annu. Rev. Fluid Mech. 13, 131156.Google Scholar
Azzalini, A., Farge, M. & Schneider, K. 2005 Nonlinear wavelet thresholding: a recursive method to determine the optimal denoising threshold. Appl. Comput. Harmon. Anal. 18, 177185.Google Scholar
Berkooz, G., Holmes, P. & Lumley, J. L. 1993 The proper orthogonal decomposition in the analysis of turbulent flows. Annu. Rev. Fluid Mech. 25, 539575.Google Scholar
Blackwelder, R. F. 1977 On the role of phase information in conditional sampling. Phys. Fluids 20, S232.Google Scholar
Cantwell, B. J. 1981 Organized motion in turbulent flow. Annu. Rev. Fluid Mech. 13, 457515.CrossRefGoogle Scholar
Cuypers, Y., Maurel, S. & Petitjeans, P. 2003 Vortex burst as a source of turbulence. Phys. Rev. Lett. 91 (19), 194502.CrossRefGoogle ScholarPubMed
Cuypers, Y., Maurel, A. & Petitjeans, P. 2004 Comparison between an experimental turbulent vortex and the Lundgren vortex. J. Turbulence 5, 30.CrossRefGoogle Scholar
Cuypers, Y., Maurel, A. & Petitjeans, P. 2006 Characterization of an experimental turbulent vortex in the physical and spectral spaces. J. Turbulence 7, 7.CrossRefGoogle Scholar
Donoho, D. & Johnstone, I. 1994 Ideal spatial adaptation via wavelet shrinkage. Biometrika 81, 425455.CrossRefGoogle Scholar
Farge, M. 1992 Wavelet transforms and their applications to turbulence. Annu. Rev. Fluid Mech. 24, 395457.CrossRefGoogle Scholar
Farge, M., Goirand, E., Meyer, Y., Pascal, F. & Wickerhauser, M. V. 1992 Improved predictability of two-dimensional turbulent flows using wavelet packet compression. Fluid Dyn. Res. 10, 229250.Google Scholar
Farge, M., Guezennec, Y. Ho, C. M. & Meneveau, C. 1990 Continuous wavelet analysis of coherent structures. In Stanford University, Studying Turbulence Using Numerical Simulation Databases. 3: Proceedings of the 1990 Summer Program (pp. 331–348).Google Scholar
Farge, M. & Rabreau, G. 1988 Transformée en ondelettes pour détecter et analyser les structures cohérentes dans les écoulements turbulents bidimensionnels. C. R. Acad. Sci. Paris 307 (série II), 14791486.Google Scholar
Farge, M., Pellegrino, G. & Schneider, K. 2001 Coherent Vortex Extraction in three-dimensional Turbulent Flows Using Orthogonal Wavelets. Phys. Rev. Lett. 87 (5), 054501.CrossRefGoogle ScholarPubMed
Farge, M., Pellegrino, G. & Schneider, K. 2003 Coherent Vortex Extraction in three-dimensional homogeneous turbulence: comparison between CVS-filtered and POF-Fourier decompositions. Phys. Fluids 15 (10), 28862896.CrossRefGoogle Scholar
Farge, M., Schneider, K. & Kevlahan, N. 1999 Non-Gaussianity and coherent vortex simulation for two-dimensional turbulence using an adaptive orthogonal wavelet basis. Phys. Fluids 11, 21872201.CrossRefGoogle Scholar
Fiedler, H. E. 1988 Coherent structures in turbulent flows. Prog. Aerosp. Sci. 25, 231269.CrossRefGoogle Scholar
Hussain, A. K. M. F. 1986 Coherent structures and turbulence. J. Fluid Mech. 173, 303356.CrossRefGoogle Scholar
Joyce, G. & Montgomery, D. 1973 Negative temperature states for the two-dimensional guiding-centre plasma. J. Plasma Phys. 10, 107.CrossRefGoogle Scholar
Kolmogorov, A. N. 1941 The local structure of turbulence in an incompressible fluid at very high Reynolds numbers. Dokl. AN SSSR 30, 299303.Google Scholar
Lundgren, T. S. 1982 Strained spiral vortex model for turbulent fine structure. Phys. Fluids 25, 21932203.CrossRefGoogle Scholar
Mallat, S. G. 1989 Multiresolution approximation and wavelet orthogonal bases of L 2. Trans. Am. Math. Soc. 315.CrossRefGoogle Scholar
Meneveau, C. 1991 Analysis of turbulence in the orthonormal wavelet representation. J. Fluid Mech. 232, 469520.Google Scholar
Robinson, S. K. 1991 Coherent motions in the turbulent boundary layer. Annu. Rev. Fluid Mech. 23, 601639.CrossRefGoogle Scholar
Ruppert-Felsot, J. E., Praud, O., Sharon, E. & Swinney, H. L. 2005 Extraction of coherent structures in a rotating turbulent flow experiment. Phys. Rev. E 72 (1), 016311.CrossRefGoogle Scholar

Ruppert-Felsot et al. supplementary movie

Movie 1a. Water channel in which a stretched vortex is visualized by injecting dyes of 4 different colors. The vortex evolution is recorded in real time at 25 images per second. The vortex is initiated by a step at the bottom of the test section and stretched by pumping water out throught two oulets located on each lateral wall. For this flow realization the flow rates have been chosen such as the vortex is permanent and does not burst.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 5.2 MB

Ruppert-Felsot et al. supplementary movie

Movie 1b. Time evolution during three successive vortex life cycles of a stretched vortex immersed in a steady laminar channel flow visualized by injecting dyes of 4 different colors and recorded in real time at 25 frames per second from above the experimental setting. For this flow realization the flow rates have been chosen such as the vortex becomes unstable and bursts.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 4.9 MB

Ruppert-Felsot et al. supplementary movie

Movie 2a. Time evolution during three successive vortex life cycles of the velocity modulus and streamlines. These diagnostics correspond to the velocity components in the observation plane measured by particle image velocimetry (PIV).

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 4.9 MB

Ruppert-Felsot et al. supplementary movie

Movie 2b. Time evolution during three successive vortex life cycles of the vorticity and streamlines. These diagnostics correspond to the velocities in the observation plane measured by particle image velocimetry (PIV) and to the vorticity component perpendicular to this plane.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 5 MB

Ruppert-Felsot et al. supplementary movie

Movie 2c. Time evolution during three successive vortex life cycles of the velocity modulus (left), of the vorticity modulus (right), of the time trace of the velocity modulus (top) measured in the location indicated by a cross on the velocity modulus inset, together with the energy spectrum computed either in the whole domain (solid line) or in a square subdomain (dashed line). During each vortex life cycle the slope of the energy spectrum varies between -2 (thick dashed line) and -1 (thin dashed line), while the slope of its time average is -5/3 (dotted line). These diagnostics correspond to the velocities in the observation plane measured by PIV and to the vorticity component perpendicular to this plane.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 4.9 MB

Ruppert-Felsot et al. supplementary movie

Movie 3a. Time evolution during three successive vortex life cycles of the vorticity modulus. In order to compare with the coherent vorticity shown on movie 3b, the color scale is adapted to the extrema of the coherent vorticity whose values are given in figure 3b. This diagnostic corresponds to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 1.6 MB

Ruppert-Felsot et al. supplementary movie

Movie 3b. Time evolution during three successive vortex life cycles of the coherent vorticity modulus, obtained by orthogonal wavelet decomposition. The color scale is adapted to the extrema of the coherent vorticity, whose values are given in figure 3b. This diagnostic corresponds to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 1.6 MB

Ruppert-Felsot et al. supplementary movie

Movie 3c. Time evolution during three successive vortex life cycles of the incoherent vorticity modulus obtained by orthogonal wavelet decomposition. The color scale is adapted to the extrema of the incoherent vorticity, whose values are given in figure 3c. This diagnostic corresponds to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 3.9 MB

Ruppert-Felsot et al. supplementary movie

Movie 4a. Spatial distribution during three successive vortex life cycles of the variance (top right), skewness (bottom left) and flatness (bottom right) of the total (green), coherent (red) and incoherent (blue) vorticities. The star moving on each curve indicates the same instant as the time evolution of the total, coherent and incoherent vorticity (top left). These diagnostics correspond to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 1.5 MB

Ruppert-Felsot et al. supplementary movie

Movie 4b. Time evolution during three successive vortex life cycles of the percent of wavelet coefficients kept to describe the coherent vorticity (top right), of the fraction of the total vorticity variance retained by the coherent (red) and incoherent (blue) vorticities (bottom left), of the coherency measure (CM) given by the signal to noise ratio (black) between the coherent and incoherent vorticity variances (bottom right). The star moving on each curve indicates the same instant as the time evolution of the total, coherent and incoherent vorticity (top left). These diagnostics correspond to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 1.5 MB

Ruppert-Felsot et al. supplementary movie

Movie 5a. Time evolution during three successive vortex life cycles of the scatter plot of the vorticity versus the stream function for the total (green), the coherent (red) and the incoherent (blue) flows, together with the corresponding time evolution of the total, coherent and incoherent vorticities. These diagnostics correspond to the vorticity component perpendicular to the observation plane computed from the PIV velocity measurements.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 2.9 MB

Ruppert-Felsot et al. supplementary movie

Movie 5b. Time evolution during three successive vortex life cycles of the probablility density function (PDF) of vorticity for the total (green), the coherent (red) and the incoherent (blue) flows, together with the corresponding time evolution of the total, coherent and incoherent vorticities. These diagnostics correspond to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 2.9 MB

Ruppert-Felsot et al. supplementary movie

Movie_5c. Time evolution during three successive vortex life cycles of the enstrophy spectrum, for the total (green), the coherent (red) and the incoherent (blue) flows, together with the corresponding time evolution of the total, coherent and incoherent vorticities. These diagnostics correspond to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 2.9 MB

Ruppert-Felsot et al. supplementary movie

Movie 6a. Time evolution during one vortex life cycle of the continuous wavelet transform (CWT) modulus of vorticity visualized in three-dimensional wavelet space with a vertical cut at abscissa x=63. The isosurfaces values and the corresponding color scale are given in the caption of figure 6. This diagnostic corresponds to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 558.1 KB

Ruppert-Felsot et al. supplementary movie

Movie 6b. Time evolution during three successive vortex life cycles of a two-dimensional cut, in space (horizontal axis) and scale (vertical axis with the smallest scales upwards), of the continuous wavelet transform (CWT) modulus of vorticity. The location of the cut varies in time to track the maximum of the vorticity modulus (top left). These diagnostics correspond to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 5 MB

Ruppert-Felsot et al. supplementary movie

Movie 7a. Time evolution during one vortex life cycle of the local intermittency measure (LIM) of vorticity visualized in three-dimensional wavelet space with a vertical cut at abscissa x=63. The isosurfaces values and the corresponding color scale are given in the caption of figure 7. This diagnostic corresponds to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 670.8 KB

Ruppert-Felsot et al. supplementary movie

Movie 7b. Time evolution during three successive vortex life cycles of the spatial distribution of the scale averaged local intermittency measure (LIM) of vorticity. This diagnostic corresponds to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 3.6 MB

Ruppert-Felsot et al. supplementary movie

Movie 8a. Time evolution during one vortex life cycle of a stretched vortex immersed in a steady laminar channel flow visualized by fluorescein and recorded at 1000 frames per second in a vertical plane located in the middle of the channel.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 683 KB

Ruppert-Felsot et al. supplementary movie

Movie 8b. Time evolution during one vortex life cycle of the continuous wavelet transform (CWT) modulus of vorticity visualized from a side in three-dimensional wavelet space, with a rainbow color scale where the strong, intermediate and weak values are respectively in read, green and blue. This diagnostic corresponds to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 1.7 MB

Ruppert-Felsot et al. supplementary movie

Movie 8c. Time evolution during one vortex life cycle of the continuous wavelet transform (CWT) modulus of vorticity visualized from the top in three-dimensional wavelet space, with a rainbow color scale where the strong, intermediate and weak values are respectively in read, green and blue. . This diagnostic corresponds to the vorticity component perpendicular to the observation plane computed from the velocity measured by PIV.

Download Ruppert-Felsot et al. supplementary movie(Video)
Video 1.3 MB