半解析对偶棱边元及其在波导不连续性问题中的应用

给出电磁波导的对偶变量变分原理,并采用对偶棱边元对波导的横截面进行半解析离散. 将波导中沿纵向均匀的区段视为子结构,运用基于Riccati方程的精细积分算法求出其出口刚度阵,然后与不均匀区段的常规有限元网格拼装即可对波导不连续性问题进行求解. 半解析对偶棱边元的采用可以在最大程度上对有限元网格进行缩减,并且能够在不增加计算量的前提下任意增加子结构的长度,从而可以将截断求解区域的人工边界设置在距离不均匀区段充分远的地方,极大地减少了近似边界条件所带来的误差. 数值算例证明这种方法具有很高的精度与效率. 关键词： 波导的不连续性 半解析辛分析 对偶棱边元 精细积分     

电磁波导的辛分析与对偶棱边元

将电磁波导的控制方程导向了Hamilton体系、辛几何的形式.以电磁场的横向分量组成对偶向量并采用分离变量法，可以得到Hamilton算子矩阵的辛本征值问题.共轭辛正交归一关系、辛本征解展开定理等均可在此应用.对于复杂横截面和填充非均匀材料的电磁波导，提出对偶棱边元，对截面半解析离散后即可进行数值求解.对偶棱边元克服了结点基有限元求解电磁场问题的困难，与常规棱边元相比在某些方面具有一定的优势. 关键词： 电磁波导 Hamilton体系 对偶变量 棱边元     

光波导的条形传递函数无限元解

将无限元引入了条形传递函数方法并用于对光波导的分析。为了对光波导的无穷区域进行处理引入了两种无限单元,一种为解析方向单元,另一种为离散方向单元。在构造无限元的过程中,为了求解方便,利用光波导标量波动方程的特点引入了二阶传递函数解法。最后以矩形波导为例计算并与其他方法进行了比较,结果表明本方法在分析光波导在近截止频率附近的传播特性方面优势明显,并可节省大量的计算。     

水平变化波导中的全局矩阵耦合简正波方法

提出了一种新的水平变化波导中声场的耦合简正波求解方法,该方法能够处理二维点源和线源问题,提供声场的双向解。该方法利用全局矩阵(DGM)一次性求解耦合模式的系数,消除了传播矩阵递推求解中存在的误差累积问题;此外,改善了现有模型中对距离函数的归一化方法,从而避免了泄露模式指数增长导致的数值溢出问题。本文还给出了绝对软海底理想波导中耦合矩阵的闭合表达式,并分析了单个阶梯下简正波耦合现象。此外,本文还计算了理想楔形波导中的声传播问题(ASA标准问题),并与解析解及COUPLE07计算结果进行了比较,结果表明该方法是一种稳定、精确的水平变化波导中的声场计算方法。      

离散修正KdV方程的解析近似解

将同伦分析法进行了推广,使之适用于求解离散修正KdV方程.获得了由指数函数表达的亮孤子解,该解析近似解与精确解符合很好.数值模拟结果说明了同伦分析法对求解复杂非线性问题的有效性和潜力.     

辛体系下电磁波导传输波的截止频率和传播常数的求解

从变分原理出发,消去纵向分量,选取横向的电场和磁场组成对偶向量,将麦克斯韦基本方程导向了哈密顿体系,即辛几何的形式.对电磁波导的横截面采用有限元方法进行离散,得到一组线性代数方程组,该方程组可以用于求解波导传输波的截止频率以及传播常数.数值算例验证了本文方法的正确性.该方法为电磁波导的分析计算提供了一种新的研究方法.     

水平变化双层波导声传播问题的耦合简正波解

为了验证现有模型的精度,导出了全反射下边界双层波导中简正波耦合矩阵的解析表达式,并将其应用到全局矩阵耦合简正波模型(Direct Global Matrix Coupled-Mode)中,使得该模型可以提供水平变化双层波导问题的标准解。文中首先利用COUPLE的简正波及耦合矩阵数值解验证了该简正波及耦合矩阵解析表达式的正确性;其次,采用改进的DGMCM模型求解了双层波导海山声传播损失,结果表明,改进后的DGMCM模型可以非常精确地求解水平变化双层波导问题,可作为求解此类问题的标准模型使用。     

基于等几何分析方法求解任意截面波导本征问题

基于等几何分析方法具有自由度花费少、高精度、高阶连续性等特点,通过加权余量法对椭圆波导本征问题的亥姆霍兹方程等几何离散得出等几何分析方程.解决了传统数值方法的求解域与几何模型的非一致性问题,实现了将问题的分析计算构架于精确几何模型基础之上.分析任意截面波导的本征问题,对不同偏心率的椭圆波导以及三角形和五边形波导的截止波数的求解结果显示等几何分析方法求解波导本征问题的高效及高精度特性.与传统方法相比,此方法以较少的自由度消耗便会达到较高的求解精度,并且数值解的收敛率较快.     

精细积分法在含各向异性介质波导不连续性问题中的应用

用精细积分法对含各向异性介质的波导不连续性问题进行了数值模拟与分析. 从矢量波动方程相对应的单变量变分形式出发, 推导出了含有各向异性介质波导横截面离散系数矩阵的表达式, 引入对偶变量, 在Hamilton体系下, 利用精细积分法求出出口刚度矩阵, 进行有限元拼装, 求解了含各向异性介质的波导不连续性问题. 算例表明了该方法的准确性和高效性. 利用本文方法还讨论了介电系数和导磁系数张量的各个分量对波导传输特性的影响. 关键词： 波导不连续性 各向异性介质 Hamilton体系 精细积分法     

半圆形波导中的电磁场

利用矩形波导和圆形波导的结论,并结合叠加法和镜象法等工具,从理论上分析给出半圆形波导的TE波模型的电磁场分量.根据金属界面的特殊性质,在一定条件下运用叠加法的延拓求出半圆形波导的电磁场分量.最后通过过渡法对一般情况下半圆形波导的电磁场分量进行求解.     

A semi-analytical finite element formulation for modeling stress wave propagation in axisymmetric damped waveguides

A semi-analytical finite element (SAFE) method is presented for analyzing the wave propagation in viscoelastic axisymmetric waveguides. The approach extends a recent study presented by the authors, in which the general SAFE method was extended to account for material damping. The formulation presented in this paper uses the cylindrical coordinates to reduce the finite element discretization over the waveguide cross-section to a mono-dimensional mesh. The algorithm is validated by comparing the dispersion results with viscoelastic cases for which a Superposition of Partial Bulk Waves solution is known. The formulation accurately predicts dispersion properties and does not show any missing root. Applications to viscoelastic axisymmetric waveguides with varying mechanical and geometrical properties are presented.     

Analysis of the Scattering Characteristic of Several Waveguide Components Using Boundary Element Method

Several kinds of waveguide components such as curved waveguide bends, arbitrary angle waveguide bends and T-junctions have been analyzed with boundary element method in this paper. A new discretization method for the boundary element method to solve the waveguide discontinuities has been given. The numerical results obtained agree well with the experimental results and numerical results in other literature. Especially, the scattering characteristics of Forded E-, H-plane T-junctions in 3mm band have been analyzed using boundary element method and the calculation results are presented.     

光波导等参元分析

本文提出一套用于分析任意介质光波导的等参元算法,它基于严格的E_z-H_z变分公式,采用八节点曲边四边形单元。这种类型的等参元分析法尤其适用于研究各种具有曲边边界的复杂光波导和光波导器件。作为例子,文中把这一方法用于分析阶跃型光纤,渐变型光纤、各向同性以及各向异性光纤耦合器。     

On the forced response of waveguides using the wave and finite element method

The forced response of waveguides subjected to time harmonic loading is treated. The approach starts with the wave and finite element (WFE) method where a segment of the waveguide is modeled using traditional finite element methods. The mass and stiffness matrices of the segment are used to formulate an eigenvalue problem whose solution yields the wave properties of the waveguide. The WFE formulation is used to obtain the response of the waveguide to a convected harmonic pressure (CHP). Since the Fourier transform of the response to a general excitation is a linear combination of the responses to CHPs, the response to a general excitation can be obtained via an inverse Fourier transform process. This is evaluated analytically using contour integration and the residue theorem. Hence, the approach presented herein enables the response of a waveguide to general loading to be found by: (a) modeling a segment of the waveguide using finite element methods and post-processing it to obtain the wave characteristics, (b) using Fourier transform and contour integration to obtain the wave amplitudes and (c) using the wave amplitudes to find the response at any point in the waveguide. Numerical examples are presented.     

左手介质矩形波导导模和表面模的场分布

从麦克斯韦方程组出发,结合电磁场的边界条件,推导出介质矩形波导导模的一般色散方程.对普通介质矩形波导和左手介质矩形波导的导模场分布分别进行了数值模拟.通过对比两种介质矩形波导的导模场分布的模拟结果,发现左手介质矩形波导Ex22模的场分布比普通介质波导Ex00模的场更集中在波导中部.同时,根据处理普通介质矩形波导的Marcatili方法.类比得到左手介质矩形波导表面模的色散方程,并数值模拟了低阶模的场分布,结果表明,能量主要集中在波导的四个角区以及波导的边缘.     

用极坐标下的广角光束传播法计算分析弯曲波导

采用了极坐标下的广角BPM(Beam Propagation Method)对弯曲波导进行了计算模拟,并与极坐标下近轴BPM方法进行了比较,取得很好的结果.     

Efficient Analysis of Millimeter Wave Ferrite Circulators by the GMRES Iterative Algorithm

The finite element method (FEM) combined with the perfectly matched layers (PML) is given for simulation of waveguide ferrite circulators. The generalized minimal residual (GMRES) iterative method is applied to solve such sparse large non-symmetric system of linear equations resulting from the use of edge-based finite element method. The formulation of FEM and the algorithm of GMRES method are described in detail. The reflection and insertion losses of millimeter wave waveguide circulator are analyzed and the results are compared with those obtained from literature.     

Analysis of temperature profiles of thermo-optic waveguides

The temperature profiles of thermo-optic waveguides are analyzed by the finite element method. The heat generated by a thin-film heater causes the temperature of a nearby waveguide to increase. The analysis results show that thermal coupling is a function of the waveguide spacing and depth. Thermal coupling increases with the waveguide depth but decreases with the waveguide spacing. Thermal coupling could be reduced by placing a cooler on the adjacent waveguide or etching a deep trench between the waveguides. The cooler can reduce the coupling, but it is not effective for deep waveguides. For the trenched structure, the temperature of the heated waveguide increases as the trench depth increases; however, the temperature of the nearby waveguide will decrease.     

A novel wavelength dispersive device with a dispersive element based on staircase-like straight and parallel arrayed waveguides

We propose a new type of arrayed waveguide grating (AWG) multiplexer/demultiplexer based on modified group refractive index. This device is composed by an array of straight and parallel waveguides of equal length and each waveguide consist of two sections with different width. The length of the two sections are changed from a waveguide to the adjacent one following a linear dependence resulting in a wavelength dispersive waveguide array. An example of the device design for silicon-on-insulator (SOI) platform is provided and numerical simulations have been carried out for various arrayed waveguide parameters. We demonstrate that the group index modification can be used for tailoring device dispersion properties, and that it can also result in new dispersion characteristics predicted numerically not observed in conventional AWGs. Additional advantages are that the demultiplexer does not necessarily require bending waveguide sections as in a conventional AWG (de)multiplexers, and thus yields highly compact devices with potentially very low insertion loss. Channel spacing of 1 nm have been predicted for sub-micron waveguides sizes. In this paper it is also proposed a novel wavefront converter based on waveguide array lens-like element with waveguides broadened sections. Numerical results for different input/output geometries are analized.     

A linear-operator formalism for bianisotropic waveguides

A bianisotropic waveguide can be defined as a cylindrical waveguide filled with bianisotropic materials, and all the conventional waveguides are special cases of the bianisotropic waveguide. In this paper, guided wave propagation in bianisotropic waveguide is analyzed by the theory of linear operators, and two types of adjoint waveguides and inner products are introduced respectively. Based on the concept of adjoint waveguides, the functional expressions of the field equations can be obtained, and from which the eigenvalue problem of the bianisotropic waveguide can be solved. Also, bi-orthogonality relations of guided modes are derived. These biorthogonality relations reported here can be used to expand electromagnetic fields in terms of a complete set of modes in straight bianisotropic waveguide. As an example of application, mode matching formulae for a discontinuity problem are given.