Source author record

Peijun Li

Peijun Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

A spectral boundary integral method for the elastic obstacle scattering problem in three dimensions

In this paper, we consider the scattering of a plane wave by a rigid obstacle embedded in a homogeneous and isotropic elastic medium in three dimensions. Based on the Helmholtz decomposition, the elastic scattering problem is reduced to a coupled boundary value problem for the Helmholtz and Maxwell equations. A novel system of boundary integral equations is formulated and a spectral boundary integral method is developed for the coupled boundary value problem. Numerical experiments are presented to demonstrate the superior performance of the proposed method.

preprint2022arXiv

An Adaptive Finite Element DtN Method for Maxwell's Equations

This paper is concerned with a numerical solution to the scattering of a time-harmonic electromagnetic wave by a bounded and impenetrable obstacle in three dimensions. The electromagnetic wave propagation is modeled by a boundary value problem of Maxwell's equations in the exterior domain of the obstacle. Based on the Dirichlet-to-Neumann (DtN) operator, which is defined by an infinite series, an exact transparent boundary condition is introduced and the scattering problem is reduced equivalently into a bounded domain. An a posteriori error estimate based adaptive finite element DtN method is developed to solve the discrete variational problem, where the DtN operator is truncated into a sum of finitely many terms. The a posteriori error estimate takes into account both the finite element approximation error and the truncation error of the DtN operator. The latter is shown to decay exponentially with respect to the truncation parameter. Numerical experiments are presented to illustrate the effectiveness of the proposed method.

preprint2022arXiv

Stability for the multi-dimensional Borg--Levinson theorem of the biharmonic operator

In this paper, we prove a conditional Hölder stability estimate for the inverse spectral problem of the biharmonic operator. The proof employs the resolvent estimate and a Weyl-type law for the biharmonic operator which were obtained by the authors in \cite{LYZ}. This work extends nontrivially the result in \cite{stefanov} from the second order Schrödinger operator to the fourth order biharmonic operator.

preprint2021arXiv

An adaptive finite element DtN method for the elastic wave scattering problem in three dimensions

Consider the elastic scattering of an incident wave by a rigid obstacle in three dimensions, which is formulated as an exterior problem for the Navier equation. By constructing a Dirichlet-to-Neumann (DtN) operator and introducing a transparent boundary condition, the scattering problem is reduced equivalently to a boundary value problem in a bounded domain. The discrete problem with the truncated DtN operator is solved by using the a posteriori error estimate based adaptive finite element method. The estimate takes account of both the finite element approximation error and the truncation error of the DtN operator, where the latter is shown to converge exponentially with respect to the truncation parameter. Moreover, the generalized Woodbury matrix identity is utilized to solve the resulting linear system efficiently. Numerical experiments are presented to demonstrate the superior performance of the proposed method.

preprint2021arXiv

An inverse source problem for the stochastic wave equation

This paper is concerned with an inverse source problem for the stochastic wave equation driven by a fractional Brownian motion. Given the random source, the direct problem is to study the solution of the stochastic wave equation. The inverse problem is to determine the statistical properties of the source from the expectation and covariance of the final-time data. For the direct problem, it is shown to be well-posed with a unique mild solution. For the inverse problem, the uniqueness is proved for a certain class of functions and the instability is characterized. Numerical experiments are presented to illustrate the reconstructions by using a truncation-based regularization method.

preprint2021arXiv

Inverse random potential scattering for elastic waves

This paper is concerned with the inverse elastic scattering problem for a random potential in three dimensions. Interpreted as a distribution, the potential is assumed to be a microlocally isotropic Gaussian random field whose covariance operator is a classical pseudo-differential operator. Given the potential, the direct scattering problem is shown to be well-posed in the distribution sense by studying the equivalent Lippmann--Schwinger integral equation. For the inverse scattering problem, we demonstrate that the microlocal strength of the random potential can be uniquely determined with probability one by a single realization of the high frequency limit of the averaged compressional or shear backscattered far-field pattern of the scattered wave. The analysis employs the integral operator theory, the Born approximation in the high frequency regime, the microlocal analysis for the Fourier integral operators, and the ergodicity of the wave field.

preprint2021arXiv

Stability for an inverse source problem of the biharmonic operator

In this paper, we study for the first time the stability of the inverse source problem for the biharmonic operator with a compactly supported potential in $\mathbb R^3$. Firstly, to connect the boundary data with the unknown source, we shall consider an eigenvalue problem for the bi-Schr$\ddot{\rm o}$dinger operator $Δ^2 + V(x)$ on a ball which contains the support of the potential $V$. We prove a Weyl-type law for the upper bounds of spherical normal derivatives of both the eigenfunctions $ϕ$ and their Laplacian $Δϕ$ corresponding to the bi-Schr$\ddot{\rm o}$dinger operator. This type of upper bounds was proved by Hassell and Tao for the Schr$\ddot{\rm o}$dinger operator. Secondly, we investigate the meromorphic continuation of the resolvent of the bi-Schr$\ddot{\rm o}$dinger operator and prove the existence of a resonance-free region and an estimate of $L^2_{\rm comp} - L^2_{\rm loc}$ type for the resolvent. As an application, we prove a bound of the analytic continuation of the data from the given data to the higher frequency data. Finally, we derive the stability estimate which consists of the Lipschitz type data discrepancy and the high frequency tail of the source function, where the latter decreases as the upper bound of the frequency increases.

preprint2021arXiv

Stability for an inverse source problem of the damped biharmonic plate equation

This paper is concerned with the stability of the inverse source problem for the damped biharmonic plate equation in three dimensions. The stability estimate consists of the Lipschitz type data discrepancy and the high frequency tail of the source function, where the latter decreases as the upper bound of the frequency increases. The stability also shows exponential dependence on the constant damping coefficient. The analysis employs Carleman estimates and time decay estimates for the damped plate wave equation to obtain an exact observability bound and depends on the study of the resonance-free region and an upper bound of the resolvent of the biharmonic operator with respect to the complex wavenumber.

preprint2020arXiv

A highly accurate boundary integral method for the elastic obstacle scattering problem

Consider the scattering of a time-harmonic plane wave by a rigid obstacle embedded in a homogeneous and isotropic elastic medium in two dimensions. In this paper, a novel boundary integral formulation is proposed and its highly accurate numerical method is developed for the elastic obstacle scattering problem. More specifically, based on the Helmholtz decomposition, the model problem is reduced to a coupled boundary integral equation with singular kernels. A regularized system is constructed in order to handle the degenerated integral operators. The semi-discrete and full-discrete schemes are studied for the boundary integral system by using the trigonometric collocation method. Convergence is established for the numerical schemes in some appropriate Sobolev spaces. Numerical experiments are presented for both smooth and nonsmooth obstacles to demonstrate the superior performance of the proposed method.

preprint2020arXiv

An adaptive finite element DtN method for the elastic wave scattering by biperiodic structures

Consider the scattering of a time-harmonic elastic plane wave by a bi-periodic rigid surface. The displacement of elastic wave motion is modeled by the three-dimensional Navier equation in an open domain above the surface. Based on the Dirichlet-to-Neumann (DtN) operator, which is given as an infinite series, an exact transparent boundary condition is introduced and the scattering problem is formulated equivalently into a boundary value problem in a bounded domain. An a posteriori error estimate based adaptive finite element DtN method is proposed to solve the discrete variational problem where the DtN operator is truncated into a finite number of terms. The a posteriori error estimate takes account of the finite element approximation error and the truncation error of the DtN operator which is shown to decay exponentially with respect to the truncation parameter. Numerical experiments are presented to illustrate the effectiveness of the proposed method.

preprint2020arXiv

An adaptive finite element DtN method for the open cavity scattering problems

Consider the scattering of a time-harmonic electromagnetic plane wave by an open cavity which is embedded in a perfectly electrically conducting infinite ground plane. This paper concerns the numerical solutions of the open cavity scattering problems in both transverse magnetic and transverse electric polarizations. Based on the Dirichlet-to-Neumann (DtN) map for each polarization, a transparent boundary condition is imposed to reduce the scattering problem equivalently into a boundary value problem in a bounded domain. An a posteriori error estimate based adaptive finite element DtN method is proposed. The estimate consists of the finite element approximation error and the truncation error of the DtN operator, which is shown to decay exponentially with respect to the truncation parameter. Numerical experiments are presented for both polarizations to illustrate the competitive behavior of the adaptive method.

preprint2020arXiv

An adaptive finite element DtN method for the three-dimensional acoustic scattering problem

This paper is concerned with a numerical solution of the acoustic scattering by a bounded impenetrable obstacle in three dimensions. The obstacle scattering problem is formulated as a boundary value problem in a bounded domain by using a Dirichlet-to-Neumann (DtN) operator. An a posteriori error estimate is derived for the finite element method with the truncated DtN operator. The a posteriori error estimate consists of the finite element approximation error and the truncation error of the DtN operator, where the latter is shown to decay exponentially with respect to the truncation parameter. Based on the a posteriori error estimate, an adaptive finite element method is developed for the obstacle scattering problem. The truncation parameter is determined by the truncation error of the DtN operator and the mesh elements for local refinement are marked through the finite element approximation error. Numerical experiments are presented to demonstrate the effectiveness of the proposed method.

preprint2020arXiv

An inverse random source problem for Maxwell's equations

This paper is concerned with an inverse random source problem for the three-dimensional time-harmonic Maxwell equations. The source is assumed to be a centered complex-valued Gaussian vector field with correlated components, and its covariance operator is a pseudo-differential operator. The well-posedness of the direct source scattering problem is established and the regularity of the electromagnetic field is given. For the inverse source scattering problem, the micro-correlation strength matrix of the covariance operator is shown to be uniquely determined by the high frequency limit of the expectation of the electric field measured in an open bounded domain disjoint with the support of the source. In particular, we show that the diagonal entries of the strength matrix can be uniquely determined by only using the amplitude of the electric field. Moreover, this result is extended to the almost surely sense by deducing an ergodic relation for the electric field over the frequencies.

preprint2019arXiv

An inverse acoustic-elastic interaction problem with phased or phaseless far-field data

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 concerns 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 establish 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 also proposed for the phaseless inverse problem. Numerical experiments are provided to demonstrate the effectiveness and robustness of the proposed methods.

preprint2019arXiv

An inverse random source problem for the time fractional diffusion equation driven by a fractional Brownian motion

This paper is concerned with the mathematical analysis of the inverse random source problem for the time fractional diffusion equation, where the source is assumed to be 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.

preprint2019arXiv

Convergence of an adaptive finite element DtN method for the elastic wave scattering by periodic structures

Consider the scattering of a time-harmonic elastic plane wave by a periodic rigid surface. The elastic wave propagation is governed by the two-dimensional Navier equation. Based on a Dirichlet-to-Neumann (DtN) map, a transparent boundary condition (TBC) is introduced to reduce the scattering problem into a boundary value problem in a bounded domain. By using the finite element method, the discrete problem is considered, where the TBC is replaced by the truncated DtN map. A new duality argument is developed to derive the a posteriori error estimate, which contains both the finite element approximation error and the DtN truncation error. An a posteriori error estimate based adaptive finite element algorithm is developed to solve the elastic surface scattering problem. Numerical experiments are presented to demonstrate the effectiveness of the proposed method.

preprint2019arXiv

Inverse scattering for the one-dimensional Helmholtz equation with piecewise constant wave speed

This paper analyzes inverse scattering for the one-dimensional Helmholtz equation in the case where the wave speed is piecewise constant. Scattering data recorded for an arbitrarily small interval of frequencies is shown to determine the wave speed uniquely, and a direct reconstruction algorithm is presented. The algorithm is exact provided data is recorded for a sufficiently wide range of frequencies and the jump points of the wave speed are equally spaced with respect to travel time. Numerical examples show that the algorithm works also in the general case of arbitrary wave speed (either with jumps or continuously varying etc.) giving progressively more accurate approximations as the range of recorded frequencies increases. A key underlying theoretical insight is to associate scattering data to compositions of automorphisms of the unit disk, which are in turn related to orthogonal polynomials on the unit circle. The algorithm exploits the three-term recurrence of orthogonal polynomials to reduce the required computation.

preprint2018arXiv

A Fast Solver for the Elastic Scattering of Multiple Particles

Consider the elastic scattering of a time-harmonic wave by multiple well separated rigid particles in two dimensions. To avoid using the complex Green's tensor of the elastic wave equation, we utilize the Helmholtz decomposition to convert the boundary value problem of the elastic wave equation into a coupled boundary value problem of Helmholtz equations. Based on single, double, and combined layer potentials with the simpler Green's function of the Helmholtz equation, we present three different boundary integral equations for the coupled boundary value problem. The well-posedness of the new integral equations are established. Computationally, a scattering matrix based method is proposed to evaluate the elastic wave for arbitrarily shaped particles. The method uses the local expansion for the incident wave and the multipole expansion for the scattered wave. The linear system of algebraic equations is solved by GMRES with fast multipole method (FMM) acceleration. Numerical results show that the method is fast and highly accurate for solving the elastic scattering problem with multiple particles.

preprint2016arXiv

An adaptive finite element PML method for the elastic wave scattering problem in periodic structures

An adaptive finite element method is presented for the elastic scattering of a time-harmonic plane wave by a periodic surface. First, the unbounded physical domain is truncated into a bounded computational domain by introducing the perfectly matched layer (PML) technique. The well-posedness and exponential convergence of the solution are established for the truncated PML problem by developing an equivalent transparent boundary condition. Second, an a posteriori error estimate is deduced for the discrete problem and is used to determine the finite elements for refinements and to determine the PML parameters. Numerical experiments are included to demonstrate the competitive behavior of the proposed adaptive method.

preprint2016arXiv

Analysis of Time-Domain Scattering by Periodic Structures

This paper is devoted to the mathematical analysis of a time-domain electromagnetic scattering by periodic structures which are known as diffraction gratings. The scattering problem is reduced equivalently into an initial-boundary value problem in a bounded domain by using an exact transparent boundary condition. The well-posedness and stability of the solution are established for the reduced problem. Moreover, a priori energy estimates are obtained with minimum regularity requirement for the data and explicit dependence on the time.

preprint2016arXiv

Analysis of Transient Acoustic-Elastic Interaction in an Unbounded Structure

Consider the wave propagation in a two-layered medium consisting of a homogeneous compressible air or fluid on top of a homogeneous isotropic elastic solid. The interface between the two layers is assumed to be an unbounded rough surface. This paper concerns the time-domain analysis of such an acoustic-elastic interaction problem in an unbounded structure in three dimensions. Using an exact transparent boundary condition and suitable interface conditions, we study an initial-boundary value problem for the coupling of the Helmholtz equation and the Navier equation. The well-posedness and stability are established for the reduced problem. Our proof is based on the method of energy, the Lax--Milgram lemma, and the inversion theorem of the Laplace transform. Moreover, a priori estimates with explicit dependence on the time are achieved for the quantities of acoustic pressure and elastic displacement by taking special test functions for the time-domain variational problem.

preprint2016arXiv

Convergence of the PML solution for elastic wave scattering by biperiodic structures

This paper is concerned with the analysis of elastic wave scattering of a time-harmonic plane wave by a biperiodic rigid surface, where the wave propagation is governed by the three-dimensional Navier equation. An exact transparent boundary condition is developed to reduce the scattering problem equivalently into a boundary value problem in a bounded domain. The perfectly matched layer (PML) technique is adopted to truncate the unbounded physical domain into a bounded computational domain. The well-posedness and exponential convergence of the solution are established for the truncated PML problem by developing a PML equivalent transparent boundary condition. The proofs rely on a careful study of the error between the two transparent boundary operators. The work significantly extend the results from the one-dimensional periodic structures to the two-dimensional biperiodic structures. Numerical experiments are included to demonstrate the competitive behavior of the proposed method.

preprint2016arXiv

Electromagnetic Scattering for Time-Domain Maxwell's Equations in an Unbounded Structure

The goal of this work is to study the electromagnetic scattering problem of time-domain Maxwell's equations in an unbounded structure. An exact transparent boundary condition is developed to reformulate the scattering problem into an initial-boundary value problem in an infinite rectangular slab. The well-posedness and stability are established for the reduced problem. Our proof is based on the method of energy, the Lax--Milgram lemma, and the inversion theorem of the Laplace transform. Moreover, a priori estimates with explicit dependence on the time are achieved for the electric field by directly studying the time-domain Maxwell equations.

preprint2016arXiv

Increasing stability for the inverse source scattering problem with multi-frequencies

Consider the scattering of the two- or three-dimensional Helmholtz equation where the source of the electric current density is assumed to be compactly supported in a ball. This paper concerns the stability analysis of the inverse source scattering problem which is to reconstruct the source function. Our results show that increasing stability can be obtained for the inverse problem by using only the Dirichlet boundary data with multi-frequencies.

preprint2016arXiv

Inverse Random Source Scattering for Elastic Waves

This paper is concerned with the direct and inverse random source scattering problems for elastic waves where the source is assumed to be driven by an additive white noise. Given the source, the direct problem is to determine the displacement of the random wave field. The inverse problem is to reconstruct the mean and variance of the random source from the boundary measurement of the wave field at multiple frequencies. The direct problem is shown to have a unique mild solution by using a constructive proof. Based on the explicit mild solution, Fredholm integral equations of the first kind are deduced for the inverse problem. The regularized Kaczmarz method is presented to solve the ill-posed integral equations. Numerical experiments are included to demonstrate the effectiveness of the proposed method.

preprint2016arXiv

Stability on the Inverse Random Source Scattering Problem for the One-Dimensional Helmholtz Equation

Consider the one-dimensional stochastic Helmholtz equation where the source is assumed to be driven by the white noise. This paper concerns the stability analysis of the inverse random source problem which is to reconstruct the statistical properties of the source such as the mean and variance. Our results show that increasing stability can be obtained for the inverse problem by using suitable boundary data with multi-frequencies.