Passive Scalar Mixing Research Articles

Optimal stirring strategies for passive scalar mixing

We address the challenge of optimal incompressible stirring to mix an initially inhomogeneous distribution of passive tracers. As a quantitative measure of mixing we adopt the H−1 norm of the scalar fluctuation field, equivalent to the (square root of the) variance of a low-pass filtered image of the tracer concentration field. First we establish that this is a useful gauge even in the absence of molecular diffusion: its vanishing as t → ∞ is evidence of the stirring flow's mixing properties in the sense of ergodic theory. Then we derive absolute limits on the total amount of mixing, as a function of time, on a periodic spatial domain with a prescribed instantaneous stirring energy or stirring power budget. We subsequently determine the flow field that instantaneously maximizes the decay of this mixing measure – when such a flow exists. When no such ‘steepest descent’ flow exists (a possible but non-generic situation), we determine the flow that maximizes the growth rate of the H−1 norm's decay rate. This local-in-time optimal stirring strategy is implemented numerically on a benchmark problem and compared to an optimal control approach using a restricted set of flows. Some significant challenges for analysis are outlined.

On the role of vortical structures for turbulent mixing using direct numerical simulation and wavelet-based coherent vorticity extraction

The influence of vortical structures on the transport and mixing of passive scalars is investigated. Initial conditions are taken from a direct numerical simulation database of forced homogeneous isotropic turbulence, with passive scalar fluctuations, driven by a uniform mean gradient, are performed for Taylor microscale Reynolds numbers (R λ) of 140 and 240, and Schmidt numbers 1/8 and 1. For each R λ, after reaching a fully developed turbulent regime, which is statistically steady, the Coherent Vorticity Extraction is applied to the flow. It is shown that the coherent part is able to preserve the vortical structures with only less than 4% of wavelet coefficients while retaining 99.9% of energy. In contrast, the incoherent part is structureless and contains negligible energy. By taking the total, coherent and incoherent velocity fields in turn as initial conditions, new simulations were performed without forcing while the uniform mean scalar gradient is maintained. It is found that the results from simulations with total and coherent velocity fields as initial conditions are very similar. In contrast, the time integration of the incoherent flow exhibits its primarily dissipative nature. The evolutions of passive scalars at Schmidt numbers 1/8 and 1 advected by the total, coherent or incoherent velocity suggest that the vortical structures retained in the coherent part play a dominant role in turbulent transport and mixing. Indeed, the total and coherent flows give almost the same time evolution of the scalar variance, scalar flux and mean scalar dissipation, while the incoherent flow only gives rise to weak scalar diffusion.

Mixing of passive tracers in the decay Batchelor regime of a channel flow

We report detailed quantitative studies of passive scalar mixing in a curvilinear channel flow, where elastic turbulence in a dilute polymer solution of high molecular weight polyacrylamide in a high viscosity water-sugar solvent was achieved. For quantitative investigation of mixing, a detailed study of the profiles of mean longitudinal and radial components of the velocity in the channel as a function of Wi was carried out. Besides, a maximum of the average value as well as a rms of the longitudinal velocity was used to determine the threshold of the elastic instability in the channel flow. The rms of the radial derivatives of the longitudinal and radial velocity components was utilized to define the control parameters of the problem, the Weissenberg Wiloc and the Péclet Pe numbers. The main result of these studies is the quantitative test of the theoretical prediction about the value of the mixing length in the decay Batchelor regime. The experiment shows large quantitative discrepancy, more than 200 times in the value of the coefficient C, which appears in the theoretical expression for the mixing length, but with the predicted scaling relation. There are two possible reasons to this discrepancy. First is the assumption made in the theory about the δ-correlated velocity field, which is in odds with the experimental observations. Second, and probably a more relevant suggestion for the significantly increased mixing length and thus reduced mixing efficiency, is the observed jets, the rare, localized, and vigorous ejection of the scalar trapped near the wall, which protrudes into the peripheral region as well as the bulk. They are first found in the recent numerical calculations and then observed in the experiment reported. The jets definitely strongly reduce the mixing efficiency in particular in the peripheral region and so can lead to considerable increase of the mixing length. We hope that this result will initiate further numerical calculations of the mixing length. Finally, we analyze statistical properties of the mixing in the decay Batchelor regime by studying the power spectra, the decay exponents scaling, the structure functions of a tracer and moments of PDF of passive scalar increments, and the temporal and spatial correlation functions and find rather satisfactory agreement with theory.

MIXING IN SUPERSONIC TURBULENCE

In many astrophysical environments, mixing of heavy elements occurs in the presence of a supersonic turbulent velocity field. Here we carry out the first systematic numerical study of such passive scalar mixing in isothermal supersonic turbulence. Our simulations show that the ratio of the scalar mixing timescale, $\tau_{\rm c}$, to the flow dynamical time, $\tau_{\rm dyn}$ (defined as the flow driving scale divided by the rms velocity), increases with the Mach number, $M$, for $M \lsim3$, and becomes essentially constant for $M \gsim3.$ This trend suggests that compressible modes are less efficient in enhancing mixing than solenoidal modes. However, since the majority of kinetic energy is contained in solenoidal modes at all Mach numbers, the overall change in $\tau_{\rm c}/\tau_{\rm dyn}$ is less than 20\% over the range $1 \lsim M \lsim 6$. At all Mach numbers, if pollutants are injected at around the flow driving scale, $\tau_{\rm c}$ is close to $\tau_{\rm dyn}.$ This suggests that scalar mixing is driven by a cascade process similar to that of the velocity field. The dependence of $\tau_{\rm c}$ on the length scale at which pollutants are injected into flow is also consistent with this cascade picture. Similar behavior is found for the variance decay timescales for scalars without continuing sources. Extension of the scalar cascade picture to the supersonic regime predicts a relation between the scaling exponents of the velocity and the scalar structure functions, with the scalar structure function becoming flatter as the velocity scaling steepens with Mach number. Our measurements of the volume-weighted velocity and scalar structure functions confirm this relation for $M\lsim 2,$ but show discrepancies at $M \gsim 3$.

The Batchelor Spectrum for Mixing of Passive Scalars in Isotropic Turbulence

We examine the support for the Batchelor spectrum from well-resolved simulations of high-Schmidt-number mixing in isotropic turbulence, and resolve a conundrum with respect to the numerical value of the prefactor, also known as the Batchelor constant. Our conclusion is that the most probable value of the most compressive principal strain rate is more relevant than its mean, at least asymptotically.

Lagrangian views on turbulent mixing of passive scalars

The Lagrangian view of passive scalar turbulence has recently produced interesting results and interpretations. Innovations in theory, experiments, simulations and data analysis of Lagrangian turbulence are reviewed here in brief. Part of the review is closely related to the so-called Kraichnan model for the advection of the passive scalar in synthetic turbulence. Possible implications for a better understanding of the passive scalar mixing in Navier-Stokes turbulence are also discussed.

Investigation of passive scalar mixing in a confined rectangular wake using simultaneous PIV and PLIF

A combined particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) system was employed to investigate the turbulent mixing in a confined liquid-phase rectangular-wake flow with a Reynolds number based on hydraulic diameter of 37,500 and a Schmidt number of 1250. The simultaneous velocity and concentration field data were analyzed for flow statistics such as turbulent fluxes, turbulent viscosity and diffusivity, and turbulent Schimdt number. The streamwise and transverse turbulent fluxes were found to be of the same magnitude. The turbulent flux vector was not aligned with the mean concentration gradient. The turbulent Schimdt number was about 0.8. The spatial correlations of turbulent fluxes and concentration fluctuations were evaluated. In the R u ′ ϕ ′ correlation field, there were a positive and a negative vertically oriented elliptical correlation region, which were symmetric around the basis point. The R v ′ ϕ ′ correlation region was a horizontally oriented ellipse with negative values of the correlation coefficient. The correlation field of R ϕ ′ ϕ ′ was also an ellipse with a horizontal major axis. The behavior of large-scale structures in both the velocity and concentration field was studied using linear stochastic estimation with a defined event of concentration fluctuation. Vortex streets were observed in the estimated velocity fields. The streamwise growth of the structure size increased linearly initially but then grew more slowly.

Scalar energy fluctuations in Large-Eddy Simulation of turbulent flames: Statistical budgets and mesh quality criterion

Large-Eddy Simulation (LES) provides space-filtered quantities to compare with measurements, which usually have been obtained using a different filtering operation; hence, numerical and experimental results can be examined side-by-side in a statistical sense only. Instantaneous, space-filtered and statistically time-averaged signals feature different characteristic length-scales, which can be combined in dimensionless ratios. From two canonical manufactured turbulent solutions, a turbulent flame and a passive scalar turbulent mixing layer, the critical values of these ratios under which measured and computed variances (resolved plus sub-grid scale) can be compared without resorting to additional residual terms are first determined. It is shown that actual Direct Numerical Simulation can hardly accommodate a sufficiently large range of length-scales to perform statistical studies of LES filtered reactive scalar-fields energy budget based on sub-grid scale variances; an estimation of the minimum Reynolds number allowing for such DNS studies is given. From these developments, a reliability mesh criterion emerges for scalar LES and scaling for scalar sub-grid scale energy is discussed.

Large Eddy Simulation of a turbulent swirling coaxial jet

This work uses the Large Eddy Simulation (LES) technique to study velocity and passive scalar mixing along with intermittency of a spatially evolving turbulent coaxial swirl jet. The simulations captured the potential core and also predicted high level turbulence intensities in the inner mixing regions. The Probability Density Functions (PDFs) and radial intermittency plots revealed an intermittent mixing behaviour especially in the outer region of the flow where the fluctuations of velocity rapidly change from rotational to irrotational and vice versa. The PDF and radial intermittency profiles exhibit Gaussian and non-Gaussian distributions close to the jet centreline and away from the centreline, respectively.

Resolution effects and scaling in numerical simulations of passive scalar mixing in turbulence

The effects of finite grid resolution on the statistics of small scales in direct numerical simulations of turbulent mixing of passive scalars are addressed in this paper. Simulations at up to 2048 3 grid points with grid spacing Δ x varied from about 2 to 1/2 Batchelor scales ( η B ) show that most conclusions on Schmidt number ( S c ) dependence from prior work at less stringent resolution remain qualitatively correct, although simulations at resolution Δ x ≈ η B are preferred and will give adequate results for many important quantities including the scalar dissipation intermittency exponent and structure functions at moderately high orders. For S c ≥ 1 , since η B = η S c − 1 / 2 (where η is the Kolmogorov scale), the requirement Δ x ≈ η B is more stringent than the corresponding criterion Δ x ≈ η for the velocity field, which is thus well resolved in simulations aimed at high Schmidt number mixing. A simple argument is given to help interpret the effects of Schmidt and Reynolds numbers on trends towards local isotropy and saturation of intermittency at high Schmidt number. The present results also provide evidence for a trend to isotropy at high Reynolds number with fixed S c = 1.0 . This is a new observation apparently not detected in less well resolved simulations in the past, and will require further investigation in the future.

Eigenmode analysis of scalar transport in distributive mixing

In this study, we explore the spectral properties of the distribution matrices of the mapping method and its relation to the distributive mixing of passive scalars. The spectral (or eigenvector-eigenvalue) decomposition of these matrices constitutes discrete approximations to the eigenmodes of the continuous advection operator in periodic flows. The eigenvalue spectrum always lies within the unit circle and due to mass conservation, always accommodates an eigenvalue equal to one with trivial (uniform) eigenvector. The asymptotic state of a fully chaotic mixing flow is dominated by the eigenmode corresponding with the eigenvalue closest to the unit circle (“dominant eigenmode”). This eigenvalue determines the decay rate; its eigenvector determines the asymptotic mixing pattern. The closer this eigenvalue value is to the origin, the faster is the homogenization by the chaotic mixing. Hence, its magnitude can be used as a quantitative mixing measure for comparison of different mixing protocols. In nonchaotic cases, the presence of islands results in eigenvalues on the unit circle and associated eigenvectors demarcating the location of these islands. Eigenvalues on the unit circle thus are qualitative indicators of inefficient mixing; the properties of its eigenvectors enable isolation of the nonmixing zones. Thus important fundamental aspects of mixing processes can be inferred from the eigenmode analysis of the mapping matrix. This is elaborated in the present paper and demonstrated by way of two different prototypical mixing flows: the time-periodic sine flow and the spatially periodic partitioned-pipe mixer.

Dynamics of passive scalars and tracers advected by a two-dimensional tripolar vortex

The dynamics of passive scalars and tracers during the formation and subsequent persistence of a laminar tripolar vortex, obtained through an unstable monopolar vortex seeded with a k = 2 azimuthal perturbation, is investigated. Two-dimensional direct numerical simulations of passive scalars with Schmidt numbers Sc = 0.1, 1, 10 and 100 are performed. The scalar variance for the four cases is analysed, as well as the different dispersion patterns up to 10 times greater than the time for formation of the tripolar vortex. During the formation of the tripole, an accelerated scalar dissipation is observed. That dissipation is connected to the advection-dominated processes associated with the growth of the perturbation mode. During that process, the patterns of mixing of the different passive scalars are very much the same as for vorticity. This stage of accelerated dissipation is preceded and followed by stages of diffusion-dominated scalar dissipation. Passive Lagrangian tracers are used to explore the transport of fluid elements during the evolution, and to provide a detailed view of the tripolar vortex formation and behaviour for longer times. Chaotic mixing was studied by examining patterns of spatial variation of finite-time Lyapunov exponents. As the perturbation grows and the tripolar vortex is formed, two large regions of regular flow, divided by a region of chaotic flow, form for each satellite. When the tripole is fully formed, it is composed of three distinct regular regions, corresponding to the core of negative vorticity and the two satellites of positive vorticity. The comparison between the evolution for vorticity, concentration, randomly distributed particles and Lyapunov exponents shows that transport occurs mainly in the regions of chaotic flow that surround the tripolar vortex after its formation. For longer times, both the chaotic/regular flow interfaces and the vorticity gradients are responsible for the integrity of the tripolar system.

Fine-scale statistics of temperature and its derivatives in convective turbulence

We study the fine-scale statistics of temperature and its derivatives in turbulent Rayleigh–Bénard convection. Direct numerical simulations are carried out in a cylindrical cell with unit aspect ratio filled with a fluid with Prandtl number equal to 0.7 for Rayleigh numbers between 107 and 109. The probability density function of the temperature or its fluctuations is found to be always non-Gaussian. The asymmetry and strength of deviations from the Gaussian distribution are quantified as a function of the cell height. The deviations of the temperature fluctuations from the local isotropy, as measured by the skewness of the vertical derivative of the temperature fluctuations, decrease in the bulk, but increase in the thermal boundary layer for growing Rayleigh number, respectively. Similarly to the passive scalar mixing, the probability density function of the thermal dissipation rate deviates significantly from a log-normal distribution. The distribution is fitted well by a stretched exponential form. The tails become more extended with increasing Rayleigh number which displays an increasing degree of small-scale intermittency of the thermal dissipation field for both the bulk and the thermal boundary layer. We find that the thermal dissipation rate due to the temperature fluctuations is not only dominant in the bulk of the convection cell, but also yields a significant contribution to the total thermal dissipation in the thermal boundary layer. This is in contrast to the ansatz used in scaling theories and can explain the differences in the scaling of the total thermal dissipation rate with respect to the Rayleigh number.

Passive scalar structures in peripheral regions of random flows

The statistical properties of the passive scalar near walls in random flows assuming a weakness of its diffusion have been investigated. Then, at advanced stages of the passive scalar mixing, its unmixed residue is concentrated in a narrow diffusive layer near the wall. The numerical simulations have revealed the structures responsible for the passive scalar transport to the bulk; these are passive scalar tongues pulled from the diffusive boundary layer. The passive scalar integrated along the wall possesses a well-pronounced scaling behavior. An analytical scheme, giving exponents of the integral passive scalar moments has been proposed. The exponents reasonably agree with the calculations in 3d.

An improved mixing model providing joint statistics of scalar and scalar dissipation

For the calculation of nonpremixed turbulent flames with thin reaction zones the joint probability density function (PDF) of the mixture fraction and its dissipation rate plays an important role. The corresponding PDF transport equation involves a mixing model for the closure of the molecular mixing term. Here, the parameterized scalar profile (PSP) mixing model is extended to provide the required joint statistics. Model predictions are validated using direct numerical simulation (DNS) data of a passive scalar mixing in a statistically homogeneous turbulent flow. Comparisons between the DNS and the model predictions are provided, which involve different initial scalar-field lengthscales.

Conditional statistics of passive-scalar mixing in a confined wake flow

Experimental results for conditional statistics in a confined turbulent wake were compared to the linear model and the transported probability density function (PDF) model for predicting the velocity conditioned on scalar mixture fraction and to a linear model for predicting the scalar mixture fraction conditioned on velocity. For velocity conditioned on scalar mixture fraction, the linear model was found to accurately predict the conditional transverse velocity, but for the conditional streamwise velocity, agreement between the experiments and the linear model was poor. The transverse conditional velocity was accurately predicted by the transported PDF model, but for the streamwise conditional velocity, a modified transported PDF model that incorporated two components of the scalar diffusivity tensor had to be used. Finally, it was found that a linear model for predicting the fluctuating scalar conditioned on velocity performed poorly due to the strongly non-Gaussian scalar PDF at the downstream locations investigated.

Scalar-energy spectra in simulations of Sc⪢1 mixing by turbulent jets using the nonlinear large-eddy simulation method

The nonlinear large-eddy simulation (nLES) method is applied to the first numerical study of passive-scalar mixing by a turbulent shear flow at a high Schmidt number (Sc⪢1). The work is intended to address inconsistencies between previous studies concerning the formation of power-law scaling in the scalar-energy spectra at viscous-convective scales. Results are reported for LES of a round turbulent jet at Sc=1024 and ReD=2000. The nLES method is first shown to recover the large-scale jet structure, including the self-similarity of far-field scalar moments. Scalar timeseries and spatial data produce the first power spectra from a LES shear-flow study that exhibits k−1 scaling at viscous-convective scales, consistent with the analysis of Batchelor [J. Fluid Mech. 5, 113 (1959)] and recent direct numerical studies of simpler Sc⪢1 flows.

Comparative experimental study of local mixing of active and passive scalars in turbulent thermal convection

We investigate experimentally the statistical properties of active and passive scalar fields in turbulent Rayleigh-Bénard convection in water, at Ra approximately 10;{10} . Both the local concentration of fluorescence dye and the local temperature are measured near the sidewall of a rectangular cell. It is found that, although they are advected by the same turbulent flow, the two scalars distribute differently. This difference is twofold, i.e., both the quantities themselves and their small-scale increments have different distributions. Our results show that there is a certain buoyant scale based on time domain, i.e., the Bolgiano time scale t_{B} , above which buoyancy effects are significant. Above t_{B} , temperature is active and is found to be more intermittent than concentration, which is passive. This suggests that the active scalar possesses a higher level of intermittency in turbulent thermal convection. It is further found that the mixing of both scalar fields are isotropic for scales larger than t_{B} even though buoyancy acts on the fluid in the vertical direction. Below t_{B} , temperature is passive and is found to be more anisotropic than concentration. But this higher degree of anisotropy is attributed to the higher diffusivity of temperature over that of concentration. From the simultaneous measurements of temperature and concentration, it is shown that two scalars have similar autocorrelation functions and there is a strong and positive correlation between them.

Genesis and bifurcations of unstable periodic orbits in a jet flow

We study the origin and bifurcations of typical classes of unstable periodic orbits in a jet flow that was introduced before as a kinematic model of chaotic advection, transport and mixing of passive scalars in meandering oceanic and atmospheric currents. A method to detect and locate the unstable periodic orbits and classify them by the origin and bifurcations is developed. We consider in detail period-1 and period-4 orbits playing an important role in chaotic advection. We introduce five classes of period-4 orbits: western and eastern ballistic ones, whose origin is associated with ballistic resonances of the fourth-order, rotational ones, associated with rotational resonances of the second and fourth orders and rotational–ballistic ones associated with a rotational–ballistic resonance. It is a new kind of unstable periodic orbits that may appear in a chaotic flow with jets and/or circulation cells. Varying the perturbation amplitude, we track out the origin and bifurcations of the orbits for each class.

On the Non-Local Geometry of Turbulence

A multi-scale methodology for the study of the non-local geometry of eddy structures in turbulence is developed. Starting from a given three-dimensional field, this consists of three main steps: extraction, characterization and classification of structures. The extraction step is done in two stages. First, a multi-scale decomposition based on the curvelet transform is applied to the full three-dimensional field, resulting in a finite set of component three-dimensional fields, one per scale. Second, by iso-contouring each component field at one or more iso-contour levels, a set of closed iso-surfaces is obtained that represents the structures at that scale. The characterization stage is based on the joint probability density function (p.d.f.), in terms of area coverage on each individual iso-surface, of two differential-geometry properties, the shape index and curvedness, plus the stretching parameter, a dimensionless global invariant of the surface. Taken together, this defines the geometrical signature of the iso-surface. The classification step is based on the construction of a finite set of parameters, obtained from algebraic functions of moments of the joint p.d.f. of each structure, that specify its location as a point in a multi-dimensional ‘feature space’. At each scale the set of points in feature space represents all structures at that scale, for the specified iso-contour value. This then allows the application, to the set, of clustering techniques that search for groups of structures with a common geometry. Results are presented of a first application of this technique to a passive scalar field obtained from 5123 direct numerical simulation of scalar mixing by forced, isotropic turbulence (Reλ = 265). These show transition, with decreasing scale, from blob-like structures in the larger scales to blob- and tube-like structures with small or moderate stretching in the inertial range of scales, and then toward tube and, predominantly, sheet-like structures with high level of stretching in the dissipation range of scales. Implications of these results for the dynamical behaviour of passive scalar stirring and mixing by turbulence are discussed.