Coupled shape oscillations and translational motion of an incompressible gas bubble in a liquid container in response to forced vibration is studied numerically. Bond number (Bo) and the ratio of the vibration amplitude to the bubble diameter (A/D) are found to be the governing non-dimensional numbers. Bubble response is studied in both 2D and 3D. Different schemes are used for 2D and 3D simulations. In 2D, the flow solver is coupled to a Volume of Fluid (VOF) algorithm to capture the interface between the two phases while in 3D the interface is captured using a level set algorithm. The oscillation outcome ranges from small amplitude and regular oscillations for small Bo and A/D to large amplitude, nonlinear, and chaotic oscillations for large Bo and A/D. Chaotic behavior occurs due to the coupling between the nonlinear shape oscillations and large amplitude oscillatory translational motion. By further increase of the forcing, the inertia of the liquid results in the formation of a liquid jet which penetrates within the bubble core and pierces the bubble and a toroidal bubble shape is formed. The toroidal bubble shape then goes through large amplitude shape oscillations and smaller bubbles are formed. A summary of the 3D simulations provides a map which shows the bubble oscillation outcome as a function of Bo and A/D. The interaction between two bubbles is studied in 2D as well and the effect of vibration amplitude, frequency and liquid to gas density ratio on the interaction force is investigated.

Coupled shape oscillations and translational motion of an incompressible gas bubble in a liquid container in response to forced vibration is studied numerically. Bond number (Bo) and the ratio of the vibration amplitude to the bubble diameter (A/D) are found to be the governing non-dimensional numbers. Bubble response is studied in both 2D and 3D. Different schemes are used for 2D and 3D simulations. In 2D, the flow solver is coupled to a Volume of Fluid (VOF) algorithm to capture the interface between the two phases while in 3D the interface is captured using a level set algorithm. The oscillation outcome ranges from small amplitude and regular oscillations for small Bo and A/D to large amplitude, nonlinear, and chaotic oscillations for large Bo and A/D. Chaotic behavior occurs due to the coupling between the nonlinear shape oscillations and large amplitude oscillatory translational motion. By further increase of the forcing, the inertia of the liquid results in the formation of a liquid jet which penetrates within the bubble core and pierces the bubble and a toroidal bubble shape is formed. The toroidal bubble shape then goes through large amplitude shape oscillations and smaller bubbles are formed. A summary of the 3D simulations provides a map which shows the bubble oscillation outcome as a function of Bo and A/D. The interaction between two bubbles is studied in 2D as well and the effect of vibration amplitude, frequency and liquid to gas density ratio on the interaction force is investigated.