Simulation procedures for typical Johnson-Mehl crystals generated under various models for random nucleation are proposed. These procedures include algorithms for simulating spatio-time-inhomogeneous Poisson processes. Empirical results for a particular class of Johnson-Mehl tessellations in two and three dimensions show remarkably different crystals.