The long-time behaviour of numerical approximations to the solutions of a s
emilinear parabolic equation undergoing a Hopf bifurcation is studied in th
is paper. The framework includes reaction-diffusion and incompressible Navi
er-Stokes equations. It is shown that the phase portrait of a supercritical
Hopf bifurcation is correctly represented by Runge-Kutta time discretizati
on. In particular, the bifurcation point and the Hopf orbits are approximat
ed with higher order. A basic tool in the analysis is the reduction of the
dynamics to a two-dimensional center manifold. A large portion of the paper
is therefore concerned with studying center manifolds of the discretizatio
n.