By tailoring the geometry of the upper boundary in turbulent Rayleigh-Benard convection we manipulate the boundary layer-interior flow interaction, and examine the heat transport using the lattice Boltzmann method. For fixed amplitude and varying boundary wavelength., we find that the exponent beta in the Nusselt-Rayleigh scaling relation, Nu - 1 proportional to Ra-beta, is maximized at lambda =lambda(max) approximate to ( 2 pi)(-1), but decays to the planar value in both the large (lambda >> lambda(max)) and small (lambda << lambda(max)) wavelength limits.