CVS and SCALES simulation of 3-D isotropic turbulence

In this work coherent vortex simulation (CVS) and stochastic coherent adaptive large eddy simulation (SCALES) simulations of decaying incompressible isotropic turbulence are compared to DNS and large eddy simulation (LES) results. Current LES relies on, at best, a zonally adapted filter width to reduce the computational cost of simulating complex turbulent flows. While there is an improvement over a uniform filter width, this approach has two limitations. First, it does not capture the high wave number components of the coherent vortices that make up the organized part of turbulent flows, thus losing essential physical information. Secondly, the flow is over-resolved in the regions between the coherent vortices, thus wasting computational resources. The SCALES approach addresses these shortcomings of LES by using a dynamic grid adaptation strategy that is able to resolve and track the most energetic coherent structures in a turbulent flow field. This corresponds to a dynamically adaptive local filter width. Unlike CVS, which we show is able to recover low order statistics with no subgrid scale (SGS) stress model, the higher compression used in SCALES necessitates that the effect of the unresolved SGS stresses must be modeled. These SGS stresses are approximated using a new dynamic eddy viscosity model based on Germano's classical dynamic procedure redefined in terms of two wavelet thresholding filters.

[1]  O. Vasilyev,et al.  Commutative discrete filtering on unstructured grids based on least-squares techniques , 2001 .

[2]  David A. Yuen,et al.  Solving PDEs using wavelets , 1998 .

[3]  Jacques Liandrat,et al.  Resolution of the 1D regularized Burgers equation using a spatial wavelet approximation , 1990 .

[4]  David A. Yuen,et al.  Modelling of viscoelastic plume-lithosphere interaction using the adaptive multilevel wavelet collocation method , 2001 .

[5]  Oleg V. Vasilyev,et al.  An Adaptive Wavelet Collocation Method for Fluid-Structure Interaction at High Reynolds Numbers , 2005, SIAM J. Sci. Comput..

[6]  W. Sweldens The Lifting Scheme: A Custom - Design Construction of Biorthogonal Wavelets "Industrial Mathematics , 1996 .

[7]  Oleg V. Vasilyev,et al.  Second-generation wavelet collocation method for the solution of partial differential equations , 2000 .

[8]  Alison L. Marsden,et al.  Construction of commutative filters for LES on unstructured meshes , 2000 .

[9]  C. Meneveau,et al.  Scale-Invariance and Turbulence Models for Large-Eddy Simulation , 2000 .

[10]  David A. Yuen,et al.  Modeling of Viscoelastic Plume-Lithosphere Interaction Using Adaptive Multilevel Wavelet Collocation , 2000 .

[11]  P. Moin,et al.  A dynamic subgrid‐scale eddy viscosity model , 1990 .

[12]  M. Farge,et al.  Non-Gaussianity and coherent vortex simulation for two-dimensional turbulence using an adaptive orthogonal wavelet basis , 1999 .

[13]  David L. Donoho,et al.  De-noising by soft-thresholding , 1995, IEEE Trans. Inf. Theory.

[14]  Ingrid Daubechies,et al.  Ten Lectures on Wavelets , 1992 .

[15]  Wim Sweldens,et al.  The lifting scheme: a construction of second generation wavelets , 1998 .

[16]  David L. Donoho,et al.  Interpolating Wavelet Transforms , 1992 .

[17]  J. Chasnov Simulation of the Kolmogorov inertial subrange using an improved subgrid model , 1991 .

[18]  S. Pope Turbulent Flows: FUNDAMENTALS , 2000 .

[19]  D. Lilly,et al.  A proposed modification of the Germano subgrid‐scale closure method , 1992 .

[20]  Oleg V. Vasilyev,et al.  Solving Multi-dimensional Evolution Problems with Localized Structures using Second Generation Wavelets , 2003 .

[21]  Nicholas K.-R. Kevlahan,et al.  An adaptive multilevel wavelet collocation method for elliptic problems , 2005 .

[22]  Parviz Moin,et al.  ADVANCES IN LARGE EDDY SIMULATION METHODOLOGY FOR COMPLEX FLOWS , 2002, Proceeding of Second Symposium on Turbulence and Shear Flow Phenomena.

[23]  M. Farge,et al.  Coherent vortex extraction in 3D turbulent flows using orthogonal wavelets. , 2001, Physical review letters.

[24]  O. Vasilyev,et al.  Stochastic coherent adaptive large eddy simulation method , 2002 .

[25]  P. Moin,et al.  A General Class of Commutative Filters for LES in Complex Geometries , 1998 .

[26]  Javier Jiménez,et al.  The structure of intense vorticity in isotropic turbulence , 1993, Journal of Fluid Mechanics.

[27]  Ugo Piomelli,et al.  Large-eddy simulation: achievements and challenges , 1999 .

[28]  P. Laguna,et al.  Signal Processing , 2002, Yearbook of Medical Informatics.

[29]  Kai Schneider,et al.  Coherent vortex simulation of three-dimensional turbulent mixing layers using orthogonal wavelets , 2005, Journal of Fluid Mechanics.

[30]  O. Vasilyev,et al.  A Fast Adaptive Wavelet Collocation Algorithm for Multidimensional PDEs , 1997 .

[31]  David A. Yuen,et al.  Modeling of compaction driven flow in poro‐viscoelastic medium using adaptive wavelet collocation method , 1998 .

[32]  D. Donoho Unconditional Bases Are Optimal Bases for Data Compression and for Statistical Estimation , 1993 .

[33]  Parviz Moin,et al.  Zonal Embedded Grids for Numerical Simulations of Wall-Bounded Turbulent Flows , 1996 .

[34]  M. Germano,et al.  Turbulence: the filtering approach , 1992, Journal of Fluid Mechanics.

[35]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[36]  M. Farge Wavelet Transforms and their Applications to Turbulence , 1992 .

[37]  Nicholas K.-R. Kevlahan,et al.  Computation of turbulent flow past an array of cylinders using a spectral method with Brinkman penalization , 2001 .

[38]  S. Mallat A wavelet tour of signal processing , 1998 .

[39]  Jie Shen,et al.  Velocity-Correction Projection Methods for Incompressible Flows , 2003, SIAM J. Numer. Anal..