Localization, Stability, and Resolution of Topological Derivative Based Imaging Functionals in Elasticity

The focus of this work is on rigorous mathematical analysis of the topological derivative based detection algorithms for the localization of an elastic inclusion of vanishing characteristic size. A filtered quadratic misfit is considered and the performance of the topological derivative imaging functional resulting therefrom is analyzed. Our analysis reveals that the imaging functional may not attain its maximum at the location of the inclusion. Moreover, the resolution of the image is below the diffraction limit. Both phenomena are due to the coupling of pressure and shear waves propagating with different wave speeds and polarization directions. A novel imaging functional based on the weighted Helmholtz decomposition of the topological derivative is, therefore, introduced. It is thereby substantiated that the maximum of the imaging functional is attained at the location of the inclusion and the resolution is enhanced and it proves to be the diffraction limit. Finally, we investigate the stability of the proposed imaging functionals with respect to measurement and medium noises.Comment: 38 pages. A new subsection 6.4 is added where we consider the case of random Lam\'e coefficients. We thought this would corrupt the statistical stability of the imaging functional but our calculus shows that this is not the case as long as the random fluctuation is weak so that Born approximation is vali

Toplogical derivative for nonlinear magnetostatic problem

The topological derivative represents the sensitivity of a domain-dependent functional with respect to a local perturbation of the domain and is a valuable tool in topology optimization. Motivated by an application from electrical engineering, we derive the topological derivative for an optimization problem which is constrained by the quasilinear equation of two-dimensional magnetostatics. Here, the main ingredient is to establish a sufficiently fast decay of the variation of the direct state at scale 1 as $|x|\rightarrow \infty$. In order to apply the method in a bi-directional topology optimization algorithm, we derive both the sensitivity for introducing air inside ferromagnetic material and the sensitivity for introducing material inside an air region. We explicitly compute the arising polarization matrices and introduce a way to efficiently evaluate the obtained formulas. Finally, we employ the derived formulas in a level-set based topology optimization algorithm and apply it to the design optimization of an electric motor.Comment: 54 pages, 9 figure

POCHODNA TOPOLOGICZNA – TEORIA I ZASTOSOWANIA

The paper is devoted to present some mathematical aspects of the topological derivative and its applications in different fields of sciences such as shape optimization and inverse problems. First the definition of the topological derivative is given and the shape optimization problem is formulated. Next the form of the topological derivative is evaluated for a mixed boundary value problem defined in a geometrical domain. Finally, an example of an application of the topological derivative in the electric impedance tomography is presented.W pracy przedstawiono matematyczne aspekty dotyczące pochodnej topologicznej oraz jej zastosowań w różnych dziedzinach nauki, takich jak optymalizacja kształtu czy problemy odwrotne. W pierwszej części podano nieformalna˛ definicje˛ pochodnej topologicznej oraz sformułowano problem optymalizacji kształtu. Następnie wyprowadzono postać pochodnej topologicznej dla mieszanego problemu brzegowego. W ostatniej części przedstawiono przykład zastosowania pochodnej topologicznej dla problemu elektrycznej tomografii impedancyjnej

Asymptotic analysis and topological derivative for 3D quasi-linear magnetostatics

In this paper we study the asymptotic behaviour of the quasilinear $curl$-$curl$ equation of 3D magnetostatics with respect to a singular perturbation of the differential operator and prove the existence of the topological derivative using a Lagrangian approach. We follow the strategy proposed in our recent previous work (arXiv:1907.13420) where a systematic and concise way for the derivation of topological derivatives for quasi-linear elliptic problems in $H^1$ is introduced. In order to prove the asymptotics for the state equation we make use of an appropriate Helmholtz decomposition. The evaluation of the topological derivative at any spatial point requires the solution of a nonlinear transmission problem. We discuss an efficient way for the numerical evaluation of the topological derivative in the whole design domain using precomputation in an offline stage. This allows us to use the topological derivative for the design optimization of an electrical machine.Comment: 26 pages, 5 figure

Inverse homogenization using the topological derivative

Purpose: The purpose of this study is to solve the inverse homogenization problem, or so-called material design problem, using the topological derivative concept. Design/methodology/approach : The optimal topology is obtained through a relaxed formulation of the problem by replacing the characteristic function with a continuous design variable, so-called density variable. The constitutive tensor is then parametrized with the density variable through an analytical interpolation scheme that is based on the topological derivative concept. The intermediate values that may appear in the optimal topologies are removed by penalizing the perimeter functional. Findings: The optimization process benefits from the intermediate values that provide the proposed method reaching to solutions that the topological derivative had not been able to find before. In addition, the presented theory opens the path to propose a new framework of research where the topological derivative uses classical optimization algorithms. Originality/value: The proposed methodology allows us to use the topological derivative concept for solving the inverse homogenization problem and to fulfil the optimality conditions of the problem with the use of classical optimization algorithms. The authors solved several material design examples through a projected gradient algorithm to show the advantages of the proposed method.This research was partially supported by Serra Húnter Research Program (Spain), PID-UTN (Research and Development Program of the National Technological University, Argentina) and CONICET (National Council for Scientific and Technical Research, Argentina). The supports of these agencies are gratefully acknowledged.Peer ReviewedPostprint (author's final draft

Topological derivative for the inverse scattering of elastic waves

To establish an alternative analytical framework for the elastic-wave imaging of underground cavities, the focus of this study is an extension of the concept of topological derivative, rooted in elastostatics and shape optimization, to three-dimensional elastodynamics involving semi-infinite and infinite solids. The main result of the proposed boundary integral approach is a formula for topological derivative, explicit in terms of the elastodynamic fundamental solution, obtained by an asymptotic expansion of the misfit-type cost functional with respect to the creation of an infinitesimal hole in an otherwise intact (semi-infinite or infinite) elastic medium. Valid for an arbitrary shape of the infinitesimal cavity, the formula involves the solution of six canonical exterior elastostatic problems, and becomes fully explicit when the vanishing cavity is spherical. A set of numerical results is included to illustrate the potential of topological derivative as a computationally efficient tool for exposing an approximate cavity topology, location, and shape via a grid-type exploration of the host solid. For a comprehensive solution to three-dimensional inverse scattering problems involving elastic waves, the proposed approach can be used most effectively as a pre-conditioning tool for more refined, albeit computationally intensive minimization-based imaging algorithms. To the authors' knowledge, an application of topological derivative to inverse scattering problems has not been attempted before; the methodology proposed in this paper could also be extended to acoustic problems