Numerical and Soft Computing Methods for Characteristic Value Problems of ODE and ODEs Systems
View this Special IssueResearch Article  Open Access
Sudden Pressurization of a Spherical Cavity in a Poroelastic Medium
Abstract
Governing equations of poroelastodynamics in time and frequency domain are derived. The continuity equation complements the momentum balance equations. After reduction for spherical symmetry (geometry and loading), the governing equations in frequency domain are solved by introducing wave potentials. The wave propagation velocities are obtained as the real parts of the characteristic equation of the coupled ODE system. Time domain solution for Dirac type boundary pressure is obtained through numerical inversion of transformed solutions. The results are compared to the solution in classical elasticity theory found in the literature.
1. Introduction
In wave mechanics, cavity expansion problems (circular and spherical) have received considerable attention since early 1930s [1]. These problems generally are amenable to exact solution, and the analytic solutions cast new light to the nature of wave propagation in solid media. The problem has practical applications in geophysics, seismology, and tunnel and mining engineering, as earthquake sources, underground detonation, and seismic probing. The exact solutions of such simple problems serve to understand more complex wave motions; these solutions can also be used as benchmark problems to assess the accuracy of numerical methods (FEM, BEM, FDM, etc.).
The literature on cavity problems in elastodynamics abounds [2–6], while corresponding work in poroelasticity has only recently been emerging [7]. An analytical solution for a dynamically loaded poroelastic column [8, 9] is also available.
This work concerns frequency domain analytical solution of dynamic pressurization of a spherical cavity in an infinite poroelastic medium. Time domain solution for Dirac type boundary pressure (pressure suddenly applied and removed) is also investigated using FFT (fast Fourier transform) technique. The developed solutions have parallel applications as the classical elastic cavity problem and clearly show the existence of a second pressure wave, the socalled slow wave.
In the sequel, a short review of Biot’s poroelasticity theory is given. The governing equations are reduced for spherical symmetry. The analytical solution in frequency domain is obtained for a permeable boundary with arbitrary but spherically symmetric pressure inside the cavity. Time domain solution is obtained via a numerical inversion using FFT.
1.1. Biot’s Theory of Poroelasticity
Unlike the classical elasticity, Biot’s theory of poroelasticity is a coupled defromation flow theory of a porous solid matrix with interstitial fluid. Biot introduced his linear quasistatic theory in 1941 [10] and later extended it to cover the dynamic range [11, 12]. An extensive review of quasistatic poroelasticity can be found in [13, 14].
The constitutive equations of linearisotropic poroelasticity are given as (summation convention applies)where are the strains in the solid , and are components of solid displacement vector and total stress tensor, is the fluid pressure, is the variation of fluid volume per unit reference volume, and is the Kronecker delta. Here, tensile and are positive, while pore pressure is positive when it is compressive. Four elastic constants, : drained Lame’s modulus, as defined in classical elasticity (dimension ≡ F/L2); : drained shear modulus, as defined in classical elasticity (dimension ≡ F/L2); : Biot’s effective stress coefficient (dimensionless); : Biot deformation modulus, corresponds to the reciprocal of constrained storage coefficient in hydrogeology (dimension ≡ F/L2),characterize the poroelastic material.
The balance of linear momentum for the solid and fluid parts gives, respectively, Here, are relative fluid displacements, and are the components of the fluid flux vector, the three densities , , and are the bulk density, fluid density, and added fluid mass density, respectively. is the porosity of the porous matrix. is the coefficient of permeability in units L^{3}T/M. , are the body forces per unit volume of solid and fluid, respectively. The third equation is sometimes unconventionally referred to as generalized Darcy’s law. The last balance law is the continuity equation for the fluid part where is fluid generation rate per unit volume.
In the absence of body forces and fluid generation, Fourier transforms of (1a), (1b), (2), (3), and (4) with respect to time give, when combined, the 3D governing equations in frequency domain It has been assumed that the medium is initially at rest. Here, , represent the Fourier transformed variables, that is, and is a coefficient defined as In dyadic form (5) reads where is the solid displacement vector. Analytical solutions of the 3D governing equations are scarse if not impossible. There are the point load solutions for infinite poroelastic medium, the fundamental solutions, for both quasistatic and dynamic poroelasticity [15–18] and poroelastic Lamb’s problem [19] available in the literature. Also a good compendium of anlytical solutions in the quasistatic range can be found in [14].
2. Governing Equations for Spherical Symmetry
Here, a spherical cavity of radius “” in an infinite PE medium (Figure 1) will be considered. By substituting the gradient and divergence in (8), governing equations in spherical coordinates can be obtained.
Let , , be the spherical components of solid displacements, because the cavity is spherically symmetric, and we assume a time varying but spherically symmetric pressure inside; waves emanating from such a source will have spherical symmetry, that is, Hence, the governing equations of 3D poroelasticity in FTS, in this case, reduce to the following: Here, denotes the only nonzero displacement for simplicity, that is, radial component .
The stress components are related to the radial displacement as Other .
3. General Solution of the Governing Equations
If one introduces a displacement potential such that then (10) and (11) become where Now, we assume solutions of the form where is a frequency dependent parameter. Substituting (14) in (12), one obtains the eigenvalue problem: Hence, are the roots of the characteristic equation With the following abbreviations: the roots of the characteristic equation are and the solution becomes It is clear that the terms represent waves propagating inwards. In an infinite medium, no waves propagate from infinity (Sommerfeld’s radiation condition [1]); discarding these terms, (19) becomes Moreover, the amplitudes , are related through (15) as introducing (21) in (20) Furthermore, since and Finally, the general solution for radial displacement and pore pressure becomes The integration constants and are to be evaluated from boundary conditions at .
3.1. Boundary Conditions
Among various possible combinations of traction, displacement, pore pressure, and fluid flux, only the permeable boundary condition will be considered here, that is, The stress boundary condition corresponds to a loading of in time domain, where is the strength of loading and is the Dirac delta function (or distribution). When (25) are introduced in (24) one obtains where
3.2. Analytical Solution in Frequency Domain
Thus, the complete analytical solutions in frequency domain for Dirac type pulse are Stresses can be found by inserting (28) in (11). The two terms in each of (28) represent two phases with different propagation characteristics. The reciprocal of real part of and corresponds to fast and slow (second) longitudinal wave velocities. The second (slow) wave is highly dispersive and highly damped.
3.3. Solution in the Classical Theory of Elasticity
The elastodynamic solution in frequency domain for suddenly pressurized spherical cavity can be derived following a similar outline and is found in [1, 20]; here we restate the form of the solution, which complies with the Fourier transform definition in (6): where is the displacement, , are the  and wave velocities in an idealized elastic medium, and other parameters are as defined previously.
4. Results in Frequency and Time Domain
In all of the computations material constants (Table 1) for Berea Sandstone are used [21].

The material parameters for corresponding elastic medium are as in Table 2.

The real and complex parts of displacement and pore pressure solutions in frequency domain are plotted in Figures 2 and 3, respectively, note that the elastic theory doesn’t involve pore pressure.
(a)
(b)
(a)
(b)
The following (Figure 4) are the dispersion curves for the fast and slow wave. As seen in the figure, the wave front of fast wave travels with a speed equal to 3137.2 m/s where as the slow wave front velocity is 1037 m/s. While, both waves display dispersion, the slow wave is rather more dispersive than the fast wave.
(a)
(b)
4.1. Time Domain Solution Using FFT Technique
The frequency domain analytical solutions (28) are extremely complex (if not impossible) to be inverse transformed (by analytical means) into time domain. Thus, time domain solutions are obtained using numerically inverting (28) by FFT (fast Fourier transform) technique.
Since Dirac impulse is highly transient, a high frequency resolution is required, thus the following parameters are chosen for FFT: seconds (analysis interval), (exponent in FFT),
which give a frequency resolution of 0.628 rad/s and a Nyquist frequency of 41177 rad/sec. The results are plotted in Figure 5.
As seen from Figure 5, the fast primary wave in poroelastic medium travels faster (arrives earlier) than wave in elastic medium, with the amplitude of elastic wave being slightly larger. The second slow wave in poroelastic wave cannot have been captured, probably because of its small amplitude and high damping and dispersion.
5. Conclusions
Analytical solution of waves propagating from the surface of a suddenly pressurized poroelastic spherical cavity with permeable boundary is obtained in frequency domain. The problem is one of the classical problems in wave propagation in elastic solids [1, 20]. Among many other possibilities of boundary conditions, this paper, which employs permeable boundary conditions, is one of the first attempts to the solution of a variety of other problems related to spherical wave propagation from cavities in poroelastic media.
Time domain solution is obtained with inverse FFT. Since, Dirac type pulse is a highly transient excitation, a high frequency resolution is necessary. The results clearly show that waves in poroelastic media propagate faster, with slightly lower amplitude, as compared to elastic media with similar material properties. The slow wave could not have been visualized in the graphics, because of small amplitude, high damping, and dispersion characteristics of this phase.
References
 A. C. Eringen and E. uhubi, Elastodynamics: Volume II, Linear Theory, Academic Press, New York, NY, USA, 1975.
 H. L. Selberg, “Transient compression waves from spherical and cylindrical cavities,” Arkiv för Fysik, vol. 5, pp. 97–108, 1952. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 F. Blake, “Spherical wave propagation in solid media,” Journal of the Acoustical Society of America, vol. 24, no. 2, pp. 211–215, 1952. View at: Publisher Site  Google Scholar
 W. Goldsmith and W. A. Allen, “Graphical representation of the spherical propagation of explosive pulses in elastic media,,” Journal of the Acoustical Society of America, vol. 27, no. 1, pp. 47–55, 1955. View at: Publisher Site  Google Scholar
 A. C. Eringen, “Elastodynamic problem concerning the spherical cavity,” Quarterly Journal of Mechanics and Applied Mathematics, vol. 10, no. 3, pp. 257–270, 1957. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. G. Hopkins, “Dynamic expansion of spherical cavities in metals,” in Progress in Solid Mechanics, I. N. Sneddon and R. Hill, Eds., vol. 1, chapter 3, pp. 83–164, NorthHolland, Amsterdam, The Netherlands, 1960. View at: Google Scholar  MathSciNet
 T. Senjuntichai and R. K. N. D. Rajapakse, “Transient response of a circular cavity in a poroelastic medium,” International Journal for Numerical & Analytical Methods in Geomechanics, vol. 17, no. 6, pp. 357–383, 1993. View at: Google Scholar
 M. Schanz and A. H. D. Cheng, “Transient wave propagation in a onedimensional poroelastic column,” Acta Mechanica, vol. 145, no. 1–4, pp. 1–18, 2000. View at: Publisher Site  Google Scholar
 M. Schanz and A. H. D. Cheng, “Dynamic analysis of a onedimensional poroviscoelastic column,” Journal of Applied Mechanics, vol. 68, no. 2, pp. 192–198, 2001. View at: Google Scholar
 M. A. Biot, “General theory of threedimensional consolidation,” Journal of Applied Physics, vol. 12, no. 2, pp. 155–164, 1941. View at: Publisher Site  Google Scholar
 M. A. Biot, “Theory of propagation of elastic waves in a fluidsaturated porous solid. I. Lowfrequency range,” Journal of the Acoustical Society of America, vol. 28, pp. 168–178, 1956. View at: Publisher Site  Google Scholar  MathSciNet
 M. A. Biot, “Theory of propagation of elastic waves in a fluidsaturated porous solid. II. Higher frequency range,” Journal of the Acoustical Society of America, vol. 28, pp. 179–191, 1956. View at: Publisher Site  Google Scholar  MathSciNet
 E. Detournay and A. H. D. Chang, “Fundamentals of Poroelasticity,” in Comprehensive Rock Engineering: Volume II Principles, Practice and Projects, chapter 5, pp. 113–171, Pergamon Press, New York, NY, USA, 1993. View at: Google Scholar
 H. F. Wang, Theory of Linear Poroelasticity with Applications to Geomechanics and Hydrogeology, Princeton University Press, Princeton, NJ, USA, 2000.
 J. Chen, “Time domain fundamental solution to Biot's complete equations of dynamic poroelasticity—part I: twodimensional solution,” International Journal of Solids and Structures, vol. 31, no. 10, pp. 1447–1490, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Chen, “Time domain fundamental solution to biot's complete equations of dynamic poroelasticity—part II: threedimensional solution,” International Journal of Solids and Structures, vol. 31, no. 2, pp. 169–202, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. H. D. Cheng, T. Badmus, and D. E. Beskos, “Integral equation for dynamic poroelasticity in frequency domain with BEM solution,” Journal of Engineering Mechanics, vol. 117, no. 5, pp. 1136–1157, 1991. View at: Publisher Site  Google Scholar
 M. Schanz and D. Pryl, “Dynamic fundamental solutions for compressible and incompressible modeled poroelastic continua,” International Journal of Solids and Structures, vol. 41, no. 15, pp. 4047–4073, 2004. View at: Publisher Site  Google Scholar
 A. J. Philippacopoulos, “Lamb’s problem for fluid saturated porous media,” Bulletin of the Seismological Society of America, vol. 78, no. 2, pp. 908–923, 1988. View at: Google Scholar
 K. F. Graff, Wave Motion in Elastic Solids, Dover, New York, NY, USA, 1975.
 M. Schanz, Wave Propagation in Viscoelastic and Poroelastic Continua: A Boundary Element Approach, Lecture Notes in Applied and Computational Mechanics, Springer, Heidelberg, Germany, 2010.
Copyright
Copyright © 2013 Mehmet Ozyazicioglu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.