BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Peijun Li

Peijun Li contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.APmath-phmath.MP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

18 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access