Sound propagation in inhomogeneous waveguides with sound-speed profiles using the multimodal admittance method

The multimodal admittance method and its improvement are presented to deal with various aspects in underwater acoustics, mostly for the sound propagation in inhomogeneous waveguides with sound-speed profiles, arbitrary-shaped liquid-like scatterers, and range-dependent environments. In all cases, the propagation problem governed by the Helmholtz equation is transformed into initial value problems of two coupled first-order evolution equations with respect to the modal components of field quantities(sound pressure and its derivative), by projecting the Helmholtz equation on a constructed orthogonal and complete local basis. The admittance matrix, which is the modal representation of Direchlet-to-Neumann operator, is introduced to compute the first-order evolution equations with no numerical instability caused by evanescent modes. The fourth-order Magnus scheme is used for the numerical integration of differential equations in the numerical implementation. The numerical experiments of sound field in underwater inhomogeneous waveguides generated by point sources are performed. Besides, the numerical results computed by simulation software COMSOL Multiphysics are given to validate the correction of the multimodal admittance method. It is shown that the multimodal admittance method is an efficient and stable numerical method to solve the wave propagation problem in inhomogeneous underwater waveguides with sound-speed profiles, liquid-like scatterers, and range-dependent environments. The extension of the method to more complicated waveguides such as horizontally stratified waveguides is available.     

三个非平行波导间的耦合

提出一种分析非平行波导的简明方法…耦合系数推广法,把非平行波导几何图形的变化直接转化成波导耦合系数中波导间距的变化．即对耦合系数进行了巧妙的推广．把波导间距的变化以指数的形式潜人到耦合微分方程中,建立了非平行三波导系统的耦合方程。这些耦合方程包含了非平行波导传输的全部信息,即非平行波导输出光的振幅和相位,非平行传输波导的长度及其夹角．非平行波导输人端口的间距,能量转换与波导参量的关系等。对这些微分方程作巧妙的数学变换,得到非平行三波导耦合系统输出光的完整分析解,分析和讨论了对称和不对称输人初始条件下三个非平行波导的耦合情况。由此可以优化非平行波导的设计参量。     

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

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

太赫兹金属线波导耦合模理论研究(英文)

基于耦合模假设,研究了描述两太赫兹金属线波导之间耦合特性的耦合模理论.通过采用与光波导耦合模理论等价的近似,研究得到了耦合方程和耦合系数表达式.进而,研究得到了两个相同金属线波导的耦合方程解及耦合特性.所提出的理论对于描述两金属线波导之间的能量传递具有重要的意义.     

一种水平变化波导中声传播问题的耦合模态法

针对介质参数及海底边界水平变化波导中的声传播问题,本文基于多模态导纳法提出一种能量守恒且便于数值稳定求解的耦合模态方法.将声压表示为一组正交完备的本地本征函数之和,对声压满足的Helmholtz方程在本地本征函数上作投影,推导出关于声压模态系数的二阶耦合模态方程组.耦合矩阵直观描述水平变化因素对模态耦合的贡献.为避免直接求解二阶耦合模态方程组可能遇到的数值发散问题,将其重构为两个耦合的一阶演化方程组,引入导纳矩阵并使用Magnus数值积分方法获得稳定的声场解.利用该耦合模态方法数值计算水平变化波导中的声场,并与COMSOL参考解比较,结果表明该耦合模态理论能够精确求解水平变化波导中的点源及分布源传播问题.     

非平行波导耦合理论研究

利用耦合模理论对非平行波导耦合特性进行了研究，得到了非平行波导的普遍性耦合方程及耦合系数表达式.进一步采用小角度近似，给出了几种情况下的计算结果，并与目前报道的利用其他方法计算的结果做了比较.研究表明，不同方法得到的波导中光功率虽然在远距离处都趋于稳定不变，但稳定功率的大小有明显的差别；另外由新提出方法得到的近距离处光功率随着传输距离的起伏变化更为明显.由于该方法同时考虑了非平行波导中光场相位和振幅对耦合方程和耦合系数的影响，因而所得结果能更全面地反映非平行波导的耦合特性，对实际应用更有意义. 关键词： 非平行波导 耦合方程 耦合系数     

Specific features of waveguide heating due to transmission of high-power microwave signals

Waveguide heating due to transmission of microwave signals is studied. Mathematical models are developed to evaluate heat liberation, and differential equations of thermal balance are derived with allowance for different working conditions of waveguides. The results prove the necessity of the further study of the effect of heat liberation in waveguides on strength and functional characteristics.     

多波导阵列耦合研究

提出了一种求解波导阵列耦合特性的方法.采用拉普拉斯变换,将同时考虑相邻和次相邻波导耦合作用的波导阵列耦合模方程组变为一个线性方程组,利用矩阵的LU分解法(将系数矩阵分解为单位下三角阵L与上三角阵U的乘积)求出该线性方程组的解.然后根据给定的初始条件进行拉普拉斯变换反演,求出耦合模方程组的解.采用该方法研究了五波导阵列,光从中间一个波导入射和从五个波导入射时的耦合特性,给出了每个波导中光功率变化规律曲线,并将所得结果与同样条件下仅考虑相邻波导耦合的结果做了比较.采用该方法可以求出任意数目波导组成的阵列,考虑任意波导之间耦合的耦合模方程组的解.     

Analysis of non-linear multilayer waveguides with Kerr-like permittivities

This paper presents a semi-analytical method for calculating the dispersion relationships for stationary non-linear TE waves guided by general multilayer waveguides with Kerr-like permittivities. The non-linear wave equations are solved rigorously for each layer, and the mode index and field distribution are systematically determined so as to satisfy the boundary conditions at all interfaces. To verify the validity of the present method, a non-linear graded-index waveguide is approximated by a number of step layers and the resultant multilayer waveguide is analysed numerically. The convergence properties of mode index and total power flow are presented, and these are compared with direct solutions of the non-linear wave equation by the numerical integration method. The dispersion relationships for a non-linear multiple-quantum-well waveguide have also been investigated under two simple permittivity models. The present method is useful for the analysis of non-linear graded-index waveguides in addition to non-linear multilayer waveguides such as the multiple-quantum-well structure.     

Theoretical Analysis of Resonant Cavities Loaded with Lossy Material

Theoretical analysis of abrupt cavities loaded with lossy material has been described in detail. Dispersion equation for waveguides with lossy material and hybrid modes (HEM) matching equations for abrupt structure were derived. Further, the equations have been solved by numerical calculation to study waveguides and abrupt cavities loaded with lossy material. The theoretical analysis is available for designing and analysis RF devices. It was used to design the RF system for Ka band gyroklystron amplifier, the calculation results accord with the cold test.     

Spinor-electron wave guided modes in coupled quantum wells structures by solving the Dirac equation

A quantum analysis based on the Dirac equation of the propagation of spinor-electron waves in coupled quantum wells, or equivalently coupled electron waveguides, is presented. The complete optical wave equations for Spin-Up (SU) and Spin-Down (SD) spinor-electron waves in these electron guides couplers are derived from the Dirac equation. The relativistic amplitudes and dispersion equations of the spinor-electron wave-guided modes in a planar quantum coupler formed by two coupled quantum wells, or equivalently by two coupled slab electron waveguides, are exactly derived. The main outcomes related to the spinor modal structure, such as the breaking of the non-relativistic degenerate spin states, the appearance of phase shifts associated with the spin polarization and so on, are shown.     

Integral equations for photonic-crystal fibers

The integral and integro-differential equations for photonic-crystal waveguides (fibers) with both infinite and finite quasi-periodic dielectric coatings have been obtained, with previous consideration of the equations for 2D periodic photonic crystals with magnetodielectric and metallic inclusions. The corresponding numerical results are presented.     

会聚激光扫描铌酸锂晶体写入光波导时的最佳曝光间距

对利用激光精细加工技术在铌酸锂(LiNbO₃)晶体中写入光波导时的最佳曝光间距进行了详细地理论分析和实验研究.通过数值求解简化后的光折变动力学方程组，对会聚激光束沿不同方向扫描晶体时的最佳曝光间距进行了详细的数值模拟.结果表明：最佳曝光间距的选取与写入光束的扫描方向无关；当聚焦激光束平行于光轴扫描晶体时，得不到波导结构，但可以采用“三明治”辐照方式得到对称的折射率分布.采用会聚的绿激光束扫描掺铁铌酸锂(LiNbO₃：Fe)晶体进行了实验研究，晶体中的光致折射率变 关键词： 激光精细加工 最佳曝光间距 光波导 3晶体')" href="#">LiNbO₃晶体     

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

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

An efficient approach toward guided mode extraction in two-dimensional photonic crystals

A rigorous, fast and efficient method is proposed for analytical extraction of guided defect modes in two-dimensional photonic crystals, where each Bloch spatial harmonic is expanded in terms of Hermite-Gauss functions. This expansion, after being substituted in Maxwell’s equations, is analytically projected in the Hilbert space spanned by the Hermite-Gauss basis functions, and then a new set of first order coupled linear ordinary differential equations with non-constant coefficients is obtained. This set of equations is solved by employing successive differential transfer matrices, whereupon defect modes, i.e. the guided modes propagating in the straight line-defect photonic crystal waveguides and coupled resonator optical waveguides, are analytically derived. In this fashion, the governing differential equations are converted into an algebraic and easy to solve matrix eigenvalue problem. Thanks to the analyticity of the proposed approach, the eigenmodes of these structures can be extracted very quickly. The validity of the obtained results is however justified by comparing them to those derived by using the standard finite-difference time-domain method.     

Numerical study of Two Photon Absorption effect on nonlinear directional couplers on coupled semiconductor waveguides

The nonlinear interactions of coupled semiconductor waveguides are important phenomena in ultrafast all-optical data processing of integrated photonic circuits. The coupling analysis of semiconductor couplers needs a careful study, because they produce large nonlinear optical responses such as Two Photon Absorption (TPA). In this work, the coupling process is studied by the evolution equations of the mode amplitude of coupled semiconductor waveguides. The formalism takes into account the effects of TPA by starting with Maxwell's equations and using the perturbation method. It is shown that free carriers produced by TPA lead to modification of refractive index. This, in turn, influences the coupling performance of semiconductor waveguides. Fast Fourier Beam Propagation Method (FFT-BPM) is used to investigate the problem, too. The inclusion of free carrier effects result in about one order of magnitude difference in the decoupling threshold.     

Nonstationary nonlinear phenomena in optical slab-waveguides. I. General theory

Generalizing the stationary coupled-mode concept, nonstationary field equations for slowly varying field envelopes are derived. These field equations are combined with the equations of motion for the non-linear polarization and the inversion as well yielding a set of coupled differential equations applicable to a variety of nonstationary phenomena in waveguides. Pulse shaping and SIT of pulses interacting with a thin layer of two-level system are considered in detail.     

Engineering optical gradient force from coupled surface plasmon polariton modes in nanoscale plasmonic waveguides

We explore the dispersion properties and optical gradient forces from mutual coupling of surface plasmon polariton(SPP) modes at two interfaces of nanoscale plasmonic waveguides with hyperbolic metamaterial cladding.With Maxwell's equations and Maxwell stress tensor,we calculate and compare the dispersion relation and optical gradient force for symmetric and antisymmetric SPP modes in two kinds of nanoscale plasmonic waveguides.The numerical results show that the optical gradient force between two coupled hyperbolic metamaterial waveguides can be engineered flexibly by adjusting the waveguide structure parameters.Importantly,an alternative way to boost the optical gradient force is provided through engineering the hyperbolic metamaterial cladding of suitable orientation.These special optical properties will open the door for potential optomechanical applications,such as optical tweezers and actuators.     

Leaky modes of channel waveguides with uniaxial anisotropy

An exact numerical method for analyzing the propagation properties of leaky modes of inhomogeneous channel optical waveguides with a complex uniaxial diagonal permittivity tensor is developed. The method is based on solving the system of integro-differential equations formulated with respect to transversal components of the magnetic field and the longitudinal component of the electric field. Some results of investigation of leaky modes of diffused channel waveguides in LiNbO₃ and LiTaO₃ crystals are given.     

Full-vector mode solver for bending waveguides based on the finite-difference frequency-domain method in cylindrical coordinate systems

A full-vector mode solver for bending waveguides is described based on the finite-difference frequency-domain method in a local cylindrical coordinate system in which the formulas are directly derived from Maxwell's equations with the help of Yee's mesh. The perfectly matched layer absorbing boundary conditions are incorporated into the present approach in order to effectively demonstrate the leaky nature of bending waveguides. A typical bending rib waveguide is considered as a numerical example to show the effectiveness of the established method.