A simple mesoscale model was developed for discontinuous dynamic recrystallization. The material is described on a grain scale as a set of (variable) spherical grains. Each grain is characterized by two internal variables: its diameter and dislocation density (assumed homogeneous within the grain). Each grain is then considered in turn as an inclusion, embedded in a homogeneous equivalent matrix, the properties of which are obtained by averaging over all the grains. The model includes: (i) a grain boundary migration equation driving the evolution of grain size via the mobility of grain boundaries, which is coupled with (ii) a dislocation-density evolution equation, such as the Yoshie–Laasraoui–Jonas or Kocks–Mecking relationship, involving strain hardening and dynamic recovery, and (iii) an equation governing the total number of grains in the system due to the nucleation of new grains. The model can be used to predict transient and steady-state flow stresses, recrystallized fractions, and grain-size distributions. A method to fit the model coefficients is also described. The application of the model to pure Ni is presented.