skip to main content

Title: Reconstructing a space-dependent source term via the quasi-reversibility method
The aim of this paper is to solve an important inverse source problem which arises from the well-known inverse scattering problem. We propose to truncate the Fourier series of the solution to the governing equation with respect to a special basis of L2. By this, we obtain a system of linear elliptic equations. Solutions to this system are the Fourier coefficients of the solution to the governing equation. After computing these Fourier coefficients, we can directly find the desired source function. Numerical examples are presented.  more » « less
Award ID(s):
Author(s) / Creator(s):
Date Published:
Journal Name:
Contemporary mathematics
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Consider the inverse random source scattering problem for the two-dimensional time-harmonic elastic wave equation with a linear load. The source is modeled as a microlocally isotropic generalized Gaussian random function whose covariance operator is a classical pseudodifferential operator. The goal is to recover the principal symbol of the covariance operator from the displacement measured in a domain away from the source. For such a distributional source, we show that the direct problem has a unique solution by introducing an equivalent Lippmann--Schwinger integral equation. For the inverse problem, we demonstrate that, with probability one, the principal symbol of the covariance operator can be uniquely determined by the amplitude of the displacement averaged over the frequency band, generated by a single realization of the random source. The analysis employs the Born approximation, asymptotic expansions of the Green tensor, and microlocal analysis of the Fourier integral operators. 
    more » « less
  2. Nguyen, Dinh-Liem ; Nguyen, Loc ; Nguyen, Thi-Phong (Ed.)
    This paper is concerned with the numerical solution to the direct and inverse electromagnetic scattering problem for bi-anisotropic periodic structures. The direct problem can be reformulated as an integro-di erential equation. We study the existence and uniqueness of solution to the latter equation and analyze a spectral Galerkin method to solve it. This spectral method is based on a periodization technique which allows us to avoid the evaluation of the quasiperiodic Green's tensor and to use the fast Fourier transform in the numerical implementation of the method. For the inverse problem, we study the orthogonality sampling method to reconstruct the periodic structures from scattering data generated by only two incident fields. The sampling method is fast, simple to implement, regularization free, and very robust against noise in the data. Numerical examples for both direct and inverse problems are presented to examine the efficiency of the numerical solvers. 
    more » « less
  3. This paper is concerned with the mathematical analysis of an inverse random source problem for the time fractional diffusion equation, where the source is driven by a fractional Brownian motion. Given the random source, the direct problem is to study the stochastic time fractional diffusion equation. The inverse problem is to determine the statistical properties of the source from the expectation and variance of the final time data. For the direct problem, we show that it is well-posed and has a unique mild solution under a certain condition. For the inverse problem, the uniqueness is proved and the instability is characterized. The major ingredients of the analysis are based on the properties of the Mittag–Leffler function and the stochastic integrals associated with the fractional Brownian motion. 
    more » « less
  4. An analysis of the structural dynamic response under uncertainty is presented. Uncertainties in load and material are modelled as intervals exploiting the interval finite element method (IFEM). To reduce overestimation and increase the computational efficiency of the solution, we do not solve the dynamic problem by an explicit step-by-step time integration scheme. Instead, our approach solves for the structural variables in the whole time domain simultaneously by an implicit scheme using discrete Fourier transform and its inverse (DFT and IDFT). Non-trivial initial conditions are handled by modifying the right-hand side of the governing equation. To further reduce overestimation, a new decomposition strategy is applied to the IFEM matrices, and both primary and derived quantities are solved simultaneously. The final solution is obtained using an iterative enclosure method, and in our numerical examples the exact solution is enclosed at minimal computational cost. 
    more » « less
  5. null (Ed.)
    This paper is concerned with the inverse scattering problem which aims to determine the spatially distributed dielectric constant coefficient of the 2D Helmholtz equation from multifrequency backscatter data associated with a single direction of the incident plane wave. We propose a globally convergent convexification numerical algorithm to solve this nonlinear and ill-posed inverse problem. The key advantage of our method over conventional optimization approaches is that it does not require a good first guess about the solution. First, we eliminate the coefficient from the Helmholtz equation using a change of variables. Next, using a truncated expansion with respect to a special Fourier basis, we approximately reformulate the inverse problem as a system of quasilinear elliptic PDEs, which can be numerically solved by a weighted quasi-reversibility approach. The cost functional for the weighted quasi-reversibility method is constructed as a Tikhonov-like functional that involves a Carleman Weight Function. Our numerical study shows that, using a version of the gradient descent method, one can find the minimizer of this Tikhonov-like functional without any advanced a priori knowledge about it. 
    more » « less