 Award ID(s):
 1912704
 NSFPAR ID:
 10182373
 Date Published:
 Journal Name:
 SIAM journal on mathematical analysis
 Volume:
 51
 ISSN:
 00361410
 Page Range / eLocation ID:
 45704603
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

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 wellposed 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

In this paper, we consider the inverse scattering problem for recovering either an isotropic or anisotropic scatterer from the measured scattered field initiated by a point source. We propose two new imaging functionals for solving the inverse problem. The first one employs a 'farfield' transform to the data which we then use to derive and provide an explicit decay rate for the imaging functional. In order to analyze the behavior of this imaging functional we use the factorization of the near field operator as well as the FunkHecke integral identity. For the second imaging functional the Cauchy data is used to define the functional and its behavior is analyzed using the Green's identities. Numerical experiments are given in two dimensions for both isotropic and anisotropic scatterers.

Abstract We consider the mapping properties of the integral operator arising in nonlocal slender body theory (SBT) for the model geometry of a straight, periodic filament. It is well known that the classical singular SBT integral operator suffers from high wavenumber instabilities, making it unsuitable for approximating the
slender body inverse problem , where the fiber velocity is prescribed and the integral operator must be inverted to find the force density along the fiber. Regularizations of the integral operator must therefore be used instead. Here, we consider two regularization methods: spectral truncation and the‐regularization of Tornberg and Shelley (2004). We compare the mapping properties of these approximations to the underlying partial differential equation (PDE) solution, which for the inverse problem is simply the Stokes Dirichlet problem with data constrained to be constant on cross sections. For the straight‐but‐periodic fiber with constant radius , we explicitly calculate the spectrum of the operator mapping fiber velocity to force for both the PDE and the approximations. We prove that the spectrum of the original SBT operator agrees closely with the PDE operator at low wavenumbers but differs at high frequencies, allowing us to define a truncated approximation with a wavenumber cutoff . For both the truncated and ‐regularized approximations, we obtain rigorous ‐based convergence to the PDE solution as : A fiber velocity with regularity gives convergence, while a fiber velocity with at least regularity yields convergence. Moreover, we determine the dependence of the ‐regularized error estimate on the regularization parameter . 
This article presents a numerical strategy for actively manipulating electromagnetic (EM) fields in layered media. In particular, we develop a scheme to characterize an EM source that will generate some predetermined field patterns in prescribed disjoint exterior regions in layered media. The proposed question of specifying such an EM source is not an inverse source problem (ISP) since the existence of a solution is not guaranteed. Moreover, our problem allows for the possibility of prescribing different EM fields in mutually disjoint exterior regions. This question involves a linear inverse problem that requires solving a severely illposed optimization problem (i.e. suffering from possible nonexistence or nonuniqueness of a solution). The forward operator is defined by expressing the EM fields as a function of the current at the source using the layered media Green’s function (LMGF), accounting for the physical parameters of the layered media. This results to integral equations that are then discretized using the method of moments (MoM), yielding an illposed system of linear equations. Unlike in ISPs, stability with respect to data is not an issue here since no data is measured. Rather, stability with respect to input current approximation is important. To get such stable solutions, we applied two regularization methods, namely, the truncated singular value decomposition (TSVD) method and the Tikhonov regularization method with the Morozov Discrepancy Principle. We performed several numerical simulations to support the theoretical framework and analyzes, and to demonstrate the accuracy and feasibility of the proposed numerical algorithms.more » « less

The aim of this paper is to solve an important inverse source problem which arises from the wellknown 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