Bubble stabilization of spectral methods: the multidimensional case