A multigrid-based preconditioned Krylov subspace method for the Helmholtz equation with PML

In this paper, we generalize the complex shifted Laplacian preconditioner to the complex shifted Laplacian-PML preconditioner for the Helmholtz equation with perfectly matched layer (Helmholtz-PML equation). The Helmholtz-PML equation is discretized by an optimal 9-point difference scheme, and the preconditioned linear system is solved by the Krylov subspace method, especially by the biconjugate gradient stabilized method (Bi-CGSTAB). The spectral analysis of the linear system is given, and a new matrix-based interpolation operator is proposed for the multigrid method, which is used to approximately invert the preconditioner. The numerical experiments are presented to illustrate the efficiency of the preconditioned Bi-CGSTAB method with the multigrid based on the new interpolation operator, also, numerical results are given for comparing the performance of the new interpolation operator with that of classic bilinear interpolation operator and the one suggested in Erlangga et al. (2006) [10].     

一阶速度-应力Biot双相各向同性介质弹性波

提出一种等价的一阶双曲型速度一应力Biot双相各向同性介质弹性波波动方程,以实现双相介质混合波场中纯快慢纵波和纯横波波场分离的问题.应用散度和旋度理论证明双相介质等价方程波场分离的可行性,采用高阶交错网格有限差分法构建高精度正演算子,推导其PML吸收边界条件和稳定性条件,并对均匀双相介质和层状非均匀双相介质模型进行数值...     

粗糙度对大高宽比槽道内流动影响的数值模拟

基于确定表面粗糙结构形状的PML模型模拟表面粗糙度对大高宽比小槽道内流动特性的影响,并与3-D模拟结果对比,结果较为一致.基于该模型对不同的流速和不同的粗糙度进行模拟,结果表明:表面粗糙度导致槽道内出现速度线性分布的流动底层.Re数相同时,单位长度压降与相对粗糙度成二次方关系.粗糙单元高度相同时,压降随Re线性增加.     

A Fourier spectral-discontinuous Galerkin method for time-dependent 3-D Schrödinger–Poisson equations with discontinuous potentials

In this paper, we propose a high order Fourier spectral-discontinuous Galerkin method for time-dependent Schrödinger–Poisson equations in 3-D spaces. The Fourier spectral Galerkin method is used for the two periodic transverse directions and a high order discontinuous Galerkin method for the longitudinal propagation direction. Such a combination results in a diagonal form for the differential operators along the transverse directions and a flexible method to handle the discontinuous potentials present in quantum heterojunction and supperlattice structures. As the derivative matrices are required for various time integration schemes such as the exponential time differencing and Crank Nicholson methods, explicit derivative matrices of the discontinuous Galerkin method of various orders are derived. Numerical results, using the proposed method with various time integration schemes, are provided to validate the method.     

使用PML技术求解Helmholtz方程的几种区域分解算法

1引言 在科学与工程计算中,时谐声波和电磁波散射现象的模拟常归结为Helmholtz方程的数值求解,这是一个重要而具有挑战性的问题.     

A PML Method for Electromagnetic Scattering from Two-dimensional Overfilled Cavities

In this paper, we consider electromagnetic scattering problems for two-dimensional overfilled cavities. A half ringy absorbing perfectly matched layer （PML） is introduced to enclose the cavity, and the PML formulations for both TM and TE polarizations are presented. Existence, uniqueness and convergence of the PML solutions are considered. Numerical experiments demonstrate that the PML method is efficient and accurate for solving cavity scattering problems.     

Efficient Analysis of Millimeter Wave Ferrite Circulators by the GMRES Iterative Algorithm

The finite element method (FEM) combined with the perfectly matched layers (PML) is given for simulation of waveguide ferrite circulators. The generalized minimal residual (GMRES) iterative method is applied to solve such sparse large non-symmetric system of linear equations resulting from the use of edge-based finite element method. The formulation of FEM and the algorithm of GMRES method are described in detail. The reflection and insertion losses of millimeter wave waveguide circulator are analyzed and the results are compared with those obtained from literature.     

时域有限差分法在光波导分析中的应用及改进

对时域有限差分算法在圆柱坐标下弱导光纤中应用时吸收边界条件的设置及其对计算精度的影响进行了讨论，提出了改善的方法。并且将这种方法应用于一般光纤和光电子晶体光纤传输特性的计算，得到了有用结果。     

基于PML边界的有限差分法及其在光波导泄漏损耗计算中的应用

利用结合完美匹配层(PML)边界的有限差分法计算了光波导的泄漏损耗.通过采用非均匀格点差分格式和反正切坐标变换的方法,有效地减小了计算量并提高了计算精度.分析了PI(polyimide)掩埋型波导和SOI(silicon-on-insulator)脊型波导两种典型结构的泄漏损耗,给出了波导结构尺寸对泄漏损耗的影响,并对有效减小泄漏损耗的方法进行了讨论.     

Analysis of a finite PML approximation for the three dimensional time-harmonic Maxwell and acoustic scattering problems

We consider the approximation of the frequency domain three-dimensional Maxwell scattering problem using a truncated domain perfectly matched layer (PML). We also treat the time-harmonic PML approximation to the acoustic scattering problem. Following work of Lassas and Somersalo in 1998, a transitional layer based on spherical geometry is defined, which results in a constant coefficient problem outside the transition. A truncated (computational) domain is then defined, which covers the transition region. The truncated domain need only have a minimally smooth outer boundary (e.g., Lipschitz continuous). We consider the truncated PML problem which results when a perfectly conducting boundary condition is imposed on the outer boundary of the truncated domain. The existence and uniqueness of solutions to the truncated PML problem will be shown provided that the truncated domain is sufficiently large, e.g., contains a sphere of radius . We also show exponential (in the parameter ) convergence of the truncated PML solution to the solution of the original scattering problem inside the transition layer.

Our results are important in that they are the first to show that the truncated PML problem can be posed on a domain with nonsmooth outer boundary. This allows the use of approximation based on polygonal meshes. In addition, even though the transition coefficients depend on spherical geometry, they can be made arbitrarily smooth and hence the resulting problems are amenable to numerical quadrature. Approximation schemes based on our analysis are the focus of future research.