I am reporting on our team's progress in investigating fundamental properties of convective shells in the deep stellar interior during advanced stages of stellar evolution. We have performed a series of 3D hydrodynamic simulations of convection in conditions similar to those in the O-shell burning phase of massive stars. We focus on characterizing the convective boundary and the mixing of material across this boundary. Results from 7683 and 15363 grids are encouragingly similar (typically within 20%). Several global quantities, including the rate of mass entrainment at the convective boundary and the driving luminosity, are related by scaling laws. We investigate the effect of several of our assumptions, including the treatment of the nuclear burning driving the convection or that of neutrino cooling. The burning of the entrained material from above the convection zone could have important implications for pre-supernova nucleosynthesis.