No abstract

The aim of this paper is to provide a fast and efficient procedure for (real-time) target identification in imaging based on matching on a dictionary of precomputed generalized polarization tensors (GPTs). The approach is based on some important properties of the GPTs and new invariants. A new shape representation is given and numerically tested in the presence of measurement noise. The stability and resolution of the proposed identification algorithm is numerically quantified. Mathematics Subject Classification (MSC2000): 35R30, 35B30

A rigorous homogenization theory is derived to describe the effective admittivity of cell suspensions. A new formula is reported for dilute cases that gives the frequency-dependent effective admittivity with respect to the membrane polarization. Different microstructures are shown to be distinguishable via spectroscopic measurements of the overall admittivity using the spectral properties of the membrane polarization. The Debye relaxation times associated with the membrane polarization tensor are shown to be able to give the microscopic structure of the medium. A natural measure of the admittivity anisotropy is introduced and its dependence on the frequency of applied current is derived. A Maxwell-Wagner-Fricke formula is given for concentric circular cells, and the results can be extended to the random cases. A randomly deformed periodic medium is also considered and a new formula is derived for the overall admittivity of a dilute suspension of randomly deformed cells.Mathematics Subject Classification (MSC2000): 35R30, 35B30.

This paper aims to mathematically advance the field of quantitative thermo-acoustic imaging. Given several electromagnetic data sets, we establish for the first time an analytical formula for reconstructing the absorption coefficient from thermal energy measurements. Since the formula involves derivatives of the given data up to the third order, it is unstable in the sense that small measurement noises may cause large errors. However, in the presence of measurement noise, the obtained formula, together with a noise regularization technique, provides a good initial guess for the true absorption coefficient. We finally correct the errors by deriving a reconstruction formula based on the least square solution of an optimal control problem and prove that this optimization step reduces the errors occurring and enhances the resolution.

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.

No abstract

We consider the optimization approach to the acousto-electric imaging problem. Assuming that the electric conductivity distribution is a small perturbation of a constant, we investigate the least-squares estimate analytically using (multiple) Fourier series, and confirm the widely observed fact that acousto-electric imaging has high resolution and is statistically stable. We also analyze the case of partial data and the case of limited-view data, in which some singularities of the conductivity can still be imaged.

We analyze the random fluctuations of several multi-scale algorithms such as the multi-scale finite element method (MsFEM) and the finite element heterogeneous multiscale method (HMM), that have been developed to solve partial differential equations with highly heterogeneous coefficients. Such multi-scale algorithms are often shown to correctly capture the homogenization limit when the highly oscillatory random medium is stationary and ergodic. This paper is concerned with the random fluctuations of the solution about the deterministic homogenization limit. We consider the simplified setting of the one dimensional elliptic equation, where the theory of random fluctuations is well understood. We develop a fluctuation theory for the multi-scale algorithms in the presence of random environments with short-range and long-range correlations.For a given mesh size h, we show that the fluctuations converge in distribution in the space of continuous paths to Gaussian processes as the correlation length ε → 0. We next derive the limit of such Gaussian processes as h → 0 and compare this limit with the distribution of the random fluctuations of the continuous model. When such limits agree, we conclude that the multi-scale algorithm captures the random fluctuations accurately and passes the corrector test. This property serves as an interesting Benchmark to assess the behavior of the multi-scale algorithm in practical situations where the assumptions necessary for the theory of homogenization are not met.What we find is that the computationally more expensive methods MsFEM, and HMM with a choice of parameter δ = h, correctly capture the random fluctuations both for short-range and long-range oscillations in the medium. The less expensive method HMM with δ < h correctly captures the fluctuations for long-range oscillations and strongly amplifies their size in media with short-range oscillations. We present a modified scheme with an intermediate computational cost that captures the random fluctuations in all cases.