Regularized Newton iteration method for a penetrable cavity with internal measurements in inverse scattering problem

In this paper, we consider the inverse scattering problem of determining the shape of a cavity with a penetrable inhomogeneous medium of compact support from one source and a knowledge of measurements placed on a curve inside the cavity. First, the boundary value problem of the partial differential equations can be transformed into an equivalent system of nonlinear and ill-posed integral equations for the unknown boundary. Then, we apply the regularized Newton iterative method to reconstruct the boundary and prove the injectivity for the linearized system. Finally, we present some numerical examples to show the feasibility of our method.     

Characterization of the scattering data for the Sturm–Liouville operator

This work studies the scattering problem on the real axis for the Sturm–Liouville equation with discontinuous leading coefficient and the real‐valued steplike potential q(x) that has different constant asymptotes as x → ± ∞ . We investigate the properties of the scattering data, obtain the main integral equations of the inverse scattering problem, and also give necessary and sufficient conditions characterizing the scattering data. Copyright © 2013 John Wiley & Sons, Ltd.     

The point source method for inverse scattering in the time domain

Many recent inverse scattering techniques have been designed for single frequency scattered fields in the frequency domain. In practice, however, the data is collected in the time domain. Frequency domain inverse scattering algorithms obviously apply to time‐harmonic scattering, or nearly time‐harmonic scattering, through application of the Fourier transform. Fourier transform techniques can also be applied to non‐time‐harmonic scattering from pulses. Our goal here is twofold: first, to establish conditions on the time‐dependent waves that provide a correspondence between time domain and frequency domain inverse scattering via Fourier transforms without recourse to the conventional limiting amplitude principle; secondly, we apply the analysis in the first part of this work toward the extension of a particular scattering technique, namely the point source method, to scattering from the requisite pulses. Numerical examples illustrate the method and suggest that reconstructions from admissible pulses deliver superior reconstructions compared to straight averaging of multi‐frequency data. Copyright © 2006 John Wiley & Sons, Ltd.     

Application of the Factorization Method to Recover Cuts with Oblique Derivative Boundary Condition

Direct and inverse problems for the scattering of cracks with mixed oblique derivative boundary conditions from the incident plane wave are considered, which describe the scattering phenomenons such as the scattering of tidal waves by spits or reefs. The solvability of the direct scattering problem is proven by using the boundary integral equation method. In order to show the equivalent boundary integral system is Fredholm of index zero, some relationships concerning the tangential potential operator is used. Due to the mixed oblique derivative boundary conditions, we cannot employ the factorization method in a usual manner to reconstruct the cracks. An alternative technique is used in the theoretical analysis such that the far field operator can be factorized in an appropriate form and fulfills the range identity theorem. Finally, we present some numerical examples to demonstrate the feasibility and effectiveness of the factorization method.     

Continuity of the scattering function and Levinson‐type formula of Klein–Gordon equation

We considered the inverse problem of scattering theory for a boundary value problem on the half line generated by Klein–Gordon differential equation with a nonlinear spectral parameter‐dependent boundary condition. We defined the scattering data, and we proved the continuity of the scattering function S(λ); in a special case, the relation for the difference of the logarithm of the scattering function, which is called the Levinson‐type formula, was obtained. Copyright © 2012 John Wiley & Sons, Ltd.     

Huygens’ principle and iterative methods in inverse obstacle scattering

The inverse problem we consider in this paper is to determine the shape of an obstacle from the knowledge of the far field pattern for scattering of time-harmonic plane waves. In the case of scattering from a sound-soft obstacle, we will interpret Huygens’ principle as a system of two integral equations, named data and field equation, for the unknown boundary of the scatterer and the induced surface flux, i.e., the unknown normal derivative of the total field on the boundary. Reflecting the ill-posedness of the inverse obstacle scattering problem these integral equations are ill-posed. They are linear with respect to the unknown flux and nonlinear with respect to the unknown boundary and offer, in principle, three immediate possibilities for their iterative solution via linearization and regularization. In addition to presenting new results on injectivity and dense range for the linearized operators, the main purpose of this paper is to establish and illuminate relations between these three solution methods based on Huygens’ principle in inverse obstacle scattering. Furthermore, we will exhibit connections and differences to the traditional regularized Newton type iterations as applied to the boundary to far field map, including alternatives for the implementation of these Newton iterations.     

On inverse problems in electromagnetic field in classical mechanics at fixed energy

In this article, we consider inverse scattering and inverse boundary value problems at sufficiently large and fixed energy for the multidimensional relativistic and nonrelativistic Newton equations in a static external electromagnetic field (V, B), V∈C², B∈C¹ in classical mechanics. Developing the approach going back to Gerver-Nadinashvili 1983's work on an inverse problem of mechanics, we obtain, in particular, theorems of uniqueness.     

On the stability and convergence of the finite section method for integral equation formulations of rough surface scattering

We consider the Dirichlet and Robin boundary value problems for the Helmholtz equation in a non‐locally perturbed half‐plane, modelling time harmonic acoustic scattering of an incident field by, respectively, sound‐soft and impedance infinite rough surfaces.Recently proposed novel boundary integral equation formulations of these problems are discussed. It is usual in practical computations to truncate the infinite rough surface, solving a boundary integral equation on a finite section of the boundary, of length 2A, say. In the case of surfaces of small amplitude and slope we prove the stability and convergence as A→∞ of this approximation procedure. For surfaces of arbitrarily large amplitude and/or surface slope we prove stability and convergence of a modified finite section procedure in which the truncated boundary is ‘flattened’ in finite neighbourhoods of its two endpoints. Copyright © 2001 John Wiley & Sons, Ltd.     

The Impedance Boundary Value Problem for the Helmholtz Equation in a Half-Plane

We prove unique existence of solution for the impedance (or third) boundary value problem for the Helmholtz equation in a half-plane with arbitrary L_∞ boundary data. This problem is of interest as a model of outdoor sound propagation over inhomogeneous flat terrain and as a model of rough surface scattering. To formulate the problem and prove uniqueness of solution we introduce a novel radiation condition, a generalization of that used in plane wave scattering by one-dimensional diffraction gratings. To prove existence of solution and a limiting absorption principle we first reformulate the problem as an equivalent second kind boundary integral equation to which we apply a form of Fredholm alternative, utilizing recent results on the solvability of integral equations on the real line in [5]. © 1997 B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Reconstruction of cracks of different types from far‐field measurements

In this paper, we deal with the acoustic inverse scattering problem for reconstructing cracks of possibly different types from the far‐field map. The scattering problem models the diffraction of waves by thin two‐sided cylindrical screens. The cracks are characterized by their shapes, the type of boundary conditions and the boundary coefficients (surface impedance). We give explicit formulas of the indicator function of the probe method, which can be used to reconstruct the shape of the cracks, distinguish their types of boundary conditions, the two faces of each of them and reconstruct the possible material coefficients on them by using the far‐field map. To test the validity of these formulas, we present some numerical implementations for a single crack, which show the efficiency of the proposed method for suitably distributed surface impedances. The difficulties for numerically recovering the properties of the crack in the concave side as well as near the tips are presented and some explanations are given. Copyright © 2009 John Wiley & Sons, Ltd.     

Iterative orthogonal direction methods for Hermitian minimum norm solutions of two consistent matrix equations

The consistent conditions and the general expressions about the Hermitian solutions of the linear matrix equations AXB=C and (AX, XB)=(C, D) are studied in depth, where A, B, C and D are given matrices of suitable sizes. The Hermitian minimum F‐norm solutions are obtained for the matrix equations AXB=C and (AX, XB)=(C, D) by Moore–Penrose generalized inverse, respectively. For both matrix equations, we design iterative methods according to the fundamental idea of the classical conjugate direction method for the standard system of linear equations. Numerical results show that these iterative methods are feasible and effective in actual computations of the solutions of the above‐mentioned two matrix equations. Copyright © 2006 John Wiley & Sons, Ltd.     

Existence of solution in elastic wave scattering by unbounded rough surfaces

We consider the two‐dimensional problem of the scattering of a time‐harmonic wave, propagating in an homogeneous, isotropic elastic medium, by a rough surface on which the displacement is assumed to vanish. This surface is assumed to be given as the graph of a function ?∈C^1,1(?). Following up on earlier work establishing uniqueness of solution to this problem, existence of solution is studied via the boundary integral equation method. This requires a novel approach to the study of solvability of integral equations on the real line. The paper establishes the existence of a unique solution to the boundary integral equation formulation in the space of bounded and continuous functions as well as in all L^p spaces, p∈[1, ∞] and hence existence of solution to the elastic wave scattering problem. Copyright © 2002 John Wiley & Sons, Ltd.     

Determination of a control function in three‐dimensional parabolic equations by Legendre pseudospectral method

A Legendre pseudospectral method is proposed for solving approximately an inverse problem of determining an unknown control parameter p(t) which is the coefficient of the solution u(x, y, z, t) in a diffusion equation in a three‐dimensional region. The diffusion equation is to be solved subject to suitably prescribed initial‐boundary conditions. The presence of the unknown coefficient p(t) requires an extra condition. This extra condition considered as the integral overspecification over the spacial domain. For discretizing the problem, after homogenization of the boundary conditions, we apply the Legendre pseudospectral method in a matrix based manner. As a results a system of nonlinear differential algebraic equations is generated. Then by using suitable transformation, the problem will be converted to a homogeneous time varying system of linear ordinary differential equations. Also a pseudospectral method for efficient solving of the resulted system of ordinary differential equations is proposed. The solution of this system gives the approximation to values of u and p. The matrix based structure of the present method makes it easy to implement. Numerical experiments are presented to demonstrate the accuracy and the efficiency of the proposed computational procedure. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 74‐93, 2012     

An hp‐adaptive Newton‐Galerkin finite element procedure for semilinear boundary value problems

In this paper, we develop an hp‐adaptive procedure for the numerical solution of general, semilinear elliptic boundary value problems in 1d, with possible singular perturbations. Our approach combines both a prediction‐type adaptive Newton method and an hp‐version adaptive finite element discretization (based on a robust a posteriori residual analysis), thereby leading to a fully hp‐adaptive Newton–Galerkin scheme. Numerical experiments underline the robustness and reliability of the proposed approach for various examples. Copyright © 2016 John Wiley & Sons, Ltd.     

A hybrid method for sound-hard obstacle reconstruction

We are interested in solving the inverse problem of acoustic wave scattering to reconstruct the position and the shape of sound-hard obstacles from a given incident field and the corresponding far field pattern of the scattered field. The method we suggest is an extension of the hybrid method for the reconstruction of sound-soft cracks as presented in [R. Kress, P. Serranho, A hybrid method for two-dimensional crack reconstruction, Inverse Problems 21 (2005) 773–784] to the case of sound-hard obstacles. The designation of the method is justified by the fact that it can be interpreted as a hybrid between a regularized Newton method applied to a nonlinear operator equation with the operator that maps the unknown boundary onto the solution of the direct scattering problem and a decomposition method in the spirit of the potential method as described in [A. Kirsch, R. Kress, On an integral equation of the first kind in inverse acoustic scattering, in: Cannon, Hornung (Eds.), Inverse Problems, ISNM, vol. 77, 1986, pp. 93–102. Since the method does not require a forward solver for each Newton step its computational costs are reduced. By some numerical examples we illustrate the feasibility of the method.     

The method of fundamental solutions for the Oseen steady‐state viscous flow past obstacles of known or unknown shapes

In this paper, the steady‐state Oseen viscous flow equations past a known or unknown obstacle are solved numerically using the method of fundamental solutions (MFS), which is free of meshes, singularities, and numerical integrations. The direct problem is linear and well‐posed, whereas the inverse problem is nonlinear and ill‐posed. For the direct problem, the MFS computations of the fluid flow characteristics (velocity, pressure, drag, and lift coefficients) are in very good agreement with the previously published results obtained using other methods for the Oseen flow past circular and elliptic cylinders, as well as past two circular cylinders. In the inverse obstacle problem the boundary data and the internal measurement of the fluid velocity are minimized using the MATLAB^© optimization toolbox lsqnonlin routine. Regularization was found necessary in the case the measured data are contaminated with noise. Numerical results show accurate and stable reconstructions of various star‐shaped obstacles of circular, bean, or peanut cross‐section.     

A boundary integral equation for the transmission eigenvalue problem for Maxwell equation

We propose a new integral equation formulation to characterize and compute transmission eigenvalues in electromagnetic scattering. As opposed to the approach that was recently developed by Cakoni, Haddar and Meng (2015) which relies on a two‐by‐two system of boundary integral equations, our analysis is based on only one integral equation in terms of the electric‐to‐magnetic boundary trace operator that results in a simplification of the theory and in a considerable reduction of computational costs. We establish Fredholm properties of the integral operators and their analytic dependence on the wave number. Further, we use the numerical algorithm for analytic nonlinear eigenvalue problems that was recently proposed by Beyn (2012) for the numerical computation of the transmission eigenvalues via this new integral equation.     

Retrieving the time-dependent thermal conductivity of an orthotropic rectangular conductor

The aim of this paper is to determine the thermal properties of an orthotropic planar structure characterized by the thermal conductivity tensor in the coordinate system of the main directions (Oxy) being diagonal. In particular, we consider retrieving the time-dependent thermal conductivity components of an orthotropic rectangular conductor from nonlocal overspecified heat flux conditions. Since only boundary measurements are considered, this inverse formulation belongs to the desirable approach of non-destructive testing of materials. The unique solvability of this inverse coefficient problem is proved based on the Schauder fixed point theorem and the theory of Volterra integral equations of the second kind. Furthermore, the numerical reconstruction based on a nonlinear least-squares minimization is performed using the MATLAB optimization toolbox routine lsqnonlin. Numerical results are presented and discussed in order to illustrate the performance of the inversion for orthotropic parameter identification.     

Mathematical basis of the linear sampling method for partially coated obstacles

We consider the inverse scattering problem of determining the shape of a partially coated obstacle D. To this end, we solve a scattering problem for the Helmholtz equation where the scattered field satisfies mixed Dirichlet–Neumann-impedance boundary conditions on the Lipschitz boundary of the scatterer D. Based on the analysis of the boundary integral system to the direct scattering problem, we propose how to reconstruct the shape of the obstacle D by using the linear sampling method.     

Calculation of the Defect Numbers of the Generalized Hilbert and Carleman Boundary Value Problems with Linear Fractional Carleman Shift

The defect numbers of the generalized Hilbert and Carleman boundary value problems with a direct or an inverse linear fractional Carleman shift of order 2 (α (α (t)) ≡ t) on the unit circle are computed. The approach followed consists of the reduction of the mentioned problems to singular integral equations with linear fractional Carleman shift and of the factorization of Hermitian matrix functions with negative determinant.