The motion induced in a layer of Boussinesq fluid by moving periodic thermal waves is obtained by numerically solving the complete nonlinear two-dimensional momentum and temperature equations. Three sets of boundary conditions are treated: rigid upper and lower boundaries with symmetrical heating; free upper boundary and rigid lower boundary with heating only at the top; free upper and lower boundaries with symmetrical heating. The nonlinear streamline patterns show that, when the velocity fluctuations are larger than the phase speed of the thermal wave and the mean flow, the convection cells have shapes governed by fluctuating nonlinear interactions. Significant mean velocities can be created even without the characteristic tilt in the convection cells expected on the basis of linear theory. Nonlinear interactions can lead to a mean shear even in the absence of motion of the thermal source. When the viscous diffusion time across the fluid layer is less than or of the same order as the period of the thermal wave, the order of magnitude of the induced mean velocity does not exceed that of the phase speed of the wave, even for intense thermal forcing.