No Arabic abstract
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 small clearance $c$ between the bubble interface and the wall. Motivated by the fact that experimentally measured migration velocity (Takemura et al. (2002, J. Fluid Mech. {bf 461}, 277)) is higher than the velocity estimated by the available analytical solution (Magnaudet et al. (2003, J. Fluid Mech. {bf 476}, 115)) using the Fax{e}n mirror image technique for $kappa(=a/(a+c))ll 1$ (here $a$ is the bubble radius), when the clearance parameter $epsilon(=c/a)$ is comparable to or smaller than unit, the numerical analysis based on the boundary-fitted finite-difference approach by solving the Stokes equation is performed to complement the experiment. To improve the understandings of a role of the squeezing flow within the bubble-wall gap, the theoretical analysis based on a soft-lubrication approach (Skotheim & Mahadevan (2004, Phys. Rev. Lett. {bf 92}, 245509)) is also performed. The present analyses demonstrate the migration velocity scales $propto{rm Ca} epsilon^{-1}V_{B1}$ (here, $V_{B1}$ and ${rm Ca}$ denote the rising velocity and the capillary number, respectively) in the limit of $epsilonto 0$.
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 $varepsilonleq 1$, the lubrication effect makes the migration velocity asymptotically $mu V_{B1}^2/(25varepsilon gamma)$ (here, $V_{B1}$, $mu$, and $gamma$ denote the rising velocity, the dynamic viscosity of liquid, and the surface tension, respectively).
Series of experiments on turbulent bubbly channel flows observed bubble clusters near the wall which can change large-scale flow structures. To gain insights into clustering mechanisms, we study the interaction of a pair of spherical bubbles rising in a vertical channel through combined experiments and modeling. Experimental imaging identifies that pairwise bubbles of 1.0 mm diameter take two preferred configurations depending on their mutual distance: side-by-side positions for a short distance ($S<5$) and nearly inline, oblique positions for a long distance ($S>5$), where $S$ is the mutual distance normalized by the bubble radius. In the model, we formulate the motions of pairwise bubbles rising at $Re=O(100)$. Analytical drag and lift, and semi-empirical, spatio-temporal stochastic forcing are employed to represent the mean acceleration and the fluctuation due to turbulent agitation, respectively. The model is validated against the experiment through comparing Lagrangian statistics of the bubbles. Simulations using this model identify two distinct timescales of interaction dynamics which elucidate the preferred configurations. For pairs initially in-line, the trailing bubble rapidly escapes from the viscous wake of the leading bubble to take the oblique position. Outside of the wake, the trailing bubble travels on a curve-line path with a slower velocity driven by potential interaction and horizontally approaches the leading bubble to become side-by-side. Moreover, statistical analysis identifies that the combination of the wake and the agitation can significantly accelerate the side-by-side clustering of in-line pairs. These results indicate positive contributions of liquid viscosity and turbulence to the formation of bubble clusters.
Droplet migration in a Hele--Shaw cell is a fundamental multiphase flow problem which is crucial for many microfluidics applications. We focus on the regime at low capillary number and three-dimensional direct numerical simulations are performed to investigate the problem. In order to reduce the computational cost, an adaptive mesh is employed and high mesh resolution is only used near the interface. Parametric studies are performed on the droplet horizontal radius and the capillary number. For droplets with an horizontal radius larger than half the channel height the droplet overfills the channel and exhibits a pancake shape. A lubrication film is formed between the droplet and the wall and particular attention is paid to the effect of the lubrication film on the droplet velocity. The computed velocity of the pancake droplet is shown to be lower than the average inflow velocity, which is in agreement with experimental measurements. The numerical results show that both the strong shear induced by the lubrication film and the three-dimensional flow structure contribute to the low mobility of the droplet. In this low-migration-velocity scenario the interfacial flow in the droplet reference frame moves toward the rear on the top and reverses direction moving to the front from the two side edges. The velocity of the pancake droplet and the thickness of the lubrication film are observed to decrease with capillary number. The droplet velocity and its dependence on capillary number cannot be captured by the classic Hele--Shaw equations, since the depth-averaged approximation neglects the effect of the lubrication film.
Hard particle erosion and cavitation damage are two main wear problems that can affect the internal components of hydraulic machinery such as hydraulic turbines or pumps. If both problems synergistically act together, the damage can be more severe and result in high maintenance costs. In this work, a study of the interaction of hard particles and cavitation bubbles is developed to understand their interactive behavior. Experimental tests and numerical simulations using computational fluid dynamics (CFD) were performed. Experimentally, a cavitation bubble was generated with an electric spark near a solid surface, and its interaction with hard particles of different sizes and materials was observed using a high-speed camera. A simplified analytical approach was developed to model the behavior of the particles near the bubble interface during its collapse. Computationally, we simulated an air bubble that grew and collapsed near a solid wall while interacting with one particle near the bubble interface. Several simulations with different conditions were made and validated with the experimental data. The experimental data obtained from particles above the bubble were consistent with the numerical results and analytical study. The particle size, density and position of the particle with respect to the bubble interface strongly affected the maximum velocity of the particles.
We consider the collapse behavior of cavitation bubbles near walls under high ambient pressure conditions. Generic configurations with different stand-off distances are investigated by numerical simulation using a fully compressible two-phase flow solver including phase change. The results show that the stand-off distance has significant effects on collapse dynamics, micro-jet formation, rebound, and maximum wall pressure. A relation between cavitation induced material damage and corresponding collapse mechanisms is obtained from pressure-impact data at the wall. We analyze the resolution dependence of collapse and rebound and the observed maximum pressure distributions. The comparison of the results on six different grid resolutions shows that main collapse features are already captured on the coarsest resolution, while the peak pressures are strongly resolution dependent.