Unconditionally Stable Complex Envelope Wave Equation PML Algorithm for Band Limited FDTD Simulations

Unconditionally stable complex envelope (CE) perfectly matched layer (PML) absorbing boundary conditions (ABCs) are presented for truncating the scalar wave-equation finite difference time domain (WE-FDTD) grids. The formulations are based on incorporating the alternating direction implicit (ADI) scheme into the CE FDTD implementations of the scalar wave-equation derived in the PML region at the domain boundaries. Numerical example carried out in two dimensional domain shows that the proposed formulations are more accurate than the classical ADI scalar wave equation PML formulations when it is used for modelling band limited electromagnetic applications.     

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

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

基于非分裂FDTD的左手介质电磁特性的研究

为验证左手介质的电磁特性,采用非分裂时域有限差分方法对左手介质的Drude模型进行建模.这种方法不需要对电场和磁场进行分裂,也不需要对PML空间进行特殊处理,吸收边界PML和工作空间可以通过参数转换来完成,并且构造的PML层为有耗介质,进入PML层的透射波将迅速衰减.它是一种准确而有效的分析色散和各向异性介质的方法,通过此方法有效地验证了左手介质的负折射效应、汇聚效应、相位补偿效应,充分验证了左手介质及其反常特性的存在性和此方法分析左手介质的有效性. 关键词： 左手介质 非分裂 时域有限差分 Drude     

Simple and Efficient PML Algorithm for Truncating Nonlinear FDTD Domains Based on the Z -Transform Method

Based on the Z-transform method, a simple, efficient and unsplit-field implementation of the Stretched Coordinate Perfectly Matched Layer (SC-PML) is proposed for truncating the nonlinear dispersive Finite Difference Time Domain (FDTD) lattices. In addition, the nonlinear FDTD formulations using the Z-transform method are reformulated with the advantage of a simple derivative process. The proposed algorithm is validated through two numerical examples carried out in one dimensional and two dimensional domains which include Lorentz dispersion as well as Kerr and Raman nonlinearities.     

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.     

Complex Envelope Crank-Nicolson Nearly PML Algorithm for the Left-handed Material FDTD Simulations

Unconditionally stable complex envelope (CE) absorbing boundary conditions (ABCs) are presented for truncating left handed material (LHM) domains. The proposed algorithm is based on incorporating the Crank Nicolson (CN) scheme into the CE finite difference time domain (FDTD) implementations of the nearly perfectly matched layer (NPML) formulations. The validity of the formulations is shown through numerical example carried out in one dimensional Lorentzian type LHM FDTD domain.     

Unconditionally Stable Z-Transform ADI-PML Algorithm for Lorentzian Double Negative Meta-Materials

Unconditionally stable formulations of the anisotropic perfectly matched layer (APML) are presented for truncating double negative (DNG) meta-material finite difference time domain (FDTD) grids. In the proposed formulations, the Z-transform theory is employed in the alternating direction implicit FDTD (ADI-FDTD) scheme to obtain update equations for the field components in the DNG meta-material domains. Numerical examples carried out in one dimensional Lorentzian type DNG meta-material domains are included to show the validity of the proposed formulations.     

PML-FDTD法在分析负折射率材料中的应用

将时域有限差分(FDTD)法引入了对负折射率材料物理现象的仿真研究.给出了二维TM波在负折射率材料中的时域差分方程,并且在吸收边界处使用了理想匹配层(PML).为了避免在迭代过程中出现的不稳定现象,在差分方程的推导中引入了Drude模型, 并对Pendry提出的由负折射率材料构成的平板透镜具有的完美成像现象进行了数值仿真验证.由仿真结果发现,完美成像现象只在平板透镜的折射率n=－1时出现,当n≠－1时则会出现近轴聚焦效应.     

SAR EVALUATION IN THE 0.25 MM HUMAN EYE MODEL EXPOSED TO 30 GHz MILLIMETER WAVE

We have evaluated the specific absorption rate (SAR) in the eye induced by the perpendicular incidence of 30 GHz uniform plane millimeter wave (MMW). An anatomically based 0.25 mm voxel model of the eye with 10 types of tissue is constructed. Debye’s dispersion equation is used to derive the dielectric properties of the eye tissues and biological materials in the neighboring region. The three-dimensional (3-D) finite-difference time-domain (FDTD) algorithm and the perfectly matched layer (PML) absorbing boundary conditions are used in the numerical simulation. The FDTD code is validated against the Mie theory. The SAR deposition in the eye is calculated and displayed, with the analysis on the SAR distribution characteristics. The average SAR values in the cornea, the lens cortex, the lens nucleus, the humor, and the whole eye are statistically derived and examined against present irradiation safety standard.     

时域有限差分法分析蚀刻衍射光栅的后向衍射效率

采用时域有限差分法(FDTD)结合理想匹配层(PML)边界条件及周期性边界条件对两种结构槽面的后向衍射效率进行了分析,在P偏振和S偏振两种情况下,均得到了比传统差分方法和射线近似方法更精确的解.数值计算结果表明:全内反射(TIR)结构在高折射率材料(如InP)的蚀刻衍射光栅(EDG)中很有效,而在低折射率材料(如SiO₂)的EDG中效果不理想.     

旋转体时域有限差分法对轴对称亚波长衍射光学元件的分析

采用旋转体时域有限差分法对轴对称亚波长衍射光学元件进行严格的矢量分析.推导了旋转体时域有限差分法的基本计算公式;给出了入射波的设置方法;采用了完全匹配层吸收边界条件;改进了平面波谱传播算法,大大简化了计算过程并提高了计算速度.对多台阶微透镜和二元亚波长微透镜进行分析给出了它们焦平面的电场强度分布.数值计算结果表明,本文的算法可以准确且高效地分析轴对称亚波长衍射光学元件.     

Mathematical simulation of nonlinear waveguiding systems based on 2D photonic crystals

A mathematical model of nonlinear two-dimensional waveguiding systems based on photonic crystals was constructed using the finite-difference time-domain method (FDTD), the total field/scattered field (TF/SF) method, and the perfectly matched layer method (PML). The constructed numerical algorithm can also be used for investigation of other two-dimensional waveguiding systems. The results of simulating the elements of particular waveguiding systems are presented.     

Three-dimensional FDTD method for optical pulse propagation analysis in microstructured optical fibers

We propose a three-dimensional (3D) finite-difference time-domain (FDTD) method to analyze the pulse propagation characteristics in microstructured optical fibers (MOFs). The computation domain size is greatly reduced by adopting the technique of moving problem space. The propagating pulse is virtually held in the buffer cell of the problem space as simulation continues. This method is capable to investigate the temporal evolution of the propagating pulse. Spectral information can be obtained by Fourier analysis. As an example, the influence of the kerr nonlinearity on the optical pulse propagation in a Lorentz dispersive MOF is demonstrated. The model is also used to simulate the nonlinear interactions between the pump spectral broadening and third harmonic generations in a highly nonlinear fused silica nanowire with good agreement with the generalized nonlinear envelop equation (GNEE) model.     

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

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

Distinguishing correct from incorrect PML proposals and a corrected unsplit PML for anisotropic,dispersive media

We show that some previous proposals for perfectly matched layer (PML) absorbers in anisotropic media or for waveguides at oblique incidence are not, in fact true PMLs; in previous work we similarly showed a failure of several PML proposals for periodic media (photonic crystals). We therefore argue that a more careful validation scheme is required for PML proposals, in contrast to past authors who have typically checked only that reflections are small for a fixed resolution, and suggest a simple validation scheme that can be readily applied to any PML proposal regardless of derivation or implementation. We demonstrate this test for a corrected, unsplit-field PML valid for anisotropic, dispersive media, implemented in both planewave-expansion and finite-difference time-domain (FDTD) methods.     

磁化铁氧体覆盖二维导体目标FDTD分析的移位算子方法

采用移位算子方法把时域有限差分法推广应用于二维磁各向异性色散介质—磁化铁氧体中.证明了电磁波横向入射二维轴向磁化铁氧体目标情形下,电磁波可按目标的轴向分解为横电波(TE波)和横磁波(TM波),且TE波的散射特性与铁氧体介质无关,而TM波的散射特性与介质电磁参量密切相关,同时对其物理原因进行了分析.通过采用移位算子方法处理磁化铁氧体频域本构关系,得到该情形下目标电磁散射的移位算子时域有限差的迭代计算公式,同时解决了电磁波在各向异性和频率色散介质中传播的问题.计算了轴向磁化铁氧体涂敷VonKarman型导体柱的TM波双站雷达散射截面,分析了铁氧体参量对目标双站雷达散射截面的影响.结果表明:恰当地选择铁氧体参量能有效地减少目标的雷达散射截面,本文时谐因子取exp〔jωt〕.     

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

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

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.     

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.     

Analysis of Millimeter-Wave Circuits with Different Dielectric Substrate Using the FDTD Technique

In this paper, The Finite Difference Time Domain Method (FDTD) technique is applied to analysis of millimeter-wave circuit which used the different dielectric substrate. Berenger's perfectly matched layers (PML) are employed as absorbing boundary conditions. And we use a Gaussian pulse to excite the power divider. We get the frequency domain information by using the Fourier transform. Experimental results in the 50GHz band show in the paper.