Two‐dimensional Boussinesq convection is studied numerically using two different methods: a filtered pseudospectral method and a high‐order accurateenoscheme. The issue whether finite time singularity occurs for initially smooth flows is investigated. In contrast to the findings of Pumir and Siggia who reported finite time collapse of the bubble cap, the present numerical results suggest that the strain rate corresponding to the intensification of the density gradient across the front saturates at the bubble cap. Consequently, the thickness of the bubble decreases exponentially. On the other hand, the bubble experiences much stronger straining and intensification of gradients at its side. As the bubble rises, a secondary front also forms from its tail. Together with the primary front, they constitute a pair of tightly bound plus and minus double vortex sheet structure which is highly unstable and vulnerable to viscous dissipation.