Permeability of vesicular Stromboli basaltic glass: Lattice Boltzmann simulations and laboratory measurements
- McGill
The permeabilities of vesicular Stromboli basaltic glasses were determined using lattice Boltzmann (LB) simulations and laboratory measurements. Lattice Boltzmann simulations were performed to simulate flow through vesicular Stromboli basaltic glasses with porosities, {phi}, from 5% to 92%. The simulations and measurements provide a power law Darcian permeability-porosity relationship k(Phi) = c ({phi}){sup 5} with c = 2.35 x 10{sup -20} from LB simulations and 5.33 x 10{sup -21} from measurements, where k is in m{sup 2}. These permeabilities of vesiculated basalts are about 1 to 2 orders of magnitude higher than in rhyolitic and dacitic volcanic rocks with the same porosity; this difference is attributed to a higher bubble interconnectivity and larger bubble apertures in our basaltic samples. The Darcian flow permeability k{sub 1} (m{sup 2}) and non-Darcian flow permeability k{sub 2} (m) are highly dependent on bubble size, D, and porosity with k{sub 1} = 7.66 x 10{sup -17}[D{sup 2{phi}3}/(1 - {phi}){sup 2}] and k{sub 2} = 2.78 x 10{sup -9}[D{phi}{sup 3}/(1 - {phi})]. Samples with power law bubble size distributions can produce higher permeabilities than samples with exponential bubble size distributions. The Darcian and non-Darcian flow regimes are delineated, demonstrating that the Darcian flow occurs at the Forchheimer number, Fo, below 0.2-1, and the transitional flow (Forchheimer flow) occurs in the Forchheimer number range 1 to 10. The correlations between friction factor, f{sub k}, and Fo are constrained by the permeability measurements, and are in good agreement with simulations: f{sub k} = (1.11 {+-} 0.17) + [(0.66 {+-} 0.39)/Fo] (measurements) and f{sub k} = (0.59 {+-} 0.49) + [(1.0 {+-} 0.01)/Fo] (LB simulations). Our results show that f{sub k} depends on k{sub 2}, pore size, and pore geometry at small Fo and tends to be a constant at large Fo. The f{sub k} - Fo correlations imply a gradual transition from Darcian to non-Darcian flow, rather than an abrupt change. Modeling the relationship between permeability created by water exsolution and depth suggests that significant increases of permeability occur at depths of {approx}100-2000 m for melts with initial water concentrations of 1-4 wt %. At these depths, for gas flow through vesicular magma with a velocity 0.1-1 m s{sup -1}, Fo is in the range {approx}0.5-47, corresponding to the transitional flow regime. For a gas flow with a velocity over {approx}10 m s{sup -1}, Fo can attain values well above the transition flow regime. Our results imply that transitional flow or turbulent flow probably prevails in vesicular magma.
- Research Organization:
- Argonne National Lab. (ANL), Argonne, IL (United States). Advanced Photon Source (APS)
- Sponsoring Organization:
- USDOE
- OSTI ID:
- 1002666
- Journal Information:
- J. Geophys. Res., Vol. 115, Issue 2010; ISSN 0148-0227
- Country of Publication:
- United States
- Language:
- ENGLISH
Similar Records
Dynamics of degassing at Kilauea Volcano, Hawaii
Equation predicts non-Darcy flow coefficient