Numerical simulation for a rising bubble interacting with a solid wall: Impact, bounce, and thin film dynamics

Using an arbitrary Lagrangian-Eulerian method on an adaptive moving unstructured mesh, we carry out numerical simulations for a rising bubble interacting with a solid wall. Driven by the buoyancy force, the axisymmetric bubble rises in a viscous liquid toward a horizontal wall, with impact on and possible bounce from the wall. First, our simulation is quantitatively validated through a detailed comparison between numerical results and experimental data. We then investigate the bubble dynamics which exhibits four different behaviors depending on the competition among the inertial, viscous, gravitational, and capillary forces. A phase diagram for bubble dynamics has been produced using the Ohnesorge number and Bond number as the two dimensionless control parameters. Finally, we turn to the late stage of the bubble rise characterized by a small flux of liquid escaping from the thin film between the wall and the bubble. Since the thin film dynamics can be accurately described by the lubrication approximation, we carry out numerical simulations to compare the simulation results with the predictions of the lubrication approximation. Remarkable agreement is obtained to further demonstrate the accuracy of the simulations.</jats:p

Taylor bubble moving in a flowing liquid in vertical channel: transition from symmetric to asymmetric shape

The velocity and shape of Taylor bubbles moving in a vertical channel in a Poiseuille liquid flow were studied for the inertial regime, characterized by large Reynolds numbers. Numerical experiments were carried out for positive (upward) and negative (downward) liquid mean velocity. Previous investigations in tube have reported that for upward flow the bubble is symmetric and its velocity follows the law of Nicklin whereas for certain downward flow conditions the symmetry is broken and the bubble rises appreciably faster. To study the bubble motion and to identify the existence of a transition, a 2D numerical code that solves the Navier-Stokes equations (through a VoF implementation) was used to obtain the bubble shape and the rise velocity for different liquid mean velocity. A reference frame located at the bubble tip as well as an irregular grid were implemented to allow for long simulation times without an excessively large numerical domain. It was observed that whenever the mean liquid velocity exceeded some critical value, bubbles adopted a symmetric final shape even though their initial shape was asymmetric. Conversely, if the mean liquid velocity was smaller than that critical value, a transition to a non-symmetric shape occurred, along with a correspondingly faster velocity. It was also found that surface tension has a stabilizing effect on the transition

Computational analysis of single rising bubbles influenced by soluble surfactant

This paper presents novel insights about the influence of soluble surfactants on bubble flows obtained by Direct Numerical Simulation (DNS). Surfactants are amphiphilic compounds which accumulate at fluid interfaces and significantly modify the respective interfacial properties, influencing also the overall dynamics of the flow. With the aid of DNS local quantities like the surfactant distribution on the bubble surface can be accessed for a better understanding of the physical phenomena occurring close to the interface. The core part of the physical model consists in the description of the surfactant transport in the bulk and on the deformable interface. The solution procedure is based on an Arbitrary Lagrangian-Eulerian (ALE) Interface-Tracking method. The existing methodology was enhanced to describe a wider range of physical phenomena. A subgrid-scale (SGS) model is employed in the cases where a fully resolved DNS for the species transport is not feasible due to high mesh resolution requirements and, therefore, high computational costs. After an exhaustive validation of the latest numerical developments, the DNS of single rising bubbles in contaminated solutions is compared to experimental results. The full velocity transients of the rising bubbles, especially the contaminated ones, are correctly reproduced by the DNS. The simulation results are then studied to gain a better understanding of the local bubble dynamics under the effect of soluble surfactant. One of the main insights is that the quasi-steady state of the rise velocity is reached without ad- and desorption being necessarily in local equilibrium

The evolution of energy in flow driven by rising bubbles

We investigate by direct numerical simulations the flow that rising bubbles cause in an originally quiescent fluid. We employ the Eulerian-Lagrangian method with two-way coupling and periodic boundary conditions. In order to be able to treat up to 288000 bubbles, the following approximations and simplifications had to be introduced: (i) The bubbles were treated as point-particles, thus (ii) disregarding the near-field interactions among them, and (iii) effective force models for the lift and the drag forces were used. In particular, the lift coefficient was assumed to be 1/2, independent of the bubble Reynolds number and the local flow field. The results suggest that large scale motions are generated, owing to an inverse energy cascade from the small to the large scales. However, as the Taylor-Reynolds number is only in the range of 1, the corresponding scaling of the energy spectrum with an exponent of -5/3 cannot develop over a pronounced range. In the long term, the property of local energy transfer, characteristic of real turbulence, is lost and the input of energy equals the viscous dissipation at all scales. Due to the lack of strong vortices the bubbles spread rather uniformly in the flow. The mechanism for uniform spreading is as follows: Rising bubbles induce a velocity field behind them that acts on the following bubbles. Owing to the shear, those bubbles experience a lift force which make them spread to the left or right, thus preventing the formation of vertical bubble clusters and therefore of efficient forcing. Indeed, when the lift is artifically put to zero in the simulations, the flow is forced much more efficiently and a more pronounced energy accumulates at large scales is achieved.Comment: 9 pages, 7 figure

Numerical Investigation of Bubble Movement in Magnetic Nanofluids

Department of Mechanical EngineeringIn this study, the idea to generate electrical energy by using waste heat is suggested. In this idea, the electrical energy can be generated by a magnetic nanofluid and bubble movement. Thus, bubble movement in a magnetic fluid is numerically investigated using the commercial CFD package COMSOL Multiphysics for effective energy generation. The slug characteristics are also investigated because it can be generated by merging each bubble. The level-set method and phase-field method are used to simulate the bubble and slug movement, respectively. For the investigation, EFH1 and EFH3 are selected as working fluidsthey are commercial magnetic fluids manufactured by Ferrotec, and each fluid contains different amounts of magnetic particles. The solvers are validated by comparing the numerical results with previous research studies and experimental data for reliable results. The properties of a fluid can be changed by solid particles if the particles are dispersed in the fluid. These particles can affect the bubble and slug characteristics, such as shape, velocity and wake. Thus, the effect of solid particles is first studied by observing the bubble and slug movement in each magnetic fluid. In the slug investigation, the effects of some parameters are also studied, such as slug length and liquid backflow. The patterns of bubble and slug movement are investigated to predict the effective condition for the energy generation. The effective bubble and slug movements are predicted by evaluating the disturbance of the surrounding magnetic fluid for each flow pattern. A magnetic force can be created, and it can affect any phenomena when the magnetic field is applied to the system. Thus, the effect of the magnetic force is also investigated because the magnetic field should be applied to the system for the electrical energy generation. Finally, the important factor for energy generation is identified by comparing the results.ope

Direct numerical simulation of the drag force in bubble swarms

This paper studies the swarm effect on the drag force in bubbly flows. From literature it is well-known that for relatively small bubbles, the drag force increases with the bubble hold-up due to additional hindrance experienced by the bubbles caused by the modified flow field. Very large (spherical cap) bubbles on the other hand may rise cooperatively. The unique capabilities of a 3D Front Tracking model were used to investigate the influence of important parameters like the gas fraction, Reynolds number and the bubble size independently. It was found that the relative drag force increases for bubbles in the range of 2 to 5 mm when the gas fraction is increased up to 13%, while the bubbles become more spherical. Also the influence of the Reynolds number and the bubble aspect ratio on the increased drag force has been determined. It can be concluded that there is only a very weak effect over several decades of the Reynolds number, while there is a strong effect of the bubble aspect ratio.\ud This also helps explaining why the increase in drag is smaller for larger bubbles: when the gas fraction is increased deformable bubbles become more spherical, thus reducing the drag force

On the lateral migration of a slightly deformed bubble rising near a vertical plane wall

Deformation-induced lateral migration of a bubble slowly rising near a vertical plane wall in a stagnant liquid is numerically and theoretically investigated. In particular, our focus is set on a situation with a short clearance $c$ between the bubble interface and the wall. Motivated by the fact that numerically and experimentally measured migration velocities are considerably higher than the velocity estimated by the available analytical solution using the Fax\'{e}n mirror image technique for $a/(a+c)\ll 1$ (here $a$ is the bubble radius), when the clearance parameter $\varepsilon(= c/a)$ is comparable to or smaller than unity, the numerical analysis based on the boundary-fitted finite-difference approach solving the Stokes equation is performed to complement the experiment. The migration velocity is found to be more affected by the high-order deformation modes with decreasing $\varepsilon$. The numerical simulations are compared with a theoretical migration velocity obtained from a lubrication study of a nearly spherical drop, which describes the role of the squeezing flow within the bubble-wall gap. The numerical and lubrication analyses consistently demonstrate that when $\varepsilon\leq 1$, the lubrication effect makes the migration velocity asymptotically $\mu V_{B1}^2/(25\varepsilon \gamma)$ (here, $V_{B1}$, $\mu$, and $\gamma$ denote the rising velocity, the dynamic viscosity of liquid, and the surface tension, respectively).Comment: 24 pages, 9 figures, J. Fluid Mech. (accepted

Bubbly and Buoyant Particle-Laden Turbulent Flows

Fluid turbulence is commonly associated with stronger drag, greater heat transfer, and more efficient mixing than in laminar flows. In many natural and industrial settings, turbulent liquid flows contain suspensions of dispersed bubbles and light particles. Recently, much attention has been devoted to understanding the behavior and underlying physics of such flows by use of both experiments and high-resolution direct numerical simulations. This review summarizes our present understanding of various phenomenological aspects of bubbly and buoyant particle-laden turbulent flows. We begin by discussing different dynamical regimes, including those of crossing trajectories and wake-induced oscillations of rising particles, and regimes in which bubbles and particles preferentially accumulate near walls or within vortical structures. We then address how certain paradigmatic turbulent flows, such as homogeneous isotropic turbulence, channel flow, Taylor-Couette turbulence, and thermally driven turbulence, are modified by the presence of these dispersed bubbles and buoyant particles. We end with a list of summary points and future research questions.Comment: 29 pages, 14 figure

Numerical Simulation of Heat Transport in Dispersed Gas-Liquid Two-Phase Flow using a Front Tracking Approach

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of heat transport in dispersed gas-liquid two-phase flow using the Front Tracking (FT) approach. Our model extends the FT model developed by van Sint Annaland et al. (2006) to non-isothermal conditions. In FT an unstructured dynamic mesh is used to represent and track the interface explicitly by a number of interconnected marker points. The Lagrangian representation of the interface avoids the necessity to reconstruct the interface from the local distribution of the fractions of the phases and, moreover, allows a direct and accurate calculation of the surface tension force circumventing the (problematic) computation of the interface curvature. The extended model is applied to predict the heat exchange rate between the liquid and a hot wall kept at a fixed temperature. It is found that the wall-to-liquid heat transfer coefficient exhibits a maximum in the vicinity of the bubble that can be attributed to the locally decreased thickness of the thermal boundary layer