Coherent substructure of turbulence near the stagnation zone of a bluff body

The evolution of freestream turbulence in crossflow about a circular cylinder was studied in order to identify the existence of a coherent substructure which is the outcome of the amplification of freesteam turbulence by the stretching mechanism in diverging flow about a bluff body. Visualization of the flow events revealed the selective stretching of cross-vortex tubes and the emergence of an organized turbulent flow pattern near the cylinder stagnation zone. Significant amplification of the total turbulent energy of the streamwise fluctuating velocity was consistently monitored. Realization of selective amplification at scales larger than the neutral scale of the stagnation flow was indicated by the variation of the discrete streamwise turbulent energy. A most amplified scale, characteristic of the energy containing eddies within the coherent substructure and commensurate with the boundary-layer thickness, was detected. Penetration of the amplified turbulence into the cylinder boundary layer led to the retardation of separation and to a concurrent decrease in the drag coefficient at subcritical cylinder-diameter Reynolds numbers

Geometry of enstrophy and dissipation, grid resolution effects and proximity issues in turbulence

We perform a multi-scale non-local geometrical analysis of the structures extracted from the enstrophy and kinetic energy dissipation-rate, instantaneous fields of a numerical database of incompressible homogeneous isotropic turbulence decaying in time obtained by DNS in a periodic box. Three different resolutions are considered: 256^3, 512^3 and 1024^3 grid points, with k_(max)η(overbar) approximately 1, 2 and 4, respectively, the same initial conditions and Re_λ ≈ 77. This allows a comparison of the geometry of the structures obtained for different resolutions. For the highest resolution, structures of enstrophy and dissipation evolve in a continuous distribution from blob-like and moderately stretched tube-like shapes at the large scales to highly stretched sheet-like structures at the small scales. The intermediate scales show a predominance of tube-like structures for both fields, much more pronounced for the enstrophy field. The dissipation field shows a tendency towards structures with lower curvedness than those of the enstrophy, for intermediate and small scales. The 256^3 grid resolution case (k_(max)η(overbar) ≈ 1) was unable to detect the predominance of highly stretched sheet-like structures at the smaller scales in both fields. The same non-local methodology for the study of the geometry of structures, but without the multi-scale decomposition, is applied to two scalar fields used by existing local criteria for the eduction of tube- and sheet-like structures in turbulence, Q and [A_ij]_+, respectively, obtained from invariants of the velocity-gradient tensor and alike in the 1024^3 case. This adds the non-local geometrical characterization and classification to those local criteria, assessing their validity in educing particular geometries. Finally, we introduce a new methodology for the study of proximity issues among structures of different fields, based on geometrical considerations and non-local analysis, by taking into account the spatial extent of the structures. We apply it to the four fields previously studied. Tube-like structures of Q are predominantly surrounded by sheet-like structures of [A_ij]_+, which appear at closer distances. For the enstrophy, tube-like structures at an intermediate scale are primarily surrounded by sheets of smaller scales of the enstrophy and structures of dissipation at the same and smaller scales. A secondary contribution results from tubes of enstrophy at smaller scales appearing at farther distances. Different configurations of composite structures are presented

On the non-local geometry of turbulence

A multi-scale methodology for the study of the non-local geometry of eddy structures in turbulence is developed. Starting from a given three-dimensional field, this consists of three main steps: extraction, characterization and classification of structures. The extraction step is done in two stages. First, a multi-scale decomposition based on the curvelet transform is applied to the full three-dimensional field, resulting in a finite set of component three-dimensional fields, one per scale. Second, by iso-contouring each component field at one or more iso-contour levels, a set of closed iso-surfaces is obtained that represents the structures at that scale. The characterization stage is based on the joint probability density function (p.d.f.), in terms of area coverage on each individual iso-surface, of two differential-geometry properties, the shape index and curvedness, plus the stretching parameter, a dimensionless global invariant of the surface. Taken together, this defines the geometrical signature of the iso-surface. The classification step is based on the construction of a finite set of parameters, obtained from algebraic functions of moments of the joint p.d.f. of each structure, that specify its location as a point in a multi-dimensional ‘feature space’. At each scale the set of points in feature space represents all structures at that scale, for the specified iso-contour value. This then allows the application, to the set, of clustering techniques that search for groups of structures with a common geometry. Results are presented of a first application of this technique to a passive scalar field obtained from 5123 direct numerical simulation of scalar mixing by forced, isotropic turbulence (Reλ = 265). These show transition, with decreasing scale, from blob-like structures in the larger scales to blob- and tube-like structures with small or moderate stretching in the inertial range of scales, and then toward tube and, predominantly, sheet-like structures with high level of stretching in the dissipation range of scales. Implications of these results for the dynamical behaviour of passive scalar stirring and mixing by turbulence are discussed

Vortices within vortices: hierarchical nature of vortex tubes in turbulence

The JHU turbulence database [1] can be used with a state of the art visualisation tool [2] to generate high quality fluid dynamics videos. In this work we investigate the classical idea that smaller structures in turbulent flows, while engaged in their own internal dynamics, are advected by the larger structures. They are not advected undistorted, however. We see instead that the small scale structures are sheared and twisted by the larger scales. This illuminates the basic mechanisms of the turbulent cascade.Comment: 2 pages, 1 low quality video, 1 high quality vide

Large-scale Magnetic Structure Formation in 3D-MHD Turbulence

The inverse cascade of magnetic helicity in 3D-MHD turbulence is believed to be one of the processes responsible for large scale magnetic structure formation in astrophysical systems. In this work we present an exhaustive set of high resolution direct numerical simulations (DNS) of both forced and decaying 3D-MHD turbulence, to understand this structure formation process. It is first shown that an inverse cascade of magnetic helicity in small-scale driven turbulence does not necessarily generate coherent large-scale magnetic structures. The observed large-scale magnetic field, in this case, is severely perturbed by magnetic fluctuations generated by the small-scale forcing. In the decaying case, coherent large-scale structure form similar to those observed astronomically. Based on the numerical results the formation of large-scale magnetic structures in some astrophysical systems, is suggested to be the consequence of an initial forcing which imparts the necessary turbulent energy into the system, which, after the forcing shuts off, decays to form the large-scale structures. This idea is supported by representative examples e.g. cluster of galaxies.Comment: 21 pages in emulateapj format. 29 figures and 1 table. Accepted for publication in APJ on 11/09/201

Magnetohydrodynamic activity inside a sphere

We present a computational method to solve the magnetohydrodynamic equations in spherical geometry. The technique is fully nonlinear and wholly spectral, and uses an expansion basis that is adapted to the geometry: Chandrasekhar-Kendall vector eigenfunctions of the curl. The resulting lower spatial resolution is somewhat offset by being able to build all the boundary conditions into each of the orthogonal expansion functions and by the disappearance of any difficulties caused by singularities at the center of the sphere. The results reported here are for mechanically and magnetically isolated spheres, although different boundary conditions could be studied by adapting the same method. The intent is to be able to study the nonlinear dynamical evolution of those aspects that are peculiar to the spherical geometry at only moderate Reynolds numbers. The code is parallelized, and will preserve to high accuracy the ideal magnetohydrodynamic (MHD) invariants of the system (global energy, magnetic helicity, cross helicity). Examples of results for selective decay and mechanically-driven dynamo simulations are discussed. In the dynamo cases, spontaneous flips of the dipole orientation are observed.Comment: 15 pages, 19 figures. Improved figures, in press in Physics of Fluid