LAUSR

Abstract Silicic magma bodies are formed by repeated injections of mobile magma and reside as a crystal-rich mush. Numerical studies of open-system events have revealed the complexity of mixing and rheological behavior. This is associated with the dilation of the crystal network and the possible occurrence of a lubricated regime. Lubrication forces are hydrodynamic interactions occurring when neighboring crystals have relative motion. The effect of such dissipative forces has not yet been explored in the case of magmatic mush. Here, we investigate the effects of lubrication on mush dynamics and on magma transport. First, we propose scaling relationships to assess the relative importance of the forces controlling the motion of one crystal within a mush by adding lubrication terms into the Basset-Boussinesq-Oseen equation that describes crystal motion in a viscous melt. We then investigate lubrication effects at the macroscopic scale with computational fluid dynamics with discrete element modeling (CFD-DEM) simulations that include these forces. We explore two cases: crystal mush sedimentation and the injection of a crystal-free magma inside a mush. We perform all simulations twice, with and without lubrication forces, and compare the results. At the grain scale, we show that three dimensionless numbers and the crystal content can describe the competition between viscous drag, buoyancy, and lubrication. Two of these numbers (Stokes and Froude numbers) have been previously employed in the context of dilute suspensions. The third is a new form of the Sommerfeld number that measures the importance of lubrication. At the macroscopic scale, simulation pairs (with and without lubrication forces) exhibit very similar behavior when in steady state. The duration of the transient regime preceding steady state, however, is increased when lubrication forces are included. Lubrication causes an apparent bulk strain hardening followed by softening at the initiation of the mush motion. Our results show that lubrication opposes dilation and the initiation of motion within the magmatic mush during this transient phase. Our results highlight the control that the crystal network exerts on magma transport and provide a novel way to evaluate when lubrication matters.