J. Fluid Mech. (2016), vol. 802, pp. 750–759. c  Cambridge University Press 2016 doi:10.1017/jfm.2016.298 750 Compressible vapour flow in conduits and fractures Herbert E. Huppert 1, 2, 3, 4, † and R. Stephen Sparks 4 1 Institute of Theoretical Geophysics, Department of Applied Mathematics and Theoretical Physics, University of Cambridge, CB3 0WA, UK 2 Faculty of Science, University of Bristol, BS8 1RJ, UK 3 School of Mathematics and Statistics, University of New South Wales, Sydney, NSW 2052, Australia 4 Department of Earth Sciences, University of Bristol, BS8 1RJ, UK (Received 16 February 2015; revised 13 November 2015; accepted 22 April 2016; first published online 10 August 2016) We consider the steady flow of a viscous compressible gas through an axisymmetric or two-dimensional porous medium whose properties in the direction of the flow are sufficiently slowly varying. The study is partly motivated by a number of different applications in the Earth sciences, including the release of magmatic volatiles from a magma chamber beneath an active volcano and the discharge of geothermal fluids. The results are also relevant to evaluating the consequences of an accidental release of carbon dioxide from a storage reservoir within the Earth, as might happen at a carbon capture and storage (CCS) site. We consider both slow, thermally equilibrated, flows and fast, adiabatic flows. Because the flow is compressible, it is the mass (and not the volume) flux which is conserved along the flow. We determine this constant mass flux and the velocity and pressure fields, both of which vary with position along the flow, as a function of all the physical parameters. We find that the resultant pressure gradient in the medium is largest at the far, low-pressure end of the conduit because the velocity is largest at that end due to the smallest density being associated with the smallest pressures. This means that the pressure in the permeable conduit is always larger than the linear pressure distribution which joins the given pressures at depth and at the surface, as would be the situation if the flow were incompressible. The detailed pressure distribution is shown to depend on the variation with depth of the quantity μT /(ka 2 ), where μ is the dynamic viscosity of the vapour, T is the external temperature, k the permeability and πa 2 the cross-sectional area of the conduit. The resultant mass flux is determined to be proportional to the mean along the flow of μT /(ka 2 ). We present two numerical illustrations of the results. Key words: compressible flows, gas dynamics, geophysical and geological flows 1. Introduction There are many circumstances where fluids flow to the Earth’s surface from depth through regions of high permeability. This paper analyses the situation wherein † Email address for correspondence: heh1@cam.ac.uk available at https:/www.cambridge.org/core/terms. https://doi.org/10.1017/jfm.2016.298 Downloaded from https:/www.cambridge.org/core. Pendlebury Library of Music, on 18 May 2017 at 11:12:43, subject to the Cambridge Core terms of use,