TY - JOUR

T1 - Turbulence is an Ineffective Mixer when Schmidt Numbers Are Large

AU - Buaria, Dhawal

AU - Clay, Matthew P.

AU - Sreenivasan, Katepalli R.

AU - Yeung, P. K.

N1 - Publisher Copyright:
©=2D⟨|∇θ|2

PY - 2021/2/19

Y1 - 2021/2/19

N2 - We solve the advection-diffusion equation for a stochastically stationary passive scalar θ, in conjunction with forced 3D Navier-Stokes equations, using direct numerical simulations in periodic domains of various sizes, the largest being 8192^{3}. The Taylor-scale Reynolds number varies in the range 140-650 and the Schmidt number Sc≡ν/D in the range 1-512, where ν is the kinematic viscosity of the fluid and D is the molecular diffusivity of θ. Our results show that turbulence becomes an ineffective mixer when Sc is large. First, the mean scalar dissipation rate ⟨χ⟩=2D⟨|∇θ|^{2}⟩, when suitably nondimensionalized, decreases as 1/logSc. Second, 1D cuts through the scalar field indicate increasing density of sharp fronts on larger scales, oscillating with large excursions leading to reduced mixing, and additionally suggesting weakening of scalar variance flux across the scales. The scaling exponents of the scalar structure functions in the inertial-convective range appear to saturate with respect to the moment order and the saturation exponent approaches unity as Sc increases, qualitatively consistent with 1D cuts of the scalar.

AB - We solve the advection-diffusion equation for a stochastically stationary passive scalar θ, in conjunction with forced 3D Navier-Stokes equations, using direct numerical simulations in periodic domains of various sizes, the largest being 8192^{3}. The Taylor-scale Reynolds number varies in the range 140-650 and the Schmidt number Sc≡ν/D in the range 1-512, where ν is the kinematic viscosity of the fluid and D is the molecular diffusivity of θ. Our results show that turbulence becomes an ineffective mixer when Sc is large. First, the mean scalar dissipation rate ⟨χ⟩=2D⟨|∇θ|^{2}⟩, when suitably nondimensionalized, decreases as 1/logSc. Second, 1D cuts through the scalar field indicate increasing density of sharp fronts on larger scales, oscillating with large excursions leading to reduced mixing, and additionally suggesting weakening of scalar variance flux across the scales. The scaling exponents of the scalar structure functions in the inertial-convective range appear to saturate with respect to the moment order and the saturation exponent approaches unity as Sc increases, qualitatively consistent with 1D cuts of the scalar.

UR - http://www.scopus.com/inward/record.url?scp=85102406648&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85102406648&partnerID=8YFLogxK

U2 - 10.1103/PhysRevLett.126.074501

DO - 10.1103/PhysRevLett.126.074501

M3 - Article

C2 - 33666462

AN - SCOPUS:85102406648

VL - 126

JO - Physical Review Letters

JF - Physical Review Letters

SN - 0031-9007

IS - 7

M1 - 074501

ER -