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

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

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

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

基于时域有限差分算法改进卷积完全匹配层的稳定性

为了克服时域有限差分算法中卷积完全匹配层对消逝波吸收效果差的缺点,提出一种在卷积完全匹配层后添加特殊吸收层的方法.在不增加物体与吸收层内层距离的情况下,通过调节特殊吸收层中两个衰减因子,使其为常数,并令吸收因子逐层从1增加到10,来增强吸收层对消逝波的吸收性能.平面波垂直入射到单层光子晶体的算例表明,添加了特殊吸收层的吸收边界在与散射体相距5个网格的情况下仍能够保持计算结果收敛,而传统的吸收边界则需要相距80个网格才能保证结果收敛,说明该方法提高了对消逝波的吸收性能.进一步在结构中采用此吸收边界来计算多层光子晶体的传输特性曲线,并将其与常规方法计算所得结果做比较,两种结果吻合较好.数值算例验证了该方法的有效性和正确性.     

应用高阶有限元-局部共形完全匹配层算法分析电磁散射问题

提出局部共形完全匹配层在高阶有限元中的实现方法,用于分析三维电磁散射问题.基于定义在局部高阶单元上的实坐标、复坐标以及参数坐标三者之间的映射关系,局部共形完全匹配层不仅能够贴近复杂结构散射体的外廓形状设置截断区域,有效地缩减有限元分析区域,而且由于保证了较高的几何拟合精度,可以有效地提高数值计算结果的准确性.利用这种结合局部共形完全匹配层的高阶有限元方法分析不同目标的电磁散射特性,数值结果表明方法的正确性和有效性.     

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

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

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

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

共形完全匹配层按层积分算法

基于共形完全匹配层(CPML)吸收边界,构造一种共形完全匹配层矢量单元按层积分算法,将多层单元积分运算叠加到一层单元中进行,用多层单元剖分,以一层单元计算矩阵元素,即保留了多层单元的几何和材料信息,又减少了单元数量和计算量.数值算例表明,这种按层积分CPML吸收效果好,计算量小,可靠,高效.     

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

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

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

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

无源光器件的传播模场分析

以波导结构的单向耦合器为例，将复坐标变量完全匹配层应用于光器件的模场分析，采用有限差分方法构造出模场的本征方程，通过数值求解得到器件的传播模场和等效折射率。结果表明等效折射率随间距变化而变化，间距越大，等效折射率越趋向于单根光波导的值，为准确分析无源光器件内部模场分布打下了基础。     

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.     

Computation of Microwave Attenuation on Coplanar Waveguide with Conformal ADI-FDTD Method

A novel Conformal Finite-Difference Time-Domain based on Alternating-Direction Implicit Method(ADI-CFDTD) used to calculate the microwave attenuation on coplanar waveguide(CPW) with complicated cross-section up to 200GHz is presented. This ADI-CFDTD algorithm is the combination of conformal algorithm dealing with the deformed cell with ADI-FDTD, and has unconditionally stable. The difference equations are derived from the Faraday's and Ampere's law. To validate its accuracy and efficiency, as an example, a CPW fabricated on lithium niobate is proposed. Numerical results demonstrate that this new ADI-CFDTD algorithm has more accuracy than that of quasi-static and CFDTD and shows good agreement with experiment data. When frequency is up to millimeter wave band the radiation loss is not negligible in total power loss.     

Analytic-numerical approach to non-linear problems in dielectric waveguides

A combined analytic-numerical approach is presented for investigations of electromagnetic signal propagation in dielectric waveguides containing a non-linear medium. The approach is based on the resolvent method for the Volterra integral equation describing an electromagnetic process. The proposed algorithm uses the exact expression for the solution to the integral equation for the electromagnetic field in the planar waveguide with time jump changing of the medium properties. In the case of a non-linear dielectric as the waveguide core the continuous change of the medium parameters is replaced by jump changes of the latter. This is the only approximation made in the algorithm. The resolvents that are constructed give the exact solution on each jump step. The results from some example numerical calculations are presented.     

An efficient locally one-dimensional finite-difference time-domain method based on the conformal scheme

An efficient conformal locally one-dimensional finite-difference time-domain(LOD-CFDTD) method is presented for solving two-dimensional(2D) electromagnetic(EM) scattering problems. The formulation for the 2D transverse-electric(TE) case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit(ADI) CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field(TF/SF) boundary and the perfectly matched layer(PML), the radar cross section(RCS) of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.     

Evaluation of optical propagation and radiation in optical waveguide using a numerical method

We introduce a mathematical model based on a concept of intrinsic mode in order to analyse and synthesise optical wave propagation and radiation occurring in a non-uniform optical waveguide used in integrated optics as optical coupler. The model is based on numerical evaluation of electromagnetic wave by applying an intrinsic field integral to evaluate the field behaviour inside the optical waveguide. To analyse the field distribution inside the non-uniform waveguide and predict the beam propagation of optical energy involved in the propagation process, it is necessary to track the motion of any observation point along the tapered waveguide itself. Physically, the rays of the spectrum undergo reflections on the waveguide boundaries until the cut-off occurs and the phenomena of radiation begin. The numerical results show good agreement with those obtained by classical methods of evaluation used bv other works.     

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

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

Application of convolution perfectly matched layer in finite element method calculation for 2D acoustic wave

A method was presented to extend the Convolution Perfectly Matched Layer(CPML), which bases on the complex coordinates transformation and complex frequency shifted stretched-coordinate metrics,to the 2D acoustic equation calculated with the method of Finite Element Method(FEM).This non-physical layer is used at the computational edge of a FEM as an Absorbing Boundary Condition(ABC) to truncate unbounded media.In this paper,the CPML equations have been presented in frequency domain and in time domain,respectively,and the calculations have been realized in the FEM software of COMSOL.The main advantage of CPML over the classical PML layer is that it is based on the unsplit components of the wave field leading to a more stable,highly effective absorption and a more facility to realize.The results of numerical simulation demonstrate that CPML has better absorbability than PML and it absorbs the outgoing energy more effectively.     

瞬变电磁场模拟中的CPML吸收边界条件

提出用于瞬变电磁法(TEM)模拟的时域有限差分(FDTD)算法中的卷积完全匹配层(CPML)吸收边界条件.首先,在磁场散度方程显式处理条件下,计算磁场z分量时,构造出一种计算时域卷积项的方法,并详细推导计算磁场z分量的表达式.而后,对均匀半空间模型进行数值计算.结果表明,提出的CPML吸收边界条件在保证计算精度的前提下,实现了瞬变电磁法时域有限差分高效正演计算.