We investigate through Direct Numerical Simulations (DNS) the statistical properties of turbulent flows in the inertial subrange for non-Newtonian power-law fluids. The structural invariance found for the vortex size distribution is achieved through a self-organized mechanism at the microscopic scale of the turbulent motion that adjusts, according to the rheological properties of the fluid, the ratio between the viscous dissipations inside and outside the vortices. Moreover, the deviations from the K41 theory of the structure functions' exponents reveal that the anomalous scaling exhibits a systematic nonuniversal behavior with respect to the rheological properties of the fluids.
PACS numbers:In many situations ranging from blood flow [1, 2] to atomization of slurries in industrial processing [3], one encounters non-Newtonian fluids in turbulent conditions. First experiments on turbulence in non-Newtonian fluids were already performed in 1959 [4]. Since then, most theoretical studies have focused on drag reduction [5,6], and the mathematical modeling of wall stresses and boundary layers [7][8][9]. For isotropic turbulence in dilute polymer solutions, De Angelis et al.[10] found through DNS that relaxation connecting different scales significantly alters the energy cascade.Intuitively, in the inertial subrange, molecular stresses should have a negligible influence on the motion and size of the eddies, regardless of the rheological nature of the fluid [11]. More precisely, even if a more complex constitutive law than a linear one is necessary to describe the stress-strain relation of a moving fluid, one should expect the statistical results obtained for the structure of Newtonian turbulence at the inertial subrange to remain valid. A relevant question that naturally arises is how the local rheological properties of the fluid must rearrange in space and time to comply with this alleged structural invariance. Here we provide an answer for this question by investigating through DNS the statistical properties of coherent structures of Newtonian and non-Newtonian turbulent flows in terms of distributions of eddy sizes and structure functions [12].For our numerical analysis, we consider a cubic box containing a non-Newtonian fluid and subjected to periodic boundary conditions in all three directions. The mathematical formulation of the fluid mechanics is based on the assumptions that we have an incompressible fluid flowing under isothermal conditions, for which the momentum and mass conservation equations reduce to,andwhere u and p are the velocity and pressure fields, respectively, Γ is a forcing term and T is the deviatoric stress tensor given by,Here E = ∇u + ∇u T /2 is the strain rate tensor anḋ γ = √ 2E : E its second principal invariant. The function µ (γ) defines the constitutive relation, which for a crosspower-law fluid is given byThe constants µ 1 and µ 2 are the lower and upper cutoffs, K is called the consistency index and n is the rheological exponent. Fluids with n > 1 are shear-thickening, while shear-thin...