An official website of the United States government
Consider the scattering of a time-harmonic acoustic plane wave by a bounded elastic obstacle which is immersed in a homogeneous acoustic medium. This paper is concerned with an inverse acoustic-elastic interaction problem, which is to determine the location and shape of the elastic obstacle by using either the phased or phaseless far-field data. By introducing the Helmholtz decomposition, the model problem is reduced to a coupled boundary value problem of the Helmholtz equations. The jump relations are studied for the second derivatives of the single-layer potential in order to deduce the corresponding boundary integral equations. The well-posedness is discussed for the solution of the coupled boundary integral equations. An efficient and high order Nyström-type discretization method is proposed for the integral system. A numerical method of nonlinear integral equations is developed for the inverse problem. For the case of phaseless data, we show that the modulus of the far-field pattern is invariant under a translation of the obstacle. To break the translation invariance, an elastic reference ball technique is introduced. We prove that the inverse problem with phaseless far-field pattern has a unique solution under certain conditions. In addition, a numerical method of the reference ball technique based nonlinear integral equations is proposed for the phaseless inverse problem. Numerical experiments are presented to demonstrate the effectiveness and robustness of the proposed methods.
more »
« less
Reconstruction techniques for quantum trees
The inverse problem of recovery of a potential on a quantum tree graph from the Weyl matrix given at a number of points is considered. A method for its numerical solution is proposed. The overall approach is based on the leaf peeling method combined with Neumann series of Bessel functions (NSBF) representations for solutions of Sturm–Liouville equations. In each step, the solution of the arising inverse problems reduces to dealing with the NSBF coefficients. The leaf peeling method allows one to localize the general inverse problem to local problems on sheaves, while the approach based on the NSBF representations leads to splitting the local problems into two‐spectrum inverse problems on separate edges and reduces them to systems of linear algebraic equations for the NSBF coefficients. Moreover, the potential on each edge is recovered from the very first NSBF coefficient. The proposed method leads to an efficient numerical algorithm that is illustrated by numerical tests.
more »
« less
- Award ID(s):
- 2308377
- PAR ID:
- 10493130
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Mathematical Methods in the Applied Sciences
- Volume:
- 47
- Issue:
- 9
- ISSN:
- 0170-4214
- Format(s):
- Medium: X Size: p. 7182-7197
- Size(s):
- p. 7182-7197
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract In this paper, we explore the inverse dynamic problem for the Dirac system on finite metric graphs, including trees and graphs with a cycle. Our primary objective is to reconstruct the graph’s topology (connectivity), determine the lengths of its edges, and identify the matrix potential function on each edge. By using only the dynamic matrix response operator as our inverse data, we adapt the leaf peeling method to recover the unknown data on a tree graph. We then introduce a new approach to reconstruct the unknown data on a graph with a cycle. Additionally, we present a novel dynamic algorithm to address the forward problem for the Dirac system on finite metric graphs.more » « less
-
We study an inverse problem for the time-dependent Maxwell system in an inhomogeneous and anisotropic medium. The objective is to recover the initial electric field $$\mathbf{E}_0$$ in a bounded domain $$\Omega \subset \mathbb{R}^3$$, using boundary measurements of the electric field and its normal derivative over a finite time interval. Informed by practical constraints, we adopt an under-determined formulation of Maxwell's equations that avoids the need for initial magnetic field data and charge density information. To address this inverse problem, we develop a time-dimension reduction approach by projecting the electric field onto a finite-dimensional Legendre polynomial-exponential basis in time. This reformulates the original space-time problem into a sequence of spatial systems for the projection coefficients. The reconstruction is carried out using the quasi-reversibility method within a minimum-norm framework, which accommodates the inherent non-uniqueness of the under-determined setting. We prove a convergence theorem that ensures the quasi-reversibility solution approximates the true solution as the noise and regularization parameters vanish. Numerical experiments in a fully three-dimensional setting validate the method's performance. The reconstructed initial electric field remains accurate even with $$10\%$$ noise in the data, demonstrating the robustness and applicability of the proposed approach to realistic inverse electromagnetic problems.more » « less
-
Abstract We study the problem of finding the resistors in a resistor network from measurements of the power dissipated by the resistors under different loads. We give sufficient conditions for local uniqueness, i.e. conditions that guarantee that the linearization of this non-linear inverse problem admits a unique solution. Our method is inspired by a method to study local uniqueness of inverse problems with internal functionals in the continuum, where the inverse problem is reformulated as a redundant system of differential equations. We use our method to derive local uniqueness conditions for other discrete inverse problems with internal functionals including a discrete analogue of the inverse Schrödinger problem and problems where the resistors are replaced by impedances and dissipated power at the zero and a positive frequency are available. Moreover, we show that the dissipated power measurements can be obtained from measurements of thermal noise induced currents.more » « less
-
Purpose The purpose of this paper is as follows: to significantly reduce the computation time (by a factor of 1,000 and more) compared to known numerical techniques for real-world problems with complex interfaces; and to simplify the solution by using trivial unfitted Cartesian meshes (no need in complicated mesh generators for complex geometry). Design/methodology/approach This study extends the recently developed optimal local truncation error method (OLTEM) for the Poisson equation with constant coefficients to a much more general case of discontinuous coefficients that can be applied to domains with different material properties (e.g. different inclusions, multi-material structural components, etc.). This study develops OLTEM using compact 9-point and 25-point stencils that are similar to those for linear and quadratic finite elements. In contrast to finite elements and other known numerical techniques for interface problems with conformed and unfitted meshes, OLTEM with 9-point and 25-point stencils and unfitted Cartesian meshes provides the 3-rd and 11-th order of accuracy for irregular interfaces, respectively; i.e. a huge increase in accuracy by eight orders for the new 'quadratic' elements compared to known techniques at similar computational costs. There are no unknowns on interfaces between different materials; the structure of the global discrete system is the same for homogeneous and heterogeneous materials (the difference in the values of the stencil coefficients). The calculation of the unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations and yields the optimal order of accuracy of OLTEM at a given stencil width. The numerical results with irregular interfaces show that at the same number of degrees of freedom, OLTEM with the 9-points stencils is even more accurate than the 4-th order finite elements; OLTEM with the 25-points stencils is much more accurate than the 7-th order finite elements with much wider stencils and conformed meshes. Findings The significant increase in accuracy for OLTEM by one order for 'linear' elements and by 8 orders for 'quadratic' elements compared to that for known techniques. This will lead to a huge reduction in the computation time for the problems with complex irregular interfaces. The use of trivial unfitted Cartesian meshes significantly simplifies the solution and reduces the time for the data preparation (no need in complicated mesh generators for complex geometry). Originality/value It has been never seen in the literature such a huge increase in accuracy for the proposed technique compared to existing methods. Due to a high accuracy, the proposed technique will allow the direct solution of multiscale problems without the scale separation.more » « less
