Bubble Motion Research Articles

Investigation of surfactant-laden bubble migration dynamics in self-rewetting fluids using lattice Boltzmann method

Self-rewetting fluids (SRFs), such as aqueous solutions of long-chain alcohols, show anomalous nonlinear (quadratic) variations of surface tension with temperature involving a positive gradient in certain ranges, leading to different thermocapillary convection compared to normal fluids (NFs). They have recently been used for enhancing thermal transport, especially in microfluidics and microgravity applications. Moreover, surface-active materials or surfactants can significantly alter interfacial dynamics by their adsorption on fluid interfaces. The coupled effects of temperature- and surfactant-induced Marangoni stresses, which arise due to surface tension gradients, on migration bubbles in SRFs remain unexplored. We use a robust lattice Boltzmann method based on central moments to simulate the two-fluid motions, capture interfaces, and compute the transport of energy and surfactant concentration fields, and systematically study the surfactant-laden bubble dynamics in SRFs. When compared to motion of bubbles in NFs, in which they continuously migrate without a stationary behavior, our results show that they exhibit dramatically different characteristics in SRFs in many different ways. Not only is the bubble motion directed toward the minimum temperature location in SRFs, but, more importantly, the bubble attains an equilibrium position. In the absence of surfactants, such an equilibrium position arises at the minimum reference temperature occurring at the center of the domain. The addition of surfactants moves the equilibrium location further upstream, which is controlled by the magnitude of the Gibbs elasticity parameter that determines the magnitude of the surface tension variation with surfactant concentration. The parabolic dependence of surface tension in SRF is parameterized by a quadratic sensitivity coefficient, which modulates this behavior. The lower this quantity, the greater is the role of surfactants modifying the equilibrium position of the bubble in SRF. Furthermore, the streamwise gradient in the surfactant concentration field influences the transient characteristics in approaching the terminal state of the bubble. These findings provide new means to potentially manipulate the bubble dynamics, and especially to tune its equilibrium states, in microchannels and other applications by exploiting the interplay between surfactants and SRFs.

Directional enhancement of underwater impact and bubble loads in neighbor two-phase fluid domains charge

The loading characteristics of underwater explosions and the dynamic behavior of bubbles are directly related to the charge structure. This study proposes a unique charge structure in a gas–liquid two-phase fluid domain. Numerical methods are used to investigate the effects of fluid layer thickness and gas–liquid ratio on underwater explosion shock wave load, bubble dynamics, and bubble pulsation load. The results show that the two-phase fluid layer significantly enhances the directional release of shock wave energy and bubble pulsation load. During the shock wave phase, a lagging wave effect appears in the liquid layer direction, causing a secondary high-energy shock, significantly increasing the specific impulse. The gas layer direction may form a pressure relief channel effect, enhancing the shock wave peak pressure. For the bubble motion phase, differences in the physical properties of the fluid layer medium lead to irregular bubble boundary movements, promoting bubble tearing and rupture. The gaseous medium converts the accumulated shock wave energy into the internal energy of the bubble, increasing its volume potential. Although this characteristic reduces the pulsation frequency, it significantly increases the specific impulse. Altering the fluid layer medium can control explosion loads and bubble movement, offering new insights for ocean engineering applications.

Dynamic behavior and transfer characteristics of CO2-enriched bubbles within Spirulina sp. suspension under various aeration conditions using the high-speed imaging technique

The biological CO2 fixation method through microalgae photosynthesis has received considerable attention to alleviate the trend of global warming. CO2-enriched gas is generally aerated into the microalgae suspension in the form of bubbles through the gas distributors. Dynamic behavior and transfer characteristics of CO2-enriched bubbles are crucial to microalgae cells growth and CO2 bio-fixation. A visual experimental system based on the high-speed camera was constructed in this work to obtain the dynamic behavior and transfer characteristics of CO2-enriched bubbles within Spirulina sp. suspension. CO2-enriched bubbles movement and dissolution characteristics were comprehensively investigated under various CO2 concentrations, gas distributor aperture size, aeration rates, and Spirulina sp. biomass densities. Experimental results indicate that the optimal CO2 dissolution mass transfer and absorption rate were attained under the CO2 concentration of 5 %, gas distributor aperture diameter of 10 μm, and aeration rate of 0.1–0.3 vvm. Moreover, as Spirulina sp. biomass density increased, the bubble average diameter decreased, and rising velocity slowed while the volumetric mass transfer coefficient and CO2 absorption rate elevated. To summarize, this work may guide future efforts to enhance the photobioreactors (PBRs) performance from the perspective of aeration conditions optimization.

Development of integrated noise analysis systems for ship propellers considering hull wake

Recently, the requirements for underwater radiated noise (URN) reductions in ships have become significant issues, driven by concerns over the impact of naval vessel performance and noise on marine ecosystems. The majority of URN generated from ships is attributed to the propeller. This study developed the integrated noise analysis systems that predict and analyze all the significant noises generated from the propellers, considering the hull wake. Since each type of noise has different sources and mechanisms, it is essential to analyze all types of noise to evaluate the propeller noise performances comprehensively. The model was developed for cavitation noise by compression and expansion motion of bubbles to evaluate the radiated noise. For non-cavitation noise, a noise model on the frequency domain was derived to cover the noise generated by changes in the wake. The unsteady pressure on the rotating blade was obtained using flow fields and wall pressure spectrum, extending to high frequencies. The developed system provides a valuable tool for propeller optimization designs by minimizing URN while considering the shape of the hull and appendages.

3D Numerical Investigation of Bubble Upflow Condensation Behaviors During Subcooled Flow Boiling in Mini-Channel with VOSET

In the subcooled flow boiling process, bubble condensation is an inevitable basic phenomenon. This paper studies the condensation phenomenon of the single, double vertical/horizontal saturated bubbles rising in a three-dimensional mini-rectangular channel based on the interface capture method VOSET (coupled volume-of-fluid and level set method) and the phase transition model. Bubble condensation behaviors are investigated at different initial diameters, inlet velocity distributions, subcooling temperatures, bubble gaps, and arrangement for the two-bubble condensing system especially. The effects of these parametric on bubble motion trajectory, shape evolution, volume variation, and condensation rate are presented. The numerical results indicated that the initial bubble size and liquid subcooling play an important role in influencing the shape and volume variation of condensing bubble behaviors significantly, while the inlet velocity distribution only affects bubble motion trajectory. Furthermore, the interaction and coalescence between the bubbles will affect the bubble behaviors and the condensation rate. Finally, the condensation heat transfer coefficients at the bubble surfaces for different cases simulated in this paper are presented, seemingly first in the literature.

Gas-rigid-flexible compound blade coupling enhanced experimental study on chaotic mixing of multiphase flow

Efficient fluid mixing is essential for process intensification. This study proposes a new method in which gas-rigid-flexible composite blades are coupled to enhance chaotic mixing in multiphase flow systems. The rigidity and flexibility of the blades were adjusted by intermittent gas injection, which increased the effectiveness of mixing of the liquid-liquid two-phase fluid. This study investigates the influence of different process parameters on the mixing efficiency and quantifies the chaotic characteristics of fluid mixing through pressure-time series analysis of multiscale entropy and the 0–1 test. A high-speed camera recorded the bubble movement in the flow field, while particle image velocimetry (PIV) revealed the enhancement of the properties of the flow field in the system due to the suspended motion of the particles. Using suitable process parameters, gas-rigid-flexible composite blade coupling significantly enhanced the mixing effect, where the mixing time of the G-RFCP system was reduced by 1.42 times compared to that of the CP system. Bubble motion, deformation, and rupture enhanced the mechanical agitation, increasing the intensity of the turbulence and chaotic behaviour. Flow-field analysis indicated a three-fold increase in the vorticity and a 1.04-fold increase in the velocity difference for the G-RFCP system compared with those of the CP system. This study provides theoretical and experimental foundations for understanding chaotic mixing in liquid-liquid two-phase fluids.

Localization–delocalization transition for light particles in turbulence

Small bubbles in fluids rise to the surface due to Archimede's force. Remarkably, in turbulent flows this process is severely hindered by the presence of vortex filaments, which act as moving potential wells, dynamically trapping light particles and bubbles. Quantifying the statistical weights and roles of vortex filaments in turbulence is, however, still an outstanding experimental and computational challenge due to their small scale, fast chaotic motion, and transient nature. Here we show that, under the influence of a modulated oscillatory forcing, the collective bubble behavior switches from a dynamically localized to a delocalized state. Additionally, we find that by varying the forcing frequency and amplitude, a remarkable resonant phenomenon between light particles and small-scale vortex filaments emerges, likening particle behavior to a forced damped oscillator. We discuss how these externally actuated bubbles can be used as a type of microscopic probe to investigate the space-time statistical properties of the smallest turbulence scales, allowing to quantitatively measure physical characteristics of vortex filaments. We develop a superposition model that is in excellent agreement with the simulation data of the particle dynamics which reveals the fraction of localized/delocalized particles as well as characteristics of the potential landscape induced by vortices in turbulence. Our approach paves the way for innovative ways to accurately measure turbulent properties and to the possibility to control light particles and bubble motions in turbulence with potential applications to oceanography, medical imaging, drug/gene delivery, chemical reactions, wastewater treatment, and industrial mixing.

Numerical Investigation of Ultrasound-Induced Bubble Motion and Coalescence

Abstract Ultrasound-induced bubble motion and coalescence can be a unique approach to neutralize cavitation clouds in many biomedical applications. However, the detailed mechanisms of bubbles’ collective behaviours have not been fully understood. This work aims to develop an Eulerian-Lagrangian computational framework to simulate the bubble motion and coalescence driven by an ultrasound field. Specifically, the bubbles are modeled as oscillating spheres governed by the Keller-Miksis equation and are tracked as Lagrangian particles moving through the background Eulerian mesh. Several forces that drive bubble dynamics are modeled, including primary and secondary Bjerknes forces, drag force. As two bubbles contacting with each other, coalescence may happen and the film drainage model is applied to simulate the process. Using the new framework, the effects of pulse amplitude, frequency and length on the motion and coalescence of bubbles with different sizes and void fraction are investigated. The new knowledge can be applied to identify the optimal pulse parameters to better control the bubble motion and coalescence, improving the effectiveness of relevant ultrasound techniques.

Study on high-order-mode oscillation and migration characteristics of single bubble in the ultrasonic field

Abstract Investigating the motion characteristics of cavitation bubbles in an ultrasonic field is important for industries such as ultrasonic cleaning and ultrasonic therapy. This paper investigated the oscillation and migration characteristics of a single cavitation bubble in a 20 kHz ultrasonic field. Peristaltic pump syringe injection tests were performed at different locations in a 20 kHz ultrasonic field, and a high-speed camera was used to observe the bubble motion characteristics. The bubble was found to resonate with the ultrasonic waves at a higher order mode, changing from spherical shape to honeycomb and golf ball shapes. Also, we found that the ultrasonic waves may prevent the buoyancy of the bubble. Through revealing the motion characteristics of the bubble in the ultrasonic field, we may control the bubble motion by ultrasonic, and strengthen the bubble cavitation effect.

Investigation of the unsteady surface pressure field under a Mach 2 compression-ramp shock/boundary-layer interaction

The shock/boundary-layer interaction induced by a 20° compression ramp in a Mach 2 flow was investigated using fast pressure-sensitive paint with a bandwidth of about 10 kHz. The mean separated flow length-scale is about two upstream boundary layer thicknesses, which indicates the interaction is weak. The primary analysis consists of cross-correlations, coherence, and time-domain filtering. Two different frequency bands were investigated: low-frequency (f < 2000 Hz; StL < 0.1) and mid-frequency (0.1 < StL < 0.26). The low-frequency band time sequences and coherence reveal the shock-foot motion is mainly correlated with the reattachment region, which is indicative of the well-established breathing motion of the separation bubble. The breathing motion is observed to occur locally and globally (spanwise-averaged). Furthermore, in the low-frequency band, fluctuations in the upstream boundary layer are moderately correlated with the reattachment region fluctuations, but show no correlation with the intermittent region fluctuations. In the mid-frequency band, the intermittent region, separation bubble and reattachment region all exhibit significant correlation with the upstream boundary layer fluctuations, with the upstream fluctuations leading. The time-sequences in this frequency band reveal broad regions of pressure fluctuations that sweep through the interaction and affect the entire interaction. There is no known turbulent source for such large-scale fluctuations and they are believed to be due to a wind tunnel phenomenon. It is concluded that the dominant low-frequency breathing motion follows an oscillator model, but there remain significant correlations to upstream fluctuations that are not tied to the dominant breathing motion and seem to follow an amplifier model.

An investigation on hydrodynamic behavior of horizontal gas–liquid intermittent flow by S-PLIF and parallel-wire conductance sensors

Horizontal gas–liquid two-phase flows widely exist in a variety of industrial fields such as the environment, chemical industry and energy. Intermittent flow, characterized by the quasi-periodic alternation motion of liquid slugs and large-scale Taylor bubbles, is the most prevalent flow pattern in horizontal gas–liquid two-phase flows. Comprehensive research on the hydrodynamic behavior of gas–liquid intermittent flow is essential for revealing the mass and heat transfer characteristics between phases, thereby establishing physical models of flow parameters. In this study, a structured planar laser-induced fluorescence (S-PLIF) system and paired parallel-wire conductance sensors (PWCSs) are designed to detect the gas–liquid interface structures and the hydrodynamic parameters. The S-PLIF system consists of a laser, a visualization box, a Ronchi ruling-plate and a pair of high-speed cameras. An optical linear prism forms a laser sheet in the transverse direction and the middle of the test section, and the cameras are positioned at 70° and 90° respectively to the plane of the laser sheet, referred to as S-PLIF70 and S-PLIF90. Under the laser excitation, bright and dark stripes generated by the Ronchi ruling-plate are deflected at the gas–liquid interface. The influence of total reflection can be effectively avoided by identifying the deflection position of the stripes, allowing for accurate detection of the gas–liquid interfacial structures. The effects of the interface fluctuations and detached bubbles on the S-PLIF70&90 results are investigated. Fluctuation and evolution characteristics of the gas–liquid interface are analyzed by probability density function (PDF) and power spectral density (PSD). Meanwhile, hydrodynamic parameters, i.e., Taylor bubble nose and tail velocity, Taylor bubble and slug length, and slug frequency, are detected by the PWCSs, and compared with previous correlations. The evolution characteristics of the interfacial structures and hydrodynamic parameters are uncovered as the flow condition changes.

Migration characteristics of bubbles moving in gas–liquid countercurrent two-phase flow in an inclined pipe

While exploring and developing oil and gas, well kick and overflow accidents are inevitable. In such an accident, if the drilling bit is not at the bottom of the well, the drilling tool is blocked, or the reservoir contains toxic gases such as hydrogen sulfide, the traditional technique for well control is no longer suitable, and the bullheading method must be adopted to kill the well. However, during bullheading, the flow patterns of gas and liquid are intricate, making gas velocity hard to predict. To solve these problems, a set of visual simulation experimental device for directional well bullheading was built to find out the effects of wellbore inclination, liquid viscosity, gas and liquid velocities, and bubble size on bubble migration characteristics in gas–liquid countercurrent (gas flowing in opposite direction to the liquid). Prediction models of bubble distribution coefficient and rising velocity considering liquid viscosity, bubble size, and wellbore inclination angle were worked out. The experimental results show that small bubbles are mainly located in the center of the wellbore and large bubbles tend to approach the wellbore wall in gas–liquid countercurrent. Increasing wellbore inclination, viscosity, and the velocity of the liquid flowing against the bubbles results in a gradual decrease in Taylor bubble migration speed, which promotes the bubbles' pressing-back, as inhibiting the upward movement of large bubbles is essential for effective bullheading operations. The prediction models of distribution coefficient and bubble migration velocity have an error of less than 10% and 9. 78%, respectively.

Theoretical study on the movements of bubbles

The radial and translational motions of multiple interacting spherical bubbles are obtained using classical Newton mechanics. It is seen that bubbles not only move in straight line, but also in circular motion. The tracks of the bubbles show that the interactions among them include attractive, repulsive and dynamic equilibrium. There are three types of straight line corresponding to attraction, coexistence of attraction and repulsion and dynamic equilibrium, and two types of circular movement corresponding to attraction and dynamic equilibrium. The results can provide an explanation for cavitation chain and profile in cavitation field.

Research of Bubble Breakup and Influencing Factors in Flotation Machine

As the core factor of flotation, bubbles have an important effect on flotation. The rotor speed, initial gas parameters and impeller structure of the flotation machine will affect the formation and movement of bubbles. It is very important to study the influence mechanism of operating parameters of flotation machine on the bubble breakup process for flotation machine design and structure optimization. This paper takes KYF-0.2m<sup>3</sup> flotation machine as the research object, establishes a single bubble analysis model, adopts the VOF (Volume of Fluid) method to analyze the influence of different initial positions of bubbles on the bubble breakup behavior, and studies the influence of impeller speed and initial position of bubbles on the bubble breakup. Result show that the breakup of bubbles mainly occurs near the stator region. With the increase of rotational speed of the impeller, the centrifugal force and the disturbance of the convection field will become greater, the time of the bubble breakup become shorter, more bubbles breakup and generate more smaller ones. With the bubble position is closer to the rotating axis of the impeller, the impact of reflow becomes stronger and the bubble breakup effect will be better, and if the bubble initial position closer to the impeller cover, the influence of impeller on the bubbles become greater and the distribution of bubbles will be more uniform.

Investigation on the oblique water entry of the flared cavity

An experimental study of the oblique water entry of a flared cavity is reported in this paper. A numerical model is established that takes gas compressibility into account. This model effectively captures the additional flow phenomena that occur due to gas expansion and contraction within the cavity. There is a significant fluctuation phenomenon during the oblique water entry of the flared cavity, which is reflected not only in the movement of the water entry bubble but also in the slamming pressure. This fluctuation phenomenon is caused by the bottom cavity of the flared cavity, which adds the additional processes of external fluid inflow and internal fluid outflow compared with the closed-bottom structure, leading to additional flow due to conventional flow separation. When the external liquid flows into the cavity, the additional flow reduces the velocity potential of the fluid near the flow separation point, causing the bubble diameter near the flow separation point to contract. At the same time, the external liquid inflow causes the gas in the cavity to compress and the slamming pressure to increase. When the liquid flows out of the cavity, the additional flow increases the velocity potential of the fluid near the flow separation point, causing the bubble diameter near the flow separation point to expand. Meanwhile, the internal fluid outflow causes the gas in the cavity to expand and the slamming pressure to decrease.

A Navier-slip model for liquid film thickness in gas–liquid slug flow in capillaries

ABSTRACT The motion of a long gas bubble in a continuous liquid phase in a small diameter channel results in a thin liquid film between the bubble and the channel wall. Generally, a free-slip (zero shear) boundary condition is applied at the gas–liquid interphase boundary which is applicable only at low gas velocities. However, at increased flow rates of gas phase, shear stress at the gas–liquid interface is non-negligible. Applying non-zero shear at the interface requires knowledge of velocity profile in the gas phase. In this work, we circumvent this problem by using a slip length like parameter to take into account the velocity gradient in the gas phase. The slip length is infinite for the free-slip condition at the interface and zero for the no-slip condition. An expression is derived for the film thickness in terms of slip length and capillary number and is compared with the experimental data and other correlations available in the literature.

Experimental and numerical study of flow boiling in a vertical upward narrow rectangular channel under low heat flux

This article experimentally and numerically investigates the flow boiling in a vertical upward narrow rectangular channel under low heat flux. VOF model for interface tracking and the phase-change Lee model are coupled to build a numerical model of a channel with a scale of 1200 mm × 250 mm × 2.75 mm. The temperature, velocity and phase distribution in the channel are obtained at different heights The bubble movement is analyzed and the two-phase flow regimes are classified. Furthermore, the local characteristic of flow boiling heat transfer along the channel and the effect of the bubble regimes on the heat transfer are studied. The results indicate that saturated nucleate boiling is the dominant boiling form under low heat flux.

Numerical Investigation on Two-Phase Flow-Electrochemical Reaction Coupling in PEM Water Electrolyzers Porous Transport Layer Under High Operating Pressures

Proton exchange membrane water electrolyzers (PEMWEs) operated at high pressures have advantages such as lower energy consumption, reduced noise, and absence of component slippage. They have a more robust potential for commercial applications than other electrolyzers. The transport characteristics of two-phase flow within the porous transport layer (PTL) and channels are not well understood, especially under high-pressure conditions. The gas bubbles generated in the catalyst layer (CL) pass through the PTL and reach the flow channel. However, oxygen evolution reaction in the anode CL may be limited if these bubbles block the pathway of liquid water, causing reactant starvation. The interaction between electrochemical reactions and gas bubbles thus directly impacts a PEMWE's efficiency. This study aims to investigate the coupling between the two-phase flow and electrochemical reactions under high-pressure conditions.A two-dimensional (2D) CFD model was developed to simulate the bubbly flow through a reconstructed PTL. The volume of fluid (VOF) method was employed to capture bubble growth, movement, and detachment, as illustrated in Figure 1. The coupling of bubble transport and electrochemical reactions was achieved through the application of a weighted averaging method. A cross-section of a stochastically reconstructed titanium felt model was generated as the computational domain; see Fig. 1(a). The computation domain consists of a liquid water channel, a PTL, and a simplified CL interface at the bottom. The CL's reaction rate is modeled as a function of local bubble coverage on the CL interface. Transient simulations of the two-phase flow of the PEMWE anode side were carried out for pressures ranging from 1 bar to 200 bar. The simulation results elucidate the flow characteristics of liquid water in the PTL under high-pressure conditions and quantitatively reveal the coupling between two-phase flow and electrochemical reactions. The behavior of the two-phase flow in PEMWE is found to be sensitive to operating pressure. As the operating pressure increases, the size of bubbles generated on the CL interface decreases, resulting in a higher permeability of liquid water than that under low pressures. When the flow rate of liquid water in the channel is sufficiently high, bubbles do not fully develop and are rapidly removed from the domain. This phenomenon significantly shifts the two-phase flow regime within the channels,. Furthermore, the activation overpotential is substantially reduced under high-pressure conditions. The present study demonstrates a numerical simulation tool for high-pressure PEMWE. A future study of a full 3D PTL model is currently underway. Figure 1

Analysis of vapor bubble diameter, departure frequency and dynamics in a single cavity

The ongoing quest for methods that enhance the efficiency of heat transfer processes, particularly those involving phase change, underscores the importance of comprehending the dynamics of vapor bubbles during boiling. This study investigates heat transfer mechanisms and the dynamics of vapor bubbles within the nucleate boiling regime. Experimental tests were conducted on a flat copper surface featuring a single cavity, focusing on analyzing vapor bubble growth and departure stages. The working fluid examined was HFE-7100 under saturated conditions. The formation and growth aspects of vapor bubbles, including their diameter (Dd) and departure frequency (f), were investigated through experimental data obtained by two different techniques: by an optical sensor of variable resistance capable of generating an analog signal from a voltage change and by a high-speed camera that captures images immediately after the instant that the bubble detached from the surface. Both techniques provide visual insights into the boiling phenomenon. An escalation in heat flux and, consequently, wall superheat resulted in an increased bubble departure frequency. Furthermore, the optical flow analysis successfully identified velocity and vorticity fields induced by micro convection, the prevailing heat transfer mode in nucleated boiling. A slight horizontal displacement trend over time was observed, which may be caused by the asymmetry in the velocity profile, confirmed by velocity peaks tending toward one direction. Vortex analysis reveals rotational motion concentrated in specific regions, with noticeable deformation near the bottom of the bubble and asymmetric movement contributing to vorticity. These findings provide insights into the vapor bubble dynamics, which is important for understanding the cavity rewetting phenomena.

Wimpled thin films via multiple motions of a bubble decorated with surface-active molecules

HypothesisThin liquid films play a crucial role in various systems and applications. Understanding the mechanisms that regulate their morphology is a scientific challenge with obvious implications for application optimization. Thin liquid films trapped between bubbles and air-liquid interface can show various configurations influenced by their deformation history and system characteristics. ExperimentsThe morphology of thin liquid films formed in the presence of surface-active molecules is here studied with interferometric techniques. Three different systems with varying interfacial properties are investigated to understand their influence on film morphology. Specific deformation histories are applied to the films to generate complex film structures. FindingsWe achieve the creation of a rather stable wimple by implementing controlled bubble motions against the air-liquid interface. We provide a criterion for wimple formation based on lubrication theory. The long-term stability of the wimple is also investigated, and more complex multi-wimple structures are experimentally produced building upon the achieved wimple stability.