全局稳定性一致的共形卷积完全匹配层

在正交网格体系中建立物理模型共形描述的基础上,针对采用扩展元胞技术的共形时域有限差分(ECT-CFDTD)方法模拟计算波导器件遇到的开放端口截断问题,给出了积分形式的共形卷积完全匹配层方法,算法具有与ECT-CFDTD相同的数值稳定性。设置不同的完全匹配层的控制参数,对波导中有消逝波存在的情况进行长时间模拟计算,分析共形卷积完全匹配层对消逝波的长效截断能力,分析卷积完全匹配层的截断误差。计算结果显示:积分形式的共形卷积完全匹配层可有效截断波导器件的开放端口。     

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

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

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

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

二维非均匀介质中大地电磁响应的计算

在二维埋均匀介质中,将大地电磁问题可看作电磁波在有耗开放波导中的传播问题,借用研究波导问题的数值模式匹配法予以解决,在计算Ex型问题时,使用Hermite基函数和非均匀无素;在计算Hx型问题时,使用变型Hermite基函数,并与有限元法进行了比较。     

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

将基于复坐标变换和复频移扩展坐标变量的卷积完全匹配层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.     

NETWORK CHARACTERIZATION OF STEP DISCONTINUITY WITH GENERALIZED CONSERVATION OF COMPLEX POWER TECHNIQUE AND ITS APPLICATION TO ANALYZING PERIODIC DIELECTRIC WAVEGUIDES

The propagation of electromagnetic waves along open periodic, dielectric waveguides is formulated in the case that surface wave is guided and propagates normally to the corrugation. Our approximate analysis with the propagation characteristics is to consider a corresponding bounded waveguide problem in which perfect electric or magnetic walls are introduced, and the periodic corrugation is regarded as consisting of step discontinuities connected by a length of uniform slab waveguide. By properly taking into account of both surface modes and only a few non-surface-modes, a novel network approach is proposed for characterizing step discontinuity based on the generalized conservation of complex power technique (GCCPT). Employing solution selection rule (SSR), we can readily derive propagation characteristics in the Bragg interaction region. A number of numerical results are shown to demonstrate the usefulness of our approach.     

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.     

Guidance of Surface Waves with Dielectric Waveguides Having Finite or Infinite Periodic Corrugations

A planar open dielectric waveguide with periodic rectangular corrugations is investigated in the case that surface wave is guided and propagates normally to the corrugation. Our approximate analysis with the propagation characteristics is to consider a corresponding bounded waveguide problem in which perfect electric or magnetic walls are introduced, and the periodic corrugation is regarded as consisting of step discontinuities connected by a length of uniform slab waveguide. By properly taking into account of both surface modes and only a few non-surface-modes, and using conservation of complex power technique (CCPT) as well as solution selection rule (SSR), we can readily derive propagation characteristics in the Bragg interaction region. The calculated results show an excellent agreement with previously published ones.     

The contribution of a lateral wave in simulating low-frequency sound fields in an irregular waveguide with a liquid bottom

A further development of a previously proposed approach to calculating the sound field in an arbitrarily irregular ocean is presented. The approach is based on solving the first-order causal mode equations, which are equivalent to the boundary-value problem for acoustic wave equations in terms of the cross-section method. For the mode functions depending on the horizontal coordinate, additional terms are introduced in the cross-section equations to allow for the multilayer structure of the medium. A numerical solution to the causal equations is sought using the fundamental matrix equation. For the modes of the discrete spectrum and two fixed low frequencies, calculations are performed for an irregular two-layer waveguide model with fluid sediments, which is close to the actual conditions of low-frequency sound propagation in the coastal zone of the oceanic shelf. The calculated propagation loss curves are used as an example for comparison with results that can be obtained for the given waveguide model with the use of adiabatic and one-way propagation approximations.     

3-D Quantum Transport Solver Based on the Perfectly Matched Layer and Spectral Element Methods for the Simulation of Semiconductor Nanodevices

A 3-D quantum transport solver based on the spectral element method (SEM) and perfectly matched layer (PML) is introduced to solve the 3-D Schr?dinger equation with a tensor effective mass. In this solver, the influence of the environment is replaced with the artificial PML open boundary extended beyond the contact regions of the device. These contact regions are treated as waveguides with known incident waves from waveguide mode solutions. As the transmitted wave function is treated as a total wave, there is no need to decompose it into waveguide modes, thus significantly simplifying the problem in comparison with conventional open boundary conditions. The spectral element method leads to an exponentially improving accuracy with the increase in the polynomial order and sampling points. The PML region can be designed such that less than -100 dB outgoing waves are reflected by this artificial material. The computational efficiency of the SEM solver is demonstrated by comparing the numerical and analytical results from waveguide and plane-wave examples, and its utility is illustrated by multiple-terminal devices and semiconductor nanotube devices.     

Quantum conformal fluctuations revisited

The conformal quantization method of Narlikar and Padmanabhan is reformulated with a view to take into account theexact propagator and to provide explicitnumerical estimates of various predictions for dust cosmologies. It is found that in spite of the divergence of quantum fluctuations at the big-bang epoch it is possible to construct wave packets which remain sharp fromt=10^–70s, say, up to the present epoch provided the present state is finely tuned to the classical one. Also, if the transition probability from the Minkowski to the FRW metric is calculated using Gaussian wave functions (with zero mean) then thet ^2/3 models withk = 0, ± 1 cannot be distinguished, i.e., a fine tuning to the flat (k=0) model does not seem to result if the conformal factor depends on time only.     

Low-frequency radio wave propagation in the earth-ionosphere waveguide disturbed by a large-scale three-dimensional irregularity

In this paper, we develop further the analytical and numerical method of solving three-dimensional problems in the theory of radio wave propagation, including three-dimensional local inhomogeneities (ionospheric disturbances or Earth’s surface irregularities). To model the Earth-ionosphere waveguide, we use the surface impedance concept, by which the irregularity extending beyond one waveguide wall has an arbitrary smooth shape, and its surface can be described by the impedance. In the scalar approximation, this problem is reduced to a two-dimensional integral equation for the irregularity surface, which, by asymptotic (kr ≫ 1) integration over the coordinate transverse to the propagation path (with allowance for terms of the order of (kr)⁻¹), is reduced to a one-dimensional integral equation, in which the integration contour is the linear contour of the irregularity. The equation is solved numerically, combining the inversion of a Volterra integral operator and successive approximations. By reducing the computer times, this method enables one to study both small-scale and large-scale irregularities. The results of numerical simulation of radio wave propagation in the presence of a powerful three-dimensional ionospheric disturbance are presented as an example. State University, St. Petersburg, Russia. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 41, No. 5, pp. 588–604, May, 1998.     

3D quantum transport solver based on the perfectly matched layer and spectral element methods for the simulation of semiconductor nanodevices

A 3D quantum transport solver based on the spectral element method (SEM) and perfectly matched layer (PML) is introduced to solve the 3D Schrödinger equation with a tensor effective mass. In this solver, the influence of the environment is replaced with the artificial PML open boundary extended beyond the contact regions of the device. These contact regions are treated as waveguides with known incident waves from waveguide mode solutions. As the transmitted wave function is treated as a total wave, there is no need to decompose it into waveguide modes, thus significantly simplifying the problem in comparison with conventional open boundary conditions. The spectral element method leads to an exponentially improving accuracy with the increase in the polynomial order and sampling points. The PML region can be designed such that less than −100 dB outgoing waves are reflected by this artificial material. The computational efficiency of the SEM solver is demonstrated by comparing the numerical and analytical results from waveguide and plane-wave examples and its utility is illustrated by multiple-terminal devices and semiconductor nanotube devices.     

Focusing and waveguiding of Lamb waves in micro-fabricated piezoelectric phononic plates

This paper presents results on the numerical and experimental studies of focusing and waveguiding of the lowest anti-symmetric Lamb wave in micro-fabricated piezoelectric phononic plates. The phononic structure was based on an AT-cut quartz plate and consisted of a gradient-index phononic crystal (GRIN PC) lens and a linear phononic plate waveguide. The band structures of the square-latticed AT-cut quartz phononic crystal plates with different filling ratios were analyzed using the finite element method. The design of a GRIN PC plate lens which is attached with a linear phononic plate waveguide is proposed. In designing the waveguide, propagation modes in square-latticed PC plates with different waveguide widths were studied and the results were served for the experimental design. In the micro-fabrication, deep reactive ion etching (Deep-RIE) process with a laboratory-made etcher was utilized to fabricate both the GRIN PC plate lens and the linear phononic waveguide on an 80 μm thick AT-cut quartz plate. Interdigital transducers were fabricated directly on the quartz plate to generate the lowest anti-symmetric Lamb waves. A vibro-meter was used to detect the wave fields and the measured results on the focusing and waveguiding of the piezoelectric GRIN PC lens and waveguide are in good accordance with the numerical predictions. The results of this study may serve as a basis for developing an active micro plate lens and related devices.     

Semi-analytical finite element analysis of elastic waveguides subjected to axial loads

Predicting the influence of axial loads on the wave propagation in structures such as rails requires numerical analysis. Conventional three-dimensional finite element analysis has previously been applied to this problem. The process is tedious as it requires that a number of different length models be solved and that the user identify the computed modes of propagation. In this paper, the more specialised semi-analytical finite element method is extended to account for the effect of axial load. The semi-analytical finite element method includes the wave propagation as a complex exponential in the element formulation and therefore only a two-dimensional mesh of the cross-section of the waveguide is required. It was found that the stiffness matrix required to describe the effect of axial load is proportional to the mass matrix, which makes the extension to existing software trivial. The method was verified by application to an aluminium rod, where after phase and group velocities of propagating waves in a rail were computed to demonstrate the method.     

基于特征基函数的球面共形微带天线阵列分析

采用全波分析法对球面共形微带天线阵列进行了分析.相比体-面积分方程,采用球并矢格林函数的面积分方程法可以大幅减少未知量的数目,进而缓解计算机内存压力.微带天线阵列表面采用曲面三角形剖分,可较精确地模拟球面特性.首先,引入边界电荷以及半Rao-Wilton-Glisson基函数,成功实现了探针馈电,并采用镜像法解决了馈电边处线积分奇异问题.然后,采用特征基函数法降低了阻抗矩阵的阶数,并采取有效措施进一步节省内存和计算时间.最后,分析计算了不同尺寸的球面共形微带天线阵列的输入阻抗及远区场特性.与文献和仿真软件结果进行比较,证明了所提出的处理方法的正确性和有效性.     

Multiobjective muffler shape optimization with hybrid acoustics modeling

This paper considers the combined use of a hybrid numerical method for the modeling of acoustic mufflers and a genetic algorithm for multiobjective optimization. The hybrid numerical method provides accurate modeling of sound propagation in uniform waveguides with non-uniform obstructions. It is based on coupling a wave based modal solution in the uniform sections of the waveguide to a finite element solution in the non-uniform component. Finite element method provides flexible modeling of complicated geometries, varying material parameters, and boundary conditions, while the wave based solution leads to accurate treatment of non-reflecting boundaries and straightforward computation of the transmission loss (TL) of the muffler. The goal of optimization is to maximize TL at multiple frequency ranges simultaneously by adjusting chosen shape parameters of the muffler. This task is formulated as a multiobjective optimization problem with the objectives depending on the solution of the simulation model. NSGA-II genetic algorithm is used for solving the multiobjective optimization problem. Genetic algorithms can be easily combined with different simulation methods, and they are not sensitive to the smoothness properties of the objective functions. Numerical experiments demonstrate the accuracy and feasibility of the model-based optimization method in muffler design.     

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.