Dehydration melting of crustal rocks may commonly occur in response to the intrusion of mafic magma in the mid- or lower crust. However, the relative importance of melt buoyancy, shear or dyking in melt generation and extraction under geologically relevant conditions is not well understood. A numerical model of the partial melting of a metapelite is presented and the model results are compared with the Ivrea-Verbano Zone in northern Italy. The numerical model uses the mixture theory approach to modelling simultaneous convection and phase change and includes special ramping and switching functions to accommodate the rheology of crystal-melt mixtures in accordance with the results of deformation experiments. The model explicitly includes both porous media flow and thermally and compositionally driven bulk convection of a restitecharged melt mass. A range of melt viscosity and critical melt fraction models is considered. General agreement was found between predicted positions of isopleths and those from the Ivrea-Verbano Zone. Maximum melt velocities in the region of porous flow are found to be 1 × 10−7 and 1 × 10−1m per year in the region of viscous flow. The results indicate that melt buoyancy alone may not be a sufficient agent for melt extraction and that extensive, vigorous convection of partially molten rocks above mafic bodies is unlikely, in accord with direct geological examples.