Permeability and melt flow in simulated and natural partially molten basaltic magmas

Abstract A Lattice–Boltzmann method is used to calculate grain-scale, low Reynolds number, single phase, interstitial porous melt flow in two independent igneous microstructures: 1) a microstructure simulated using a stochastic algorithm for progressive crystallization, and 2) a natural microstructure obtained through X-ray Computed Tomography (CT) of a partially melted basalt. In both cases, the error in the calculated flow field due to finite discretization and domain size has been estimated or demonstrated to be insignificant. Visually, the porous melt flow tends to localize into high flux channels, especially with increasing crystallinity, and this impression is quantitatively confirmed by increasing skewness of the velocity distribution in the direction parallel to the imposed pressure gradient driving the flow. A change from uniform to localized melt flow in naturally occurring situations may have a profound effect on the distribution of trace elements and, if the melt is reactive, the chemical and structural evolution of the igneous microstructure through localized phase change. Permeabilities of both microstructures are determined from the calculated steady flow field. The permeabilities are then fitted to two different correlation models which are based on the Rumpf–Gupte and Carman–Kozeny relations. In these models, we express permeability as functions of melt fraction and either the mean crystal length or specific surface area. Extrapolated permeability to higher melt fractions using both correlation models for the partially melted basalt is shown to be within a factor of four of experimentally determined permeability of a similar sample. We determine permeability estimates in the relatively unexplored melt fraction range of 30%–80% and find consistency with previous work for permeability estimates in the melt fraction range of 20%–30%.

[1]  Nikolaus von Bargen,et al.  Permeabilities, interfacial areas and curvatures of partially molten systems: Results of numerical computations of equilibrium microstructures , 1986 .

[2]  Bryant,et al.  Prediction of relative permeability in simple porous media. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[3]  Q. Zou,et al.  On pressure and velocity boundary conditions for the lattice Boltzmann BGK model , 1995, comp-gas/9611001.

[4]  Skordos,et al.  Initial and boundary conditions for the lattice Boltzmann method. , 1993, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[5]  D. Kohlstedt,et al.  Stress‐driven melt segregation in partially molten rocks , 2003 .

[6]  B. Marsh Solidification fronts and magmatic evolution , 1996, Mineralogical Magazine.

[7]  D. Rothman,et al.  Two-fluid flow in sedimentary rock: simulation, transport and complexity , 1997, Journal of Fluid Mechanics.

[8]  Matthaeus,et al.  Recovery of the Navier-Stokes equations using a lattice-gas Boltzmann method. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[9]  Martin O. Saar,et al.  Numerical models of the onset of yield strength in crystal–melt suspensions , 2001 .

[10]  O. Bachmann,et al.  On the Origin of Crystal-poor Rhyolites: Extracted from Batholithic Crystal Mushes , 2004 .

[11]  R. Maier,et al.  Pore-Scale Flow and Dispersion , 1998 .

[12]  Qinjun Kang,et al.  Simulation of dissolution and precipitation in porous media , 2003 .

[13]  William J. Bosl,et al.  A study of porosity and permeability using a lattice Boltzmann simulation , 1998 .

[14]  W. Carlson,et al.  Plagioclase-chain networks in slowly cooled basaltic magma , 1999 .

[15]  F. Dullien Porous Media: Fluid Transport and Pore Structure , 1979 .

[16]  Robert S. Bernard,et al.  Boundary conditions for the lattice Boltzmann method , 1996 .

[17]  M. Avrami Kinetics of Phase Change. II Transformation‐Time Relations for Random Distribution of Nuclei , 1940 .

[18]  Daniel H. Rothman,et al.  MACROSCOPIC MANIFESTATIONS OF MICROSCOPIC FLOWS THROUGH POROUS MEDIA: Phenomenology from Simulation , 1996 .

[19]  Hans Rumpf,et al.  Einflüsse der Porosität und Korngrößenverteilung im Widerstandsgesetz der Porenströmung , 1971 .

[20]  Daniel H. Rothman,et al.  Cellular‐automaton fluids: A model for flow in porous media , 1988 .

[21]  B. Marsh,et al.  Igneous microstructures from kinetic models of crystallization , 2006 .

[22]  D. McKenzie,et al.  Percolation threshold and permeability of crystallizing igneous rocks: The importance of textural equilibrium , 2004 .

[23]  A. Philpotts,et al.  PHYSICAL PROPERTIES OF PARTLY MELTED THOLEIITIC BASALT , 1996 .

[24]  S. Maaløe,et al.  The permeability controlled accumulation of primary magma , 1982 .

[25]  G. Taylor Dispersion of soluble matter in solvent flowing slowly through a tube , 1953, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[26]  E. Bruce Watson,et al.  Grain-scale permeabilities of texturally equilibrated, monomineralic rocks , 1998 .

[27]  P. Carman,et al.  Flow of gases through porous media , 1956 .

[28]  B. Marsh On bimodal differentiation by solidification front instability in basaltic magmas, part 1: basic mechanics , 2002 .

[29]  E. Watson,et al.  Nonlinear pressure diffusion in a porous medium: Approximate solutions with applications to permeability measurements using transient pulse decay method , 2001 .

[30]  Charles A. Williams,et al.  Reassessment of pore shapes in microstructurally equilibrated rocks, with implications for permeability of the upper mantle , 2003 .

[31]  Roland Glantz,et al.  Calibration of a Pore-Network Model by a Pore-Morphological Analysis , 2003 .

[32]  U. Faul Permeability of partially molten upper mantle rocks from experiments and percolation theory , 1997 .

[33]  Dominique d'Humières,et al.  Multireflection boundary conditions for lattice Boltzmann models. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[34]  M. Avrami,et al.  Kinetics of Phase Change 2 , 1940 .

[35]  D. Shirley Compaction of Igneous Cumulates , 1986, The Journal of Geology.

[36]  C. Pan,et al.  Pore-scale modeling of saturated permeabilities in random sphere packings. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.