[en] We present for the first time, the self-consistent solution of the two-dimensional, time-dependent equations of radiation-hydrodynamics governing the accretion of matter onto the highly magnetized polar caps of luminous x-ray pulsars. The calculations show a structure in the accretion column very different from previous one-zone uniform models. We have included all the relevant magnetic field corrections to both the hydrodynamics and the radiative transport. We include a new theory for the diffusion and advection of both radiation energy density and photon number density. For initially uniformly accreting models with super-Eddington flows, we have uncovered evidence of strong radiation-driven outflowing optically thin radiation filled regions of the accretion column embedded in optically-thick inflowing plasma. The development of these photon ''bubbles'' have growth times on the order of a millisecond and show fluctuations on sub-millisecond timescales. The photon bubbles are likely to be a consequence of convective over-stability and may result in observable fluctuations in the emitted luminosity leading to luminosity dependent changes in the pulse profile. This may provide important new diagnostics for conditions in accreting x-ray pulsars. 13 refs., 18 figs