Rayleigh-Bernard convection offers a unique flow situation in which buoyancy-generated turbulence can be studied in isolation, free of the complicating influence of shear production of turbulence. The objective of this paper is to examine and model important aspects of buoyancy-generated turbulence using numerical simulation (DNS) data of Rayleigh-Bernard convection. In particular, we examine the pressure-strain and pressure-temperature gradient correlations, turbulent transport of Reynolds stress, and assumptions pertaining to algebraic modeling of thermal flux and Reynolds stress.