Source author record

Jiguang Sun

Jiguang Sun 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

12works
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

12 published item(s)

preprint2020arXiv

A Multilevel Spectral Indicator Method for Eigenvalues of Large Non-Hermitian Matrices

Recently a novel family of eigensolvers, called spectral indicator methods (SIMs), was proposed. Given a region on the complex plane, SIMs first compute an indicator by the spectral projection. The indicator is used to test if the region contains eigenvalue(s). Then the region containing eigenvalues(s) is subdivided and tested. The procedure is repeated until the eigenvalues are identified within a specified precision. In this paper, using Cayley transformation and Krylov subspaces, a memory efficient multilevel eigensolver is proposed. The method uses less memory compared with the early versions of SIMs and is particularly suitable to compute many eigenvalues of large sparse (non-Hermitian) matrices. Several examples are presented for demonstration.

preprint2020arXiv

A new finite element approach for the Dirichlet eigenvalue problem

In this paper, we propose a new finite element approach, which is different than the classic Babuska-Osborn theory, to approximate Dirichlet eigenvalues. The Dirichlet eigenvalue problem is formulated as the eigenvalue problem of a holomorphic Fredholm operator function of index zero. Using conforming finite elements, the convergence is proved using the abstract approximation theory for holomorphic operator functions. The spectral indicator method is employed to compute the eigenvalues. A numerical example is presented to validate the theory.

preprint2020arXiv

A priori and a posteriori error estimates for the quad-curl eigenvalue problem

In this paper, we propose a new family of H(curl^2)-conforming elements for the quad-curl eigenvalue problem in 2D. The accuracy of this family is one order higher than that in [32]. We prove a priori and a posteriori error estimates. The a priori estimate of the eigenvalue with a convergence order 2(s-1) is obtained if the eigenvector u\in H^{s+1}(Ω). For the a posteriori estimate, by analyzing the associated source problem, we obtain lower and upper bounds for the eigenvector in an energy norm and an upper bound for the eigenvalues. Numerical examples are presented for validation.

preprint2020arXiv

Finite Element Approximation of the Modified Maxwell's Stekloff Eigenvalues

The modified Maxwell's Stekloff eigenvalue problem arises recently from the inverse electromagnetic scattering theory for inhomogeneous media. This paper contains a rigorous analysis of both the eigenvalue problem and the associated source problem on Lipschitz polyhedra. A new finite element method is proposed to compute Stekloff eigenvalues. By applying the Babuska-Osborn theory, we prove an error estimate without additional regularity assumptions. Numerical results are presented for validation.

preprint2020arXiv

Finite Element Approximation of Transmission Eigenvalues for Anisotropic Media

The transmission eigenvalue problem arises from the inverse scattering theory for inhomogeneous media and has important applications in many qualitative methods. The problem is posted as a system of two second order partial differential equations and is essentially nonlinear, non-selfadjoint, and of higher order. It is nontrivial to develop effective numerical methods and the proof of convergence is challenging. In this paper, we formulate the transmission eigenvalue problem for anisotropic media as an eigenvalue problem of a holomorphic Fredholm operator function of index zero. The Lagrange finite elements are used for discretization and the convergence is proved using the abstract approximation theory for holomorphic operator functions. A spectral indicator method is developed to compute the eigenvalues. Numerical examples are presented for validation.

preprint2020arXiv

Finite Element Calculation of Photonic Band Structures for Frequency Dependent Materials

We consider the calculation of the band structure of frequency dependent photonic crystals. The associated eigenvalue problem is nonlinear and it is challenging to develop effective convergent numerical methods. In this paper, the band structure problem is formulated as the eigenvalue problem of a holomorphic Fredholm operator function of index zero. Lagrange finite elements are used to discretize the operator function. Then the convergence of the eigenvalues is proved using the abstract approximation theory for holomorphic operator functions. A spectral indicator method is developed to practically compute the eigenvalues. Numerical examples are presented to validate the theory and show the effectiveness of the proposed method.

preprint2020arXiv

Quality-Bayesian approach to inverse acoustic source problems with partial data

A quality-Bayesian approach, combining the direct sampling method and the Bayesian inversion, is proposed to reconstruct the locations and intensities of the unknown acoustic sources using partial data. First, we extend the direct sampling method by constructing a new indicator function to obtain the approximate locations of the sources. The behavior of the indicator is analyzed. Second, the inverse problem is formulated as a statistical inference problem using the Bayes' formula. The well-posedness of the posterior distribution is proved. The source locations obtained in the first step are coded in the priors. Then an Metropolis-Hastings MCMC algorithm is used to explore the posterior density. Both steps use the same physical model and measured data. Numerical experiments show that the proposed method is effective even with partial data.

preprint2019arXiv

Revisiting the Jones eigenproblem in fluid-structure interaction

The Jones eigenvalue problem first described by D.S. Jones in 1983 concerns unusual modes in bounded elastic bodies: time-harmonic displacements whose tractions and normal components are both identically zero on the boundary. This problem is usually associated with a lack of unique solvability for certain models of fluid-structure interaction. The boundary conditions in this problem appear, at first glance, to rule out {\it any} non-trivial modes unless the domain possesses significant geometric symmetries. Indeed, Jones modes were shown to not be possible in most $C^\infty$ domains (see article by T. Hargé 1990). However, we should in this paper that while the existence of Jones modes sensitively depends on the domain geometry, such modes {\it do} exist in a broad class of domains. This paper presents the first detailed theoretical and computational investigation of this eigenvalue problem in Lipschitz domains. We also analytically demonstrate Jones modes on some simple geometries.

preprint2016arXiv

A spectral projection method for transmission eigenvalues

In this paper, we consider a nonlinear integral eigenvalue problem, which is a reformulation of the transmission eigenvalue problem arising in the inverse scattering theory. The boundary element method is employed for discretization, which leads to a generalized matrix eigenvalue problem. We propose a novel method based on the spectral projection. The method probes a given region on the complex plane using contour integrals and decides if the region contains eigenvalue(s) or not. It is particularly suitable to test if zero is an eigenvalue of the generalized eigenvalue problem, which in turn implies that the associated wavenumber is a transmission eigenvalue. Effectiveness and efficiency of the new method are demonstrated by numerical examples.

preprint2016arXiv

Direct and Inverse Elastic Scattering From Anisotropic Media

Assume a time-harmonic elastic wave is incident onto a penetrable anisotropic body embedded into a homogeneous isotropic background medium. We propose an equivalent variational formulation in a truncated bounded domain and show the uniqueness and existence of weak solutions by applying the Fredholm alternative and using properties of the Dirichlet-to-Neumann map in both two and three dimensions. The Fréchet derivative of the near-field solution operator with respect to the boundary of the scatterer is derived. As an application, we design a descent algorithm for recovering the interface from the near-field data of one or several incident directions and frequencies. Numerical examples in 2D are presented to show the validity and accuracy of our methods.

preprint2015arXiv

Recursive integral method for transmission eigenvalues

Recently, a new eigenvalue problem, called the transmission eigenvalue problem, has attracted many researchers. The problem arose in inverse scattering theory for inhomogeneous media and has important applications in a variety of inverse problems for target identification and nondestructive testing. The problem is numerically challenging because it is non-selfadjoint and nonlinear. In this paper, we propose a recursive integral method for computing transmission eigenvalues from a finite element discretization of the continuous problem. The method, which overcomes some difficulties of existing methods, is based on eigenprojectors of compact operators. It is self-correcting, can separate nearby eigenvalues, and does not require an initial approximation based on some a priori spectral information. These features make the method well suited for the transmission eigenvalue problem whose spectrum is complicated. Numerical examples show that the method is effective and robust.

preprint2013arXiv

A mixed FEM for the quad-curl eigenvalue problem

The quad-curl problem arises in the study of the electromagnetic interior transmission problem and magnetohydrodynamics (MHD). In this paper, we study the quad-curl eigenvalue problem and propose a mixed method using edge elements for the computation of the eigenvalues. To the author's knowledge, it is the first numerical treatment for the quad-curl eigenvalue problem. Under suitable assumptions on the domain and mesh, we prove the optimal convergence. In addition, we show that the divergence-free condition can be bypassed. Numerical results are provided to show the viability of the method.