On the rise of an ellipsoidal bubble in water: oscillatory paths and liquid-induced velocity

2001 ◽  
Vol 440 ◽  
pp. 235-268 ◽  
Author(s):  
KJETIL ELLINGSEN ◽  
FRÉDÉRIC RISSO

This work is an experimental study of the rise of an air bubble in still water. For the bubble diameter considered, path oscillations develop in the absence of shape oscillations and the effect of surfactants is shown to be negligible. Both the three-dimensional motion of the bubble and the velocity induced in the liquid are investigated. After the initial acceleration stage, the bubble shape remains constant and similar to an oblate ellipsoid with its symmetry axis parallel to the bubble-centre velocity, and with constant velocity magnitude. The bubble motion combines path oscillations with slow trajectory displacements. (These displacements, which consist of horizontal drift and rotation about a vertical axis, are shown to have no influence on the oscillations). The bubble dynamics involve two unstable modes which have the same frequency and are π/2 out of phase. The primary mode develops first, leading to a plane zigzag trajectory. The secondary mode then grows, causing the trajectory to progressively change into a circular helix. Liquid-velocity measurements are taken up to 150 radii behind the bubble. The nature of the liquid flow field is analysed from systematic comparisons with potential theory and direct numerical simulations. The flow is potential in front of the bubble and a long wake develops behind. The wake structure is controlled by two mechanisms: the development of a quasi-steady wake that spreads around the non-rectilinear bubble trajectory; and the wake instability that generates unsteady vortices at the bubble rear. The velocities induced by the wake vortices are small compared to the bubble velocity and, except in the near wake, the flow is controlled by the quasi-steady wake.

2010 ◽  
Vol 650 ◽  
pp. 287-306 ◽  
Author(s):  
JULIÁN MARTÍNEZ MERCADO ◽  
DANIEL CHEHATA GÓMEZ ◽  
DENNIS VAN GILS ◽  
CHAO SUN ◽  
DETLEF LOHSE

Three-dimensional particle tracking velocimetry (PTV) and phase-sensitive constant temperature anemometry in pseudo-turbulence – i.e. flow solely driven by rising bubbles – were performed to investigate bubble clustering and to obtain the mean bubble rise velocity, distributions of bubble velocities and energy spectra at dilute gas concentrations (α ≤ 2.2 %). To characterize the clustering the pair correlation function G(r, θ) was calculated. The deformable bubbles with equivalent bubble diameter db = 4–5 mm were found to cluster within a radial distance of a few bubble radii with a preferred vertical orientation. This vertical alignment was present at both small and large scales. For small distances also some horizontal clustering was found. The large number of data points and the non-intrusiveness of PTV allowed well-converged probability density functions (PDFs) of the bubble velocity to be obtained. The PDFs had a non-Gaussian form for all velocity components and intermittency effects could be observed. The energy spectrum of the liquid velocity fluctuations decayed with a power law of −3.2, different from the ≈ −5/3 found for homogeneous isotropic turbulence, but close to the prediction −3 by Lance & Bataille (J. Fluid Mech., vol. 222, 1991, p. 95) for pseudo-turbulence.


Author(s):  
Deify Law ◽  
Thomas G. Shepard

The present work relates to the dynamics of single bubbles accelerating through a converging nozzle. There are two main aspects to this study. First, this expands upon a previously used two-dimensional model [1] by providing three-dimensional volume of fluid (VOF) simulations that show better agreement with experiments. The VOF model is employed to perform simulations using the commercial computational fluid dynamics (CFD) code ANSYS FLUENT. Second, the present work uses experimental high-speed camera results in conjunction with simulation results to demonstrate bubble time trace and velocity information. Time series of the average liquid velocity at the atomizer exit orifice when the bubble exits as determined via simulation are reported. The passing of a bubble through the nozzle is found to cause a significant fluctuation in the exit velocity that is coupled to the liquid and gas dynamics upstream of the exit.


Author(s):  
Kamran Siddiqui ◽  
Wajid A. Chishty

The paper reports on an experimental study conducted to investigate the dynamics of gas bubbles when injected from an orifice that is subjected to a liquid cross-flow. The experiments were conducted over a range of gas and liquid flow rates and at various orientations of the liquid channel. An image processing algorithm was developed for the detection and characterization of the bubbles in both temporal and spatial coordinates. The transient behaviour of bubble dynamics at different channel orientations under different liquid and gas flow conditions are presented in the paper. Results show that the equivalent bubble diameter in the vertical plane decreased with increase in time as well as with increase in the gas-to-liquid velocity ratio. The channel inclination has no appreciable effect on the bubble diameter. The streamwise bubble velocity showed significant transient fluctuations, which diminished with an increase in the channel inclination angle. The fluctuations of the bubble vertical velocity were found to be unaffected by the channel inclination angle. The magnitudes of both streamwise and vertical velocities were influenced by the channel inclination, however, the trends were found to be opposite.


2011 ◽  
Vol 133 (2) ◽  
Author(s):  
Deify Law ◽  
Samuel T. Jones ◽  
Theodore J. Heindel ◽  
Francine Battaglia

The external loop airlift reactor (ELALR) is a modified bubble column reactor that is composed of two vertical columns that are interconnected with two horizontal tubes and is often preferred over traditional bubble column reactors because they can operate over a wider range of conditions. In the present work, the gas-liquid flow dynamics in an ELALR was simulated using an Eulerian–Eulerian ensemble-averaging method in two-dimensional (2D) and three-dimensional (3D) coordinate systems. The computational fluid dynamics (CFD) simulations were compared to experimental measurements from a 10.2 cm diameter ELALR for superficial gas velocities ranging from 1 cm/s to 20 cm/s. The effect of specifying a mean bubble diameter to represent the gas phase in the CFD modeling was investigated, and 2D and 3D simulations were found to be in good agreement with the experimental data. The ELALR flow regimes were compared for the reactor operating in bubble column, closed vent, and open vent modes, and the 2D simulations qualitatively predicted the behavior of bubble growth in the downcomer. However, it was found that 3D simulations were necessary to capture the physics of the ELALR for gas holdup, bulk density differences, and riser superficial liquid velocity.


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.


2017 ◽  
Vol 140 (4) ◽  
Author(s):  
Aaron S. Alexander ◽  
Arvind Santhanakrishnan

Savonius vertical axis wind turbines (VAWTs) typically suffer from low efficiency due to detrimental drag production during one half of the rotational cycle. The present study examines a stator assembly created with the objective of trapping cylindrical flow for application in a Savonius VAWT. While stator assemblies have been studied in situ around Savonius rotors in the past, they have never been isolated from the rotor to determine the physics of the flow field, raising the likelihood that a moving rotor could cover up deficiencies attributable to the stator design. The flow field created by a stator assembly, sans rotor, is studied computationally using three-dimensional (3D) numerical simulations in the commercial computational fluid dynamics (CFD) package Star-CCM+. Examination of the velocity and pressure contours at the central stator plane shows that the maximum induced velocity exceeded the freestream velocity by 65%. However, flow is not sufficiently trapped in the stator assembly, with excess leakage occurring between the stator blades due to adverse pressure gradients and momentum loss from induced vorticity. A parametric study was conducted on the effect of the number of stator blades with simulations conducted with 6, 12, and 24 blades. Reducing the blade number resulted in a reduction in the cohesiveness of the internal swirling flow structure and increased the leakage of flow through the stator. Two unique energy loss mechanisms have been identified with both caused by adverse pressure gradients induced by the stator.


2006 ◽  
Vol 65 (6) ◽  
pp. 429-439 ◽  
Author(s):  
Keisuke Kushiro ◽  
Jun Maruta

Author(s):  
Ahmed M Nagib Elmekawy ◽  
Hassan A Hassan Saeed ◽  
Sadek Z Kassab

Three-dimensional CFD simulations are carried out to study the increase of power generated from Savonius vertical axis wind turbines by modifying the blade shape and blade angel of twist. Twisting angle of the classical blade are varied and several proposed novel blade shapes are introduced to enhance the performance of the wind turbine. CFD simulations have been performed using sliding mesh technique of ANSYS software. Four turbulence models; realizable k -[Formula: see text], standard k - [Formula: see text], SST transition and SST k -[Formula: see text] are utilized in the simulations. The blade twisting angle has been modified for the proposed dimensions and wind speed. The introduced novel blade increased the power generated compared to the classical shapes. The two proposed novel blades achieved better power coefficients. One of the proposed models achieved an increase of 31% and the other one achieved 32.2% when compared to the classical rotor shape. The optimum twist angel for the two proposed models achieved 5.66% and 5.69% when compared with zero angle of twist.


Flow ◽  
2021 ◽  
Vol 1 ◽  
Author(s):  
Pablo Ouro ◽  
Maxime Lazennec

Graphical Abstract


Energies ◽  
2020 ◽  
Vol 14 (1) ◽  
pp. 49
Author(s):  
Zheng Yuan ◽  
Jin Jiang ◽  
Jun Zang ◽  
Qihu Sheng ◽  
Ke Sun ◽  
...  

In the array design of the vertical axis wind turbines (VAWT), the wake effect of the upstream VAWT on the downstream VAWT needs to be considered. In order to simulate the velocity distribution of a VAWT wake rapidly, a new two-dimensional numerical method is proposed, which can make the array design easier and faster. In this new approach, the finite vortex method and vortex particle method are combined to simulate the generation and evolution of the vortex, respectively, the fast multipole method (FMM) is used to accelerate the calculation. Based on a characteristic of the VAWT wake, that is, the velocity distribution can be fitted into a power-law function, a new correction model is introduced to correct the three-dimensional effect of the VAWT wake. Finally, the simulation results can be approximated to the published experimental results in the first-order. As a new numerical method to simulate the complex VAWT wake, this paper proves the feasibility of the method and makes a preliminary validation. This method is not used to simulate the complex three-dimensional turbulent evolution but to simulate the velocity distribution quickly and relatively accurately, which meets the requirement for rapid simulation in the preliminary array design.


Sign in / Sign up

Export Citation Format

Share Document