scholarly journals Energy spectra in turbulent bubbly flows

2016 ◽  
Vol 791 ◽  
pp. 174-190 ◽  
Author(s):  
Vivek N. Prakash ◽  
J. Martínez Mercado ◽  
Leen van Wijngaarden ◽  
E. Mancilla ◽  
Y. Tagawa ◽  
...  

We conduct experiments in a turbulent bubbly flow to study the nature of the transition between the classical $-5/3$ energy spectrum scaling for a single-phase turbulent flow and the $-3$ scaling for a swarm of bubbles rising in a quiescent liquid and of bubble-dominated turbulence. The bubblance parameter (Lance & Bataille J. Fluid Mech., vol. 222, 1991, pp. 95–118; Rensen et al., J. Fluid Mech., vol. 538, 2005, pp. 153–187), which measures the ratio of the bubble-induced kinetic energy to the kinetic energy induced by the turbulent liquid fluctuations before bubble injection, is often used to characterise bubbly flow. We vary the bubblance parameter from $b=\infty$ (pseudoturbulence) to $b=0$ (single-phase flow) over 2–3 orders of magnitude (0.01–5) to study its effect on the turbulent energy spectrum and fluctuations in liquid velocity. The probability density functions (PDFs) of the fluctuations in liquid velocity show deviations from the Gaussian profile for $b>0$, i.e. when bubbles are present in the system. The PDFs are asymmetric with higher probability in the positive tails. The energy spectra are found to follow the $-3$ scaling at length scales smaller than the size of the bubbles for bubbly flows. This $-3$ spectrum scaling holds not only in the well-established case of pseudoturbulence, but surprisingly in all cases where bubbles are present in the system ($b>0$). Therefore, it is a generic feature of turbulent bubbly flows, and the bubblance parameter is probably not a suitable parameter to characterise the energy spectrum in bubbly turbulent flows. The physical reason is that the energy input by the bubbles passes over only to higher wavenumbers, and the energy production due to the bubbles can be directly balanced by the viscous dissipation in the bubble wakes as suggested by Lance & Bataille (1991). In addition, we provide an alternative explanation by balancing the energy production of the bubbles with viscous dissipation in the Fourier space.

1996 ◽  
Vol 118 (2) ◽  
pp. 377-382 ◽  
Author(s):  
V. E. Nakoryakov ◽  
O. N. Kashinsky ◽  
V. V. Randin ◽  
L. S. Timkin

Gas-liquid bubbly flow was investigated in vertical pipes for different flow conditions: fully developed turbulent downward flow in a 42.3 mm diameter pipe and upward flow in a 14.8 mm diameter pipe with liquid of elevated viscosity. Wall shear stress, local void fraction, and liquid velocity profiles, shear stress, and velocity fluctuations were measured using an electrodiffusional method. Results obtained demonstrate the existence of “universal” near-wall velocity distribution in a downward bubbly flow. The reduction of turbulent fluctuations is observed in downward flow as compared to a single-phase turbulent flow. The development of bubble-induced liquid velocity fluctuations in a “laminar” bubbly flow was studied.


2020 ◽  
Vol 10 (17) ◽  
pp. 6029 ◽  
Author(s):  
Hepeng Zheng ◽  
Yun Zhang ◽  
Yuan Wang ◽  
Lifeng Zhang ◽  
Jun Peng ◽  
...  

The intensification of Typhoon Lekima (2019) is simulated with the Weather Research and Forecasting model to study the atmospheric horizontal kinetic energy (HKE) spectra and corresponding spectral HKE budgets under the control of real tropical cyclone (TC). The results show that the TC has the ability to modify the canonical atmospheric energy spectrum during its evolution, which is dominated by its rotational mode. With the intensification of Lekima, the HKE spectrum in the troposphere swells over the central mesoscale and develops an arc-like shape. The stronger the TC, the more pronounced the arc-like shape is and the smaller scale it extends to. The roles various physical processes play at different heights and horizontal scales during the intensification of Lekima are investigated and the dependence of the effect of physical processes on scale and height is revealed. Meanwhile, the potential relationship between the intensification of TC, the activation of energy activity at smaller scales, and the downscale extension of the arc-like spectral shape is found.


2011 ◽  
Vol 133 (9) ◽  
Author(s):  
M. E. Shawkat ◽  
C. Y. Ching

The liquid turbulence kinetic energy transfer between the liquid and gas phases was investigated for upward air-water bubbly flow in a 200 mm diameter pipe. The liquid and gas axial momentum equations were analyzed to estimate the interfacial drag from experimental measurements, and hence the liquid turbulence production due to the relative velocity of the bubbles. The liquid turbulence production due to the bubbles was significantly higher than that due to the liquid shear. The liquid turbulence kinetic energy budget indicates that the turbulence production due to the bubbles is approximately balanced by the viscous dissipation, estimated assuming an isotropic turbulence structure, with negligible dissipation due to the bubbles. The liquid turbulence kinetic energy spectra showed an addition of energy at length scales in the range corresponding to the bubble diameter. A model for the turbulence energy production spectra due to the bubbles is proposed and used to investigate the spectral turbulence energy budget. The model indicates that when there is a liquid turbulence augmentation, most of the production occurs in the low wave number range with only a small overlap with the viscous dissipation region. In the case of a turbulence suppression, most of the bubble production occurs in the same wave number range as the viscous dissipation.


2008 ◽  
Vol 15 ◽  
pp. 17-22 ◽  
Author(s):  
M. V. Kurgansky

Abstract. In a forced-dissipative barotropic model of the atmosphere on a spherical planet, by following mathematical techniques in (Thompson, P. D.: The equilibrium energy spectrum of randomly forced two-dimensional turbulence, Journal of the Atmospheric Sciences, 30, 1593–1598, 1973) but applying them in a novel context of the discrete spectrum on a rotating sphere, the "minus 2" energy spectrum for wavenumbers much greater than a characteristic wavenumber of the baroclinic forcing has been obtained if the forcing is taken in the simplest and most fundamental form. Some observation-based atmospheric kinetic energy spectra, with their slopes lying between "minus 2" and "minus 3" laws, are discussed from the perspective of the deduced "minus 2" energy spectrum.


2001 ◽  
Vol 437 ◽  
pp. 203-228 ◽  
Author(s):  
ROBERT F. MUDDE ◽  
TAKAYUKI SAITO

The hydrodynamical similarities between the bubbly flow in a bubble column and in a pipe with vertical upward liquid flow are investigated. The system concerns air/water bubbly flow in a vertical cylinder of 14.9 cm inner diameter. Measurements of the radial distribution of the liquid velocity, gas fraction and the bubble velocity and size are performed using laser Doppler anemometry for the liquid velocity and a four-point optical fibre probe for the gas fraction, bubble velocity and size. The averaged gas fraction was 5.2% for the bubble column (with a superficial liquid velocity of zero) and 5.5% for the bubbly pipe flow at a superficial liquid velocity of 0.175 m s−1. From a hydrodynamical point of view, the two modes of operation are very similar. It is found that in many respects the bubbly pipe flow is the superposition of the flow in the bubble column mode and single-phase flow at the same superficial liquid velocity.The radial gas fraction profiles are the same and the velocity profiles differ only by a constant offset: the superficial liquid velocity. This means that the well-known large-scale liquid circulation (in a time-averaged sense) of the bubble column is also present in the bubbly pipe flow. For the turbulence intensities it is found that the bubbly pipe flow is like the superposition of the bubble column and the single-phase flow at the superficial liquid velocity of the pipe flow, the former being at least an order of magnitude higher than the latter. The large vortical structures that have been found in the bubble columns are also present in the bubbly pipe flow case, partly explaining the much higher ‘turbulence’ levels observed.


2012 ◽  
Vol 2012 ◽  
pp. 1-20 ◽  
Author(s):  
K. Ekambara ◽  
R. Sean Sanders ◽  
K. Nandakumar ◽  
J. H. Masliyah

Modelling of gas-liquid bubbly flows is achieved by coupling a population balance equation with the three-dimensional, two-fluid, hydrodynamic model. For gas-liquid bubbly flows, an average bubble number density transport equation has been incorporated in the CFD code CFX 5.7 to describe the temporal and spatial evolution of the gas bubbles population. The coalescence and breakage effects of the gas bubbles are modeled. The coalescence by the random collision driven by turbulence and wake entrainment is considered, while for bubble breakage, the impact of turbulent eddies is considered. Local spatial variations of the gas volume fraction, interfacial area concentration, Sauter mean bubble diameter, and liquid velocity are compared against experimental data in a horizontal pipe, covering a range of gas (0.25 to 1.34 m/s) and liquid (3.74 to 5.1 m/s) superficial velocities and average volume fractions (4% to 21%). The predicted local variations are in good agreement with the experimental measurements reported in the literature. Furthermore, the development of the flow pattern was examined at three different axial locations ofL/D= 25, 148, and 253. The first location is close to the entrance region where the flow is still developing, while the second and the third represent nearly fully developed bubbly flow patterns.


Author(s):  
Takamichi Hiroi ◽  
Tatsuya Hamada ◽  
Chiharu Kawakita

Abstract Friction drag and characteristic of bubbles in horizontal water channel are investigated at bulk liquid velocity Um = 1 ∼ 5 m/s (Reynolds number Rem = 16,000 ∼ 120,000 (based on the channel height)) and mean void fraction α = 0.5, 1, 2 %. Firstly, shear stress sensor is applied to investigate the relation between friction drag with bubbles and bulk liquid velocity. Friction drag in the bubbly flow is larger than it in the single-phase flow at Um = 1 ∼ 2 m/s. It in the bubbly flow, however, decreases with the mean liquid velocity. Furthermore, it in bubbly flow is smaller than it in the single-phase flow at Um ≥ 3 m/s. Secondly, fiber optic sensor is applied to investigate the void fraction distribution, bubble diameter distribution and streamwise velocity of bubbles. The peak value of the void fraction decreases with increasing of the bulk liquid velocity at Um = 2 ∼ 5 m/s. Bubbles exist in only y/δ = 0 ∼ 0.5 at Um = 1, 2 m/s. The mean bubble velocity increases with the bulk liquid velocity. The mean Sauter diameter decreases with increasing of the bulk liquid velocity. It appears that the high void ratio near the wall causes the increasing of drag and friction drag decreases when bubbles exist in whole of upper half of channel.


Tellus ◽  
1981 ◽  
Vol 33 (1) ◽  
pp. 102-104 ◽  
Author(s):  
Tsing-Chang Chen ◽  
Joseph J. Tribbia

Sign in / Sign up

Export Citation Format

Share Document