The relaxation of two-dimensional rolls in Rayleigh-Benard convection

Large aspect ratio, two-dimensional, periodic convection layers containing a Boussinesq fluid of finite Prandtl number bounded by rigid or free horizontal surfaces are investigated numerically. The fluid equations are solved using both a standard pseudospectral and a Fourier integral method for the time evolution of finite initial perturbations, both random thermal perturbations and localized roll disturbances, into a final equilibrium state. The suggestion that a Fourier integral solution method is required to yield roll relaxation, the two-dimensional process increasing the convection wavelength to values larger than critical, is investigated. Roll relaxation is found for both free-slip and no-slip surfaces using either solution method as long as the initial state is chosen to be of the form of a localized roll disturbance. A wide variety of simulations are performed and roll relaxation is found to be independent of the periodic domain length, weakly dependent on the Rayleigh number and dependent upon the magnitude of the initial localized roll disturbances.