有耗介质空间完全匹配层吸收边界条件及其应用

 针对Gedney提出的完全匹配层(PML)无法用于有耗各向同性计算域的截断及其对倏逝波的衰减不理想等问题,提出了一种扩展方法。扩展的PML的主要思想是在各向异性的PML中引入与有耗介质空间相一致的复介电常数和复磁导率,使之可以与有耗介质计算域相匹配。通过给PML的张量介电常数、张量磁导率增加衰减因子以加速倏逝波的衰减。构造了PML吸收效果验证模型,数值结果证明了扩展的PML在处理有耗介质计算域截断问题中的有效性。利用该吸收边界条件,采用时域有限差分法计算了电磁脉冲作用下地面铺设电缆的电磁脉冲响应,计算结果和试验结果取得了较好的一致。     

粒子模拟软件吸收边界的研究

 实现了时域有限差分法的Gedney完全匹配层吸收边界条件，并进行了数值验证和参数优化。从结果可以看出其有较低的反射，性能优良；对吸收媒质厚度和电导率分布阶数的参数优化结果与Gedney经验值基本一致。将此吸收边界模块加入到国产粒子模拟软件CHIPIC中，同美国商业粒子模拟软件MAGIC相比有更好的吸收效果。通过对2维极坐标下磁控管的模拟，证明了此吸收边界具有实际应用价值。     

孔隙介质与液体界面上波动方程数值模拟的改进 方案——MPML边界条件的应用*

由于波场数值模拟的计算区域是有限的,为了压制截断边界造成的人工边界反射,通常采用完全匹配层作为吸收边界条件。但是在含液-固界面的弹性介质中进行正演模拟时,完全匹配层边界条件容易造成稳定性方面的问题。因此,该文将多轴完全匹配层应用于该类介质的正演模拟,以改善完全匹配层边界条件引起的不稳定问题。通过在具有弹性海底的海洋环境模型以及充填液体的井孔模型中进行弹性波方程的正演模拟,分别采用传统的完全匹配层和多轴完全匹配层作为吸收边界条件,对比验证了多轴完全匹配层在含液-固界面弹性介质中数值模拟的有效性和稳定性。     

双相各向同性介质BISQ模型弹性波方程数值模拟方法

本文从Taylor级数展开式出发,推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式;从本构方程和运动方程推导出了BISQ模型双相介质一阶双曲型应力-速度弹性波方程交错网格任意偶数阶精度差分格式以及推导出二维双相各向同性介质完全吸收层边界条件公式和相应的高阶交错网格差分格式。通过数值模拟研究表明,该方法边界吸收效果好,稳定性好,能够高精度模拟双相介质中地震弹性波场,且计算效率也高。     

二阶声场波动方程的非分裂完全匹配层算法

在声场仿真中,完全匹配层(Perfectly Matched Layer,PML)是一种十分有效的吸收边界并得到广泛应用。为了解决基于二阶声场波动方程数值仿真中的吸收边界问题,提出了一种非分裂PML算法。首先,基于伸缩坐标变换,推导了PML算法的频域表达式。然后,通过构造辅助微分方程,得到了非分裂PML的时域表达式。最后,进行了相关理论分析和数值仿真,结果表明:相对于已有的声场分裂PML算法,该算法在保持相同的吸收效率的同时,能较大地节约存储空间,提高计算效率,且更易于实现。      

二维单步交替方向隐式时域有限差分法吸收边界性能分析

给出了两种适用于二维单步交替方向隐式时域有限差分(2-D Leapfrog ADI-FDTD)方法的吸收边界:Mur边界和卷积完全匹配层(CPML)边界。单步交替方向隐式时域有限差分(Leapfrog ADI-FDTD)方法是一种无条件稳定的全隐式差分算法,由于二维空间Leapfrog ADI-FDTD的迭代同时存在显式和隐式方程,故而不同电磁分量的边界条件也存在差异。从原理出发,推导了适用于2-D Leapfrog ADI-FDTD方法的CPML边界条件,并与一阶Mur边界进行比较,利用自由空间的反射误差来表征两种边界的吸收性能,简要总结了两种吸收边界的优缺点。     

An FDTD analysis of photonic crystal waveguides comprising third-order nonlinear materials

A finite difference time domain (FDTD) study of two-dimensional photonic crystals containing nonlinear materials is presented in this paper. An appropriate Z-transform oriented formulation of the FDTD method for the simulation of third-order nonlinear Kerr- and Raman-type media is analyzed and applied to model nonlinear photonic crystal waveguide structures. For their reflectionless termination a novel perfectly matched layer (PML) is proposed and evaluated comparatively to other periodic and inhomogeneous absorbers. Furthermore, the absorbing efficiency of the proposed PML is investigated varying its parameters.     

Absorbing boundary conditions for scalar waves in anisotropic media. Part 1: Time harmonic modeling

With the ultimate goal of devising effective absorbing boundary conditions (ABCs) for general anisotropic media, we investigate the accuracy aspects of local ABCs designed for the scalar anisotropic wave equation in the frequency domain (time harmonic case). The ABC analyzed in this paper is the perfectly matched discrete layers (PMDL). PMDL is a simple variant of perfectly matched layers (PML) and is equivalent to rational approximation-based local ABCs. Specifically, we derive a sufficient condition for PMDL to accurately absorb wave modes with outgoing group velocities and this condition turns out to be a simple bound on the PMDL parameters. The reflection coefficient derived in this paper clearly reveals that the PMDL absorption is based on group velocities, and not phase velocities, and hence a PMDL can be designed to correctly identify and accurately absorb all outgoing wave modes (even those with opposing signs of phase and group velocities). The validity of the sufficient condition is demonstrated through a series of frequency domain simulations. In part 2 of this paper [S. Savadatti, M.N. Guddati, Absorbing boundary conditions for scalar waves in anisotropic media. Part 2: Time-dependent modeling, J. Comput. Phys. (2010), http://dx.doi.org/10.1016/j.jcp.2010.05.017], the accuracy condition presented here is shown to govern both the well-posedness and accuracy aspects of PMDL designed for transient (time-dependent) modeling of scalar waves in anisotropic media.     

完全匹配层吸收边界在2D3V等离子体粒子模拟中的实现

基于Berenger完全匹配层的基本思想,导出了高斯单位制下2维完全匹配层吸收介质的控制方程及参数设置方法,并进行了数值验证和参数优化。依据所得结论编制了相应的吸收边界模块加入到2D3V等离子体粒子模拟程序PLASIM中,取得了良好的吸收效果。基于完全匹配层吸收边界的实现,在小窗口内对超强激光预脉冲与低密度等离子体的长时间相互作用过程进行了模拟研究,结果表明:激光强度和作用时间是影响等离子体密度改变的重要因素,超强激光其长时间的预脉冲可以通过有质动力的累积效应在低密度等离子体中引起较大的密度变化。     

常规分裂和非分裂完全匹配层吸收边界比较研究

分别对常规分裂和非分裂完全匹配层吸收边界在完全弹性介质和孔隙介质的应用效果进行了比较研究。比较了两种吸收边界应用于完全弹性介质和孔隙介质所需计算量、存储量大小以及吸收效果。为比较吸收效果,首先给出两种吸收边界下声场快照图的比较,然后利用吸收边界附近的接收点波形,计算和分析了两种吸收边界在固定入射角度处随时间变化的边界反射幅度,以及对不同入射角度声波的边界反射系数。研究结果表明,无论是完全弹性介质还是孔隙介质,常规非分裂完全匹配层吸收边界比分裂完全匹配层吸收边界对垂直入射乃至一般入射角度的声波都有更好的吸收效果,但需要更大的存储量和计算量。      

色散媒质中的理想匹配层吸收边界及其误差分析

根据色散媒质的性质及电磁波在不同媒质交界面上的反射的机理，证明并分析了在理想匹配层吸收边界中引入色散可实现有效吸收入射到吸收边界上的电磁波，从而达到减少时域有限差分法计算区域和计算误差的目的。     

Frequency domain finite-element and spectral-element acoustic wave modeling using absorbing boundaries and perfectly matched layer

Wave propagation modeling as a vital tool in seismology can be done via several different numerical methods among them are finite-difference, finite-element, and spectral-element methods (FDM, FEM and SEM). Some advanced applications in seismic exploration benefit the frequency domain modeling. Regarding flexibility in complex geological models and dealing with the free surface boundary condition, we studied the frequency domain acoustic wave equation using FEM and SEM. The results demonstrated that the frequency domain FEM and SEM have a good accuracy and numerical efficiency with the second order interpolation polynomials. Furthermore, we developed the second order Clayton and Engquist absorbing boundary condition (CE-ABC2) and compared it with the perfectly matched layer (PML) for the frequency domain FEM and SEM. In spite of PML method, CE-ABC2 does not add any additional computational cost to the modeling except assembling boundary matrices. As a result, considering CE-ABC2 is more efficient than PML for the frequency domain acoustic wave propagation modeling especially when computational cost is high and high-level absorbing performance is unnecessary.     

3D-Scalar Fourier Eigenvector Expansion Method (Fourier-EEM) for analyzing optical waveguide discontinuities

In this paper we propose a Fourier Eigenvector Expansion Method (Fourier-EEM) with coordinate-stretching type PML absorbing boundary conditions. This technique has been used to analyze, under the scalar assumption, 3D optical devices with abrupt discontinuities in the longitudinal direction (i.e. Multimode Interference devices (MMI)). The obtained results clearly confirm the accuracy and efficiency of the proposed approach to calculate the transmission and reflection characteristics of photonic devices.     

The moment method for boundary layer problems in Brownian motion theory

We apply Grad's moment method, with Hermite moments and Marshak-type boundary conditions, to several boundary layer problems for the Klein-Kramers equation, the kinetic equation for noninteracting Brownian particles, and study its convergence properties as the number of moments is increased. The errors in various quantities of physical interest decrease asymptotically as inverse powers of this number; the exponent is roughly three times as large as in an earlier variational method, based on an expansion in the exact boundary layer eigenfunctions. For the case of a fully absorbing wall (the Milne problem) we obtain full agreement with the recent exact solution of Marshall and Watson; the relevant slip coefficient, the Milne length, is reproduced with an accuracy better than 10^–6. We also consider partially absorbing walls, with specular or diffuse reflection of nonabsorbed particles. In the latter case we allow for a temperature difference between the wall and the medium in which the particles move. There is noa priori reason why our method should work only for Brownian dynamics; one may hope to extend it to a broad class of linear transport equations. As a first test, we looked at the Milne problem for the BGK equation. In spite of the completely different analytic structure of the boundary layer eigenfunctions, the agreement with the exact solution is almost as good as for the Klein-Kramers equation.     

Novel optical directional coupler based on surface plasmon polaritons

In this paper, finite difference time domain (FDTD) method and perfect matching layer (PML) absorbing boundary condition are adopted to simulate and analyze a novel optical directional coupler (ODC) based on surface plasmon polaritons (SPPs). Transmittance at each output port of the novel ODC with different coupling region lengths shows it follows the general regulations of a conventional ODC. Especially, its transverse size is of nanoscale. The extreme power position offset between the two output ports is proved to be connected with the real part of Ag's complex refractive index. The excess loss and isolation of the ODC are, respectively, 0.57 and 25.9 dB for 1550 nm telecommunication wavelength, when the length of the coupling region equals half of its coupling length.     

Near-field performance analysis of locally-conformal perfectly matched absorbers via Monte Carlo simulations

In the numerical solution of some boundary value problems by the finite element method (FEM), the unbounded domain must be truncated by an artificial absorbing boundary or layer to have a bounded computational domain. The perfectly matched layer (PML) approach is based on the truncation of the computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence. In this paper, we present the near-field numerical performance analysis of our new PML approach, which we call as locally-conformal PML, using Monte Carlo simulations. The locally-conformal PML method is an easily implementable conformal PML implementation, to the problem of mesh truncation in the FEM. The most distinguished feature of the method is its simplicity and flexibility to design conformal PMLs over challenging geometries, especially those with curvature discontinuities, in a straightforward way without using artificial absorbers. The method is based on a special complex coordinate transformation which is ‘locally-defined’ for each point inside the PML region. The method can be implemented in an existing FEM software by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. We first introduce the analytical derivation of the locally-conformal PML method for the FEM solution of the two-dimensional scalar Helmholtz equation arising in the mathematical modeling of various steady-state (or, time-harmonic) wave phenomena. Then, we carry out its numerical performance analysis by means of some Monte Carlo simulations which consider both the problem of constructing the two-dimensional Green’s function, and some specific cases of electromagnetic scattering.     

非线性欧拉方程组的完美匹配层吸收边界条件

对于非线性Euler方程,提出一类基于完美匹配层(PML)技术的吸收边界条件。首先对线性化的Euler方程设计出PML公式,然后将线性化Euler方程中的通量函数替换成相对应的非线性通量函数,得到非线性的PML方程。考虑到PML方程中包含有一个刚性的源项,文中采用一种隐显Runge-Kutta方法来求解空间半离散后得到的ODE系统。数值实验表明设计的非线性PML吸收边界条件优于传统的特征边界条件。     

An FDTD absorbing boundary condition for anisotropic media based on NPML

Based on the nearly perfectly matched layer (NPML) theory, a finite-difference time-domain (FDTD) absorbing boundary condition (ABC) is presented for truncating three-dimensional (3-D) anisotropic medium. In the proposed technique, the complex coordinate stretching in the NPML scheme and the spatial interpolation method are employed. The associated ABC formulations have the advantage of simplicity in the FDTD implementations. The radiation fields of an electric dipole in anisotropic media are calculated using the presented ABC. The results are numerically verified by the comparison with the reference solutions. Furthermore, in order to clearly show the effective absorbing performance of the proposed method, the reflection coefficient and time-dependent relative error for different layers NPML absorbing boundary are also simulated.     

2.5维圆波导开放边界截断的修正完全匹配层方法

 在时域有限差分算法中,完全匹配层是一种高效的吸收边界条件。通过在标准完全匹配层(PML)方程中添加具有物理意义的电导率和导磁率附加项,提出了修正完全匹配层(MPML)方法,给出了圆柱坐标下3维MPML方程,并应用于2.5维圆波导开放边界的截断计算。计算结果表明,在不同频率和模式的激励下,加入平滑电导率和导磁率附加项后,与标准PML方法相比,MPML方法的相对误差降低8~10 dB。     

A staggered-grid finite-difference method with perfectly matched layers for poroelastic wave equations.

A particle velocity-strain, finite-difference (FD) method with a perfectly matched layer (PML) absorbing boundary condition is developed for the simulation of elastic wave propagation in multidimensional heterogeneous poroelastic media. Instead of the widely used second-order differential equations, a first-order hyperbolic leap-frog system is obtained from Biot's equations. To achieve a high accuracy, the first-order hyperbolic system is discretized on a staggered grid both in time and space. The perfectly matched layer is used at the computational edge to absorb the outgoing waves. The performance of the PML is investigated by calculating the reflection from the boundary. The numerical method is validated by analytical solutions. This FD algorithm is used to study the interaction of elastic waves with a buried land mine. Three cases are simulated for a mine-like object buried in "sand," in purely dry "sand" and in "mud." The results show that the wave responses are significantly different in these cases. The target can be detected by using acoustic measurements after processing.