STABILITY ANALYSIS OF YEE SCHEMES TO PML AND UPML FOR MAXWELL EQUATIONS

We utilize Fourier methods to analyze the stability of the Yee difference schemes for Berenger PML （perfectly matched layer） as well as the UPML （uniaxial perfectly matched layer） systems of two-dimensional Maxwell equations. Using a practical spectrum stability concept, we find that the two schemes are spectrum stable under the same conditions for mesh sizes. Besides, we prove that the UPML schemes with the same damping in both directions are stable. Numerical examples are given to confirm the stability analysis for the PML method.     

A Source Transfer Domain Decomposition Method for Helmholtz Equations in Unbounded Domain Part II: Extensions

In this paper we extend the source transfer domain decomposition method (STDDM) introduced by the authors to solve the Helmholtz problems in two-layered media, the Helmholtz scattering problems with bounded scatterer, and Helmholtz problems in 3D unbounded domains. The STDDM is based on the decomposition of the domain into non-overlapping layers and the idea of source transfer which transfers the sources equivalently layer by layer so that the solution in the final layer can be solved using a PML method defined locally outside the last two layers. The details of STDDM is given for each extension. Numerical results are presented to demonstrate the efficiency of STDDM as a preconditioner for solving the discretization problem of the Helmholtz problems considered in the paper.     

双重弹性波波场分离数值模拟及相关理论证明

提出一种等效的双重弹性波波场分离数值模拟方法,用于模拟纯纵波和纯横波分离模式的质点振动速度、位移以及散度场和旋度场,并将该方法应用于全弹性波波动方程数值模拟中.同时,详细推导双重弹性波波场分离波动方程的高阶交错网格有限差分数值计算公式及其稳定性条件、数值频散关系和完全匹配层(PML)吸收边界条件.理论分析和数值计算均表明,该方法可以实现高精度双重弹性波波场分离数值模拟,且纯纵波和纯横波得到完全分离,边界吸收效果较好.与前人工作相比,存储量和计算时间均得到有效改善,数值计算结果进一步验证了该方法的优越性.     

FDTD Analysis of EW Wave Circulating by a Magnetized Ferrite Body in Free Space

An efficient numerical method has been devised for the study of wave circulating by a magnetised ferrite body of an arbitrary shape. It is a finite-difference time-domain formulation that incorporates Mur's absorbing boundary conditions and a perfectly matched layer. Several shapes of interest have been studied, including spheres, circular cylinders. The electromagnetic fields inside ferrite and the power-density distribution on the ferrite's surface normal to the bias external magnetic field are obtained. It is observed that an incident plane wave would circulate around the magnetised ferrite body in an open space as if the ferrite were placed inside a waveguide junction circulator.