Hybrid numerical technique combining the finite-difference beam-propagation method and the finite-difference time-domain method

A technique combining a scalar wide-angle finite-difference beam-propagation method with a scalar finite-difference time-domain method is proposed in cylindrical coordinates. The propagating field in a lensed coreless fiber is simulated, and the effect of a lens on the far-field pattern is evaluated.     

层状介质时域有限差分方法斜入射平面波引入新方式

应用二维时域有限差分方法分析层状介质中的目标散射时,在总场-散射场边界斜入射平面波源用常规方法难以引入,因为在总场-散射场边界处设置的入射波实际上包含了入射脉冲以及各分层界面的反射和多次反射.为解决这个问题,提出了斜入射平面波的混合引入方式,即对总场-散射场的四个边界面采取不同的处理方式.对于总场-散射场的纵向侧边界,用含有斜入射角度的修正一维时域有限差分方法,只要在自由空间位置加入入射脉冲就会自行产生由各分层界面形成的反射波,包括多次反射.同时,把纵向总场-散射场侧边界向下延伸,使得总场-散射场下边界位于完全匹配层内,这样透射波和散射波均为外向行波而被吸收.对于总场-散射场的上边界,由于完全位于自由空间中,边界上各点的入射波将是总场-散射场纵向边界角点处入射波的带有时间延迟的复制.数值模拟结果表明了本文所提出方法的正确性和有效性. 关键词： 时域有限差分 层状介质 斜入射平面波 修正一维麦克斯韦方程     

斜入射分层线性各向异性等离子体电磁散射时域有限差分方法分析

提出了斜入射分层线性各向异性等离子体电磁散射的时域有限差分(FDTD)方法,通过将二维麦克斯韦方程等价地转换为一维麦克斯韦方程,避免了用二维时域有限差分方法分析该散射问题,极大地提高了计算效率.分析推导了TEz和TMz波斜入射线性分层各向异性等离子体电磁散射的FDTD方法,然后通过该方法计算不同入射角的各向异性等离子板的电磁波反射系数,并与其解析解进行比较,结果表明该方法的准确性和有效性.最后,将该算法应用于计算涂覆分层各向异性等离子体金属板在不同入射角下的反射系数,分析了不同入射角对反射系数的影响.     

Mode-coupling analysis of three-dimensional microdisk resonators by the finite-difference time-domain technique

Quality factor enhancement due to mode coupling is observed in a three-dimensional microdisk resonator. The microdisk, which is vertically sandwiched between air and a substrate, with a radius of 1 microm, a thickness of 0.2 microm, and a refractive index of 3.4, is considered in a finite-difference time-domain (FDTD) numerical simulation. The mode quality factor of the fundamental mode HE71 decreases with an increase of the refractive index of the substrate, n(sub), from 2.0 to 3.17. However, the mode quality factor of the first-order mode HE72 reaches a peak value at n(sub) = 2.7 because of the mode coupling between the fundamental and the first-order modes. The variation of mode field distributions due to the mode coupling is also observed. This mechanism may be used to realize high-quality-factor modes in microdisks with high-refractive-index substrates.     

球面电磁波斜入射建筑物内场的数值分析

为了分析窄带电磁脉冲源近、远场辐照大型建筑物内部空间电磁场分布特点,采用时域有限差分方法,对比分析平面波、球面波场源正面斜入射多层多单元建筑物的计算结果,并较为全面地分析了各房间各水平面中心电场强度幅值、各层各水平面电场强度最大值分布情况。建筑物各层相对应高度面上场强分布近似,在电磁波传播方向上窗户房间内部区域场强较强,两种波场源入射结果基本一致;在其余区域场强较弱,特别是在较大空间内的区域,球面波场源入射下场强相对更弱。其中通过与待模拟建筑物外形相紧凑性匹配的球面波场源构建,模拟仿真了窄带球面电磁波在真空空间的辐射传播,计算结果验证所构建球面波近场源准确可行。     

球面电磁波斜入射建筑物内场的数值分析

为了分析窄带电磁脉冲源近、远场辐照大型建筑物内部空间电磁场分布特点,采用时域有限差分方法,对比分析平面波、球面波场源正面斜入射多层多单元建筑物的计算结果,并较为全面地分析了各房间各水平面中心电场强度幅值、各层各水平面电场强度最大值分布情况。建筑物各层相对应高度面上场强分布近似,在电磁波传播方向上窗户房间内部区域场强较强,两种波场源入射结果基本一致;在其余区域场强较弱,特别是在较大空间内的区域,球面波场源入射下场强相对更弱。其中通过与待模拟建筑物外形相紧凑性匹配的球面波场源构建,模拟仿真了窄带球面电磁波在真空空间的辐射传播,计算结果验证所构建球面波近场源准确可行。     

TE-TE mode coupling at oblique incidence in a periodic dielectric waveguide

The characteristics of reflection of a TE mode into another TE mode in a planar dielectric waveguide with a sinusoidally modulated index are investigated for the general case in which the wavevectors are not aligned with the grating vector. A singular perturbation procedure is used to deduce the asymptotically exact coupled mode equations which govern the interaction. An analysis of the coupled-mode equations has shown the existence of the Brewster phenomenon at the angle of incidence θ_iB=45°. For angles of incidence sufficiently greater than the Brewster angle, the reflection coefficient can be considerably larger than that for the normal incidence. In the spectral response curves of the reflection coefficient, a higher maximum and a smaller base bandwidth are possible for the oblique incidence as compared to that for the normal incidence. The angular sensitivity of the reflection characteristics increases as the angle of incidence changes from the normal incidence to the grazing incidence.     

Uniform stable conformal convolutional perfectly matched layer for enlarged cell technique conformal finite-difference time-domain method

Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer(CPML) is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain(ECT-CFDTD) method is used to simulate the wave propagation inside a perfect electric conductor(PEC) waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.     

Acoustic analysis of the vocal tract during vowel production by finite-difference time-domain method

The vocal tract shape is three-dimensionally complex. For accurate acoustic analysis, a finite-difference time-domain method was introduced in the present study. By this method, transfer functions of the vocal tract for the five Japanese vowels were calculated from three-dimensionally reconstructed magnetic resonance imaging (MRI) data. The calculated transfer functions were compared with those obtained from acoustic measurements of vocal tract physical models precisely constructed from the same MRI data. Calculated transfer functions agreed well with measured ones up to 10 kHz. Acoustic effects of the piriform fossae, epiglottic valleculae, and inter-dental spaces were also examined. They caused spectral changes by generating dips. The amount of change was significant for the piriform fossae, while it was almost negligible for the other two. The piriform fossae and valleculae generated spectral dips for all the vowels. The dip frequencies of the piriform fossae were almost stable, while those of the valleculae varied among vowels. The inter-dental spaces generated very small spectral dips below 2.5 kHz for the high and middle vowels. In addition, transverse resonances within the oral cavity generated small spectral dips above 4 kHz for the low vowels.     

Subpixel smoothing for conductive and dispersive media in the finite-difference time-domain method

Staircasing of media properties is one of the intrinsic problems of the finite-difference time-domain method, which reduces its accuracy. There are different approaches for solving this problem, and the most successful of them are based on correct approximation of inverse permittivity tensor epsilon(-1) at the material interface. We report an application of this tensor method for conductive and dispersive media. For validation, comparisons with analytical solutions and various other subpixel smoothing methods are performed for the Mie scattering from a small sphere.     

时域有限差分法在建筑声学中的应用及前景

介绍了声波方程的基本差分格式及稳定条件、数值色散、吸收边界条件等数值计算理论，例举了前人用时域有限差分法对噪声传播过程的模拟和室内声学中座椅吸声低谷效应模拟的模拟结果。本文指出，由于时域有限差分法的特点使其具有在模拟脉冲响应方而的特别优势，因而，应用这一技术研究厅堂的声学特性，尤其是低频特性，将有重要意义。     

色散介质电磁特性时域有限差分分析的Newmark方法

根据Debye模型、Drude模型和Lorentz模型3种常见色散介质模型频域极化率的特点,利用频域到时域的转换关系jω→?/?t,将极化矢量P与电场强度E的频域关系转换成时域内关于P的二阶微分方程,其对3种色散介质模型皆适用,具有统一的形式.然后采用相比于中心差分具有更高精度的Newmark两步算法(Newmark-β-γ法)求解该方程,进而得到E→P的递推公式,再结合本构关系得到D→E的时域递推式.实现了色散介质电磁场量的时域有限差分迭代计算.数值计算结果表明该方法是适用于3种色散介质模型的通用算法,并且相比于移位算子时域有限差分方法等以中心差分为基础的离散方案具有更高的计算精度.     

Locally conformal method for acoustic finite-difference time-domain modeling of rigid surfaces

The finite-difference time-domain method is a simple but powerful numerical method for simulating full-wave acoustic propagation and scattering. However, the method can demand a large amount of computational resources. Traditionally, continuously curved boundaries are represented in a stair-step fashion and thus accurately modeling scattering from a boundary will require a finer discretization than would otherwise be necessary for modeling propagation in a homogeneous medium. However, a fine discretization might not be practical due to limited computational resources. A locally conformal technique is presented here for modeling acoustic scattering from continuously curved rigid boundaries. This technique is low cost, simple to implement, and gives better results for the same grid discretization than the traditional stair-step representation. These improvements can be traded for a coarser discretization which reduces the computational burden. The improved accuracy of this technique is demonstrated for a spherical scatterer.     

Full-wave parallel dispersive finite-difference time-domain modeling of three-dimensional electromagnetic cloaking structures

A parallel dispersive finite-difference time-domain (FDTD) method for the modeling of three-dimensional (3-D) electromagnetic cloaking structures is presented in this paper. The permittivity and permeability of the cloak are mapped to the Drude dispersion model and taken into account in FDTD simulations using an auxiliary differential equation (ADE) method. It is shown that the correction of numerical material parameters and the slow switching-on of source are necessary to ensure stable and convergent single-frequency simulations. Numerical results from wideband simulations demonstrate that waves passing through a three-dimensional cloak experience considerable delay comparing with the free space propagations, as well as pulse broadening and blue-shift effects.     

Auxiliary differential equation: efficient implementation in the finite-difference time-domain method

An efficient algorithm for modeling dispersive media in finite-difference time-domain methods is presented. It is based on the auxiliary differential equation method for treatment of Lorentz media with an arbitrary number of relaxations. The algorithm shows excellent accuracy of second order in time and space and is efficient in both memory requirements and computational effort.     

共形单步交替方向隐式时域有限差分方法及其改进

提出了一种基于共形网格技术的共形单步交替方向隐式时域有限差分(CLeapfrog ADI-FDTD)方法。与常规FDTD方法相比,此方法能够减小由于目标边界不契合网格划分而引入的阶梯近似误差,提高算法计算不规则目标时的精度;同时算法稳定性更强,计算效率更高。由于引入共形技术后显著降低了原差分法的无条件稳定性,本文利用增长矩阵本征值方法理论分析了算法的稳定性,然后采用了一种改进的共形面积计算方法,在此基础上提出了一种稳定性更高的改进的共形单步交替方向隐式时域有限差分(ICLeapfrog ADI-FDTD)方法。数值算例验证了ICLeapfrog ADI-FDTD是一种具有高稳定性和高精度的高效算法。     

一种基于数字信号处理技术的改进通用色散介质移位算子时域有限差分方法

借鉴数字信号处理技术中的无限脉冲响应滤波器实现的思想,提出一种适用于通用色散介质的改进移位算子时域有限差分（shift operator finite-difference time-domain,SO-FDTD）方法.与原SO-FDTD方法相比,改进的SO-FDTD方法内存占用减少33%以上,同时计算时间也有所减少.最后通过一阶、二阶色散介质雷达散射截面计算验证了方法的通用性和有效性. 关键词： 时域有限差分 色散介质 移位算子     

Effective optical response of silicon to sunlight in the finite-difference time-domain method

The frequency dependent dielectric permittivity of dispersive materials is commonly modeled as a rational polynomial based on multiple Debye, Drude, or Lorentz terms in the finite-difference time-domain (FDTD) method. We identify a simple effective model in which dielectric polarization depends both on the electric field and its first time derivative. This enables nearly exact FDTD simulation of light propagation and absorption in silicon in the spectral range of 300-1000 nm. Numerical precision of our model is demonstrated for Mie scattering from a silicon sphere and solar absorption in a silicon nanowire photonic crystal.     

坐席低谷效应的时域有限差分多项模拟

介绍了交错网格下,以声压和质点振动速度为声场参量的室内声学时域有限差分递推格式,给出了在建筑边界截断计算区域的边界条件差分格式。在水平地面房间模型下,应用时域有限差分法模拟了坐席低谷效应对于空间和时间的多种特性;此外,模拟了地面起坡时的低谷效应。     

Simpler and more intuitive interpretation of stability condition for finite-difference time-domain method

In using a popular numerical electromagnetic solver, finite-different time-domain (FDTD) method, temporal and spatial discretization must obey a so-called stability condition. This short paper provides straightforward and easy-to-understand interpretation for the condition than a frequently referred traditional one.