Source author record

Ruming Zhang

Ruming Zhang 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

9works
4topics
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

9 published item(s)

preprint2024arXiv

Fast convergent PML method for scattering with periodic surfaces: the exceptional case

In the author's previous paper (Zhang et al. 2022), exponential convergence was proved for the perfectly matched layers (PML) approximation of scattering problems with periodic surfaces in 2D. However, due to the overlapping of singularities, an exceptional case, i.e., when the wave number is a half integer, has to be excluded in the proof. However, numerical results for these cases still have fast convergence rate and this motivates us to go deeper into these cases. In this paper, we focus on these cases and prove that the fast convergence result for the discretized form. Numerical examples are also presented to support our theoretical results.

preprint2024arXiv

The PML-Method for a Scattering Problem for a Local Perturbation of an Open Periodic Waveguide

The perfectly matched layers method is a well known truncation technique for its efficiency and convenience in numerical implementations of wave scattering problems in unbounded domains. In this paper, we study the convergence of the perfectly matched layers (PML) for wave scattering from a local perturbation of an open waveguide in the half space above the real line, where the refractive index is a function which is periodic along the axis of the waveguide and equals to one above a finite height. The problem is challenging due to the existence of guided waves, and a typical way to deal with the difficulty is to apply the limiting absorption principle. Based on the Floquet-Bloch transform and a curve deformation theory, the solution from the limiting absorption principle is rewritten as the integral of a coupled family of quasi-periodic problems with respect to the quasi-periodicity parameter on a particularly designed curve. By comparing the Dirichlet-to-Neumann maps on a straight line above the locally perturbed periodic layer, we finally show that the PML method converges exponentially with respect to the PML parameter. Finally, the numerical examples are shown to illustrate the theoretical results.

preprint2022arXiv

A nonuniform mesh method for wave scattered by periodic surfaces

In this paper, we propose a new nonuniform mesh method to simulate acoustic scattering problems in two dimensional periodic structures with non-periodic incident fields numerically. As existing methods are difficult to extend to higher dimensions, we have designed the new method with such extensions in mind. With the help of the Floquet-Bloch transform, the solution to the original scattering problem is written as an integral of a family of quasi-periodic problems. These are defined in bounded domains for each value of the Floquet parameter which varies in a bounded interval. The key step in our method is the numerical approximation of the integral by a quadrature rule adapted to the regularity of the family of quasi-periodic solutions. We design a nonuniform mesh with a Gaussian quadrature rule applied on each subinterval. We prove that the numerical method converges exponentially with respect to both the number of subintervals and the number of Gaussian quadrature points. Some numerical experiments are provided to illustrate the results.

preprint2022arXiv

A spectral decomposition method to approximate DtN maps in complicated waveguides

In this paper, we propose a new spectral decomposition method to simulate waves propagating in complicated waveguides. For the numerical solutions of waveguide scattering problems, an important task is to approximate the Dirichlet-to-Neumann map efficiently. From previous results, the physical solution can be decomposed into a family of generalized eigenfunctions, thus we can write the Dirichlet-to-Neumann map explicitly by these functions. From the exponential decay of the generalized eigenfunctions, we approximate the Dirichlet-to-Neumann (DtN) map by a finite truncation and the approximation is proved to converge exponentially. With the help of the truncated DtN map, the unbounded domain is truncated into a bounded one, and a variational formulation for the problem is set up in this bounded domain. The truncated problem is then solved by a finite element method. The error estimation is also provided for the numerical algorithm and numerical examples are shown to illustrate the efficiency of the algorithm.

preprint2022arXiv

High order complex contour discretization methods to simulate scattering problems in locally perturbed periodic waveguides

In this paper, two high order complex contour discretization methods are proposed to simulate wave propagation in locally perturbed periodic closed waveguides. As is well known the problem is not always uniquely solvable due to the existence of guided modes. The limiting absorption principle is a standard way to get the unique physical solution. Both methods are based on the Floquet-Bloch transform which transforms the original problem to an equivalent family of cell problems. The first method, which is designed based on a complex contour integral of the inverse Floquet-Bloch transform, is called the CCI method. The second method, which comes from an explicit definition of the radiation condition, is called the decomposition method. Due to the local perturbation, the family of cell problems are coupled with respect to the Floquet parameter and the computational complexity becomes much larger. To this end, high order methods to discretize the complex contours are developed to have better performances. Finally we give the convergence results which we confirm with numerical examples.

preprint2020arXiv

Spectrum decomposition of translation operators in periodic waveguide

Scattering problems in periodic waveguides are interesting but also challenging topics in mathematics, both theoretically and numerically. Due to the existence of eigenvalues, the unique solvability of these problems is not always guaranteed. To obtain a unique solution that is "physically meaningful", the limiting absorption principle (LAP) is a commonly used method. LAP assumes that the limit of a family of solutions with absorbing media converges, as the absorption parameter tends to 0, and the limit is the "physically meaningful solution". It is also called the LAP solution in this paper. It has been proved that the LAP holds for periodic waveguides in [Hoa11]. In this paper, we consider the spectrum decomposition of periodic translation operators. With the curve integral formulation and a generalized Residue theorem, the operator is explicitly described by its eigenvalues and generalized eigenfunctions, which are closely related to Bloch wave solutions. Then the LAP solution is decomposed into generalized eigenfunctions. This gives a better understanding of structure of the scattered fields.

preprint2016arXiv

A Convergent Numerical Scheme for Scattering of Aperiodic Waves from Periodic Surfaces Based on the Floquet-Bloch Transform

Periodic surface structures are nowadays standard building blocks of optical devices. If such structures are illuminated by aperiodic time-harmonic incident waves as, e.g., Gaussian beams, the resulting surface scattering problem must be formulated in an unbounded layer including the periodic surface structure. An obvious recipe to avoid the need to discretize this problem in an unbounded domain is to set up an equivalent system of quasiperiodic scattering problems in a single (bounded) periodicity cell via the Floquet-Bloch transform. The solution to the original surface scattering problem then equals the inverse Floquet-Bloch transform applied to the family of solutions to the quasiperiodic problems, which simply requires to integrate these solutions in the quasiperiodicity parameter. A numerical scheme derived from this representation hence completely avoids the need to tackle differential equations on unbounded domains. In this paper, we provide rigorous convergence analysis and error bounds for such a scheme when applied to a two-dimensional model problem, relying upon a quadrature-based approximation to the inverse Floquet-Bloch transform and finite element approximations to quasiperiodic scattering problems. Our analysis essentially relies upon regularity results for the family of solutions to the quasiperiodic scattering problems in suitable mixed Sobolev spaces. We illustrate our error bounds as well as efficiency of the numerical scheme via several numerical examples.

preprint2016arXiv

A Floquet-Bloch transform based numerical method for scattering from locally perturbed periodic surfaces

Scattering problems for periodic structures have been studied a lot in the past few years. A main idea for numerical solution methods is to reduce such problems to one periodicity cell. In contrast to periodic settings, scattering from locally perturbed periodic surfaces is way more challenging. In this paper, we introduce and analyze a new numerical method to simulate scattering from locally perturbed periodic structures based on the Bloch transform. As this transform is applied only in periodic domains, we firstly rewrite the scattering problem artificially in a periodic domain. With the help of the Bloch transform, we secondly transform this problem into a coupled family of quasiperiodic problems posed in the periodicity cell. A numerical scheme then approximates the family of quasiperiodic solutions (we rely on the finite element method) and backtransformation provides the solution to the original scattering problem. In this paper, we give convergence analysis and error bounds for a Galerkin discretization in the spatial and the quasiperiodicity's unit cells. We also provide a simple and efficient way for implementation that does not require numerical integration in the quasiperiodicity, together with numerical examples for scattering from locally perturbed periodic surfaces computed by this scheme.

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.