Péclet-number dependence of small-scale anisotropy of passive scalar fluctuations under a uniform mean gradient in isotropic turbulence

2020 ◽  
Vol 898 ◽  
Author(s):  
Tatsuya Yasuda ◽  
Toshiyuki Gotoh ◽  
Takeshi Watanabe ◽  
Izumi Saito

2003 ◽  
Vol 56 (6) ◽  
pp. 615-632 ◽  
Author(s):  
RA Antonia ◽  
P Orlandi

Previous reviews of the behavior of passive scalars which are convected and mixed by turbulent flows have focused primarily on the case when the Prandtl number Pr, or more generally, the Schmidt number Sc is around 1. The present review considers the extra effects which arise when Sc differs from 1. It focuses mainly on information obtained from direct numerical simulations of homogeneous isotropic turbulence which either decays or is maintained in steady state. The first case is of interest since it has attracted significant theoretical attention and can be related to decaying turbulence downstream of a grid. Topics covered in the review include spectra and structure functions of the scalar, the topology and isotropy of the small-scale scalar field, as well as the correlation between the fluctuating rate of strain and the scalar dissipation rate. In each case, the emphasis is on the dependence with respect to Sc. There are as yet unexplained differences between results on forced and unforced simulations of homogeneous isotropic turbulence. There are 144 references cited in this review article.


1998 ◽  
Vol 358 ◽  
pp. 135-175 ◽  
Author(s):  
L. MYDLARSKI ◽  
Z. WARHAFT

The statistics of a turbulent passive scalar (temperature) and their Reynolds number dependence are studied in decaying grid turbulence for the Taylor-microscale Reynolds number, Rλ, varying from 30 to 731 (21[les ]Peλ[les ]512). A principal objective is, using a single (and simple) flow, to bridge the gap between the existing passive grid-generated low-Péclet-number laboratory experiments and those done at high Péclet number in the atmosphere and oceans. The turbulence is generated by means of an active grid and the passive temperature fluctuations are generated by a mean transverse temperature gradient, formed at the entrance to the wind tunnel plenum chamber by an array of differentially heated elements. A well-defined inertial–convective scaling range for the scalar with a slope, nθ, close to the Obukhov–Corrsin value of 5/3, is observed for all Reynolds numbers. This is in sharp contrast with the velocity field, in which a 5/3 slope is only approached at high Rλ. The Obukhov–Corrsin constant, Cθ, is estimated to be 0.45–0.55. Unlike the velocity spectrum, a bump occurs in the spectrum of the scalar at the dissipation scales, with increasing prominence as the Reynolds number is increased. A scaling range for the heat flux cospectrum was also observed, but with a slope around 2, less than the 7/3 expected from scaling theory. Transverse structure functions of temperature exist at the third and fifth orders, and, as for even-order structure functions, the width of their inertial subranges dilates with Reynolds number in a systematic way. As previously shown for shear flows, the existence of these odd-order structure functions is a violation of local isotropy for the scalar differences, as is the existence of non-zero values of the transverse temperature derivative skewness (of order unity) and hyperskewness (of order 100). The ratio of the temperature derivative standard deviation along and normal to the gradient is 1.2±0.1, and is independent of Reynolds number. The refined similarity hypothesis for the passive scalar was found to hold for all Rλ, which was not the case for the velocity field. The intermittency exponent for the scalar, μθ, was found to be 0.25±0.05 with a possible weak Rλ dependence, unlike the velocity field, where μ was a strong function of Reynolds number. New, higher-Reynolds-number results for the velocity field, which smoothly follow the trends of Mydlarski & Warhaft (1996), are also presented.


2018 ◽  
Vol 840 ◽  
pp. 215-237 ◽  
Author(s):  
Aurore Loisy ◽  
Aurore Naso ◽  
Peter D. M. Spelt

We investigate the dispersion of a passive scalar such as the concentration of a chemical species, or temperature, in homogeneous bubbly suspensions, by determining an effective diffusivity tensor. Defining the longitudinal and transverse components of this tensor with respect to the direction of averaged bubble rise velocity in a zero mixture velocity frame of reference, we focus on the convective contribution thereof, this being expected to be dominant in commonly encountered bubbly flows. We first extend the theory of Kochet al.(J. Fluid Mech., vol. 200, 1989, pp. 173–188) (which is for dispersion in fixed beds of solid particles under Stokes flow) to account for weak inertial effects in the case of ordered suspensions. In the limits of low and of high Péclet number, including the inertial effect of the flow does not affect the scaling of the effective diffusivity with respect to the Péclet number. These results are confirmed by direct numerical simulations performed in different flow regimes, for spherical or very deformed bubbles and from vanishingly small to moderate values of the Reynolds number. Scalar transport in arrays of freely rising bubbles is considered by us subsequently, using numerical simulations. In this case, the dispersion is found to be convectively enhanced at low Péclet number, like in ordered arrays. At high Péclet number, the Taylor dispersion scaling obtained for ordered configurations is replaced by one characterizing a purely mechanical dispersion, as in random media, even if the level of disorder is very low.


2014 ◽  
Vol 748 ◽  
pp. 241-277 ◽  
Author(s):  
D. P. G. Foures ◽  
C. P. Caulfield ◽  
P. J. Schmid

AbstractWe consider the nonlinear optimisation of the mixing of a passive scalar, initially arranged in two layers, in a two-dimensional plane Poiseuille flow at finite Reynolds and Péclet numbers, below the linear instability threshold. We use a nonlinear-adjoint-looping approach to identify optimal perturbations leading to maximum time-averaged energy as well as maximum mixing in a freely evolving flow, measured through the minimisation of either the passive scalar variance or the so-called mix-norm, as defined by Mathew, Mezić & Petzold (Physica D, vol. 211, 2005, pp. 23–46). We show that energy optimisation appears to lead to very weak mixing of the scalar field whereas the optimal mixing initial perturbations, despite being less energetic, are able to homogenise the scalar field very effectively. For sufficiently long time horizons, minimising the mix-norm identifies optimal initial perturbations which are very similar to those which minimise scalar variance, demonstrating that minimisation of the mix-norm is an excellent proxy for effective mixing in this finite-Péclet-number bounded flow. By analysing the time evolution from initial perturbations of several optimal mixing solutions, we demonstrate that our optimisation method can identify the dominant underlying mixing mechanism, which appears to be classical Taylor dispersion, i.e. shear-augmented diffusion. The optimal mixing proceeds in three stages. First, the optimal mixing perturbation, energised through transient amplitude growth, transports the scalar field across the channel width. In a second stage, the mean flow shear acts to disperse the scalar distribution leading to enhanced diffusion. In a final third stage, linear relaxation diffusion is observed. We also demonstrate the usefulness of the developed variational framework in a more realistic control case: mixing optimisation by prescribed streamwise velocity boundary conditions.


2021 ◽  
Vol 929 ◽  
Author(s):  
John M. Lawson ◽  
Bharathram Ganapathisubramani

By coupling direct numerical simulation of homogeneous isotropic turbulence with a localised solution of the convection–diffusion equation, we model the rate of transfer of a solute (mass transfer) from the surface of small, neutrally buoyant, axisymmetric, ellipsoidal particles (spheroids) in dilute suspension within a turbulent fluid at large Péclet number, $\textit {Pe}$ . We observe that, at $\textit {Pe} = O(10)$ , the average transfer rate for prolate spheroids is larger than that of spheres with equivalent surface area, whereas oblate spheroids experience a lower average transfer rate. However, as the Péclet number is increased, oblate spheroids can experience an enhancement in mass transfer relative to spheres near an optimal aspect ratio $\lambda \approx 1/4$ . Furthermore, we observe that, for spherical particles, the Sherwood number $\textit {Sh}$ scales approximately as $\textit {Pe}^{0.26}$ over $\textit {Pe} = 1.4\times 10^{1}$ to $1.4\times 10^{4}$ , which is below the $\textit {Pe}^{1/3}$ scaling observed for inertial particles but consistent with available experimental data for tracer-like particles. The discrepancy is attributed to the diffusion-limited temporal response of the concentration boundary layer to turbulent strain fluctuations. A simple model, the quasi-steady flux model, captures both of these phenomena and shows good quantitative agreement with our numerical simulations.


2000 ◽  
Vol 42 (1-2) ◽  
pp. 319-324 ◽  
Author(s):  
H. Rubin ◽  
A. Rabideau

This study presents an approximate analytical model, which can be useful for the prediction and requirement of vertical barrier efficiencies. A previous study by the authors has indicated that a single dimensionless parameter determines the performance of a vertical barrier. This parameter is termed the barrier Peclet number. The evaluation of barrier performance concerns operation under steady state conditions, as well as estimates of unsteady state conditions and calculation of the time period requires arriving at steady state conditions. This study refers to high values of the barrier Peclet number. The modeling approach refers to the development of several types of boundary layers. Comparisons were made between simulation results of the present study and some analytical and numerical results. These comparisons indicate that the models developed in this study could be useful in the design and prediction of the performance of vertical barriers operating under conditions of high values of the barrier Peclet number.


Sign in / Sign up

Export Citation Format

Share Document