• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

Flickr Twitter UniPHY Group iResearch App Facebook

Author Select

FULL-TEXT OPTIONS:

Phys. Fluids 24, 025102 (2012); http://dx.doi.org/10.1063/1.3683556 (14 pages)

Wavelet decomposition of forced turbulence: Applicability of the iterative Donoho-Johnstone threshold

Jesse W. Lord1, Mark P. Rast1, Christopher Mckinlay2, John Clyne2, and Pablo D. Mininni3

1Laboratory for Atmospheric and Space Physics, Department of Astrophysical and Planetary Sciences, University of Colorado, Boulder, Colorado 80309-0391, USA
2Computational and Information Systems Laboratory, National Center for Atmospheric Research (NCAR), Boulder, Colorado 80307-3000, USA
3Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Argentina

View MapView Map

(Received 9 May 2011; accepted 18 January 2012; published online 14 February 2012)

We examine the decomposition of forced Taylor-Green and Arn'old-Beltrami-Childress (ABC) flows into coherent and incoherent components using an orthonormal wavelet decomposition. We ask whether wavelet coefficient thresholding based on the Donoho-Johnstone criterion can extract a coherent vortex signal while leaving behind Gaussian random noise. We find that no threshold yields a strictly Gaussian incoherent component, and that the most Gaussian incoherent flow is found for data compression lower than that achieved with the fully iterated Donoho-Johnstone threshold. Moreover, even at such low compression, the incoherent component shows clear signs of large-scale spatial correlations that are signatures of the forcings used to drive the flows.

© 2012 American Institute of Physics

Article Outline

  1. INTRODUCTION
  2. FORCED TURBULENCE SIMULATIONS
  3. COHERENT VORTEX EXTRACTION
  4. CVE ANALYSIS OF TAYLOR-GREEN AND ABC FLOWS
  5. A CLOSER LOOK AT THE ITERATIVE DONOHO-JOHNSTONE THRESHOLD
    1. Is the incoherent vorticity Gaussianly distributed?
    2. Does the incoherent flow helicity indicate random vector orientations?
    3. Does the incoherent kinetic energy scale as k 2 ?
  6. SPATIAL AND TEMPORAL CORRELATIONS
  7. CONCLUSION

KEYWORDS and PACS

PACS

ARTICLE DATA

PUBLICATION DATA

ISSN

1070-6631 (print)  
1089-7666 (online)

  1. M. Farge, K. Schneider, and N. Kevlahan, “Non-gaussianity and coherent vortex simulation for two-dimensional turbulence using an adaptive orthogonal wavelet basis,” Phys. Fluids 11(8), 2187 (1999)PHFLE6000011000008002187000001.
  2. M. Farge, G. Pellegrino, and K. Schneider, “Coherent vortex extraction in 3D turbulent flows using orthogonal wavelets,” Phys. Rev. Lett. 87(5), 054501 (2001). [MEDLINE]
  3. B. Kadoch, M. Oliveira Domingues, I. Broemstrup, L. Larchevêque, K. Schneider, and M. Farge, “Coherent vorticity extraction in 3D homogeneous isotropic turbulence: Influence of the Reynolds number and geometrical statistics,” Braz. J. Phys. 39(2), 531 (2009).
  4. D. L. Donoho and I. M. Johnstone, “Ideal spatial adaptation by wavelet shrinkage,” Biometrika 81(3), 425 (1994). [ISI]
  5. A. Azzalini, M. Farge, and K. Schneider, “Nonlinear wavelet thresholding: A recursive method to determine the optimal denoising threshold,” Appl. Comput. Harmon. Anal. 18, 177 (2005).
  6. P. D. Mininni, A. Alexakis, and A. Pouquet, “Large-scale flow effects, energy transfer, and self-similarity on turbulence,” Phys. Rev. E 74(1), 016303 (2006). [MEDLINE]
  7. S. A. Orszag, “On the elimination of aliasing in finite-difference schemes by filtering high-wavenumber components,” J. Atmos. Sci. 28(6), 1074 (1971). [ISI]
  8. F. Jacobitz, L. Liechtenstein, K. Schneider, and M. Farge, “On the structure and dynamics of sheared and rotating turbulence: Direct numerical simulation and wavelet-based coherent vortex extraction,” Phys. Fluids 20(4), 045103 (2008)PHFLE6000020000004045103000001.
  9. N. Okamoto, K. Yoshimatsu, K. Schneider, M. Farge, and Y. Kaneda, “Coherent vortices in high resolution direct numerical simulation of homogeneous isotropic turbulence: A wavelet viewpoint,” Phys. Fluids 19(11), 115109 (2007)PHFLE6000019000011115109000001.
  10. N. Okamoto, K. Yoshimatsu, K. Schneider, M. Farge, and Y. Kaneda, “Coherent vorticity simulation of three-dimensional forced homogeneous isotropic turbulence,” Multiscale Model. Simul. 9(3), 1144 (2011)MMSUBT000009000003001144000001.
  11. M. Farge, K. Schneider, G. Pellegrino, A. A. Wray, and R. S. Rogallo, “Coherent vortex extraction in three-dimensional homogeneous turbulence: Comparison between CVS-wavelet and POD-Fourier decompositions,” Phys. Fluids 15(10), 2886 (2003)PHFLE6000015000010002886000001. [ISI]
  12. O. Roussel, K. Schneider, and M. Farge, “Coherent vortex extraction in 3D homogeneous turbulence: Comparison between orthogonal and biorthogonal wavelet decomposition,” J. Turbul. 6, 11 (2005).
  13. M. M. Rogers and P. Moin, “Helicity fluctuations in incompressible turbulent flows,” Phys. Fluids 30(9), 2662 (1987)PFLDAS000030000009002662000001.
  14. R. M. Kerr, “Histograms of helicity and strain in numerical turbulence,” Phys. Rev. Lett. 59(7), 783 (1987). [ISI] [MEDLINE]
  15. A. Tsinober, E. Kit, and T. Dracos, “Experimental investigation of the field of velocity gradients in turbulent flows,” J. Fluid Mech. 242, 169 (1992). [Inspec] [ISI]
  16. K. Schneider, M. Farge, G. Pellegrino, and M. M. Rogers, “Coherent vortex simulation of three-dimensional turbulent mixing layers using orthogonal wavelets,” J. Fluid Mech. 534, 39 (2005).
  17. K. Pearson, “On the criterion that a given system of deviations from the probable in the case of a correlated system of variables is such that it can be reasonably supposed to have arisen from random sampling,” Philos. Mag. 50(302), 157 (1900).
  18. T. W. Anderson and D. A. Darling, “Asymptotic theory of certain `goodness of fit' criteria based on stochastic processes,” Ann. Math. Stat. 23(2), 193 (1952). [ISI]
  19. T. W. Anderson and D. A. Darling, “A test of goodness of fit,” J. Am. Stat. Assoc. 49(268), 765 (1954).

Figures (5) Tables (1)

Figures (click on thumbnails to view enlargements)

FIG.1
Logarithm of the probability density of the velocity P(ux) in (a) and vorticity Px) in (b). The coherent (dashed) and incoherent (dotted-dashed) components of the original (solid) Taylor-Green flow were computed using the Coifman12 wavelet in the CVE analysis (iterative Donoho-Johnstone threshold applied to vorticity magnitude). Dotted curves indicate Gaussian probability density functions (PDFs) with the same mean and variance. Insets show incoherent flow velocity and vorticity PDFs and the Gaussian fits over a narrower range of values. The flow is quite isotropic with nearly identical PDFs for the y and z velocity and vorticity, though the z velocity is more Gaussian with somewhat smaller variance.

FIG.1 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.2
Relative kinetic helicity distribution in the (a) Taylor-Green and (b) ABC flows. Black curves plot values for the original flow (solid) and the coherent (dashed) and incoherent (dotted-dashed) components, the later found using the Coifman12 wavelet in the CVE analysis (iterative Donoho-Johnstone threshold applied to the vorticity vector magnitude). Red and blue curves show the same components obtained using the most Gaussian vorticity and 60% compression thresholds, respectively. Green curves plot values obtained when applying the iterative Donoho-Johnstone threshold to individual vorticity component amplitudes. Insets display incoherent results over a limited range, and dotted curves in inset (a) indicate helicity distribution of randomized vector components (Sec. 5B).

FIG.2 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.3
Kinetic energy spectra of the (a) Taylor-Green and (b) ABC flows. Black curves for the original flow (solid), and coherent (dashed), and incoherent (dotted-dashed) components, the later determined using the Coifman12 wavelet in the CVE analysis (iterative Donoho-Johnstone threshold applied to the vorticity vector magnitude). Red and blue curves show the same components obtained using the most Gaussian vorticity and 60% compression thresholds, respectively. Green curves for values obtained by applying the iterative Donoho-Johnstone threshold to individual vorticity component amplitudes. Note that the wavelet decomposition introduces power above the Orszag frequency kmax into both the coherent and incoherent components (Sec. 3), but that these sum to zero in amplitude. Fiducial k−5/3 and k2 lines shown for reference.

FIG.3 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.4
In (a), the Anderson-Darling test statistic A2 measuring the Gaussianity of the incoherent component of the flow as a function of compression, and in (b) A2 as a function of threshold value ADJ. Filled and empty symbols indicate analysis on the Taylor-Green and ABC solutions, respectively. The incoherent component was extracted using the Haar (circles), Coifman12 (squares), and Coifman30 (triangles) wavelets. Black curves highlight Coifman12 and Haar analysis of the Taylor-Green flow. The incoherent A2 values obtained for thresholds (applied to vorticity vector magnitude of the Coifman12 wavelet decomposition) yielding the most Gaussian incoherent component are shown in red and those yielding 60% compression in blue (60% compression symbols for Taylor-Green and ABC flows directly overlap). The incoherent A2 values obtained using the iterative Donoho-Johnstone threshold applied to the individual vorticity amplitudes are shown by the six disconnected points with highest compression on the right of (a) without crosses. Symbols with superimposed crosses indicate results for incoherent flows determined by application of the single iteration threshold to the vorticity vector magnitude, as done in previous studies.

FIG.4 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

FIG.5
In (a), spatial autocorrelation ρz of the Taylor-Green enstrophy as a function of the spatial lag Δz, and in (b), the temporal correlation ρτ of Taylor-Green vorticity ωx as a function of temporal lag τ (vorticity chosen here only to reduce the computational cost of time series analysis). Black curves for the original flow (solid) and coherent (dashed) and incoherent (dotted-dashed) components, the later determined using the Coifman12 wavelet in the CVE analysis (iterative Donoho-Johnstone threshold applied to the vorticity vector magnitude). Red and blue curves show results for the same components obtained using the most Gaussian vorticity and 60% compression thresholds, respectively, and in green, those for components identified using the iterative Donoho-Johnstone threshold applied to the individual vorticity amplitudes. The coherent components (dashed lines) in (a) underlie the solid curve and are thus hidden from view. The double curve for each of the incoherent components (dotted-dashed lines) in (a) reflects grid oscillation of the incoherent component at the Nyquist frequency, not seen in the original flow, due to enhanced power above the Orszag frequency kmax as a result of the wavelet filtering (Sec. 3).

FIG.5 Download High Resolution Image (.zip file) | Export Figure to PowerPoint

Tables

Table I. Compression achieved and kinetic energy and enstrophy retained in the coherent component of the flow after employing the iterative Donoho-Johnstone filter threshold on two-flows with three wavelets. Results differ when application of the filter is based on the vorticity magnitude |ω| vs. component amplitudes |ωx, y, z|.

View Table


Close
Google Calendar
ADVERTISEMENT

close