Hostname: page-component-78c5997874-dh8gc Total loading time: 0 Render date: 2024-11-05T04:58:46.175Z Has data issue: false hasContentIssue false

Flow in deformable porous media. Part 2 Numerical analysis – the relationship between shock waves and solitary waves

Published online by Cambridge University Press:  26 April 2006

Marc Spiegelman
Affiliation:
Lamont-Doherty Geological Observatory of Columbia University, Palisades, NY 10964, USA

Abstract

Using numerical schemes, this paper demonstrates how viscous resistance to volume changes modifies the simplest shock wave solutions presented in Part 1. For an initial condition chosen to form a step-function shock, viscous resistance causes the shock to disperse into a rank-ordered wavetrain of solitary waves. Large obstructions in flux produce large-amplitude, slow-moving wavetrains while smaller shocks shed small-amplitude waves. While the viscous resistance term is initially important over a narrow boundary layer, information about obstructions in the flux can propagate over many compaction lengths through the formation of non-zero wavelength porosity waves. For large-amplitude shocks, information can actually propagate backwards relative to the matrix. The physics of dispersion is discussed and a physical argument is presented to parameterize the amplitude of the wavetrain as a function of the amplitude of the predicted shock. This quantitative relationship between the prediction of shocks and the development of solitary waves also holds when mass transfer between solid and liquid is included. Melting causes solitary waves to decrease in amplitude but the process is reversible and freezing can cause small perturbations in the fluid flux to amplify into large-amplitude waves. These model problems show that the equations governing volume changes of the matrix are inherently time dependent. Perturbations to steady-state solutions propagate as nonlinear waves and these problems demonstrate several initial conditions that do not relax to steady state. If these equations describe processes such as magma migration in the Earth, then these processes should be inherently episodic in space and time.

Type
Research Article
Copyright
© 1993 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Barcilon, V. & Lovera, O. 1989 Solitary waves in magma dynamics. J. Fluid Mech. 204, 121133.Google Scholar
Barcilon, V. & Richter, F. M. 1986 Non-linear waves in compacting media. J. Fluid Mech. 164, 429448.Google Scholar
Dodd, R., Eilbeck, J., Gibbon, J. & Morris, H. 1982 Solitons and Non-linear Wave Equations. Academic.
Drazin, P. & Johnson, R. 1989 Solitons: an Introduction. Cambridge University Press.
McKenzie, D. 1984 The generation and compaction of partially molten rock. J. Petrol. 25, 713765.Google Scholar
Olson, P. & Christensen, U. 1986 Solitary wave propagation in a fluid conduit within a viscous matrix. J. Geophys. Res. 91, 63676374.Google Scholar
Press, W., Flannery, B., Teukolsky, S. & Vertterling, W. 1986 Numerical Recipes. Cambridge University Press.
Richter, F. M. 1986 Simple models of trace element fractionation during melt segregation. Earth Planet. Sci. Lett. 77, 333344.Google Scholar
Richter, F. M. & McKenzie, D. 1984 Dynamical models for melt segregation from a deformable matrix. J. Geol. 92, 729740.Google Scholar
Scott, D. & Stevenson, D. 1984 Magma solitons. Geophys. Res. Lett. 11, 11611164.Google Scholar
Scott, D. & Stevenson, D. 1986 Magma ascent by porous flow. J. Geophys. Res. 91, 92839296.Google Scholar
Sparks, D. W. & Parmentier, E. M. 1991 Melt extraction from the mantle beneath spreading centers. Earth Planet. Sci. Lett. 105, 368377.Google Scholar
Spiegelman, M. 1989 Melting and melt migration: The physics of flow in deformable porous media. Ph.D. thesis. University of Cambridge.
Spiegelman, M. 1990 Focusing on Freezing: A new mechanism for lateral melt migration at mid-ocean ridges. EOS Trans. Am. Geophys. Union, 71, 1829.Google Scholar
Spiegelman, M. 1991 2-D or not 2-D: Understanding melt migration near a sloping, freezing boundary. EOS, Trans Am. Geophys. Union 72, 265.Google Scholar
Spiegelman, M. 1993 Flow in deformable porous media. Part 1. Simple analysis. J. Fluid Mech. 247, 1738.Google Scholar
Whitham, G. 1974 Linear and Non-linear Waves. Wiley.