  The Rayleigh (Ra) and Prandtl (Pr) number scaling of the Nusselt number Nu,
the Reynolds number Re, the temperature fluctuations, and the kinetic and
thermal dissipation rates is studied for (numerical) homogeneous
Rayleigh-Benard turbulence, i.e., Rayleigh-Benard turbulence with periodic
boundary conditions in all directions and a volume forcing of the temperature
field by a mean gradient. This system serves as model system for the bulk of
Rayleigh-Benard flow and therefore as model for the so called ``ultimate regime
of thermal convection''. With respect to the Ra dependence of Nu and Re we
confirm our earlier results \cite{loh03} which are consistent with the
Kraichnan theory \cite{kra62} and the Grossmann-Lohse (GL) theory
\cite{gro00,gro01,gro02,gro04}, which both predict $Nu \sim Ra^{1/2}$ and $Re
\sim Ra^{1/2}$. However the Pr dependence within these two theories is
different. Here we show that the numerical data are consistent with the GL
theory $Nu \sim Pr^{1/2}$, $Re \sim Pr^{-1/2}$. For the thermal and kinetic
dissipation rates we find $\eps_\theta/(\kappa \Delta^{2}L^{-2}) \sim (Re
Pr)^{0.87}$ and $\eps_u/(\nu^3 L^{-4}) \sim Re^{2.77}$, also both consistent
with the GL theory, whereas the temperature fluctuations do not depend on Ra
and Pr. Finally, the dynamics of the heat transport is studied and put into the
context of a recent theoretical finding by Doering et al. \cite{doe05}.
