A boundary integral method was developed for simulating the motion and deformation of a viscous drop in an axisymmetric ambient Stokes flow near a rigid wall and for direct calculating the stress on the wall. Numerical experiments by the method were performed for different initial stand-off distances of the drop to the wall, viscosity ratios, combined surface tension and buoyancy parameters and ambient flow parameters. Numerical results show that due to the action of ambient flow and buoyancy the drop is compressed and stretched respectively in axial and radial directions when time goes. When the ambient flow action is weaker than that of the buoyancy the drop raises and bends upward and the stress on the wall induced by drop motion decreases when time advances. When the ambient flow action is stronger than that of the buoyancy the drop descends and becomes flatter and flatter as time goes. In this case when the initial stand-off distance is large the stress on the wall increases as the drop evolutes but when the stand-off distance is small the stress on the wall decreases as a result of combined effects of ambient flow, buoyancy and the stronger wall action to the flow. The action of the stress on the wall induced by drop motion is restricted in an area near the symmetric axis, which increases when the initial stand-off distance increases. When the initial stand-off distance increases the stress induced by drop motion decreases substantially. The surface tension effects resist the deformation and smooth the profile of the drop surfaces. The drop viscosity will reduce the deformation and migration of the drop.
A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier-Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving inter-face is captured by the level set function, and the interface velocity is resolved by ＂one-side＂ velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the ＂shock wave＂-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating.
In the inviscid and incompressible fluid flow regime,surface tension effects on the behaviour of an initially spherical buoyancy-driven bubble rising in an infinite and initially stationary liquid are investigated numerically by a volume of fluid （VOF） method. The ratio of the gas density to the liquid density is 0.001, which is close to the case of an air bubble rising in water. It is found by numerical experiment that there exist four critical Weber numbers We1,~We2,~We3 and We4, which distinguish five different kinds of bubble behaviours. It is also found that when 1≤We2, the bubble will finally reach a steady shape, and in this case after it rises acceleratedly for a moment, it will rise with an almost constant speed, and the lower the Weber number is, the higher the speed is. When We 〉We2, the bubble will not reach a steady shape, and in this case it will not rise with a constant speed. The mechanism of the above phenomena has been analysed theoretically and numerically.
In the inviscid and incompressible fluid flow regime, surface tension effects on the behavior of two initially spherical bubbles with same size rising axisymmetrically in an infinite and initially stationary liquid are investigated numerically with the VOF method. The numerical experiments are performed for two bubbles with two different bubble distances. The ratio of gas density to liquid density is 0.001, which is close to the case of air bubbles rising in water. In the case of Dis = 2.5, where Dis is defined as the ratio of the distance between the bubble centroids to the radius of the bubble, it is found from numerical experiments that there exist four critical Weber numbers We1 , We2 , We3 and We4 , which are in between 10 and 100, 3 and 4, 1.5 and 1.8, and 0.2 and 0.3, respectively. In the case of Dis = 2.3, similar phenomena also appear but the corresponding four critical Weber numbers are lower than those in the case of Dis = 2.5. The mechanism of the above phenomena is analyzed theoretically and numerically.
WANG HanZHANG Zhen-yuYANG Yong-mingZHANG Hui-sheng
In the incompressible fluid flow regime, without taking consideration of surface tension effects, the viscosity effects on the behavior of an initially spherical buoyancy-driven bubble rising in an infinite and initially stationary liquid are investigated numerically by the Volume Of Fluid （VOF） method. The ratio of the gas density to the liquid density is taken as 0.001, and the gas viscosity to the liquid viscosity is 0.01, which is close to the case of an air bubble rising in water. It is found by numerical experiments that there exist two critical Reynolds numbers Re1 and Re2 , which are in between 30 and 50 and in between 10 and 20, respectively. As Re 〉 Re1 the bubble will have the transition to toroidal form, and the toroidal bubble will break down into two toroidal bubbles. In this case viscosity will damp the development of the liquid jet and delay the formation of the toroidal bubble. As Re〈Re1 the transition will not happen. As Re2 〈 Re 〈 Re1, the bubble will split from its rim into a toroidal bubble and a spherical cap-like bubble, and as Re〈Re2 the splitting will not occur and the bubble can finally reach a stationary shape. With the decrease of the Reynolds number, the stationary shape changes from spherical-cap bubble with skirt to dimpled peach-like bubble. Before the bubble reaches its stationary shape the vortex structure in the flow field varies with time. The vortex structure corresponding to bubble stationary shape varies with the Reynolds number. It is also found that there exists another critical Reynolds number Re^＊ which is in between Re1 and Re2 , and as Re 〈 Re^＊, after the bubble rises in an accelerating manner for a moment, it will rise with an almost constant speed, and the speed increases with increasing Reynolds number. As Re 〉 Re^＊, it will not rise with a constant speed. The mechanism of the above phenomena has been analyzed theoretically and numerically.
WANG HanZHANG Zhen-yuYANG Yong-mingZHANG Hui-sheng
A numerical method for simulating the motion and deformation of an axisymmetric bubble or drop rising or falling in another infinite and initially stationary fluid is developed based on the volume of fluid （VOF） method in the frame of two incompressible and immiscible viscous fluids under the action of gravity, taking into consideration of surface tension effects. A comparison of the numerical results by this method with those by other works indicates the validity of the method. In the frame of inviseid and incompressible fluids without taking into consideration of surface tension effects, the mechanisms of the generation of the liquid jet and the transition from spherical shape to toroidal shape during the bubble or drop deformation, the increase of the ring diameter of the toroidal bubble or drop and the decrease of its cross-section area during its motion, and the effects of the density ratio of the two fluids on the deformation of the bubble or drop are analysed both theoretically and numerically.
