Chinese Journal of Theoretical and Applied Mechanics ›› 2019, Vol. 51 ›› Issue (5): 1285-1295.DOI: 10.6052/0459-1879-19-071

• Fluid Mechanics • Previous Articles     Next Articles


Zhang Yang***,Chen Ke***2)(),You Yunxiang***,Sheng Li††   

  1. * State Key Laboratory of Ocean Engineering$,$ Shanghai Jiao Tong University$,$ Shanghai 200240, China
    ? School of Mechanical Engineering$,$ Purdue University$,$ West Lafayette 47907$,$ United States
    ** Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration$,$ Shanghai 200240, China
    ?? No. 92537 Unit of PLA$,$ Beijing 100161, China
  • Received:2019-03-26 Online:2019-09-18 Published:2019-09-30
  • Contact: Chen Ke


It is not only interesting but also complex that buoyancy-driven bubbles rising in viscous liquids. In particular, the interactions between bubbles and boundaries (e.g., solid walls) is relevant in practical applications and these interactions may have a significant effect in the global behaviors of the multi-phase fluids. In this work, the rising, collision and bouncing to a horizontal solid wall for a single bubble are studied by axisymmetric computations. The incompressible, density-variable Navier-Stokes equations with surface tension are used to describe the gas-liquid flow and are solved by a tree-based finite volume method (FVM). The evolution of bubble shape is implemented by using a volume of fluid (VOF) approach that combines a balanced surface tension force calculation and a height-function curvature estimation. To finely resolve the local but fast topological evolutions of bubble, the technique of adaptive mesh refinement (AMR) is used. Starting with the basic phenomenon of bubble impacting and bouncing, we explore the effects of Galilei number Ga and approach velocity $U_{\rm a}$. To study the bubble behaviors under different conditions both qualitatively and quantitatively, the evolution of the velocity vector field and a lot of parameters such as bouncing height $H$, bouncing period $T,$rising velocity $U$, axial coordinate $z$and coefficient of restitution $C_{\rm r }$are analyzed. Based on the results, we find that the bubble bouncing behaviors are pretty sensitive to the Galilei number. The increase of Ga promotes the bouncing and signifies the deformation of bubbles, increasing the collisions of bubbles, bouncing parameters and the coefficient of restitution. However, the value of $C_{\rm r}$is nearly unaffected by the variation of the approach velocity, indicating that $U_{\rm a}$is not a governing parameter for the bubble bouncing motion.

Key words: bubble bouncing dynamics, AMR-VOF, Galilei number, approach velocity

CLC Number: