THE BUBMAC METHOD-A NUMERICAL METHOD BUBBLE DYNAMICS IN VORTEX CORE

This paper presents a numerical method called BUBMAC method,whichis based on the Marker And Cell(MAC) technique,to numerically simulate the problemof bubble dynamics in vortex core. With the incorporation of the azimuthal velocity intothe momentum equation, the model takes into account the complete interaction betweenthe bubble dynamics and the vortex flow field.The three momentum equations (N-Seqs.) are solved numerically by the finite-difference methed, and the motion of bubblesurface is described by tracing massless marker particles distributed only on the surface ofthe bubble. With some important medifications to the original MAC method, the numeri-cal accuracy of the method is conisderably increased. The validation of the proposed BUB-MAC method is checked by comparing the numerical results with some available analyticalsolutions in the cases of spherical bubble evolution and with numerical results in the casesof nonspherical collapse of bubble and bubble evolution in tip-vortex core. The compar-iso

THE BUBMAC METHOD-A NUMERICAL METHOD BUBBLE DYNAMICS IN VORTEX CORE

This paper presents a numerical method called BUBMAC method. which is based on the Marker And Cell (MAC) technique, to numerically simulate the problem of bubble dynamics in vortex core. With the incorporation of the azimuthal velocity into the momentum equation, the model takes into account the complete interaction between the bubble dynamics and the vortex flow field. The three momentum equations (N-S eqs. ) are solved numerically by the finite-difference method, and the motion of bubble surface is described by tracing massless marker particles distributed only on the surface of the bubble. With some important modifications to the original MAC method, the numerical accuracy of the method is conisderably increased. The validation of the proposed BUB-MAC method is checked by comparing the numerical results with some available analytical solutions in the cases of spherical bubble evolution and with numerical results in the cases of nonspherical collapse of bubble and bubble evolution in tip-vortex core. The comparisons show the high capability and reasonable numerical accuracy of the proposed method.