The effective diffusivity of ordered and freely evolving bubbly suspensions

2018 ◽  
Vol 840 ◽  
pp. 215-237 ◽  
Author(s):  
Aurore Loisy ◽  
Aurore Naso ◽  
Peter D. M. Spelt
Keyword(s):  
Numerical Simulations ◽  
Peclet Number ◽  
Random Media ◽  
Chemical Species ◽  
Effective Diffusivity ◽  
Passive Scalar ◽  
Solid Particles ◽  
Péclet Number ◽  
Mechanical Dispersion ◽  
Fixed Beds

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.

Dispersion in fixed beds

1985 ◽  
Vol 154 ◽  
pp. 399-427 ◽  
Author(s):  
Donald L. Koch ◽  
John F. Brady
Keyword(s):  
Peclet Number ◽  
Molecular Diffusion ◽  
Volume Fraction ◽  
Effective Diffusivity ◽  
Ensemble Averaging ◽  
Péclet Number ◽  
Physical Mechanisms ◽  
Molecular Diffusivity ◽  
Fixed Beds ◽  
Effective Diffusivities

A macroscopic equation of mass conservation is obtained by ensemble-averaging the basic conservation laws in a porous medium. In the long-time limit this ‘macro-transport’ equation takes the form of a macroscopic Fick's law with a constant effective diffusivity tensor. An asymptotic analysis in low volume fraction of the effective diffusivity in a bed of fixed spheres is carried out for all values of the Péclet number ℙ = Ua/Df, where U is the average velocity through the bed. a is the particle radius and Df is the molecular diffusivity of the solute in the fluid. Several physical mechanisms causing dispersion are revealed by this analysis. The stochastic velocity fluctuations induced in the fluid by the randomly positioned bed particles give rise to a convectively driven contribution to dispersion. At high Péclet numbers, this convective dispersion mechanism is purely mechanical, and the resulting effective diffusivities are independent of molecular diffusion and grow linearly with ℙ. The region of zero velocity in and near the bed particles gives rise to non-mechanical dispersion mechanisms that dominate the longitudinal diffusivity at very high Péclet numbers. One such mechanism involves the retention of the diffusing species in permeable particles, from which it can escape only by molecular diffusion, leading to a diffusion coefficient that grows as ℙ2. Even if the bed particles are impermeable, non-mechanical contributions that grow as ℙ ln ℙ and ℙ2 at high ℙ arise from a diffusive boundary layer near the solid surfaces and from regions of closed streamlines respectively. The results for the longitudinal and transverse effective diffusivities as functions of the Péclet number are summarized in tabular form in §6. Because the same physical mechanisms promote dispersion in dilute and dense fixed beds, the predicted Péclet-number dependences of the effective diffusivities are applicable to all porous media. The theoretical predictions are compared with experiments in densely packed beds of impermeable particles, and the agreement is shown to be remarkably good.

scholarly journals The dynamics of miscible viscous fingering from onset to shutdown

2018 ◽  
Vol 837 ◽  
pp. 520-545 ◽  
Author(s):  
Japinder S. Nijjer ◽  
Duncan R. Hewitt ◽  
Jerome A. Neufeld
Keyword(s):  
Numerical Simulations ◽  
Flow Regime ◽  
Peclet Number ◽  
Viscous Fingering ◽  
Linear Stability Theory ◽  
Late Time ◽  
Exchange Flow ◽  
Mixing Rate ◽  
Péclet Number ◽  
Single Finger

We examine the full ‘life cycle’ of miscible viscous fingering from onset to shutdown with the aid of high-resolution numerical simulations. We study the injection of one fluid into a planar two-dimensional porous medium containing another, more viscous fluid. We find that the dynamics are distinguished by three regimes: an early-time linearly unstable regime, an intermediate-time nonlinear regime and a late-time single-finger exchange-flow regime. In the first regime, the flow can be linearly unstable to perturbations that grow exponentially. We identify, using linear stability theory and numerical simulations, a critical Péclet number below which the flow remains stable for all times. In the second regime, the flow is dominated by the nonlinear coalescence of fingers which form a mixing zone in which we observe that the convective mixing rate, characterized by a convective Nusselt number, exhibits power-law growth. In this second regime we derive a model for the transversely averaged concentration which shows good agreement with our numerical experiments and extends previous empirical models. Finally, we identify a new final exchange-flow regime in which a pair of counter-propagating diffusive fingers slow exponentially. We derive an analytic solution for this single-finger state which agrees well with numerical simulations. We demonstrate that the flow always evolves to this regime, irrespective of the viscosity ratio and Péclet number, in contrast to previous suggestions.

Passive scalar statistics in high-Péclet-number grid turbulence

1998 ◽  
Vol 358 ◽  
pp. 135-175 ◽  
Author(s):  
L. MYDLARSKI ◽  
Z. WARHAFT
Keyword(s):  
Reynolds Number ◽  
Velocity Field ◽  
Peclet Number ◽  
Passive Scalar ◽  
Structure Functions ◽  
Order Structure ◽  
Grid Turbulence ◽  
Temperature Derivative ◽  
Péclet Number ◽  
Transverse Temperature

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.

scholarly journals Dispersion in the large-deviation regime. Part 2. Cellular flow at large Péclet number

2014 ◽  
Vol 745 ◽  
pp. 351-377 ◽  
Author(s):  
P. H. Haynes ◽  
J. Vanneste
Keyword(s):  
Peclet Number ◽  
Eigenvalue Problems ◽  
Large Deviation ◽  
Effective Diffusivity ◽  
Rate Function ◽  
Homogenization Theory ◽  
Péclet Number ◽  
Stagnation Points ◽  
Scalar Dispersion ◽  
Scalar Concentration

AbstractA standard model for the study of scalar dispersion through the combined effect of advection and molecular diffusion is a two-dimensional periodic flow with closed streamlines inside periodic cells. Over long time scales, the dispersion of a scalar released in this flow can be characterized by an effective diffusivity that is a factor$\mathit{Pe}^{1/2}$larger than molecular diffusivity when the Péclet number$\mathit{Pe}$is large. Here we provide a more complete description of dispersion in this regime by applying the large-deviation theory developed in Part 1 of this paper. Specifically, we derive approximations to the rate function governing the scalar concentration at large time$t$by carrying out an asymptotic analysis of the relevant family of eigenvalue problems. We identify two asymptotic regimes and, for each, make predictions for the rate function and spatial structure of the scalar. Regime I applies to distances$|\boldsymbol {x}|$from the scalar release point that satisfy$|\boldsymbol {x}|= O(\mathit{Pe}^{1/4} t)$. The concentration in this regime is isotropic at large scales, is uniform along streamlines within each cell, and varies rapidly in boundary layers surrounding the separatrices between adjacent cells. The results of homogenization theory, yielding the$O(\mathit{Pe}^{1/2})$effective diffusivity, are recovered from our analysis in the limit$|\boldsymbol {x}|\ll \mathit{Pe}^{1/4} t$. Regime II applies when$|\boldsymbol {x}|=O(\mathit{Pe}\, t/{\rm log}\, \mathit{Pe})$and is characterized by an anisotropic concentration distribution that is localized around the separatrices. A novel feature of this regime is the crucial role played by the dynamics near the hyperbolic stagnation points. A consequence is that in part of the regime the dispersion can be interpreted as resulting from a random walk on the lattice of stagnation points. The two regimes overlap so that our asymptotic results describe the scalar concentration over a large range of distances$|\boldsymbol {x}|$. They are verified against numerical solutions of the family of eigenvalue problems yielding the rate function.

scholarly journals Optimal mixing in two-dimensional plane Poiseuille flow at finite Péclet number

2014 ◽  
Vol 748 ◽  
pp. 241-277 ◽  
Author(s):  
D. P. G. Foures ◽  
C. P. Caulfield ◽  
P. J. Schmid
Keyword(s):  
Scalar Field ◽  
Poiseuille Flow ◽  
Peclet Number ◽  
Passive Scalar ◽  
Two Dimensional ◽  
Plane Poiseuille Flow ◽  
Péclet Number ◽  
Dimensional Plane ◽  
Scalar Variance ◽  
Optimal Mixing

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.

scholarly journals Chaotic advection at large Péclet number: Electromagnetically driven experiments, numerical simulations, and theoretical predictions

2014 ◽  
Vol 26 (1) ◽  
pp. 013601 ◽  
Author(s):  
Aldo Figueroa ◽  
Patrice Meunier ◽  
Sergio Cuevas ◽  
Emmanuel Villermaux ◽  
Eduardo Ramos
Keyword(s):  
Numerical Simulations ◽  
Peclet Number ◽  
Chaotic Advection ◽  
Péclet Number ◽  
Theoretical Predictions ◽  
Large Peclet Number

Approximate analysis of the containment of contaminated sites prior to remediation

2000 ◽  
Vol 42 (1-2) ◽  
pp. 319-324 ◽  
Author(s):  
H. Rubin ◽  
A. Rabideau
Keyword(s):  
Steady State ◽  
Peclet Number ◽  
Dimensionless Parameter ◽  
Contaminated Sites ◽  
Unsteady State ◽  
Péclet Number ◽  
Vertical Barrier ◽  
Time Period ◽  
Barrier Performance ◽  
Vertical Barriers

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