We present the first post core collapse models of initially rotating star clusters, using the numerical solution of an orbit-averaged 2D Fokker-Planck equation. Based on the code developed by Einsel & Spurzem (1999), we have improved the speed and the stability and included the steady three-body binary heating source. We have confirmed that rotating clusters, whether they are in a tidal field or not, evolve significantly faster than non-rotating ones. Consequences for observed shapes, density distribution, and kinematic properties of young and old star clusters are discussed. The results are compared with gaseous and 1D Fokker-Planck models in the non-rotating case.

1- Simulations of star clusters and galactic nuclei: rotation and other challenges (R. Spurzem & C. Einsel) 2- Relaxation of elliptical galaxies (C. Boily) 3- Simulations with an Hybrid Nbody-Scf code (M. Hemsendorf) These are three two-page contributions to the Rutgers galaxy dynamics held in august, 1998. Each will be included in the proceedings to be published by the Astronomical Society of the Pacific.

The influence of rotation on the dynamical evolution of collisional stellar systems is investigated by solving the orbit-averaged Fokker-Planck equation in (E,J_z)-space. We find that large amounts of initial rotation drive the system into a phase of strong mass loss while it is moderately contracting. The core is rotating even faster than before although angular momentum is transported outwards. At the same time the core is heating. Given these features this phase can be associated with the gravo-gyro 'catastrophe' found by Hachisu (1979). The increase in central angular momentum levels off after about 2-3 initial half-mass relaxation times indicating that the source of this `catastrophe' is depleted. Finally, the central angular velocity increases again, but with a rather small power of the central density -- the same power as for the central velocity dispersion during self-similar contraction towards the collapse. The rotation curves flatten and the ellipticity variations decrease with time, but their shapes are very similar. These results suggest the existence of a self-similar solution for a rotating cluster as well. The maximum values of rotational velocity and ellipticity occur at about the half mass radius.