The past few years have seen dramatic improvements in the scope and realism of star cluster simulations. Accurate treatments of stellar evolution, coupled with robust descriptions of all phases of binary evolution, have been incorporated self-consistently into several dynamical codes, allowing for the first time detailed study of the interplay between stellar dynamics and stellar physics. The coupling between evolution, dynamics, and the observational appearance of the cluster is particularly strong in young systems and those containing large numbers of primordial binary systems, and important inroads have been made in these areas, particularly in N-body simulations. I discuss some technical aspects of the current generation of N-body integrators, and describe some recent results obtained using these codes.