Single-scattering parabolic equation solutions for elastic media propagation, including Rayleigh waves

The parabolic equation method with a single-scattering correction allows for accurate modeling of range-dependent environments in elastic layered media. For problems with large contrasts, accuracy and efficiency are gained by subdividing vertical interfaces into a series of two or more single-scattering problems. This approach generates several computational parameters, such as the number of interface slices, an iteration convergence parameter τ, and the number of iterations n for convergence. Using a narrow-angle approximation, the choices of n=1 and τ=2 give accurate solutions. Analogous results from the narrow-angle approximation extend to environments with larger variations when slices are used as needed at vertical interfaces. The approach is applied to a generic ocean waveguide that includes the generation of a Rayleigh interface wave. Results are presented in both frequency and time domains.     

Leaky Stoneley waves at a moving interphase boundary

The possibility of the existence of a leaky Stoneley wave at a uniformly moving interphase boundary defined as a jump in the acoustic parameters of an isotropic elastic medium is discussed. It is shown that the motion of the interphase boundary exerts dissimilar effects on the orientation of the wave normals of partial waves forming the Stoneley wave, and this results in different Doppler shifts of the frequencies of the partial waves in the laboratory frame of reference.     

Excitation of Rayleigh and Stoneley surface acoustic waves by distributed seismic sources

Using the integral Fourier-transform technique, we obtain a solution in integral form to the problem of excitation of elastic waves in a homogeneous isotropic solid half-space and the bordering homogeneous gas by the time-dependent forces which are arbitrarily distributed in a solid over the plane parallel to the interface of the media. Different configurations of the force sources are analyzed from the viewpoint of excitation of different types of seismoacoustic waves. Expressions for the time-averaged radiated powers of the Stoneley wave at the gas–solid interface and the Rayleigh wave at the solid–vacuum interface as well as analytical expressions for the Rayleigh wave displacements, which are valid for large distances from the source, are obtained for the harmonic dependence of forces on time. Excitation of a Rayleigh wave by the point sources oriented vertically, i.e., along the normal to the surface of elastic half-space, and horizontally, i.e., parallel to this surface, is analyzed in detail. Analytical expressions for the Rayleigh-wave radiated power are obtained. The dependences of these powers on the source orientation and depth are derived. It is shown that the Rayleigh-wave radiated power decreases with distance between the point of the force application and the boundary and turns to zero for a source depth of about 17.5% of the wavelength of the transverse wave in the case of a horizontally oriented subsurface source and a medium with identical Lamé parameters λ and μ. This power increases and reaches a relative maximum when the source depth becomes equal to about 42.4% of the wavelength of the transverse wave and then exponentially falls off as the source depth increases. This maximum is about 5.5% of the surface-source radiated power.     

Radiation of surface Stoneley waves by distributed force sources acting at the solid earth-atmosphere interface

Using the Fourier transform, we find an integral solution describing the excitation of seismoacoustic waves in the solid Earth and the atmosphere by time-dependent forces arbitrarily distributed over the interface between the media. The solid Earth and the atmosphere are modeled by an isotropic solid half-space and a homogeneous gaseous half-space, respectively. Depending on the types of the excited surface and bulk waves, classification of the corresponding force distributions is performed. In the case of harmonic sources, an expression for the period-averaged radiated power of the Stoneley wave is obtained. For arbitrary time dependence of the forces, we find an expression describing the the Stoneley-wave energy radiated during the entire time of the source operation. __________ Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 50, No. 7, pp. 624–637, July 2007.     

Thermodynamic compatible model of microfractured porous media and Stoneley waves

Nonstationary theory of two-velocity continuum describing the propagation of acoustic waves inmicrofractured porousmedia is based on general physical principles: the first law of thermodynamics, the conservation laws, the kinematic relationships in the metric tensor and the Galilean principle of relativity. As a physical application, the theory of the Stoneley wave in microfractured porous media is developed. The simulation results are compared with the results of physical measurement of the Stoneley wave parameters in the boreholes. It is shown that an additional fluid transport through fractures makes it possible to satisfactorily correlate the experimental and theoretical data. In general, the developed theory is a nonlinear physical model of fluid dynamics in fractured porous media.     

Interfacial waves with free-surface boundary conditions: an approach via a model equation

In a two-fluid system where the lower fluid is bounded below by a rigid bottom and the upper fluid is bounded above by a free surface, two kinds of solitary waves can propagate along the interface and the free surface: classical solitary waves characterized by a solitary pulse or generalized solitary waves with nondecaying oscillations in their tails in addition to the solitary pulse. The classical solitary waves move faster than the generalized solitary waves. The origin of the nonlocal solitary waves can be understood from a physical point of view. The dispersion relation for the above system shows that short waves can propagate at the same speed as a “slow” solitary wave. The interaction between the solitary wave and the short waves creates a nonlocal solitary wave. In this paper, the interfacial-wave problem is reduced to a system of ordinary differential equations by using a classical perturbation method, which takes into consideration the possible resonance between short waves and “slow” solitary waves. In the past, classical Korteweg–de Vries type models have been derived but cannot deal with the resonance. All solutions of the new system of model equations, including classical as well as generalized solitary waves, are constructed. The domain of validity of the model is discussed as well. It is also shown that fronts connecting two conjugate states cannot occur for “fast” waves. For “slow” waves, fronts exist but they have ripples in their tails.     

High order matched interface and boundary methods for the Helmholtz equation in media with arbitrarily curved interfaces

This work overcomes the difficulty of the previous matched interface and boundary (MIB) method in dealing with interfaces with non-constant curvatures for optical waveguide analysis. This difficulty is essentially bypassed by avoiding the use of local cylindrical coordinates in the improved MIB method. Instead, novel jump conditions are derived along global Cartesian directions for the transverse magnetic field components. Effective interface treatments are proposed to rigorously impose jump conditions across arbitrarily curved interfaces based on a simple Cartesian grid. Even though each field component satisfies the scalar Helmholtz equation, the enforcement of jump conditions couples two transverse magnetic field components, so that the resulting MIB method is a full-vectorial approach for the modal analysis of optical waveguides. The numerical performance of the proposed MIB method is investigated by considering interface problems with both constant and general curvatures. The MIB method is shown to be able to deliver a fourth order of accuracy in all cases, even when a high frequency solution is involved.     

Interaction of re-radiated Rayleigh waves with fatigue cracks

The usefulness of ultrasonic re-radiated Rayleigh waves for the detection of small surface-connected fatigue cracks in titanium alloy has been demonstrated. The variation in the frequency spectrum of re-radiated Rayleigh waves in the presence of a fatigue crack has been observed, and the sensitivity of this acoustic propagation mode in a scanning operation has been studied.     

Optimal nonlocal boundary control of the wide-angle parabolic equation for inversion of a waveguide acoustic field

This paper applies the concept of optimal boundary control for solving inverse problems in shallow water acoustics. To treat the controllability problem, a continuous analytic adjoint model is derived for the Claerbout wide-angle parabolic equation (PE) using a generalized nonlocal impedance boundary condition at the water-bottom interface. While the potential of adjoint methodology has been recently demonstrated for ocean acoustic tomography, this approach combines the advantages of exact transparent boundary conditions for the wide-angle PE with the concept of adjoint-based optimal control. In contrast to meta-heuristic approaches the inversion procedure itself is directly controlled by the waveguide physics and, in a numerical implementation based on conjugate gradient optimization, many fewer iterations are required for assessment of an environment that is supported by the underlying subbottom model. Furthermore, since regularization schemes are particularly important to enhance the performance of full-field acoustic inversion, special attention is devoted to the application of penalization methods to the adjoint optimization formalism. Regularization incorporates additional information about the desired solution in order to stabilize ill-posed inverse problems and identify useful solutions, a feature that is of particular importance for inversion of field data sampled on a vertical receiver array in the presence of measurement noise and modeling uncertainty. Results with test data show that the acoustic field and the bottom properties embedded in the control parameters can be efficiently retrieved.     

Frequency shift of Rayleigh waves excited by a source moving perpendicularly to the interface between two media

Characteristic features of the Doppler frequency shift of Rayleigh and bulk waves excited in an elastic halfspace by a source vertically moving in a contacting gaseous halfspace (atmosphere) are studied.     

Modeling the propagation of ultrasonic waves in the interface region between two bonded elements

An important task in nondestructive materials evaluation is the development of techniques to characterize the bond quality of adherent joints. Binding forces are nonlinear and cause a nonlinear modulation of transmitted and reflected ultrasonic waves. As a consequence, the higher harmonics generated by an insonified monochromatic wave give information about the adhesive bonds. The local binding forces in thin bonded interfaces can be obtained by the amplitudes of the ultrasonic waves of the insonified frequency and its higher harmonics as transmitted through the interface. Additional phase measurements may enable one to obtain the evaluation of the full hysteretic cycle of the interaction force. In order to gain a deeper understanding of the interface region and to improve the technique, numerical simulations of the ultrasonic wave propagation through specimens of two bonded elements can be used. A simple model based on the local interaction simulation approach (LISA) is described in this contribution, and a comparison between the results of the simulations and the experimental data is presented. Besides its intrinsic relevance for NDE, the problem considered in this paper may be very useful to analyze and test models for the simulation of ultrasonic wave propagation in nonclassical nonlinear mesoscopic elastic materials.     

Anisotropy measured with shear and Rayleigh waves in rolled plates

The Rayleigh wave velocity is derived for propagation on polycrystalline metal plates with weakly orthorhombic anisotropy modelling the rolling texture. Use is made of the texture weakness to solve the velocity in a perturbation scheme. For the fully orthorhombic case, there was found no simple relation between the velocity difference of Rayleigh waves propagating along the principal axes of texture and the birefringence of shear waves propagated in the thickness direction. However, for the cubic metals the result is reduced to their linear relation to the substitution of crystallite orientation distribution function expanded with spherical harmonics. This relation has been obtained by Sayers. The application of texture-independent measurement of principal stress difference is discussed.     

The kinetic boundary layer for the Fokker-Planck equation with absorbing boundary

The Fokker-Planck equation for the distribution of position and velocity of a Brownian particle is a particularly simple linear transport equation. Its normal solutions and an apparently complete set of stationary boundary layer solutions can be determined explicitly. By a numerical algorithm we select linear combinations of them that approximately fulfill the boundary condition for a completely absorbing plane wall, and that approach a linearly increasing position space density far from the wall. Various aspects of these approximate solutions are discussed. In particular we find that the extrapolated asymptotic density reaches zero at a distance x_M beyond the wall. We find x_M=1.46 in units of the velocity persistence length of the Brownian particle. This study was motivated by certain problems in the theory of diffusion-controlled reactions, and the results might be used to test approximate theories employed in that field.     

Periodic solution of the generalized Rayleigh equation

The periodic solutions of a strongly cubic nonlinear oscillator whose motion is described with the generalized Rayleigh equation are studied. Approximate analytic solving methods are introduced. A new method based on homotopy and averaging is developed to determine the limit cycle motion. The obtained analytical solutions are compared with those calculated by the elliptic harmonic balance method with generalized Fourier series and Jacobian elliptic functions. Three types of cubic nonlinearity are considered: the coefficients of the linear and cubic terms are positive, the coefficient of the linear term is positive and that of the cubic term is negative and the opposite case. Comparisons of the analytical solution and numerical solution, obtained by using the Runge-Kutta method, are illustrated with examples.     

Rayleigh and Lamb waves at liquid-solid boundaries

The equations governing Rayleigh and Lamb mode propagation are examined for free and for liquid-loaded solids. Examples are given to show under what conditions the free-solid approach yields acceptable solutions for the velocities and under what conditions the more involved liquid-loaded solid formulism must be used.