Bubbly Flow Research Articles

An entropy measure-based study on flow pattern of gas–liquid two-phase flow in a U-Tube

U-tubes are widely applied in gas–liquid two-phase transportation in chemical engineering. The diverse flow patterns within these tubes significantly affect the pressure loss, heat transfer efficiency, and even the fluid-induced vibration amplitude of the tubes. This study explores the complex flow pattern features in a U-tube in a vertical plane and focuses on recognizing them. For the acquisition and classification of flow patterns, a Computational Fluid Dynamics (CFD) model for gas–liquid two-phase flow is first established, and its quantitative calculation error is ensured to be less than 5%. Then, the spatiotemporal evolution characteristics of flow patterns is analyzed. The real-time pressure drop response is chosen as the representation signal, and its nonlinear features in the time and frequency domain under different flow patterns are explored. A nonlinear time series is constructed by extracting a segment from the real-time pressure drop data, and six entropy measures are applied to analyze and identify them. Finally, the sensitivity of entropy measures to both the time series lengths and the tested sections are evaluated. Results show that there are six typical flow patterns in a U-tube. According to most entropy measures, the bubble flow has the highest complexity; however, the plug flow presents the lowest complexity. In the U-bend, pressure drop signals for the bubble and annular flows show random fluctuations within a specific range, in contrast to the marked periodicity in plug flow signals, while wavy and slug flows exhibit intermittent peak values. Including the upstream and downstream straight pipes in the analysis, rather than focusing solely on the U-bend, significantly increases the complexity of the stratified, plug, and slug flows. Fuzzy entropy is an effective tool for identifying the six flow patterns, demonstrating good resilience to variations in the length of the data series. This characteristic makes it highly useful for real-time identification of flow patterns in the U-bend sections of non-transparent U-tubes, offering considerable potential in chemical equipment.

Effect of solid loading and particle size on bubble behavior and flow field structure in slurry bubble column

The complexity of fluid dynamics in a slurry bubble column reactor introduces significant uncertainty in reactor design and scale-up. This paper investigates the hydrodynamic performance of the gas–liquid–solid system within the reactor by employing computational fluid dynamics-population balance modeling numerical simulations alongside particle image velocimetry (PIV) experiments. The effect of superficial gas velocity and particle conditions on the overall gas holdup were analyzed, focusing on the effects of particle size and solid concentration on bubble size, bubble behavior, flow field structure, and local gas holdup distribution at high superficial gas velocities. Bubble size was evaluated using calibrated image measurements, and the impact of varying solid conditions was thoroughly explored. The results revealed that an increase in solid size correlated with higher gas holdup and smaller bubble sizes, whereas a greater solid concentration resulted in decreased gas holdup and larger bubble sizes. PIV experiments indicated that bubbles exhibited a tendency to migrate toward the central region of the reactor, leading to the formation of larger bubbles that accelerated the rise of surrounding bubbles, while smaller bubbles near the wall moved downward. As the slurry bed height increased, the range of local gas holdup distribution expanded, resulting in a symmetrical distribution of radial local gas holdup in the fully developed stage at a height of 0.16 m.

An improved computational fluid dynamics-discrete bubble model for multiphase bubble flow simulation with a wide grid/bubble size ratio

An improved computational fluid dynamics-discrete bubble model for multiphase bubble flow simulation with a wide grid/bubble size ratio

Bubble boundary R-CNN: A multitask model for segmenting and reconstructing overlapping bubbles

When measuring the dynamic characteristics of bubbles using image-based methods, the images including complex overlapping of bubbles with irregular shape pose a huge challenge to current bubble reconstruction algorithms. Based on the Mask R-CNN architecture, a multitask model “Bubble Boundary R-CNN” is proposed for detection, segmentation and shape reconstruction of overlapping bubbles with high deformation at high gas velocities. By leveraging shared deep features from detection and segmentation, model achieves flexible bubble reconstruction and accurate size distribution extraction from high gas velocity bubble flows. To overcome the limitations of Mask R-CNN and address the challenges posed by high-deformation and overlapping bubbles, a semantic segmentation auxiliary head and the Boundary-preserving Mask Head (BMask Head) are introduced to enhance edge information and extract bubble boundary features. Additionally, an Overlapping Boundary Prediction Network (OBPN) is proposed to reconstruct boundary features that are lost due to bubble occlusion and predict the complete bubble shapes. The optimized model shows its ability to flexibly and accurately reconstruct overlapped and deformed bubbles by varying occlusion rates and circularity. The new model has also proven its reliability when applied to bubble images from a bubble column at gas velocities of up to 0.067 m3/s. In future work, more efforts are needed to further enhance the detection capability of the model and reduce missed and multiple detections at high gas velocity.

Experimental study on pressure drop fluctuations for boiling flow at various gravity levels

Experiments were conducted on the pressure drop fluctuations of R245fa boiling through a horizontal microchannel with a diameter of 0.94 mm and a heated length of 170 mm at various gravity levels. The visualization results show that under inlet saturation conditions, three typical periodic oscillation modes are observed: continuous two-phase flow, bubbly/annular alternating flow, and churn/annular alternating flow. The flow oscillation characteristics within the channel are closely related to bubble dynamics, the vapor–liquid interface evolution, and the upstream compressible vapor volume. It has also been found that increasing the mass flux or decreasing the heat flux can enhance the suppression effect on pressure drop oscillation. In addition, under inlet subcooling conditions, although the total pressure drop increases with the increase in gravity level (1–2.78 g), and the amplitude of pressure drop fluctuation first increases and then decreases with the increase in gravity level, with the transition point occurring approximately at an = 1.10 g. The impact of gravity level on the amplitude of pressure drop fluctuation is primarily determined by the relative dominance of bubbly flow and annular flow within the channel.

Flow boiling heat transfer in multi-stage enhanced open microchannels with micro/nano structures

In the realm of flow boiling research, open microchannels heat sink is garnering more attention due to its potential in enhancing heat transfer and mitigating flow instability. Nevertheless, growing heat dissipation needs necessitate further enhancement of the heat removing capability of flow boiling within open microchannels. In this paper, experimental research on flow boiling in open microchannels utilized a new environmentally friendly refrigerant SF-33 and presented a novel heat transfer enhancement strategy, specifically, a multi-region enhancement strategy for diverse flow pattern stages. A copper heat sink with 25 open microchannels was fabricated, and a layer of copper powder particles was sintered on the channel bottom to create a microporous structured surface in about 1/3 of the inlet region. Several layers of copper wire mesh with micropores were sintered in about 1/3 of the outlet region, and the nanostructures were grown on the wire mesh by chemical treatment, resulting in a micro-nanocomposite structure. The final product formed the multi-stage enhanced open microchannels (MSEOM). The heat transfer performance and pressure drop of SF-33 saturated flow boiling in MSEOM under different operating conditions were determined, and the transitions of flow patterns were visualized to analyze the mechanism of two-phase heat transfer. The study demonstrated that despite SF-33 having a low latent heat of vaporization, its flow boiling heat transfer coefficient in MSEOM was high. This was attributed to the enhanced effect of multi-region surface modifications on heat transfer. Additionally, a vapor plug breakup-dispersed bubble flow pattern was observed for the first time, which resulted from the combined influence of the low surface tension property of SF-33 and the boosted nucleate boiling triggered by the micro-nanostructures in the outlet region.

Bubble deformation characteristics and the effect on particle removal efficiency in wet scrubbers

Single-bubble formation and rise are fundamental to the prediction of decontamination factors (DF) under the complex hydrodynamics of bubble flow and the gas–liquid flow pattern. However, bubbles in the water for particle control are not always perfectly spherical or ellipsoidal. To date, studies of bubble washing for particle removal have mainly focused on models with simple spherical and ellipsoidal shapes or the impact of operating parameters on removal efficiency. No attempt has been reported in the literature to determine the effect of bubble shape on particle removal efficiency. An investigation of the influence of bubble deformation during the bubble rising period was carried out using a high-speed camera and image analysis techniques to study the bubble shape. The experimental results show that the deformation of a bubble is directly correlated with its rising velocity, rising tail, equivalent radius and residence time. These parameters are generally estimated using ideal values, which are crucial for the accuracy of efficiency calculation model results. SPARC-90 code is an important scrubbing modelling tool developed for prediction of bubble hydrodynamics and the corresponding aerosol retention behavior. The potential for particle removal with a single rising deformed bubble was recently determined using a prototype of the SPARC-90 code with the experimental results. The results of SPARC-90 calculations with various equivalent diameters obtained in this experiment have provided more accurate decontamination factors compared to those obtained with a round bubble and a default bubble diameter.

Cavitating flow through a Venturi using CFD: effects of inlet and outlet pressures

Hydrodynamic cavitation is a complex physical phenomenon that appeared in hydraulic systems (pumps, turbines, valves, Venturi tubes …etc.) when the fluid pressure decreases below the saturated vapor pressure. The works carried out in this study aimed to get a better understanding of the cavitating flow phenomenon. For this, we have numerically studied a cavitating bubbly flow through a venturi nozzle. The cavitation model is selected and solved using commercial code computational fluid dynamics (CFD). The numerical results are compared to the previous experimental and numerical data.The obtained results showed that the effect of the inlet pressure (10, 7, 5 and 2 bars) and outlet pressures (6.5, 6, 5, 4, 3.5, 1.5, and 1 bars) of the Venturi on pressure, velocity of the fluid flow and the vapor fraction formation. We found that the inlet and outlet pressures of the venturi are strongly affected on the evolution of the pressure, velocity and vapor fraction formation in the cavitating flow. The results presented in this study provide insight and useful guidance to the designer of hydrodynamic cavitation devices, identifying key geometrical and physical parameters of cavitation onset and intensification that can be manipulated to achieve the desired level of cavitation flow.

Dynamic Sliding Mode Control of Spherical Bubble for Cavitation Suppression

Cavitation is a disadvantageous phenomenon that occurs when fluid pressure drops below its vapor pressure. Under these conditions, bubbles form in the fluid. When these bubbles flow into a high-pressure area or tube, they erupt, causing harm to mechanical parts such as centrifugal pumps. The difference in pressure in a fluid is the result of varying temperatures. One way to eliminate cavitation is to reduce the radius of the bubbles to zero before they reach high-pressure areas, using a robust approach. In this paper, sliding mode control is used for this purpose due to its invariance property. To force the radius of the bubbles toward zero and prevent chattering, a new dynamic sliding mode control approach is used. In dynamic sliding mode control, chattering is removed by passing the input control through a low-pass filter, such as an integrator. A general model of the spherical bubble is used, transferred to the state space, and then a state proportional-integral feedback is applied to obtain a linear system with a new input control signal. A comparison is also made with traditional sliding mode control using state feedback, providing a trusted comparison.

Experimental study on condensing flow and heat transfer characteristics in minichannels under mechanical vibration conditions

Minichannel heat exchangers are often used in the condensers of vehicle air conditioners, and vibration is inevitable in vehicle driving. However, the influence of mechanical vibration on the characteristics of minichannel condensing flow and heat transfer has been rarely studied. Therefore, the experimental study of condensing flow and heat transfer in minichannels under mechanical vibration was carried out in this paper. With R134a as the working medium, the condensation heat transfer coefficient, frictional pressure drop and flow pattern images were measured in the saturation temperature of 40 °C, mass flux density of 200 kg/(m2·s), vapor quality range of 0–1, vibration frequency of 15–50 Hz and amplitude of 0–2.4 mm. For annular flow, intermittent flow and bubble flow, the vibration can enhance the wave of gas-liquid interface, and promote the breakup of long bubbles and the polymerization of small bubbles. Vibration enhances heat transfer in most cases, but weakens heat transfer in some specific cases. The frictional pressure drops under vibration conditions are greater than those under static conditions, which indicates that the influence mechanism of vibration on heat transfer and pressure drop is different. This paper attempts to analyze the complex mechanism of vibration affecting minichannel condensing heat transfer and pressure drop, which provides theoretical support for further study of condenser optimization design under vibration conditions.

Investigating of vibration and noise characteristics in two-phase gas-liquid flow through a spiral capillary tube

The vaporization of refrigerant in the liquid-phase in spiral capillary tube induces complex two-phase flow that causes vibrations and noise, which affect the performance and stability of refrigeration systems. Therefore, the flow patterns, vibrations, flow-borne and structure-borne noises in spiral capillary tube are explored, and the effects of coil diameters, pitches and temperatures are analyzed. The results showed that the flow upstream the vaporization point was dominant by liquid-phase, then changed to bubbly flow downstream the vaporization point under various structures, and further changed to mist flow near the outlet, while it changed to gas-phase flow near the outlet under high temperature. The flow-borne noise was mainly affected by the flow turbulence, and the total sound pressure level (TSPL) increased by 17.1 % on average as the temperature increased from 309.6 to 317.6 K, and rose with the increase in pitch. As the coil diameter increased, the TSPL first increased and then decreased. Moreover, the maximum deformation of structure increased by 58.33 % as the diameter increased from 35 to 50 mm, but changed little with various pitches. The changing trends of vibration intensity caused the similar variations of structure-borne noise. The TSPL of structure-borne noise increased by 9.14 % with the diameter increasing from 35 to 50 mm, but it changed little at various pitches. The TSPL of structure-borne noise was much higher than flow-borne noise, which meant that improving the structural vibration can control the noise of refrigerant flow in spiral capillary tube.

Self-similar solution for laminar bubbly flow evolving from a vertical plate

The development of a bubble plume from a vertical gas-evolving electrode is driven by buoyancy and hydrodynamic bubble dispersion. This canonical fluid mechanics problem is relevant for both thermal and electrochemical processes. We adopt a mixture model formulation for the two-phase flow, considering variable density (beyond Boussinesq), viscosity and hydrodynamic bubble dispersion. Introducing a new change of coordinates, inspired by the Lees–Dorodnitsyn transformation, we obtain a new self-similar solution for the laminar boundary layer equations. The results predict a wall gas fraction and gas plume thickness that increase with height to the power of 1/5 before asymptotically reaching unity and scaling with height to the power 2/5, respectively. The vertical velocity scales with height to the power of 3/5. Our analysis shows that self-similarity is only possible if gas conservation is entirely formulated in terms of the gas specific volume instead of the gas fraction.

Relation between bubble behaviors in spray sheet and bubbly flow inside an air-induction fan nozzle: A visualization study

The present study aims to elucidate the relationship between bubble behaviors in the spray sheet and internal bubbly flow of the air-induction nozzle. An experimental work was performed using the visualization technique. Effects of the air inlet position and spray pressure were investigated. The results show that compared with the bubbles inside the air-induction nozzle, bubbles in the spray sheet have smaller volume but larger average diameter. Disturbance propagates in the horizontal air-inlet segment. When the air inlet position shifts toward the nozzle outlet, overall bubble volume inside the nozzle decreases by about 56%, while in the spray sheet, the bubble volume decreases by about 77%. Bubble breakup causes a decrease in overall bubble volume as bubbles travel from the inner flow passage of the nozzle to the environment. The coalescence and compression of bubbles induce the increase in average bubble diameter. When the spray pressure increases from 0.1 to 0.3 MPa, both the total bubble volume and average bubble diameter increase.

Experimental and numerical analysis of the fluid flow behavior of a tank corner impacting a water surface

The interaction of a tank impacting a water surface is an extremely complex nonlinear multiphase flow phenomenon. In this study, experiments and numerical simulations are used to systematically investigate the flow physics and load characteristics of a tank corner impacting a water surface. Free surface flow at different fall heights (200–800 mm) and inclination angles (0°–15°) was obtained through free fall experiments. The volume of fluids method and overset grid technology were used to simulate the water impact process of a three-dimensional structure accurately. For typical bubble flows, the numerical and experimental results agree well. On the basis of the three-dimensional flow characteristics and pressure distribution, flow behaviors, such as fluid climbing, corrugation disturbances, and air cavity effects, are analyzed. Bubble flow has a significant effect on the behavior mode of the impact load. In particular, the bubbles at the upper wall play a key role in the load characteristics at different locations. In addition, the influences of corrugations inside the tank's corner and the impact velocity on fluid flow were investigated. These results provide beneficial references for an in-depth understanding of the fluid flow and load characteristics between a tank and fluid.

Void drift in narrow rectangular channels for highly advective bubbly flows: A macroscopic drift-flux model derived from two-fluid local simulations

This article addresses the issue of reduced models to describe turbulent two-phase flows in industrial applications. This work is an application based on Bois (2021) theoretical derivations to connect the spatially-averaged drift-flux model (DFM) to the local Reynolds-averaged two-fluid model (TFM). It presents new opportunities to calibrate closure laws or propose new models based on finer-scale descriptions. Highly convective steady bubbly flows with wall-peaking void fraction profiles are considered. Dispersed bubbles are weakly deformable. A new model dedicated to void fraction dispersion in flat rectangular channels is developed for application into the DFM. The model is calibrated based on reference local Euler-Euler two-fluid simulations of pressurised water and steam mono-dispersed bubbles, in adiabatic conditions. 64 conditions are considered to cover regularly the 3D parameter space in void fraction, Reynolds and Eötvös numbers. The new formulation proposed shows a very good fit with the post-processed predictions of local CFD simulations along with a substantial reduction in grid requirement. The model contains sub-filter correlations between vapour concentration and mixture or relative velocities. This model derived from local simulations is validated a posteriori in an industrial component scale code where the full space-filtered DFM is implemented and resolved. This work proves the benefits and feasibility of one-by-one model development and calibration based on two-fluid simulations.

Non-negligible buoyancy effect on bubbles travelling in horizontal microchannels of comparable size at small Bond numbers

When a gas bubble is transported by a liquid flow confined within a horizontal channel of comparable size, buoyancy effects are usually assumed to be negligible if the Bond number of the flow is less than unity. However, recent experimental studies showed that buoyancy may still significantly impact the bubble dynamics even when Bo≪1, provided that the flow speed is sufficiently small, such that the viscous and inertial forces are weak. To derive a new criterion to assess the significance of buoyancy on the flow of small bubbles in horizontal microchannels, we have performed systematic numerical simulations using the free software Basilisk, covering a wide range of Bond, capillary and Reynolds numbers, Bo=0.004−0.4, Ca=10−4−0.5 and Re=0−100, and exploring the bubble-to-channel diameter ratios db/D=0.2−0.9. We demonstrate that the nondimensional group Bo(db/D)2/Ca is effective in assessing the importance of buoyancy in flows with negligible inertial effects. When Bo(db/D)2/Ca<0.1, buoyancy effects are negligible and the bubble travels along the channel axis. When Bo(db/D)2/Ca>10, buoyancy effects dominate and the bubble travels in the vicinity of the upper wall. For intermediate values, bubbles take equilibrium positions between the channel centre and the wall, and the threshold Bo(db/D)2/Ca=1 is effective in predicting whether bubbles will travel closer to the channel centre or to the wall. The capillary number at the denominator describes the lift force generated by the deformation of the bubble due to viscous shear, which acts towards the channel centre thus opposing buoyancy. Inertial forces significantly alter the shape and equilibrium position of the bubble when the Weber number of the flow exceeds unity and Bo/We<1. Under the action of inertia, bubbles of size comparable to the channel diameter become elongated but remain centred. Smaller bubbles migrate towards the channel wall and do not exhibit a preferential near-wall circumferential position.

Hydrodynamics of molten media bubble columns for hydrogen production through methane pyrolysis

Methane pyrolysis using a molten media bubble column reactor is a promising technique for hydrogen production with low carbon dioxide emissions at a feasible price. Understanding the bubble dynamics in molten media is essential to elucidate the reaction mechanisms and establish design requirements for efficient reactors. Computational fluid dynamics provides an effective means to understand the hydrodynamics in opaque molten media. This research used the volume of fluid method to study the effects of gas injection rate as well as variations in gas and molten media (iron, aluminum, and a salt mixture of sodium bromide and potassium bromide in a 48.7:51.3 molar ratio) properties on bubble dynamics. The computational model was first validated using existing experimental and empirical observations. This study makes fundamental contributions to the understanding of bubble dynamics in molten media. First, it was confirmed that gas properties had a small effect on bubble dynamics. The difference in bubble diameters between argon at ambient temperature and 1600 °C was less than 10%. Second, it was found that the volumetric gas injection rate and molten media properties significantly impacted the bubble dynamics, including the bubble diameter and flow regime. Future work will build on these findings to recommend appropriate operating conditions and molten media for specific pyrolysis reactor designs.

An experimental study on ice melting processes under point-source bubble flows at different flow rates

The anti-/de-icing capability of ships and offshore structures in the polar regions is of importance to ensure the safety of operation. The bubble anti-/de-icing method has great application potential. Here, a point-source bubbler system is developed to study the ice melting processes under point-source bubble flows, especially ice melting stage, bubble distribution, ice melting rate, and final ice morphology. The ice melting process is divided into flat, concave, and holed ice stages. With the increase of the flow rate, the duration of the ice melting process gradually decreases while that of the flat ice stage increases and that of the concave ice stage decreases. The number density of bubbles at the concave ice stage is the smallest and the average contact area of bubbles at the concave ice stag is the largest of the three stages. The average contact areas of bubbles at 1.0 L/min and 1.5 L/min are significantly larger than those at 0.5 L/min and 2.0 L/min at the concave ice stage. When the flow rate increases from 0.5 L/min to 2.0 L/min, the melting rate in the height direction increases by 95.4% while the melting rate in the radial direction increases by 61.8%. The cross-sectional profile of the final ice morphology gradually becomes steeper as the flow rate rises. The findings of this work provide insights into the ice melting mechanism under bubble flows and are helpful to the optimization of related applications.

Sensitized Reynolds stress modeling of a bubbly jet emerging into a water cross-flow

The present work is concerned with the computational study of an air–water mixture stream emanating from a nozzle submerged in a water cross-flow inside a rectangular open channel to form a curved, intensively dispersed, bubbly jet. The work focuses on evaluating the predictive performance of an eddy-resolving turbulence model applied to the described case of a multiphase gas-liquid flow system characterized by a variety of closely coupled phenomena, including: turbulence anisotropy-induced secondary motion, free surface flow, bubbly jet propagation, and the varying interaction dynamics of the carrier water flow with the air bubble dispersion. Correspondingly, an appropriately extended differential near-wall Reynolds stress model, which describes the dynamics of unresolved subscale structures within the Sensitized Reynolds-Averaged Navier–Stokes (RANS) computational framework, is currently used in conjunction with a two-way coupled Euler–Lagrange approach to simulate this gas-water bubbly flow, with the operating and boundary conditions following the experimental reference of Zhang and Zhu (2013). Accordingly, the conditions at the free surface are adequately derived for all components of the residual turbulence stress tensor and the corresponding length scale determining variable. It is shown that the statistics of the bubble phase can be accurately captured, and the proper-orthogonal-decomposition analysis of the dynamic properties of the flow reveals at least two large-scale transient effects associated with the bubbly jet. Furthermore, in the preliminary part of this work it is shown that the currently applied turbulence model can be successfully used to correctly capture the effects of Reynolds stress anisotropy in the single-phase open-channel configuration, representing the incoming flow field of the main two-phase flow configuration.

Numerical investigation on the slamming loads of a truncated trimaran hull entering regular waves

For ship slamming, waves can disturb the water entry progress and thus affect the slamming load characteristics. This paper studied the slamming loads of a trimaran entering regular waves based on the CFD method. The results of a trimaran model impacting calm water and wedge impacting waves are used to verify the present method. The free surface evolution and load characteristics of the trimaran model considering wave effects are discussed. The results showed that the position at which the structure touches the wave surface depends on the formation, deformation and escape of bubble flows, which cause obvious asymmetrical flow. The enclosure of the air cavity induces a step increase in pressure, and the pressure peak appears when the bubble fully escapes. Three typical positions, wave crest, wave cross and wave through, are chosen to explore the influence of waves on slamming loads. The impact load shows a 98% difference at different wave locations. In ship load forecasting, the position of the waves relative to the hull is noteworthy. The influence of the water entry velocity on the slamming load is further discussed. Bubble retention related to the water entry velocity is the key factor affecting the peak load.