Here comes the patch
Line 64 as well: + alpha/rhob*(liftCoeff + Cvm*rhob*DDtUa) should be + alpha/rhob*(liftCoeff...
Bug in bubbleFoam solver: UEqns.H, line 46