Shock compression of polycrystalline copper was numerically investigated by use of a discrete-element model to highlight underlying heterogeneous and nonequilibrium processes at the grain level. The average diameter of model grains was 12 \ensuremath{\mu}m. Results show highly transient vortical flow fields and strong particle velocity dispersion that are consistent with the experimental results of Mescheryakov and his associates. Characteristic times for these phenomena were on the order of acoustic propagation times across the grains. The number of vortices increased with shock strength, but their size decreased almost inversely. Ejection of copper particles from the back free surface of the specimen was also observed. The cause of ejection is grain boundary cracking.