Three-dimensional time-harmonic fundamental solutions for a fluid-saturated poroelastic half-space with partially permeable free surface

2016 ◽  
Vol 40 (14) ◽  
pp. 1903-1935 ◽  
Author(s):  
Ahmad Pooladi ◽  
Mohammad Rahimian ◽  
Ronald Y. S. Pak
1998 ◽  
Vol 22 (4B) ◽  
pp. 485-499
Author(s):  
K.L. Chowdhury

The solution of the axisymmetric boundary value problem of an isotropic elastic dielectric half space subjected to charge distribution on its rigid polarization free surface is constructed by Hankel transforms. For the problem of an electric point dipole applied at origin, exact expressions for the components of displacement and polarization vectors and the potential fields are obtained in terms of Bessel function and fundamental solutions 1/R and e-mR/R, R being the distance from the source point. The electric field is determined both inside and outside the polarized region. In the particular case of a continuous electric charge distribution with density of the form l/(r2+h2)1/2, the mechanical and electric stresses on the surface of the semi-space are derived. The MathematicsTM software is used to present the numerical results on graphs depicting the variation of surface stresses for the particular charge distributions.


This paper is concerned with the study of transient response of a transversely isotropic elastic half-space under internal loadings and displacement discontinuities. Governing equations corresponding to two-dimensional and three-dimensional transient wave propagation problems are solved by using Laplace–Fourier integral transforms and Laplace−Hankel integral transforms, respectively. Explicit general solutions for displacements and stresses are presented. Thereafter boundary-value problems corresponding to internal transient loadings and transient displacement discontinuities are solved for both two-dimensional and three-dimensional problems. Explicit analytical solutions for displacements and stresses corresponding to internal loadings and displacement discontinuities are presented. Solutions corresponding to arbitrary loadings and displacement discontinuities can be obtained through the application of standard analytical procedures such as integration and Fourier expansion to the fundamental solutions presented in this article. It is shown that the transient response of a medium can be accurately computed by using a combination of numerical quadrature and a numerical Laplace inversion technique for the evaluation of integrals appearing in the analytical solutions. Comparisons with existing transient solutions for isotropic materials are presented to confirm the accuracy of the present solutions. Selected numerical results for displacements and stresses due to a buried circular patch load are presented to portray some features of the response of a transversely isotropic elastic half-space. The fundamental solutions presented in this paper can be used in the analysis of a variety of transient problems encountered in disciplines such as seismology, earthquake engineering, etc. In addition these fundamental solutions appear as the kernel functions in the boundary integral equation method and in the displacement discontinuity method.


Geophysics ◽  
1987 ◽  
Vol 52 (8) ◽  
pp. 1106-1126 ◽  
Author(s):  
Perry A. Eaton ◽  
Gerald W. Hohmann

An important element of electromagnetic (EM) prospecting is survey design; numerical modeling algorithms may be used to calculate signal‐to‐geologic‐noise (S/N) ratios to compare different survey configurations and measured responses quantitatively. Our models consist of a prismatic three‐dimensional (3-D) target in a conductive half‐space which may contain an overburden conductor; the models are energized by a time‐varying current transmitted in a loop of wire. The signal is the scattered or anomalous response of the target, while the geologic noise is either the response of the half‐space or the anomalous response of the overburden conductor. For typical loop sizes in exploration, the coincident‐loop configuration has a relatively high S/N ratio and thus a relatively high capability to resolve the target in the case of half‐space noise. Measurements made with the horizontal‐loop, moving‐coil configuration can be just as effective if the coil separation is one and one‐half to two times the depth of burial of the target and the transmitting and receiving coils are on opposite sides of the target. For coil positions on one side of the target, the S/N ratio decreases with increasing separation. The advantage in resolving power provided by the coincident loop’s superior S/N ratio diminishes as the size of the loop increases. For the case of noise due to the overburden conductor, the horizontal‐loop configuration with a large coil separation is optimal. If the depth of the target is unknown, the fixed‐loop, roving‐receiver configuration is useful for detecting the target but poor in resolving its depth because its S/N ratio is the least sensitive to the depth. With the fixed‐loop configuration, galvanic effects enhance the detectability of the target in a conductive half‐space, but inhibit detection if an overburden conductor is present. Regarding the S/N ratio, there does not appear to be any advantage in measuring the step response of a 3-D target in a conductive environment versus measuring the impulse response. The shapes of their respective S/N anomalies are essentially the same and the maximum impulse S/N ratio is 10 to 30 percent larger than the maximum step S/N ratio, though it occurs later in time by a factor of about 1.7. Although transient S/N ratios for a 3-D target in a conductive host reach a maximum value and then decrease with increasing time, harmonic S/N ratios do not necessarily reach a maximum value at an intermediate frequency. For all three survey configurations and both types of noise, target depths, and half‐space conductivities studied here, maximum transient S/N ratios are larger than harmonic S/N ratios. Peak step S/N ratios are 30 to 50 percent larger than corresponding in‐phase ratios in the case of half‐space noise, and several times larger in the case of the overburden conductor. A phase rotation of the target’s response due to the conductive host appears to amplify the quadrature S/N ratio relative to the in‐phase S/N ratio. However, in‐phase S/N ratios are always much larger than quadrature S/N ratios over the range of host resistivities used in this study.


Author(s):  
Loukas F. Kallivokas ◽  
Jacobo Bielak

Abstract This paper is concerned with the numerical solution by the finite element method of transient and time-harmonic three-dimensional acoustic scattering problems in infinite and semi-infinite domains. Its main objective is to illustrate how a local second-order surface-only infinite element — either transient or time-harmonic — developed recently for the three-dimensional wave equation in a full-space can be applied readily to scattering problems with penetrable objects near a planar free surface. Taking a problem in structural acoustics as a prototype, the combined infinite element-finite element method is used here to determine the total and scattered pressure patterns generated when a traveling plane wave impinges upon a structure of general geometry submerged in an acoustic fluid in half-space. One key feature of this methodology is that the ordinary differential equations that result from the spatial discretization maintain the symmetry and sparsity associated with problems defined only over interior domains; the resulting equations can then be solved by standard step-by-step time integration techniques. Thus, the combination of low bandwidth matrices with the ease of use of the infinite elements places the method in an ideal position to meet the large computational demands typically associated with large-scale underwater acoustics problems.


Geophysics ◽  
1984 ◽  
Vol 49 (6) ◽  
pp. 822-825 ◽  
Author(s):  
Walter L. Pilant

Plane‐wave directivity patterns for both P- and S-waves approaching a free surface are well known (Knopoff et al., 1957, Figure 3–5). These have been shown to apply in a reciprocal manner to time‐harmonic S-waves emanating from vertical and horizontal sources (Miller and Pursey, 1954; Cherry, 1962) in both two‐dimensional (2-D) and three‐dimensional (3-D) cases. Knopoff and Gilbert (1959) showed that the plane‐wave directivity patterns also apply to the first motions seen in the impulsive‐source case (3-D) and Pilant (1979, sec. 9–6) showed that they held in the equivalent 2-D problem. Theoretical expressions for these patterns are given by Pilant (ibid) as [Formula: see text] and [Formula: see text] where [Formula: see text] is measured from the vertical and the positive z-axis is into the medium. The x-axis lies along the free surface and the quantity [Formula: see text]. For angles greater than critical [Formula: see text], the proper expression for the square root is given by [Formula: see text] Thus for angles of incidence (or take‐off) greater than [Formula: see text], both [Formula: see text] and [Formula: see text] become complex numbers and lead to phase‐shift induced waveform changes as the S-waves interact with the free surface. The functions [Formula: see text] and [Formula: see text] are shown in Figure 1 for the angular range 34–37 degrees which includes the angle [Formula: see text] degrees. For this example, [Formula: see text] corresponding to a Poisson’s ratio equal to one‐quarter. The null in [Formula: see text] and the maximum in [Formula: see text] are clearly seen.


2012 ◽  
Vol 594-597 ◽  
pp. 2757-2762 ◽  
Author(s):  
Xi Luo ◽  
Xian Wei Zeng ◽  
Li Qun Tang

Based on Biot’s elastodynamic theory for poroelastic media, the dynamic response of a poroelastic half-space due to a time-harmonic concentrated vertical load applied at the free surface is investigated. Different from previous treatments of the free surface as either fully permeable or fully impermeable, the free surface of a pororelastic half-space is treated in this study as a more realistic semi-permeable boundary condition, i.e. the permeability of the free surface is considered. The governing equation for axisymmetric motion of a poroelastic half-space is solved by applying the Hankel integral transform. Numerical results are presented to show the effects of semi-permeable boundary condition on the dynamic response of poroelastic half-space.


1967 ◽  
Vol 63 (4) ◽  
pp. 1341-1367 ◽  
Author(s):  
R. D. Gregory

AbstractThis paper is principally concerned with the two-dimensional time harmonic vibrations of an elastic half-space y ≥ 0, containing a submerged cavity in the form of an infinite circular cylinder. Two sequences of line source potentials are obtained which are singular along the axis of the cylinder, satisfy the free surface conditions on y = 0, and represent outgoing waves at infinity. (The radiation condition.) It is proved that any solution of the governing equations which satisfies the free surface conditions and consists of outgoing waves at infinity is expansible as a sum of these fundamental source potentials, with coefficients to be determined from the boundary conditions on the cylinder only.The requirement of outgoing waves is carefully discussed and it is shown that the conditions taken give rise to, and are satisfied by, potential fields which would be regarded intuitively as representing outgoing waves.


Sign in / Sign up

Export Citation Format

Share Document