Published online by Cambridge University Press: 12 March 2014
An isolated buoyant thermal in very viscous fluid has been shown to attain a self-similar form at large times which grows as   $t^{1/2}$  (Whittaker & Lister J. Fluid Mech., vol. 606, 2008, pp. 295–324). For large values of the Rayleigh number
 $t^{1/2}$  (Whittaker & Lister J. Fluid Mech., vol. 606, 2008, pp. 295–324). For large values of the Rayleigh number   $\mathit{Ra}$  (based on the conserved total buoyancy), the similarity solution is slender with a roughly spherical head at the top and a long tail that contains most of the buoyancy and extends down to the origin. We investigate the time-dependent behaviour of the thermal numerically; both the long-time behaviour in terms of perturbations to the similarity solution and the short-time evolution from a spherical initial condition. Using a spectral method, we find the growth rates of the linear perturbations and their spatial structure in similarity space. All eigenmodes decay monotonically for
 $\mathit{Ra}$  (based on the conserved total buoyancy), the similarity solution is slender with a roughly spherical head at the top and a long tail that contains most of the buoyancy and extends down to the origin. We investigate the time-dependent behaviour of the thermal numerically; both the long-time behaviour in terms of perturbations to the similarity solution and the short-time evolution from a spherical initial condition. Using a spectral method, we find the growth rates of the linear perturbations and their spatial structure in similarity space. All eigenmodes decay monotonically for   $\mathit{Ra}\lesssim 360$ , while for larger
 $\mathit{Ra}\lesssim 360$ , while for larger   $\mathit{Ra}$  the dominant (slowest decaying or fastest growing) eigenmodes are oscillatory with waves propagating up the tail. Above a critical value
 $\mathit{Ra}$  the dominant (slowest decaying or fastest growing) eigenmodes are oscillatory with waves propagating up the tail. Above a critical value   $\mathit{Ra}_c \approx 10\, 000$ , the steady solution becomes unstable to a limit cycle. A one-dimensional reduction to horizontally integrated quantities hints at a theoretical explanation for the oscillatory behaviour, but does not reproduce the loss of stability. Investigation of the initial transient at large
 $\mathit{Ra}_c \approx 10\, 000$ , the steady solution becomes unstable to a limit cycle. A one-dimensional reduction to horizontally integrated quantities hints at a theoretical explanation for the oscillatory behaviour, but does not reproduce the loss of stability. Investigation of the initial transient at large   $\mathit{Ra}$  reveals that an initially spherical thermal can rise
 $\mathit{Ra}$  reveals that an initially spherical thermal can rise   $O(100)$  times its initial diameter before approaching its final self-similar shape. The presence of a rigid horizontal floor below the thermal makes a quantitative difference of around 10 % to the rate of rise at large
 $O(100)$  times its initial diameter before approaching its final self-similar shape. The presence of a rigid horizontal floor below the thermal makes a quantitative difference of around 10 % to the rate of rise at large   $\mathit{Ra}$ .
 $\mathit{Ra}$ .