No Arabic abstract
The proliferation of turbulence in subcritical wall-bounded shear flows involves spatially localised coherent structures. Turbulent spots correspond to finite-time nonlinear responses to pointwise disturbances and are regarded as seeds of turbulence during transition. The rapid spatial decay of the turbulent fluctuations away from a spot is accompanied by large-scale flows with a robust structuration. The far field velocity field of these spots is investigated numerically using spectral methods in large domains in four different flow scenarios (plane Couette, plane Poiseuille, Couette-Poiseuille and a sinusoidal shear flow). At odds with former expectations, the planar components of the velocity field decay algebraically. These decay exponents depend only on the symmetries of the system, which here depend on the presence of an applied gradient, and not on the Reynolds number. This suggests an effective two-dimensional multipolar expansion for the far field, dominated by a quadrupolar flow component or, for asymmetric flow fields, by a dipolar flow component.
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.
Compressibility effects in a turbulent transport of temperature field are investigated applying the quasi-linear approach for small Peclet numbers and the spectral $tau$ approach for large Peclet numbers. Compressibility of a fluid flow reduces the turbulent diffusivity of the mean temperature field similarly to that for particle number density and magnetic field. However, expressions for the turbulent diffusion coefficient for the mean temperature field in a compressible turbulence are different from those for the mean particle number density and the mean magnetic field. Combined effect of compressibility and inhomogeneity of turbulence causes an increase of the mean temperature in the regions with more intense velocity fluctuations due to a turbulent pumping. Formally, this effect is similar to a phenomenon of compressible turbophoresis found previously [J. Plasma Phys. {bf 84}, 735840502 (2018)] for non-inertial particles or gaseous admixtures. Gradient of the mean fluid pressure results in an additional turbulent pumping of the mean temperature field. The latter effect is similar to turbulent barodiffusion of particles and gaseous admixtures. Compressibility of a fluid flow also causes a turbulent cooling of the surrounding fluid due to an additional sink term in the equation for the mean temperature field. There is no analog of this effect for particles.
With the aim of efficiently simulating three-dimensional multiphase turbulent flows with a phase-field method, we propose a new discretization scheme for the biharmonic term (the 4th-order derivative term) of the Cahn-Hilliard equation. This novel scheme can significantly reduce the computational cost while retaining the same accuracy as the original procedure. Our phase-field method is built on top of a direct numerical simulation solver, named AFiD (www.afid.eu) and open-sourced by our research group. It relies on a pencil distributed parallel strategy and a FFT-based Poisson solver. To deal with large density ratios between the two phases, a pressure split method [1] has been applied to the Poisson solver. To further reduce computational costs, we implement a multiple-resolution algorithm which decouples the discretizations for the Navier-Stokes equations and the scalar equation: while a stretched wall-resolving grid is used for the Navier-Stokes equations, for the Cahn-Hilliard equation we use a fine uniform mesh. The present method shows excellent computational performance for large-scale computation: on meshes up to 8 billion nodes and 3072 CPU cores, a multiphase flow needs only slightly less than 1.5 times the CPU time of the single-phase flow solver on the same grid. The present method is validated by comparing the results to previous studies for the cases of drop deformation in shear flow, including the convergence test with mesh refinement, and breakup of a rising buoyant bubble with density ratio up to 1000. Finally, we simulate the breakup of a big drop and the coalescence of O(10^3) drops in turbulent Rayleigh-Benard convection at a Rayleigh number of $10^8$, observing good agreement with theoretical results.
A mean-field theory of differential rotation in a density stratified turbulent convection has been developed. This theory is based on a combined effect of the turbulent heat flux and anisotropy of turbulent convection on the Reynolds stress. A coupled system of dynamical budget equations consisting in the equations for the Reynolds stress, the entropy fluctuations and the turbulent heat flux has been solved. To close the system of these equations, the spectral tau approach which is valid for large Reynolds and Peclet numbers, has been applied. The adopted model of the background turbulent convection takes into account an increase of the turbulence anisotropy and a decrease of the turbulent correlation time with the rotation rate. This theory yields the radial profile of the differential rotation which is in agreement with that for the solar differential rotation.