Effective dynamic properties of 3D composite materials containing rigid penny-shaped inclusions

The propagation of time-harmonic plane elastic waves in infinite elastic composite materials consisting of linear elastic matrix and rigid penny-shaped inclusions is investigated in this paper. The inclusions are allowed to translate and rotate in the matrix. First, the three-dimensional (3D) wave scattering problem by a single inclusion is reduced to a system of boundary integral equations for the stress jumps across the inclusion surfaces. A boundary element method (BEM) is developed for solving the boundary integral equations numerically. Far-field scattering amplitudes and complex wavenumbers are computed by using the stress jumps. Then the solution of the single scattering problem is applied to estimate the effective dynamic parameters of the composite materials containing randomly distributed inclusions of dilute concentration. Numerical results for the attenuation coefficient and the effective velocity of longitudinal and transverse waves in infinite elastic composites containing parallel and randomly oriented rigid penny-shaped inclusions of equal size and equal mass are presented and discussed. The effects of the wave frequency, the inclusion mass, the inclusion density, and the inclusion orientation or the direction of the wave incidence on the attenuation coefficient and the effective wave velocities are analysed. The results presented in this paper are compared with the available analytical results in the low-frequency range.     

The scattering of steady-state SH waves in a bi-material half space with multiple cylindrical elastic inclusions

The scattering of steady-state SH waves in a bi-material half space with multiple cylindrical elastic inclusions is presented analytically. Mirror method and multi-polar coordinate systems are developed to solve the complex boundary value problem. Considering the displacement and stress continuity conditions, a series of integral equations for unknown coefficients are obtained and solved by truncation. The solution is used to calculate the dynamic stress concentration factor around the edge of the inclusion, and the analysis and numerical results are discussed. The results show that degree of dynamic stress concentration around the circular inclusion is influenced by the incident angle, the incident wave number, and the parameters of medium.     

Employment of the scattering amplitude in solving the diffraction problems for waves in a halfspace

The problem on the diffraction of an acoustic wave by a finite-size scatterer (inclusion) located in a halfspace is considered. The method of solving this problem is based on the use of the scattering amplitude of the inclusion. A formula analogous to the Green formula is presented. It allows one to determine the scattering amplitude of the inclusion for an arbitrary incident wave (determined by the directional pattern of the source of primary waves) from the scattering amplitude corresponding to plane incident waves. The algorithm is presented for solving the problem on the operation of an acoustically opaque radiator in a halfspace whose boundary is characterized by an arbitrary reflection coefficient. As an example, the problem is solved on the generation of low-frequency oscillations by a sphere with an acoustically soft boundary near an acoustically hard or soft boundary of the halfspace.     

Scattering of SH guided wave by a circular inclusion in an infinite piezoelectric material strip

By the means of the guided wave theory, the scattering of SH guided wave by a circular inclusion in an infinite piezoelectric material strip is investigated. With the aid of repeated image superposition, the analytical expression of scattering wave is conducted, which satisfies the stress free and electric insulation conditions on the upper and lower horizontal boundaries of the strip. According to the boundary condition, integral equation is set up and analytical expression of dynamic stress concentration factor and electric field intensity concentration factor are obtained. The influence of the order of guided waves, the physical parameters of the medium and position of the circular inclusion on the dynamic stress concentration factor and electric field intensity concentration factor are analyzed and compared with the existing literature in calculating example.     

A Weak Form of the Conjugate Gradient FFT Method for Two-Dimensional Elastodynamics

The problem of two-dimensional scattering of elastic waves by an elastic inclusion can be formulated in terms of a domain integral equation, in which the grad-div operator acts on a vector potential. The vector potential is the spatial convolution of a Green's function with the product of the density and the displacement over the domain of interest. A weak form of the integral equation for the unknown displacement is obtained by testing it with rooftop functions. This method shows excellent numerical performance.     

Diffraction of guided waves by the end face of a dielectric slab waveguide terminated with a finite conductive strip

The diffraction of guided waves by the end face of a dielectric slab waveguide short circuited with a finite conductive strip is analyzed. An integral equation technique is employed to formulate the corresponding boundary problem. The unknown term in this integral equations is the electric field E(x) on the terminal plane of the waveguide. The homogeneous term is determined from the incident guided wave. A method of moments technique is employed to compute approximately the electric field E(x) by using Laguerre functions as describing and testing functions. The reflection coefficients of the guided waves are computed by using the approximate expression of the E(x) field. Numerical results are given for several guide and conductor plate dimensions.     

On formulation of a transition matrix for electroporoelastic medium and application to analysis of scattered electroseismic wave

On the basis of Pride's theory (1994) which couples Biot's theory for poroelastic medium (1956) and Maxwell equations via flux/force transport equations, we extend Yeh et al. (2004) approach for poroelastic medium to develop a transition matrix for electroporoelastic medium. The transition matrix, which relates the coefficients of scattered waves to those of incident waves, is then derived through the application of Betti's third identity and the associated orthogonality conditions for the electroporoelastic medium. To illustrate the application, a simple case of the scattering problem of a spherical electroporoelastic inclusion, embedded within the surrounding electroporoelastic medium subjected to an incident plane compressional wave is considered.     

Simulation of the vector wave field of a low-frequency sound source in the ocean: Some important features

An approach to the simulation of low frequency vector wave fields in stratified media (mainly in the ocean) is considered. The approach is characterized by an improved stability with respect to dividing the medium into many layers of arbitrary thickness. The model for the sound field of a point source is based on an integral representation of two-dimensional, cylindrically symmetric vector wave fields in inhomogeneous media, so that the contributions of all types of waves are included automatically. The model medium is subdivided into N horizontally homogeneous layers for which 4(N?1) equations are formulated to satisfy the boundary conditions between adjacent layers. The method of the generalized Schmidt matrix is used to obtain the coefficients of the equations; these coefficients are substituted into the expressions (of the Fourier-Bessel integral type) for the local parameters of the field. The latter are calculated according to the numerical procedure, and the results are used to model the distributions of the acoustic pressure and the horizontal and vertical components of the particle velocity in liquid and elastic media. The instability of the calculation procedure may result in a disagreement between the model and the exact solution. However, the disagreement is shown to occur mainly in models containing excessively thick layers. A way for improving the stability of the numerical model is suggested. The simulation results are compared with the exact analytical solution for the simplest example and with the results obtained according to the commonly used generalized matrix procedure (the benchmark problem). The examples of the practical application of the model for investigating more complex seismoacoustic wave fields in the ocean are presented.     

Multi-step processes in high-energy elastic and inelastic nuclear scattering

Multi-step processes in elastic and inelastic nuclear scattering at intermediate and high energies are investigated using a formulation whereby a finite number of channels are explicitly treated while all the other channels are approximately accounted for through a “second-order potential matrix”. Within the framework of the eikonal approximation the problem reduces to a finite system of first-order coupled integro-differential equations with non-local potentials which depend on the two-body density matrix of the target nucleus. The relationship of the above formulation to the DWIA, the close-coupling method, and the Glauber multiple scattering model is examined. This approach is applied to the elastic and inelastic (2⁺, 4.43 MeV) scattering of 1 GeV nucleons by ¹²C. The corrections to the DWIA are sizeable, and the inelastic scattering appears to be very sensitive to the multi-step contributions and the nuclear structure.     

Rigorous and approximate methods for modeling wave scattering from a locally perturbed perfectly conducting surface

Rigorous and approximate methods are considered for solving the problem of harmonic plane wave scattering from a plane surface arbitrarily perturbed along one dimension on a finite interval. This problem is treated using the Fredholm integral equations of the second kind and the Kirchhoff and Rayleigh approximations. The estimates of the computational efficiency of the integral equation method and the Rayleigh approximation are compared by calculating fields scattered from random rough surfaces in the resonance region (i.e., when the roughness height is comparable to or smaller than the incident wavelength) for an arbitrary incidence of a plane wave. Scattering patterns calculated using the integral equations and the Kirchhoff approximation are discussed in the case of large-scale random rough surface scattering. Particular attention is paid to scattering at near-grazing incidence.     

Dispersion and Attenuation of Surface Acoustic Waves of Various Polarizations on a Stress-Free Randomly Rough Surface of Solid

An approach to obtaining the dispersion equation of surface acoustic waves (SAWs) on a stress-free, randomly rough surface of an anisotropic elastic medium is suggested. The problem is solved in the approximation of a weakly rough surface using Green′s function technique. The dispersion and attenuation of sagittally and shear horizontally (SH) polarized SAWs are investigated both analytically and numerically for a three-dimensionally (3D) and a two-dimensionally (2D) rough surface of an isotropic medium. The results for 2D roughness are shown to be contained in the more general expressions for the 3D case, and the connection between the results for the 3D and the 2D cases is pointed out. Dispersion relations are derived for SAWs of both polarizations propagating in an arbitrary direction along a 2D rough surface. The SAW attenuation mechanisms are investigated at various incidence angles. It is concluded that all three mechanisms (viz. scattering into bulk transverse, longitudinal, and Rayleigh surface acoustic waves) are involved for the Rayleigh and SH polarized SAWs at certain incidence angles, whereas at the other angles only some of the mechanisms are. The criterion for the existence of SH polarized SAWs on a rough surface is considered. A possible increase of the SAW phase velocity on a rough surface compared with that for a flat boundary is discussed. In the limit λ a (where a is the roughness correlation length) simple explicit expressions for the phase velocities of Rayleigh and SH polarized SAWs are derived. A comparison of the results obtained herein with those of other workers is presented.     

Simulation of electromagnetic scattering by random discrete particles using a hybrid FE-BI-CBFM technique

A hybrid finite element–boundary integral–characteristic basis function method (FE-BI-CBFM) is proposed for an efficient simulation of electromagnetic scattering by random discrete particles. Specifically, the finite element method (FEM) is used to obtain the solution of the vector wave equation inside each particle and the boundary integral equation (BIE) using Green's functions is applied on the surfaces of all the particles as a global boundary condition. The coupling system of equations is solved by employing the characteristic basis function method (CBFM) based on the use of macro-basis functions constructed according to the Foldy–Lax multiple scattering equations. Due to the flexibility of FEM, the proposed hybrid technique can easily deal with the problems of multiple scattering by randomly distributed inhomogeneous particles that are often beyond the scope of traditional numerical methods. Some numerical examples are presented to demonstrate the validity and capability of the proposed method.     

Model problem of wave refraction by a periodically nonuniform boundary of the solar plasma

The problem of the diffraction of Alfvén and magnetosonic waves by a boundary between two media in the form of a plane interface modulated by a running sinusoidal wave is considered in the linear approximation on the basis of the mathematical apparatus of integral equations of solar magnetohydrodynamics. The results obtained are analyzed and discussed.     

Scattering by cavities of arbitrary shape in an infinite plate and associated vibration problems

This paper presents a solution for the displacement of a uniform elastic thin plate with an arbitrary cavity, modelled using the biharmonic plate equation. The problem is formulated as a system of boundary integral equations by factorizing the biharmonic equation, with the unknown boundary values expanded in terms of a Fourier series. At the edge of the cavity we consider free-edge, simply-supported and clamped boundary conditions. Methods to suppress ill-conditioning which occurs at certain frequencies are discussed, and the combined boundary integral equation method is implemented to control this problem. A connection is made between the problem of an infinite plate with an arbitrary cavity and the vibration problem of an arbitrarily shaped finite plate, using the jump discontinuity present in single-layer distributions at the boundary. The first few frequencies and modes of displacement are computed for circular and elliptic cavities, which provide a check on our numerics, and results for the displacement of an infinite plate are given for four specific cavity geometries and various boundary conditions.     

Scattering of monochromatic longitudinal waves on a planar crack of arbitrary shape in a fluid-saturated poroelastic medium

Scattering of monochromatic longitudinal waves on a planar crack of arbitrary shape in a saturated poroelastic medium is considered. The medium is described by Biot’s constitutive equations, the crack sides are fluid permeable. The problem is reduced to a two-dimensional integral equation for the crack opening vector. Gaussian approximating functions are used for discretization of this equation. For such functions, the elements of the matrix of discretized problem are combinations of four standard one-dimensional integrals that can be tabulated. As a result, numerical integration is not needed. For regular grids of approximating nodes, this matrix has Toeplitz’s structure, and matrix-vector products can be calculated by the fast Fourier transform technique. The latter accelerates substantially the process of iterative solution of the discretized problem. Calculation of crack opening vectors, differential, and total cross-sections of circular and elliptic cracks are performed for longitudinal incident waves orthogonal to the crack surfaces. Dependencies of these characteristics on the medium permeability and wavefrequency are studied. Comparison of a crack in the poroelastic medium and in a dry elastic medium with the same porosity and skeleton elastic properties is presented.     

Modified null field method for elastic wave scattering from a partially debonded fiber-SH-waves

Debonded fibers influence the macro-mechanical behavior of fiber-reinforced composite materials. Debonded fibers contribute to the initiation and growth of cracks at the fiber/matrix interface. To examine such problems, the scattering of elastic SH-waves (problem of anti-plane strain) from debonded fibers is studied with a numerical method which can handle the mixed boundary conditions on the partially bonded fiber. A modification of the null field of T-matrix method is developed for this purpose. The modification is achieved by the introduction of a mathematical surface. The simultaneous solution of the integral representations of the field scattered by the mathematical surface and the actual fiber surface give rise to sufficient equations that permit solution of the debonded fiber problem. The scattering cross-section and the far field amplitude are calculated as a function of frequency, fiber properties and debonding area. This will be of interest in structural applications where such cross-sections can be used to compute the dynamical effective properties of damaged composites.     

Low-energy elastic scattemng of pions by the deuteron

The problem of the elastic scattering of pions by a deuteron is considered using the separable representation of the two-body t-matrix. The Faddeev equations are reduced to a set of one-dimensional integral equations by separating the angular variables. The dependence of the π-d scattering length on the form of two-body interaction and on the values of the π-N scattering lengths is studied in the case of a one-term nonlocal potential with separable variables. The π-d scattering length proves to be practically independent of the two-body interaction form, and is essentially dependent on the values of the π-N scattering lengths.     

Plane wave interpolation in direct collocation boundary element method for radiation and wave scattering: numerical aspects and applications

The classical boundary element formulation for the Helmholtz equation is rehearsed, and its limitations with respect to the number of variables needed to model a wavelength are explained. A new type of interpolation for the potential is then described in which the usual boundary element shape functions are modified by the inclusion of a set of plane waves, propagating in a range of directions. This is termed the plane wave basis boundary element method. The modifications needed to the classical procedures, in terms of integration of the element matrices, and location of collocation points are described. The well-known Singular Value Decomposition solution technique, which is adopted here for the solution of the system matrix equation in its complex form, is briefly outlined. The conditioning of the system matrix is analysed for a simple radiation problem. The corresponding diffraction problem is also analysed and results are compared with analytical and classical boundary element solutions. The CHIEF method is adopted to enhance the quality of the solution, particularly in the vicinity of irregular frequencies. The plane wave basis boundary element method is then applied to two problems: scattering of plane waves by an elliptical cylinder and the multiple circular cylinder plane wave scattering problem. In both cases results are compared with analytical solutions. The results clearly demonstrate that the new method is considerably more efficient than the classical approach. For a given number of degrees of freedom, the frequency for which accurate results can be obtained, using the new technique, can be up to three or four times higher than that of the classical method. This makes the method a powerful new addition to our tools for tackling high-frequency radiation and scattering problems.     

Dynamic polarization of 6,7Li nuclei in the elastic scattering of α particles

The elastic scattering of α particles on weakly bound ^6,7Li cluster nuclei is considered with allowance for their dynamic polarization within the three-particle model. The considered states of the αd and αt continuums are projected on a finite basis of stationary wave packets, which allows the total three-body problem to be reduced to a matrix problem. The results from calculating α + ^6,7Li-elastic scattering differential cross sections are considered as illustration of the three-particle approach and is compared to the results of other authors.     

Influence of a thin surface layer on light scattering by spherical particles

The Mie problem with modified boundary conditions that take into account the influence of a thin surface layer on the scattering of an electromagnetic wave by a spherical particle is considered. Analytical equations are derived for the partial amplitudes of scattered waves and forced oscillations. These equations are applicable in the case of anisotropy and gyrotropy of an optical response from the surface layer.