If you are exporting VTU files, you can script Paraview to extract the bubble properties. A pipeline such as:
- contour Epg = 0.9
- marching cubes to segment the contours
- vtk mass properties to extract bubble size, etc.
Bubble rise velocity is a little more difficult. Need to match bubbles between time steps and calculate the displacement.