The effects of viscoelasticity on the expansion of gas bubbles arranged in a hexagonal array in a polymeric fluid are investigated. The expansion is driven by the diffusion of a soluble gas from the liquid phase, and the rate of expansion is controlled by a combination of gas diffusion, fluid rheology and surface tension. In the diffusion limited case, the initial growth rate is slow due to small surface area, whereas at high diffusivity initial growth is rapid and resisted only by background solvent viscosity. In this high Deborah number limit, we see a two stage expansion in which there is an initial rapid expansion up to the size at which the elastic stresses balance the pressure difference. Beyond this time the bubble expansion is controlled by the relaxation of the polymer. We also illustrate how viscoelasticity affects the shape of the bubble. In addition to a full finite element calculation of the two-dimensional flow, two one-dimensional approximations valid in the limits of small and large gas area fractions are presented. We show that these approximations give accurate predictions of the evolution of the bubble area, but give less accurate predictions of the bubble shape. (c) 2006 Elsevier B.V. All rights reserved.