AbstractMagma mixing is widely recognized as a means of producing compositional diversity and preconditioning magmas for eruption. However, the processes and associated time scales that produce the commonly observed expressions of magma mixing are poorly understood, especially under crystal‐rich conditions. Here we introduce and exemplify a parameterized method to predict the characteristic mixing time of crystals in a crystal‐rich magma mush that is subject to open‐system reintrusion events. Our approach includes novel numerical simulations that resolve multiphase particle‐fluid interactions. It also quantifies the crystal mixing by calculating both the local and system‐wide progressive loss of the spatial correlation of individual crystals throughout the mixing region. Both inertial and viscous time scales for bulk mixing are introduced. Estimated mixing times are compared to natural examples and the time for basaltic mush systems to become well mixed can be on the order of 10 days.