Buoyancy-driven crack propagation: A mechanism for magma migration. (English) Zbl 0602.73112
A solution has been obtained for steady propagation of a two-dimensional fluid fracture driven by buoyancy in an elastic medium. The problem is formulated in terms of an integro-differential equation governing the elastic deformation, coupled with the differential equation of lubrication theory for viscous flow in the crack.
The numerical treatment of this system is carried out in terms of an eigenfunction expansion of the cavity shape, in which the coefficients are found by use of a nonlinear constrained optimization technique. When suitably non-dimensionalized, the solution appears to be unique. It exhibits a semi-infinite crack of constant width following the propagating fracture. For each value of the stress intensity factor of the medium, the width and propagation speed are determined. The results are applied to the problem of the vertical ascent of magma through the earth’s mantle and crust. Values obtained for the crack width and ascent velocity are in accord with observations. This mechanism can explain the high ascent velocities required to quench diamonds during a Kimberlite eruption. The mechanism can also explain how basaltic eruptions can carry large mantle rocks (xenoliths) to the surface.
The numerical treatment of this system is carried out in terms of an eigenfunction expansion of the cavity shape, in which the coefficients are found by use of a nonlinear constrained optimization technique. When suitably non-dimensionalized, the solution appears to be unique. It exhibits a semi-infinite crack of constant width following the propagating fracture. For each value of the stress intensity factor of the medium, the width and propagation speed are determined. The results are applied to the problem of the vertical ascent of magma through the earth’s mantle and crust. Values obtained for the crack width and ascent velocity are in accord with observations. This mechanism can explain the high ascent velocities required to quench diamonds during a Kimberlite eruption. The mechanism can also explain how basaltic eruptions can carry large mantle rocks (xenoliths) to the surface.
MSC:
74L05 | Geophysical solid mechanics |
74R05 | Brittle damage |
76R99 | Diffusion and convection |
74S30 | Other numerical methods in solid mechanics (MSC2010) |
76T99 | Multiphase and multicomponent flows |
74G70 | Stress concentrations, singularities in solid mechanics |
65R20 | Numerical methods for integral equations |
86A99 | Geophysics |
Keywords:
buoyancy-driven; Chebychev polynomials; two-dimensional fluid fracture; integro-differential equation; elastic deformation; differential equation of lubrication theory; viscous flow; eigenfunction expansion of the cavity shape; nonlinear constrained optimization technique; semi-infinite crack of constant width; stress intensity factor; propagation speed; vertical ascent of magmaCitations:
Zbl 0568.76100References:
[1] | Carmichael, Phil. Trans. R. Soc. Lond. 286 pp 373– (1977) |
[2] | Atkinson, J. Geophys. Res. 89 pp 4077– (1984) |
[3] | Anderson, Rev. Geophys. Space Phys. 15 pp 77– (1977) |
[4] | Weertman, J. Geophys. Res. 76 pp 1171– (1971) |
[5] | Spence, J. Geophys. Res. 90 pp 575– (1985) |
[6] | Spence, Proc. R. Soc. Lond. 400 pp 289– (1985) |
[7] | Secor, Geophys. Res. Lett. 2 pp 510– (1975) |
[8] | DOI: 10.1016/0148-9062(77)90740-9 · doi:10.1016/0148-9062(77)90740-9 |
[9] | Pollard, J. Geophys. Res. 81 pp 975– (1976) |
[10] | Pollard, Geophys. Res. Lett. 3 pp 513– (1976) |
[11] | DOI: 10.1016/0040-1951(73)90021-8 · doi:10.1016/0040-1951(73)90021-8 |
[12] | Persikov, Dokl. Akad. Nauk SSSR 260 pp 47– (1981) |
[13] | Muller, Proc. Lunar Planet. Sci. Conf. 11 pp 1979– (1980) |
[14] | DOI: 10.1146/annurev.ea.12.050184.001025 · doi:10.1146/annurev.ea.12.050184.001025 |
[15] | DOI: 10.1016/0031-9201(86)90004-X · doi:10.1016/0031-9201(86)90004-X |
[16] | Delaney, Am. J. Sci. 282 pp 856– (1982) |
[17] | Delaney, U.S. Geol. Survey Prof. Paper 286 pp 373– (1981) |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.