Topological Sensitivity Based Far-Field Detection of Elastic Inclusions

The aim of this article is to present and rigorously analyze topological sensitivity based algorithms for detection of diametrically small inclusions in an isotropic homogeneous elastic formation using single and multiple measurements of the far-field scattering amplitudes. A $L^2-$cost functional is considered and a location indicator is constructed from its topological derivative. The performance of the indicator is analyzed in terms of the topological sensitivity for location detection and stability with respect to measurement and medium noises. It is established that the location indicator does not guarantee inclusion detection and achieves only a low resolution when there is mode-conversion in an elastic formation. Accordingly, a weighted location indicator is designed to tackle the mode-conversion phenomenon. It is substantiated that the weighted function renders the location of an inclusion stably with resolution as per Rayleigh criterion.Comment: 31 pages, 1 figur

Elastodynamics of radially inhomogeneous spherically anisotropic elastic materials in the Stroh formalism

A method is presented for solving elastodynamic problems in radially inhomogeneous elastic materials with spherical anisotropy, i.e.\ materials such that $c_{ijkl}= c_{ijkl}(r)$ in a spherical coordinate system ${r,\theta,\phi}$. The time harmonic displacement field $\mathbf{u}(r,\theta ,\phi)$ is expanded in a separation of variables form with dependence on $\theta,\phi$ described by vector spherical harmonics with $r$-dependent amplitudes. It is proved that such separation of variables solution is generally possible only if the spherical anisotropy is restricted to transverse isotropy with the principal axis in the radial direction, in which case the amplitudes are determined by a first-order ordinary differential system. Restricted forms of the displacement field, such as $\mathbf{u}(r,\theta)$, admit this type of separation of variables solutions for certain lower material symmetries. These results extend the Stroh formalism of elastodynamics in rectangular and cylindrical systems to spherical coordinates.Comment: 15 page

A vector fast multipole algorithm for low frequency problems

Instead of the traditional factorization of the scalar Green's function by using scalar addition theorem in the lowfrequency fast multipole algorithm (LF-FMA), we adopt the vector addition theorem (VAT) for the factorization of the dyadic Green's function to realize memory savings for large scale problems. We validate this factorization and use it to develop a low-frequency vector fast multipole algorithm (LF-VFMA) for low-frequency problems. © 2010 IEEE.published_or_final_versionThe URSI International Symposium on Electromagnetic Theory (EMTS 2010), Berlin, Germany, 16-19 August 2010. In Proceedings of the URSI International Symposium on Electromagnetic Theory, 2010, p. 620-62

Generalized linear sampling method for elastic-wave sensing of heterogeneous fractures

A theoretical foundation is developed for active seismic reconstruction of fractures endowed with spatially-varying interfacial condition (e.g.~partially-closed fractures, hydraulic fractures). The proposed indicator functional carries a superior localization property with no significant sensitivity to the fracture's contact condition, measurement errors, and illumination frequency. This is accomplished through the paradigm of the $F_\sharp$-factorization technique and the recently developed Generalized Linear Sampling Method (GLSM) applied to elastodynamics. The direct scattering problem is formulated in the frequency domain where the fracture surface is illuminated by a set of incident plane waves, while monitoring the induced scattered field in the form of (elastic) far-field patterns. The analysis of the well-posedness of the forward problem leads to an admissibility condition on the fracture's (linearized) contact parameters. This in turn contributes toward establishing the applicability of the $F_\sharp$-factorization method, and consequently aids the formulation of a convex GLSM cost functional whose minimizer can be computed without iterations. Such minimizer is then used to construct a robust fracture indicator function, whose performance is illustrated through a set of numerical experiments. For completeness, the results of the GLSM reconstruction are compared to those obtained by the classical linear sampling method (LSM)

Elastic scattering coefficients and enhancement of nearly elastic cloaking

The concept of scattering coefficients has played a pivotal role in a broad range of inverse scattering and imaging problems in acoustic and electromagnetic media. In view of their promising applications, we introduce the notion of scattering coefficients of an elastic inclusion in this article. First, we define elastic scattering coefficients and substantiate that they naturally appear in the expansions of elastic scattered field and far field scattering amplitudes corresponding to a plane wave incidence. Then an algorithm is developed and analyzed for extracting the elastic scattering coefficients from multi-static response measurements of the scattered field. Moreover, the estimate of the maximal resolving order is provided in terms of the signal-to-noise ratio. The decay rate and symmetry of the elastic scattering coefficients are also discussed. Finally, we design scattering-coefficients-vanishing structures and elucidate their utility for enhancement of nearly elastic cloaking

Elastic scattering coefficients and enhancement of nearly elastic cloaking

The concept of scattering coefficients has played a pivotal role in a broad range of inverse scattering and imaging problems in acoustic and electromagnetic media. In view of their promising applications, we introduce the notion of scattering coefficients of an elastic inclusion in this article. First, we define elastic scattering coefficients and substantiate that they naturally appear in the expansions of elastic scattered field and far field scattering amplitudes corresponding to a plane wave incidence. Then an algorithm is developed and analyzed for extracting the elastic scattering coefficients from multi-static response measurements of the scattered field. Moreover, the estimate of the maximal resolving order is provided in terms of the signal-to-noise ratio. The decay rate and symmetry of the elastic scattering coefficients are also discussed. Finally, we design scattering-coefficients-vanishing structures and elucidate their utility for enhancement of nearly elastic cloaking