In this work, we present a variational multiscale model for grain growth in face-centered cubic nanocrystalline (nc) metals. In particular, grain-growth-induced stress softening and the resulting relaxation phenomena are addressed. The behavior of the polycrystal is described by a conventional Taylor-type averaging scheme in which the grains are treated as two-phase composites consisting of a grain interior phase and a grain boundary-affected zone. Furthermore, a grain-growth law that captures the experimentally observed characteristics of the grain coarsening phenomena is proposed. To this end, the grain size is not taken as constant and varies according to the proposed stress-driven growth law. Several parametric studies are conducted to emphasize the influence of the grain-growth rule on the overall macroscopic response. Finally, the model is shown to provide a good description of the experimentally observed grain-growth-induced relaxation in nc-copper.