High-order FDTD methods for transverse electromagnetic systems in dispersive inhomogeneous media

This Letter introduces a novel finite-difference time-domain (FDTD) formulation for solving transverse electromagnetic systems in dispersive media. Based on the auxiliary differential equation approach, the Debye dispersion model is coupled with Maxwell's equations to derive a supplementary ordinary differential equation for describing the regularity changes in electromagnetic fields at the dispersive interface. The resulting time-dependent jump conditions are rigorously enforced in the FDTD discretization by means of the matched interface and boundary scheme. High-order convergences are numerically achieved for the first time in the literature in the FDTD simulations of dispersive inhomogeneous media.     

基于移位算子时域有限差分的色散薄层节点修正算法

提出了一种时域有限差分（FDTD）计算中色散介质薄层问题处理的新算法.对于厚度小于一个元胞尺度的电小尺寸色散介质薄层问题,采用将元胞内电位移矢量和磁感应强度加权平均的方法,求得薄层所在元胞内修正点处的等效介质参数.然后根据常见色散介质模型,包括Debye模型、Lorenz模型、Drude模型等,介电常数和磁导率可以表示为jω分式多项式的特点,结合频域到时域的转换关系（即用/t代替jω）和移位算子方法得到了修正点处的时域本构关系,进而获得时域递推计算式.数值结果表明,该方法具有通用、节省计算时间、节省内存和计算精度良好等优点. 关键词： 色散介质薄层 节点修正 移位算子 时域有限差分     

Transient acoustic wave propagation in rigid porous media: a time-domain approach

Wave propagation of acoustic waves in porous media is considered. The medium is assumed to have a rigid frame, so that the propagation takes place in the air which fills the material. The Euler equation and the constitutive relation are generalized to take into account the dispersive nature of these media. It is shown that the connection between the fractional calculus and the behavior of materials with memory allows time-domain wave equations, the coefficients of which are no longer frequency dependent, to be worked out. These equations are suited for direct and inverse scattering problems, and lead to the complete determination of the porous medium parameters.     

频率相关材料的计算模型分析

系统地介绍了Debye、Drude、Lorentz模型相对介电常量的推导过程,并给出了其相应结果;再用有损耗的Drude模型讨论了当前引起广泛关注的左手介质的性质;最后通过推导和举例,进一步揭示了各种模型在实际研究中的重要作用.     

Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions

In this paper, we will present a unified formulation of discontinuous Galerkin method (DGM) for Maxwell's equations in linear dispersive and lossy materials of Debye type and in the artificial perfectly matched layer (PML) regions. An auxiliary differential equation (ADE) method is used to handle the frequency-dependent constitutive relations with the help of auxiliary polarization currents in the computational and PML regions. The numerical flux for the dispersive lossy Maxwell's equations with the auxiliary polarization current variables is derived. Various numerical results are provided to validate the proposed formulation.     

On the stability of nonlinear waves in integrable models

A new method of stability investigation is presented for solutions of nonlinear equations integrable with the help of the inverse scattering transform (IST). The stability problem for periodic nonlinear waves in weakly dispersive media is solved with respect to transverse perturbations. It is shown that for positive dispersion media one-dimensional waves are unstable, and for negative dispersion such waves are stable.     

高阶色散介质的改进移位算子时域有限差分方法

提出一种适用于高阶Debye,Drude,Lorentz及其混合模型的改进移位算子的时域有限差分(SO-FDTD)方法。从介质极化率函数出发,将其写成一阶或二阶有理分式求和的形式,并在随时间步推进计算的过程中,通过引入中间变量和设置临时变量,克服了常规SO-FDTD将高阶模型直接转化为有理分式所导致的计算复杂性和内存占用量大的问题。同时,改进SO-FDTD方法的时域推进计算步骤具有通用性,克服了常规递归卷积(RC-FDTD)方法对各种高阶模型具有不同计算公式,因而不能形成通用计算程序的问题。最后,通过空气-高阶色散介质界面的反射系数计算验证了算法的有效性和通用性。     

含色散与非线性媒质的光子晶体的简正耦合模

利用传输矩阵方法及对色散媒质采用洛伦兹振子模型，研究了一维含色散媒质的光子晶体的简正耦合模，并计算了该结构的复有效折射率。在光学双稳态中，与仅含非线性媒质的光子晶体相比，得到了低得多的阈值入射光强。     

Effect of the intrinsic nonlinearity on the propagation of ultrashort optical pulses in carbon nanotubes in dispersive nonmagnetic dielectric media

The Maxwell equations for the electromagnetic field that propagates in carbon nanotubes (CNTs) placed in dispersive nonmagnetic dielectric media are analyzed with allowance for the intrinsic nonlinearity of the medium. The dependences of the pulse on the initial pulse amplitude, dispersion constants, and nonlinearity are revealed.     

一种处理Cole-Cole色散媒质的FDTD改进方案

对柯尔-柯尔(Cole-Cole)色散媒质的时域有限差分(FDTD)方案尝试给出改进:(1)保留了Cole-Cole经验公式中的静态电导率项;(2)扩展到三维情形;(3)人体乳房组织的模型参数源于实验测量.数值算例结果初步证实改进方案的可行性和有效性.改进方案适用于模拟更一般柯尔-柯尔色散媒质中的电波传播.     

Propagation of a self-induced transparency pulse in a spatially dispersive medium

Propagation of a self-induced transparency pulse in spatially dispersive media is analyzed. A generalized model of two-level systems allowing for excitonic energy transfer is proposed. Periodic and soliton solutions to the governing equations are found. Estimates show that the spatial dispersion of a pulse propagating in a slow-light medium can be substantially enhanced and become important in the case of resonant transition and sufficiently long pulse duration.     

On the modeling of effective constitutive parameters of bianisotropic media by a periodic unfolding method in time and frequency domains

In this paper, we compute the equivalent electromagnetic parameters of a composite periodic material by using the unfolding homogenization method. We are therefore able to solve the Maxwell equations in the unit cell instead of the whole periodic structure. We choose bianisotropic 3D media as they constitute the more general case and we consider dispersive materials to introduce the frequency dependency. The simulations are based on various scenarios studied in the literature. We compare our results with several analytical formalisms and prove the effectiveness of our numerical method in the particular case of the time domain homogenization.     

Interior penalty DG methods for Maxwell’s equations in dispersive media

In this paper, we develop a fully-discrete interior penalty discontinuous Galerkin method for solving the time-dependent Maxwell’s equations in dispersive media. The model is described by a vector integral–differential equation. Our scheme is proved to be unconditionally stable and achieve optimal error estimates in both L² norm and energy norm. The scheme is implemented and numerical results supporting our analysis are presented.     

Posteriori Error Estimation for an Interior Penalty Discontinuous Galerkin Method for Maxwell's Equations in Cold Plasma

In this paper, we develop a residual-based a posteriori error estimator for the time-dependent Maxwell's equations in the cold plasma. Here we consider a semi-discrete interior penalty discontinuous Galerkin (DG) method for solving the governing equations. We provide both the upper bound and lower bound analysis for the error estimator. This is the first posteriori error analysis carried out for the Maxwell's equations in dispersive media.     

An efficient multi-domain spectral algorithm for the simulation of dispersive metallic structures in the time domain

In this article, we present a multi-domain formulation of the spectral time domain algorithm for the simulation of dispersive materials. We propose a leap-frog time-stepping scheme similar to the finite-difference time domain method in order to minimize memory usage. Dispersive material behavior is modelled in the frequency domain and used in our time-domain algorithm by introducing auxiliary differential equations for the macroscopic polarization. Absorbing boundary conditions are presented that can be used with dispersive materials. The numerical investigation of structures with material parameters fitted to experimental data shows excellent agreement with theoretical calculations.     

一种处理色散介质问题的通用时域有限差分方法

色散介质的介电系数是频率的函数,使本构关系在时域成为卷积关系.这就给用时域有限差分方法计算色散介质中波的散射和传播带来了困难.现有算法往往要针对不同色散介质模型推导相应的递推公式,算法的通用性较差.本文完善和发展了移位算子-时域有限差分方法,使之成为一种处理色散介质电磁问题的通用方法.首先,证明了常见的三种色散介质模型(德拜模型、洛伦兹模型和德鲁模型)的介电系数均可以写成适于移位算子法计算的有理分式函数形式.然后,用/t代替jω,过渡到时域,再引入时域移位算子z_t代替时间微分算子来处理有理分式函数形式的介电系数,给出离散时域本构关系的表示式,进而导出时域有限差分方法当中电位移矢量和电场强度之间的关系.最后,计算了几种色散介质的电磁散射,数值结果表明了本文方法和程序的通用性和正确有效性. 关键词： 时域有限差分方法 色散介质 移位算子     

基于半解析递归卷积的通用色散介质FDTD方法

基于数字信号处理中的半解析递归卷积（SARC）算法,提出了一种用于色散介质电磁特性分析的半解析递归卷积时域有限差分方法（SARC FDTD）. 该方法既保持了FDTD方法处理复杂目标的灵活性,又吸取了线性系统SARC算法的绝对稳定性和高精度、低内存、高效率等优点,只需给出色散介质模型的极点和对应系数,即可运用SARC FDTD递推公式和通用程序进行计算,具有较好的通用性. 通过Debye,Drude和Lorentz三种常用色散介质模型对算法进行了验证. 关键词： 时域有限差分 色散介质 半解析递归卷积     

一种色散介质FDTD通用吸收边界

由单轴各向异性介质所满足的场方程出发,并根据相位匹配原理以得到介质中的色散关系和无反射条件.结合频域到时域的转换关系（即用/t代替jω）和移位算子时域有限差分（SO-FDTD）方法提出了一种适用于各向同性常见色散介质模型,包括德拜模型、洛伦兹模型、德鲁德模型等的通用UPML吸收边界.数值结果表明了该算法的通用性和高效性. 关键词： 时域有限差分 吸收边界 移位算子 色散介质     

色散金属光学特性分析通用时域有限差分方法

结合数字信号处理中的半解析递归卷积(SARC)算法,提出一种用于金属光学特性分析的改进半解析递归时域有限差分方法(SARC FDTD)。该方法保持了SARC算法的绝对稳定性和高精度、低内存等优点,又利用梯形近似使FDTD递推公式更加简洁,计算效率进一步提高,对于各种常用金属色散模型,均只需给出模型的极点和对应的系数,即可运用该算法的程序进行计算,具有较好的通用性。最后,通过金属的高阶Lorentz,Drude-Lorentz和Drude-CP三种常用色散模型对算法进行了验证。