Kr. Helfrich, THERMO-VISCOUS FINGERING OF FLOW IN A THIN GAP - A MODEL OF MAGMA FLOW IN DIKES AND FISSURES, Journal of Fluid Mechanics, 305, 1995, pp. 219-238
Flow of a fluid with a strongly temperature-dependent viscosity in a f
inite-length slot is analysed as a model of magma flow in dikes. The s
lot walls are held at a fixed temperature, thus cooling and increasing
the viscosity of the fluid as it moves along the gap. Poiseuille flow
and temperature advection, averaged across the slot, are used to stud
y the stability of this basic one-dimensional flow to lateral perturba
tions. A linear stability analysis shows that for sufficiently strong
cooling and viscosity increase with decreasing temperature, the flow i
s unstable to a fingering instability. Warm fluid is focused into rela
tively fast flowing zones and suffers only modest cooling, while cold,
slow flowing regions experience more cooling and an increase in visco
sity, which acts to locally clog the slot. The necessary condition for
instability is the presence of multiple solutions for velocity (fast,
intermediate and slow branches) in the basic one-dimensional flow. Th
e intermediate branch, where the thermal adjustment lengthscale is com
parable to the slot length, is unstable and the analysis indicates tha
t the instability continues onto the slow branch. The parametric regio
ns of instability and the growth rates are dependent on the choice of
boundary conditions at the slot entrance (i.e. the magma source): eith
er uniform flux, or uniform pressure. The latter case is the more geop
hysically realistic and has the larger unstable region and growth rate
s. Numerical solutions of the nonlinear equations show that at finite-
amplitude the hot, low-viscosity, fast-flowing fingers continue to spe
ed up, while the slow, cold regions continue to cool and slow down. At
the slot exit fluid issues from the gap in isolated hot, low-viscosit
y spouts separated by zones of cold, nearly still fluid. Application o
f the model to geophysical settings indicates that the instability is
expected for realistic parameter values. The model may help explain th
e observed focusing of fissure eruptions.