Thebubblesriseupandburstatthefreesurfaceisacomplextwo-phaseprocess.AfreeenergylatticeBoltzmannmethod(LBM)modelisadoptedinthispapertostudythisphenomenon.Theinterfacecapturingtechnique[Zhengetal.,2006]isusedtodealwiththehighdensityratioproblem.TheLaplacelawandtheair-waterinterfacecapturingabilityarevalidatedforthemultiphasemodel.Theinteractionbetweenthesinglebubbleormultiplebubblesandthefreesurfacearestudiedbythemultiphasemodel.Theforceactingonthebubbleandtheevolutionofthefreesurfaceisstudied.Meanwhile,effectoftheinitialdistancebetweentwoadjacentbubblesoninteractioneffectsofmultiplebubblesisinvestigatedaswell.