No Arabic abstract
The similarity of the two-point correlation tensor along the streamwise direction in the axi-symmetric jet far-field is analyzed, herein its utility in spectral theory. A separable two-point correlation coefficient has been the basis for the argument that the energy-optimized basis functions along the streamwise direction are Fourier modes (from the approach of equilibrium similarity theory). This would naturally be highly desirable both from a computational and an analytical perspective. The present work, however, shows that the two-point correlation tensor multiplied by the Jacobian is not displacement invariant even in logarithmically stretched coordinates. This result directly impacts the motivation for a Fourier-based representation of the correlation function in spectral space in relation to the Proper Orthogonal Decomposition (POD) of the field. It is demonstrated that a displacement invariant form of the kernel is impossible to achieve using the suggested coordinate transformations from earlier works. This inability is shown to be related to the fundamental differences between the turbulent flow at hand and the ideal case of homogeneous turbulence.
In the current work the reconstruction of the far-field region of the turbulent axi-symmetric jet is performed in order to investigate the modal turbulence kinetic energy production contributions. The reconstruction of the field statistics is based on a semi-analytical Lumley Decomposition (LD) of the PIV sampled field using stretched amplitude decaying Fourier modes (SADFM), derived in Hodv{z}ic et al. 2019, along the streamwise coordinate. It is shown that, a wide range of modes obtain a significant amount of energy directly from the mean flow, and are therefore not exclusively dependent on a Richardson-like energy cascade even in the $kappa$-range in which the energy spectra exhibit the $-5/3$-slope. It is observed that the $-7/3$-range in the cross-spectra is fully reconstructed using a single mode in regions of high mean shear, and that shear-stresses are nearly fully reconstructed using the first two modes. These results indicate that most of the energy production related to shear-stresses is related to the first LD mode.
The current work presents a tensor formulation of the Lumley Decomposition (LD), introduced in its original form by Lumley (1967b), allowing decompositions of turbulent flow fields in curvilinear coordinates. The LD in his form is shown to enable semi-analytical decompositions of self-similar turbulent flows in general coordinate systems. The decomposition is applied to the far-field region of the fully developed turbulent axi-symmetric jet, which is expressed in stretched spherical coordinates in order to exploit the self-similar nature of the flow while ensuring the self-adjointness of the LD integral. From the LD integral it is deduced that the optimal eigenfunctions in the streamwise direction are stretched amplitude-decaying Fourier modes (SADFM). The SADFM are obtained from the LD integral upon the introduction of a streamwise-decaying weight function in the vector space definition. The wavelength of the Fourier modes is linearly increasing in the streamwise direction with an amplitude which decays with the -3/2 power of distance from the virtual origin. The streamwise evolution of the SADFM re-sembles reversed wave shoaling known from surface waves. The energy- and cross-spectra obtained from these SADFM exhibit a -5/3- and a -7/3-slope region, respectively, as would be expected for regular Fourier modes in homogeneous and constant shear flows. The approach introduced in this work can be extended to other flows which admit to equilibrium similarity, such that a Fourier-based decomposition along inhomogeneous flow directions can be performed.
The seminal Batchelor-Greens (BG) theory on the hydrodynamic interaction of two spherical particles of radii a suspended in a viscous shear flow neglects the effect of the boundaries. In the present paper we study how a plane wall modifies this interaction. Using an integral equation for the surface traction we derive the expression for the particles relative velocity as a sum of the BGs velocity and the term due to the presence of a wall at finite distance, z_0. Our calculation is not the perturbation theory of the BG solution, so the contribution due to the wall is not necessarily small. The distance at which the wall significantly alters the particles interaction scales as z_0^{3/5}. The phase portrait of the particles relative motion is different from the BG theory, where there are two singly-connected regions of open and closed trajectories both of infinite volume. For finite z_0, there is a new domain of closed (dancing) and open (swapping) trajectories. The width of this region behaves as 1/z_0. Along the swapping trajectories, that have been previously observed numerically, the incoming particle is turning back after the encounter with the reference particle, rather than passing it by, as in the BG theory. The region of dancing trajectories has infinite volume. We found a one-parameter family of equilibrium states, overlooked previously, whereas the pair of spheres flows as a whole without changing its configuration. These states are marginally stable and their perturbation yields a two-parameter family of the dancing trajectories, where the particle is orbiting around a fixed point in a frame co-moving with the reference particle. We suggest that the phase portrait obtained at z_0>>a is topologically stable and can be extended down to rather small z_0 of several particle diameters. We confirm this by direct numerical simulations of the Navier-Stokes equations with z_0=5a.
For wall-bounded turbulent flows, Townsends attached eddy hypothesis proposes that the logarithmic layer is populated by a set of energetic and geometrically self-similar eddies. These eddies scale with a single length scale, their distance to the wall, while their velocity scale remains constant across their size range. To investigate the existence of such structures in fully developed turbulent pipe flow, stereoscopic particle image velocimetry measurements were performed in two parallel cross-sectional planes, spaced apart by a varying distance from 0 to 9.97$R$, for $Re_tau = 1310$, 2430 and 3810. The instantaneous turbulence structures are sorted by width using an azimuthal Fourier decomposition, allowing us to create a set of average eddy velocity profiles by performing an azimuthal alignment process. The resulting eddy profiles exhibit geometric self-similar behavior in the $(r,theta)$-plane for eddies with spanwise length scales ($lambda_theta/R$) spanning from 1.03 to 0.175. The streamwise similarity is further investigated using two-point correlations between the two planes, where the structures exhibit a self-similar behaviour with length scales ($lambda_theta/R$) ranging from approximately $0.88$ to $0.203$. The candidate structures thereby establish full three-dimensional geometrically self-similarity for structures with a volumetric ratio of $1:80$. Beside the geometric similarity, the velocity magnitude also exhibit self-similarity within these ranges. However, the velocity scale depends on eddy size, and follow the trends based on the scaling arguments proposed by cite{Perry1986}.
In an earlier paper (Wan et al. 2012), the authors showed that a similarity solution for anisotropic incompressible 3D magnetohydrodynamic (MHD) turbulence, in the presence of a uniform mean magnetic field $vB_0$, exists if the ratio of parallel to perpendicular (with respect to $vB_0$) similarity length scales remains constant in time. This conjecture appears to be a rather stringent constraint on the dynamics of decay of the energy-containing eddies in MHD turbulence. However, we show here, using direct numerical simulations, that this hypothesis is indeed satisfied in incompressible MHD turbulence. After an initial transient period, the ratio of parallel to perpendicular length scales fluctuates around a steady value during the decay of the eddies. We show further that a Taylor--Karman-like similarity decay holds for MHD turbulence in the presence of a mean magnetic field. The effect of different parameters, including Reynolds number, DC field strength, and cross-helicity, on the nature of similarity decay is discussed.