Phase Distribution in Buoyancy-Driven Bubbly Flows

Author(s):  
Asghar Esmaeeli ◽  
Chan Ching ◽  
Mamdouh Shoukri

This study aims to investigate the effect of topology change on the rise velocity of bubbly flows and the phase distribution in a channel at a moderate Reynolds number. A front tracking/finite difference method is used to solve the momentum equation inside and outside deformable bubbles. It is found that bubble/bubble coalescence enhances the average rise velocity of the bubbles dramatically and also increases the fluctuations of the liquid velocity. Examination of the pair distribution function shows that the flow becomes more non-homogeneous as a result of topology change.

2002 ◽  
Vol 466 ◽  
pp. 53-84 ◽  
Author(s):  
BERNARD BUNNER ◽  
GRÉTAR TRYGGVASON

Direct numerical simulations of the motion of up to 216 three-dimensional buoyant bubbles in periodic domains are presented. The bubbles are nearly spherical and have a rise Reynolds number of about 20. The void fraction ranges from 2% to 24%. Part 1 analysed the rise velocity and the microstructure of the bubbles. This paper examines the fluctuation velocities and the dispersion of the bubbles and the ‘pseudo-turbulence’ of the liquid phase induced by the motion of the bubbles. It is found that the turbulent kinetic energy increases with void fraction and scales with the void fraction multiplied by the square of the average rise velocity of the bubbles. The vertical Reynolds stress is greater than the horizontal Reynolds stress, but the anisotropy decreases when the void fraction increases. The kinetic energy spectrum follows a power law with a slope of approximately −3.6 at high wavenumbers.


2017 ◽  
Vol 816 ◽  
pp. 94-141 ◽  
Author(s):  
Aurore Loisy ◽  
Aurore Naso ◽  
Peter D. M. Spelt

Various expressions have been proposed previously for the rise velocity of gas bubbles in homogeneous steady bubbly flows, generally a monotonically decreasing function of the bubble volume fraction. For suspensions of freely moving bubbles, some of these are of the form expected for ordered arrays of bubbles, and vice versa, as they do not reduce to the behaviour expected theoretically in the dilute limit. The microstructure of weakly inhomogeneous bubbly flows not being known generally, the effect of microstructure is an important consideration. We revisit this problem here for bubbly flows at small to moderate Reynolds number values for deformable bubbles, using direct numerical simulation and analysis. For ordered suspensions, the rise velocity is demonstrated not to be monotonically decreasing with volume fraction due to cooperative wake interactions. The fore-and-aft asymmetry of an isolated ellipsoidal bubble is reversed upon increasing the volume fraction, and the bubble aspect ratio approaches unity. Recent work on rising bubble pairs is used to explain most of these results; the present work therefore forms a platform of extending the former to suspensions of many bubbles. We adopt this new strategy also to support the existence of the oblique rise of ordered suspensions, the possibility of which is also demonstrated analytically. Finally, we demonstrate that most of the trends observed in ordered systems also appear in freely evolving suspensions. These similarities are supported by prior experimental measurements and attributed to the fact that free bubbles keep the same neighbours for extended periods of time.


2021 ◽  
Vol 932 ◽  
Author(s):  
Vikash Pandey ◽  
Dhrubaditya Mitra ◽  
Prasad Perlekar

We present a direct numerical simulation (DNS) study of buoyancy-driven bubbly flows in the presence of large-scale driving that generates turbulence. On increasing the turbulence intensity: (a) the bubble trajectories become more curved and (b) the average rise velocity of the bubbles decreases. We find that the energy spectrum of the flow shows a pseudo-turbulence scaling for length scales smaller than the bubble diameter and a Kolmogorov scaling for scales larger than the bubble diameter. We conduct a scale-by-scale energy budget analysis to understand the scaling behaviour observed in the spectrum. Although our bubbles are weakly buoyant, the statistical properties of our DNS are consistent with the experiments that investigate turbulence modulation by air bubbles in water.


Author(s):  
Alain Cartellier

Systems involving swarms of bubbles in an otherwise laminar continuous phase are common in industrial processes. In some cases, the gas is injected to ensure a given chemical reaction (bubble columns in oil industry) or to sustain a biochemical process (aeration tanks in waste water treatment plants). Gas inclusions can also appear due to the reaction itself (electrolysis cells, anaerobic digestion). In others circumstances, the gas phase is chemically passive and it is introduced mainly to favor mixing and/or separation (flotation devices). In these processes, it is desirable to access parameters such as the pressure drop, the mean void fraction, the bubble size distributions. In addition, their optimum functioning often depends on the transverse distribution of phasic quantities. Even if break-up/coalescence mechanisms are discarded, it happens that refined descriptions of such laminar dispersed flows has not yet reached a truly predictive status. On one hand, the Reynolds stresses reduce to the so-called bubble-induced agitation (or pseudo-turbulence) so that the interactions between inclusions and shear-induced turbulence need not to be accounted for. Yet, another complexity emerges because of strong and non-trivial couplings between phases. In particular, bubble-bubble interactions have a crucial effect on the induced agitation and consequently on the phase distribution. How to properly account for these interactions in an average description is still a matter of controversy. This presentation will highlight the importance of coupling mechanisms arising in laminar bubbly flows. Available experiments will be presented that illustrate the variety of phase organizations observed in stable Poiseuille bubbly flows [1–8]. It will be shown that some characteristics such as the mean void fraction and the wall shear stress are accessible through simplified models based on axial momentum balances [9,10]. On another hand, predictions of the phase distribution require solving transverse mechanical equilibria: the later are sensitive to many parameters, and in addition, they involve various coupling modes between phases. To overcome the corresponding modeling difficulties, a hybrid model has been developed in the spirit of approaches combining kinetic theory and classical continuum mechanics [see for example 11–13]. Compared with classical Eulerian two-fluid model, this framework provides, at least in the limit of dilute systems, a mean to derive closure laws [14–17]. These improvements will be illustrated for the interfacial momentum exchanges and the extra deformation tensor. In particular, the behavior of these coupling terms near walls will be shown to have important consequences on the phase distribution by the mediation of the continuous phase velocity profile. Concerning dispersion mechanisms, experimental information available on the bubble-induced agitation and on the dispersed phase microstructure in uniform flows will be summarized [18–20]. These observations will be connected with some characteristic features of the equations governing the perturbed liquid velocity field and the pair density distribution, and derived in the framework of the hybrid model. For finite particulate Reynolds numbers, estimates of the agitation tensors will be shown to be feasible using numerical simulations of two-body interactions [17]. Finally, the relevance of local closures for the induced-agitation for predicting phase distributions in confined systems will be debated, and the corresponding modeling issues will be underlined.


Author(s):  
Tomio Okawa ◽  
Kazuhiro Torimoto ◽  
Masanori Nishiura ◽  
Isao Kataoka ◽  
Michitsugu Mori

Experiments were conducted to clarify the single bubble rise characteristics in turbulent flows in vertical flow channels. It was revealed that the rise velocity of a bubble relative to the time-averaged local liquid velocity could be much smaller in turbulent upflows than in stagnant liquids. The reduction of relative velocity was estimated to be caused by the two factors: turbulence in continuous phase and steep velocity gradient near wall; new correlations for describing these two effects were proposed. The relative velocity between the phases significantly affects the lateral phase distribution in multidimensional simulation of bubbly two-phase flow and the present correlations can give reasonable predictions for the relative velocity in turbulent flow. It is hence expected that the new correlations can contribute to the further improvement of the simulation models of bubbly two-phase flows.


Microscopy ◽  
2019 ◽  
Vol 68 (4) ◽  
pp. 301-309 ◽  
Author(s):  
Xiaoke Mu ◽  
Andrey Mazilkin ◽  
Christian Sprau ◽  
Alexander Colsmann ◽  
Christian Kübel

Abstract Imaging the phase distribution of amorphous or partially crystalline organic materials at the nanoscale and analyzing the local atomic structure of individual phases has been a long-time challenge. We propose a new approach for imaging the phase distribution and for analyzing the local structure of organic materials based on scanning transmission electron diffraction (4D-STEM) pair distribution function analysis (PDF). We show that electron diffraction based PDF analysis can be used to characterize the short- and medium-range order in aperiodically packed organic molecules. Moreover, we show that 4D-STEM-PDF does not only provide local structural information with a resolution of a few nanometers, but can also be used to image the phase distribution of organic composites. The distinct and thickness independent contrast of the phase image is generated by utilizing the structural difference between the different types of molecules and taking advantage of the dose efficiency due to use of the full scattering signal. Therefore, this approach is particularly interesting for imaging unstained organic or polymer composites without distinct valence states for electron energy loss spectroscopy. We explore the possibilities of this new approach using [6,6]-phenyl-C61- butyric acid methyl ester (PC61BM) and poly(3-hexylthiophene-2,5-diyl) (P3HT) as the archetypical and best-investigated semiconductor blend used in organic solar cells, compare our phase distribution with virtual dark-field analysis and validate our approach by electron energy loss spectroscopy.


Author(s):  
T-C Kuo ◽  
A-S Yang ◽  
C-C Chieng

The coupled Eulerian-Lagrangian approach was used to study the effects of bubble size and high-pressure transport behaviour on the phase distribution mechanisms in vertically upward air-water two-phase bubbly flows. The approach solves the conservation equations of liquid phase in Eulerian space and equations of motion in conjunction with the random walk method for dispersed air bubbles in Lagrangian space. Numerical calculations were performed under conditions of three bubble diameters (2.8, 4.0 and 5.0 mm) and two different pressure levels (0.1 and 7.17 MPa) to explore the flow and void fraction development phenomena. Simulation results indicate the tendency of higher slip ratios and the movement of the void fraction peak towards the flow core for larger gas bubbles. In the pressure range 0.1-7.17 MPa, predictions reveal that the effect of high-pressure transport behaviour on the phase distribution is insignificant.


Author(s):  
G. Tryggvason ◽  
J. Lu ◽  
S. Biswas

Recent DNS studies of bubbly flows in channels are discussed. Simulations of nearly spherical bubbly flows in vertical channels show that the bubbles move towards the wall for upflow and away from the wall for downflow in such a way that the core is in hydrostatic equilibrium. For down flow the wall layer is free of bubbles but for upflow there is an excess of bubbles in the wall layer. The liquid velocity in the core is uniform. For laminar downflow the velocity in the wall layer can be computed analytically but for upflow the velocity is strongly influenced by the presence of the bubbles. Results for turbulent flow show similar behavior and for downflow the velocity is given (almost) by the law of the wall. Several simulations are used to examine the effect of void fraction and bubble size for turbulent downflow.


Sign in / Sign up

Export Citation Format

Share Document