skip to main content


Title: Analysis of topological derivative as a tool for qualitative identification
The concept of topological derivative has proved effective as a qualitative inversion tool for a wave-based identification of finite-sized objects. Although for the most part, this approach remains based on a heuristic interpretation of the topological derivative, a first attempt toward its mathematical justification was done in Bellis et al (2013 Inverse Problems 29 075012) for the case of isotropic media with far field data and inhomogeneous refraction index. Our paper extends the analysis there to the case of anisotropic scatterers and background with near field data. Topological derivative-based imaging functional is analyzed using a suitable factorization of the near fields, which became achievable thanks to a new volume integral formulation recently obtained in Bonnet (2017 J. Integral Equ. Appl. 29 271–95). Our results include justification of sign heuristics for the topological derivative in the isotropic case with jump in the main operator and for some cases of anisotropic media, as well as verifying its decaying property in the isotropic case with near field spherical measurements configuration situated far enough from the probing region.  more » « less
Award ID(s):
1813492
NSF-PAR ID:
10158766
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Inverse problems
Volume:
35
Issue:
10
ISSN:
0266-5611
Page Range / eLocation ID:
104007
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In this paper, we consider the inverse scattering problem for recovering either an isotropic or anisotropic scatterer from the measured scattered field initiated by a point source. We propose two new imaging functionals for solving the inverse problem. The first one employs a 'far-field' transform to the data which we then use to derive and provide an explicit decay rate for the imaging functional. In order to analyze the behavior of this imaging functional we use the factorization of the near field operator as well as the Funk-Hecke integral identity. For the second imaging functional the Cauchy data is used to define the functional and its behavior is analyzed using the Green's identities. Numerical experiments are given in two dimensions for both isotropic and anisotropic scatterers.

     
    more » « less
  2. ABSTRACT The circular-crack model has been widely used in seismology to infer earthquake stress drop. A common assumption is that the background medium is isotropic, although many earthquakes occur in geologically anisotropic settings. In this article, we study the effect of anisotropy on stress drop for a circular crack model and present explicit formalism in both static and kinematic cases. In the static case, we obtain the relationship between stress drop and slip for a circular crack model in an arbitrarily anisotropic medium. Special attention is given to the transversely isotropic (TI) medium. The static formalism is useful in understanding stress drop, but not all quantities are observables. Therefore, we resort to the kinematic case, from which we can infer stress drop using recorded far-field body waves. In the kinematic case, we assume that the crack ruptures circularly and reaches the final displacement determined by the static solutions. The far-field waveforms show that the corner frequency will change with different anisotropic parameters. Finally, we calculate the stress drops for cracks in isotropic and anisotropic media using the far-field waveforms. We find that in an isotropic medium, only shear stress acting on the crack surface contributes to shear slip. However, in a TI medium, if the anisotropy symmetry axis is not perpendicular or parallel to the crack surface, a normal stress (normal to the crack surface) can produce a shear slip. In calculating stress drop for an earthquake in an anisotropic medium using far-field body waves, a large error may be introduced if we ignore the possible anisotropy in the inversion. For a TI medium with about 18% anisotropy, the misfit of inferred stress drop could be up to 41%. Considering the anisotropic information, we can further improve the accuracy of stress-drop inversion. 
    more » « less
  3. Frequency-domain wave scattering problems that arise in acoustics and electromagnetism can be often described by the Helmholtz equation. This work presents a boundary integral equation method for the Helmholtz equation in 3-D multilayered media with many doubly periodic smooth layers. Compared with conventional quasi-periodic Green’s function method, the new method is robust at all scattering parameters. A periodizing scheme is used to decompose the solution into near-and far-field contributions. The near-field contribution uses the free-space Green’s function in an integral equation on the interface in the unit cell and its immediate eight neighbors; the far-field contribution uses proxy point sources that enclose the unit cell. A specialized high-order quadrature is developed to discretize the underlying surface integral operators to keep the number of unknowns per layer small. We achieve overall linear computational complexity in the number of layers by reducing the linear system into block tridiagonal form and then solving the system directly via block LU decomposition. The new solver is capable of handling a 100-interface structure with 961.3k unknowns to 10−5 accuracy in less than 2 hours on a desktop workstation. 
    more » « less
  4. null (Ed.)
    Abstract This paper is concerned with the inverse scattering problem for the three-dimensional Maxwell equations in bi-anisotropic periodic structures.The inverse scattering problem aims to determine the shape of bi-anisotropic periodic scatterers from electromagnetic near-field data at a fixed frequency.The factorization method is studied as an analytical and numerical tool for solving the inverse problem.We provide a rigorous justification of the factorization method which results in the unique determination and a fast imaging algorithm for the periodic scatterer.Numerical examples for imaging three-dimensional periodic structures are presented to examine the efficiency of the method. 
    more » « less
  5. Anisotropy in additive manufacturing (AM), particularly in the material extrusion process, plays a crucial role in determining the actual structural performance, including the stiffness and strength of the printed parts. Unless accounted for, anisotropy can compromise the objective performance of topology-optimized structures and allow premature failures for stress-sensitive design domains. This study harnesses process-induced anisotropy in material extrusion-based 3D printing to design and fabricate stiff, strong, and lightweight structures using a two-step framework. First, an AM-oriented anisotropic strength-based topology optimization formulation optimizes the structural geometry and infill orientations, while assuming both anisotropic (i.e., transversely isotropic) and isotropic infill types as candidate material phases. The dissimilar stiffness and strength interpolation schemes in the formulation allow for the optimized allocation of anisotropic and isotropic material phases in the design domain while satisfying their respective Tsai–Wu and von Mises stress constraints. Second, a suitable fabrication methodology realizes anisotropic and isotropic material phases with appropriate infill density, controlled print path (i.e., infill directions), and strong interfaces of dissimilar material phases. Experimental investigations show up to 37% improved stiffness and 100% improved strength per mass for the optimized and fabricated structures. The anisotropic strength-based optimization improves load-carrying capacity by simultaneous infill alignment along the stress paths and topological adaptation in response to high stress concentration. The adopted interface fabrication methodology strengthens comparatively weaker anisotropic joints with minimal additional material usage and multi-axial infill patterns. Furthermore, numerically predicted failure locations agree with experimental observations. The demonstrated framework is general and can potentially be adopted for other additive manufacturing processes that exhibit anisotropy, such as fiber composites. 
    more » « less