A dynamical model has been developed to study the evolution of a quark
-gluon system towards equilibrium with the gluons equilibrating prior
to the quarks and the quarks executing random Brownian motion in the g
luonic heat bath. We estimate the thermalization times for various qua
rk flavors and heavy flavor production rates within this model.