基于Z变换的高阶完全匹配层实现

基于拉伸坐标完全匹配层(SC-PML)公式和Z变换方法,提出以非分裂场形式来实现具有多极点的高阶完全匹配层的高效算法,来截断时域有限差分(FDTD)网格.在吸收性能方面,高阶PML同时具有传统PML和复频率偏移完全匹配层(CFS-PML)二者的优点.提供的数值算例是二维TE极化电磁波与无限长且有限宽度的理想电导体(PEC)薄片的相互作用.仿真结果显示,高阶PML公式在衰减低频行波和隐失波及减少后期反射方面效果好,比传统SC-PML和复频率偏移的卷积完全匹配层(CPML)算法有更好的吸收性能.     

卷积完全匹配层截断3维金属矩形波导的应用研究

 讨论了高功率微波源模拟中波导开放边界截断的需求，分析了不同类型完全匹配层(PML)的特点，选用卷积形式PML截断色散波导器件的开放边界。从自由空间电磁波的平面波解和分裂形式的PML出发，构造了未分裂形式的PML，用傅里叶变换的卷积定理，导出了直角坐标系下卷积完全匹配层(CPML)介质中电磁场的迭代形式的离散方程。在不同频率和模式激励源作用下，模拟计算了CPML截断矩形波导开放边界的性能，数值结果表明最大相对误差都小于-70 dB，远好于Mur吸收边界的截断效果。     

基于辅助差分方程的完全匹配层在时域多分辨率分析算法中的应用与性能分析

将基于辅助微分方程的完全匹配层(ADE-PML)吸收边界条件引入到基于Daubechies尺度函数的时域多分辨率分析算法中. 与目前广泛应用的Berenger完全匹配层(PML)和各向异性介质完全匹配层(APML) 相比, 该吸收边界条件的实现更加容易且更节省内存. 数值结果表明, ADE-PML在吸收传播模和低频凋落模方面均优于PML和APML.     

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

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

卷积完全匹配层在两维声波有限元计算中的应用

将基于复坐标变换和复频移扩展坐标变量的卷积完全匹配层Convolution Perfectly Machted Layer(CPML)引入到两维声波方程的有限元(FEM)计算中,该匹配层作为一种吸收边界条件Absorbing Boundary Condition(ABC)应用在有限元计算的边界截断上。文中分别给出了频域和时域的CPML方程的表达形式,并在有限元计算软件COMSOL中完成数值计算。相对于经典的PML,CPML最大的优势在于它不需要把场分裂开,这使其具有更好的稳定性和更高的吸收性能,且更易于实现。数值计算结果表明,CPML边界层有着比PML更好的吸收效果,其更有效的吸收了进入其中的声场能量。      

CPML在空中核爆电磁脉冲三维数值模拟中的应用

在分析不同类型完全匹配层(PML,perfectl ymatched layer)特点的基础上,选用了卷积形式PML(CPML,convolutional perfectly matched layer)截断空中核爆电磁脉冲(NEMP,nuclear electromagnetic pulse)数值模拟的开放边界.从自由空间中电磁波的平面波解和分裂形式的PML出发,构造了未分裂形式的PML,应用付里叶变换的卷积定理,推导出三维旋转椭球-双曲正交坐标系下CPML介质中电磁场的迭代形式的离散方程.获得了很好的截断效果,展示了CPML在NEMP数值模拟中的应用前景.     

卷积完全匹配层在空爆电磁脉冲2维模拟中的应用

 在分析不同类型完全匹配层特点的基础上，选用了卷积形式完全匹配层（CPML）截断空中核爆电磁脉冲数值模拟的开放边界。从自由空间中电磁波的平面波解和分裂形式的完全匹配层出发，构造了未分裂形式的完全匹配层，应用傅里叶变换的卷积定理，推导出2维旋转椭球-双曲正交坐标系下卷积形式完全匹配层介质中电磁场的迭代形式的离散方程。计算表明，采用CPML吸收边界方法使得截断的外边界处的场的计算误差大大减少。     

一种计算和分析二维光子晶体缺陷模式的方法

通过改进时域有限差分(FDTD)法,计算和分析了二维光子晶体的缺陷模式。运用一维时域有限差分算法和线性插值法在总场散射场(TF-SF)连接边界引入入射平面波,采用完全匹配层(PML)技术对外行波进行了有效吸收。计算和分析结果表明,在光子晶体非对称方向入射的平面波能激发所有的缺陷模式,选取合适的探测点位置收集电场值,经快速傅里叶变换(FFT)能得到所有的共振峰值。另外,采用该方法研究了二维正方介质柱光子晶体缺陷模的共振频率与缺陷介质柱半径和介电常量之间的关系。结果表明通过改变缺陷的半径和介电常量大小可以在光子晶体禁带中一定的范围内调节缺陷模式的共振频率大小。     

基于广角FD-BPM的PML边界处理方法

本文提出了基于三阶Padé近似广角有限差分光束传输法(FD-BPM)的简单有效的完美匹配层(PML)边界处理方法.分别给出了该方法应用于斜波导和多模干涉(MMI)数值模拟的结果,并与用完全透明边界条件处理的模拟计算结果进行了比较.本文最后对基于广角FD-BPM的PML边界处理方法进行了优化.     

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

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

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

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

Absorbing boundary conditions for nonlinear Euler and Navier–Stokes equations based on the perfectly matched layer technique

Absorbing boundary conditions for the nonlinear Euler and Navier–Stokes equations in three space dimensions are presented based on the perfectly matched layer (PML) technique. The derivation of equations follows a three-step method recently developed for the PML of linearized Euler equations. To increase the efficiency of the PML, a pseudo mean flow is introduced in the formulation of absorption equations. The proposed PML equations will absorb exponentially the difference between the nonlinear fluctuation and the prescribed pseudo mean flow. With the nonlinearity in flux vectors, the proposed nonlinear absorbing equations are not formally perfectly matched to the governing equations as their linear counter-parts are. However, numerical examples show satisfactory results. Furthermore, the nonlinear PML reduces automatically to the linear PML upon linearization about the pseudo mean flow. The validity and efficiency of proposed equations as absorbing boundary conditions for nonlinear Euler and Navier–Stokes equations are demonstrated by numerical examples.     

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

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

一种新时域交替隐式差分算法在散射问题中的应用

提出一种新时域交替隐式有限差分(ADI-FDTD)算法格式. 传统ADI-FDTD算法的 场量步进方程涉及周围若干网格的较多场量，导致两个区域的步进方程处理较困难：一个是 邻近完全匹配层(PML)和散射场交界区，另一个是邻近连接边界区. 特别是后者，考虑入射 波影响需对场量所在区域判断，根据不同情况对原有方程进行修正，一维和二维散射问题 相对简单，可三维问题修正式有数十种之多而几乎无法完成. 本方法基于分裂场形式的ADI- FDTD技术，使得散射场区和PML吸收层区的表达形式完全一致，从而忽略两者差别.另 关键词： 时域交替隐式有限差分算法 电磁散射     

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.     

A novel Fourier based 3D full-vectorial beam propagation method

In this paper a Fourier based full-vectorial beam propagation method (BPM) with perfectly matched layer (PML) absorbing boundary conditions is presented. The method starts with the slowly varying envelope coupled wave equations, including coordinate-stretching type PMLs, which are discretised using the Galerkin method with Fourier basis functions in the transversal direction. The proposed method is then assessed in a wide variety of situations, including 2D propagation of TE and TM waves and 3D full-vectorial situations. The obtained results show that good accuracy can be obtained even for complex strongly-guiding waveguides with a small number of basis functions.     

The Analysis of Post-Mount in W-Band Circular Groove Guide

Finite-difference time-domain (FDTD) method is used to calculate the driving-point impedance of the circular groove guide diode mount. A simple way is proposed to approximate the curvilinear metal boundary of circular groove guide. Discrete equations are obtained from the Ampere's law and Faraday's law. Perfectly matched layer (PML) is applied on the open boundary of circular groove guide. At last program is made to obtain the field components. The driving-point impedance can be achieved from the ratio of the Fourier transform of the voltage and the Fourier transform of the current. The cutoff frequency and field pattern of the dominant mode are also obtained with well agreement with the theoretical results, this verifies the validity of the FDTD program.     

A novel unsplit perfectly matched layer for the second-order acoustic wave equation

When solving acoustic field equations by using numerical approximation technique, absorbing boundary conditions (ABCs) are widely used to truncate the simulation to a finite space. The perfectly matched layer (PML) technique has exhibited excellent absorbing efficiency as an ABC for the acoustic wave equation formulated as a first-order system. However, as the PML was originally designed for the first-order equation system, it cannot be applied to the second-order equation system directly. In this article, we aim to extend the unsplit PML to the second-order equation system. We developed an efficient unsplit implementation of PML for the second-order acoustic wave equation based on an auxiliary-differential-equation (ADE) scheme. The proposed method can benefit to the use of PML in simulations based on second-order equations. Compared with the existing PMLs, it has simpler implementation and requires less extra storage. Numerical results from finite-difference time-domain models are provided to illustrate the validity of the approach.     

A perfectly matched layer approach to the nonlinear Schrödinger wave equations

Absorbing boundary conditions (ABCs) are generally required for simulating waves in unbounded domains. As one of those approaches for designing ABCs, perfectly matched layer (PML) has achieved great success for both linear and nonlinear wave equations. In this paper we apply PML to the nonlinear Schrödinger wave equations. The idea involved is stimulated by the good performance of PML for the linear Schrödinger equation with constant potentials, together with the time-transverse invariant property held by the nonlinear Schrödinger wave equations. Numerical tests demonstrate the effectiveness of our PML approach for both nonlinear Schrödinger equations and some Schrödinger-coupled systems in each spatial dimension.