基于高阶叠层矢量基函数的E-H时域有限元方法分析谐振腔和波导结构

将高阶叠层矢量基函数用于E-H时域有限元方法,电场和磁场用相同的基函数展开并同时求解,时间离散采用Crank-Nicolson差分格式使得时间步长的选取摆脱稳定性条件的限制,同时采用完美匹配层来截断计算区域.对三维谐振腔及波导结构进行数值模拟与分析,结果表明,相较于低阶基函数,高阶叠层矢量基函数可以有效提高E-H时域有限元方法的计算精度.     

矢量有限元方法精确计算轴对称?谐振腔高阶模

采用矢量有限元方法计算了任意轴对称谐振腔高阶模的本征频率.这种方法能够将数学模型和物理理论很好地结合起来,避免了一般有限元方法计算谐振腔本征模式中可能出现伪根的问题,并且可以通过采用二阶标量和矢量基函数,在较少的网格数的情况下得到很高的计算精度.基于该方法的程序Cafe为进一步精确研究下一代直线对撞机失谐结构的尾场奠定了基础.     

基于曲四面体单元剖分的CN E-H TDFEM方法

从采用Crank-Nicolson差分格式的基于麦克斯韦旋度方程的E-H时域有限元方法出发,将展开电场和磁场的叠层矢量基函数与曲四面体单元相结合.对金属球谐振腔及介质填充的圆柱形谐振腔的数值模拟表明：相较于规则四面体单元,在剖分单元数目相同的情况下,曲四面体单元离散表面弯曲结构可以获得更高的计算精度.同时,与0.5阶基函数结合曲四面体单元相比,1.0阶基函数与曲四面体单元结合可以用更少的单元数及未知量数目来获得更高的计算精度.     

曲边六面体矢量单元的改进与应用

改进了曲边六面体矢量单元的矢量基函数.简要回顾了原六面体矢量单元的定义之后,介绍了新的矢量基函数的构造方法,矢量基函数间的正交性因此有所改善,并有利于处理有限元Dirichlet边界条件.结合完全匹配层,将新单元应用于三维电磁散射问题,并进行了数值计算.     

三阶矢量有限元方法精确计算轴对称谐振腔高阶模

为了进一步提高数值求解谐振腔高阶模的精度，本文提出了三阶矢量有限元方法，并针对二阶矢量有限元轴对称谐振腔高阶模计算程序Cafe对曲线边界的计算能力较差和计算速度较慢的缺点做了改进. 在这些改进的基础上编制了三阶矢量有限元轴对称谐振腔高阶模计算程序meshmatrix3，得到了很好的结果.     

高阶辛算法的稳定性与数值色散性分析

利用Maxwell方程的哈密尔顿函数,导出对应的欧拉-哈密尔顿方程.利用辛积分技术与高阶交错差分技术,建立求解三维时域Maxwell方程的高阶辛算法;结合电磁场中的物理概念,借助矩阵分析和张量分析理论,获得高阶时域方法及高阶辛算法的稳定性和数值色散性的统一处理新方法.用数值结果证实方法的正确性,与FDTD算法和其它时域高阶方法相比,高阶辛算法具有较大的计算优势,为电磁计算提供了新的途径.     

求解Kuramoto-Sivashinsky方程的平移基无单元Galerkin方法

Kuramoto-Sivashinsky方程是一种可以描述复杂混沌现象的高阶非线性演化方程.方程中高阶导数项的存在,使得传统无单元Galerkin方法采用高次多项式基函数构造形函数时,形函数违背了一致性条件.因此,本文提出了一种采用平移多项式基函数的无单元Galerkin方法.与传统无单元Galerkin方法相比,该方法在方程离散时依然采用Galerkin进行离散,但形函数的构造采用了基于平移多项式基函数的移动最小二乘近似.通过对具有行波解和混沌现象的Kuramoto-Sivashinsky方程的数值模拟,验证了本文方法的有效性.     

并行时域电场积分方程方法在动态海面二维电磁散射中的应用

采用并行时域电场积分方程方法对动态海面的二维瞬态散射特性进行研究。为了保证该方法的后期稳定性,时间基函数和空间基函数采用二阶B样条基函数和三角基函数,矩阵元素采用时间维度精确解析、空间奇异部分精确解析进行计算;为了减少对无限海面进行截断带来的边缘效应,入射波采用锥形调制高斯脉冲;结合信息传递接口(MPI)技术和稀疏矩阵压缩存储技术,对不同时刻的海面进行瞬态散射分析。大量的数值算例证明了该方法在计算动态海面的二维瞬态散射问题时的正确性,还可以保证后期的稳定性,提高计算效率,减少对计算机内存需求。     

一类正交增强的阶谱六面体矢量单元构造

利用阶谱正交多项式的优越性,从单元系统矩阵形成的角度分析阶谱六面体矢量单元的基函数,提出具有正交增强的阶谱六面体矢量单元(ORHHVFE)构造目的;并利用该目的构造的阶谱六面体矢量单元与其它阶谱六面体矢量单元(HHVFE)进行了金属腔本征模问题的数值计算对比实验.结果表明ORHHVFE具有与其它HHVFE同等的数值计算精度,且由ORHHVFE形成的有限元系统矩阵条件数获得了极大的改善.     

基于辛Runge-Kutta-Nystrom方法的雷达散射截面计算

从基本的差分概念和Maxwell方程出发, 引入电磁场方程的Hamilton函数.提出一种基于Runge-Kutta-Nystrom辛算法的高阶时域有限差分方法，该方法保持了系统的相空间体积不变和总能量不变，并导出了迭代公式.在此基础上计算了一种金属圆柱的雷达散射截面.计算结果表明该方法的正确性及快速、精确的特性. 关键词： 雷达散射截面 高阶算法 辛Runge-Kutta-Nystrom方法 时域有限差分     

高阶CIP数值方法及其在相关物理问题中的应用

利用函数的高阶空间导数值构建其高次插值,得到高阶CIP(Constrained Interpolation Profile)数值算法,并在此基础上模拟研究等离子体物理中著名的伏拉索夫-泊松(Vlasov-Poisson)方程相关物理问题.高阶CIP数值方法具有更高数值精度,从而可以在同等精度的情况下减少计算格点数,加速数值计算速度.     

��ά�߽�CIP��ֵ��������Ӧ��

利用函数及其高阶导数值构造五次插值函数近似网格单元内的真实解,改进数值求解双曲类偏微分方程的CIP数值算法。基于之前的一维高阶CIP数值算法思想,不同于利用时间分裂技术,发展了二维高阶CIP数值算法。改进后的算法具有五阶数值精度和显示格式的优点。     

A new volume-of-fluid method with a constructed distance function on general structured grids

A second-order volume-of-fluid method (VOF) is presented for interface tracking and sharp interface treatment on general structured grids. Central to the new method is a second-order distance function construction scheme on a general structured grid based on the reconstructed interface. A novel technique is developed for evaluating the interface normal vector using the distance function. With the normal vector, the interface is reconstructed from the volume fraction function via a piecewise linear interface calculation (PLIC) scheme on the computational domain. Several numerical tests are conducted to demonstrate the accuracy and efficiency of the present method. In general, the new VOF method is more efficient than both the high-order level set and the coupled level set and volume-of-fluid (CLSVOF) methods. The results from the new method are better than those from the benchmark VOF method, particularly in the under-resolved regions, and are comparable to those from the CLSVOF method. Breaking waves over a submerged bump and around a wedge-shaped bow are simulated to demonstrate the application of the new method and sharp interface treatment in a two-phase flow solver on curvilinear grids. The computational results are in good agreement with the available experimental measurements.     

Extension of the spectral volume method to high-order boundary representation

In this paper, the spectral volume method is extended to the two-dimensional Euler equations with curved boundaries. It is well-known that high-order methods can achieve higher accuracy on coarser meshes than low-order methods. In order to realize the advantage of the high-order spectral volume method over the low order finite volume method, it is critical that solid wall boundaries be represented with high-order polynomials compatible with the order of the interpolation for the state variables. Otherwise, numerical errors generated by the low-order boundary representation may overwhelm any potential accuracy gains offered by high-order methods. Therefore, more general types of spectral volumes (or elements) with curved edges are used near solid walls to approximate the boundaries with high fidelity. The importance of this high-order boundary representation is demonstrated with several well-know inviscid flow test cases, and through comparisons with a second-order finite volume method.     

Efficient quadrature-free high-order spectral volume method on unstructured grids: Theory and 2D implementation

An efficient implementation of the high-order spectral volume (SV) method is presented for multi-dimensional conservation laws on unstructured grids. In the SV method, each simplex cell is called a spectral volume (SV), and the SV is further subdivided into polygonal (2D), or polyhedral (3D) control volumes (CVs) to support high-order data reconstructions. In the traditional implementation, Gauss quadrature formulas are used to approximate the flux integrals on all faces. In the new approach, a nodal set is selected and used to reconstruct a high-order polynomial approximation for the flux vector, and then the flux integrals on the internal faces are computed analytically, without the need for Gauss quadrature formulas. This gives a significant advantage over the traditional SV method in efficiency and ease of implementation. For SV interfaces, a quadrature-free approach is compared with the Gauss quadrature approach to further evaluate the accuracy and efficiency. A simplified treatment of curved boundaries is also presented that avoids the need to store a separate reconstruction for each boundary cell. Fundamental properties of the new SV implementation are studied and high-order accuracy is demonstrated for linear and non-linear advection equations, and the Euler equations. Several well known inviscid flow test cases are utilized to show the effectiveness of the simplified curved boundary representation.     

Effects of high-order dispersions on dark-bright vector soliton propagation and interaction

The dynamics of dark-bright vector solitons is investigated in a birefringent fiber with the high-order dispersions,and their effects on vector soliton propagation and interaction are analyzed using the numerical method.The combined role of the high-order dispersions,such as the third-order dispersion (TOD) and the fourth-order dispersion (FOD),may cause various deformation of the vector soliton and enhance interaction.These effects depend strictly on the sign of the high-order dispersions.Results indicate that the disadvantageous effects can be reduced effectively via proper mapping of the high-order dispersions.     

基于递归T矩阵的离散随机散射体散射特性研究

本文根据电磁场矢量球波函数多极点展开原理及矢量叠加定理提出了递归T矩阵算法的矢量形式,并且基于矢量递归T矩阵算法建立了多散射球模拟离散随机散射体散射的三维电磁散射模型.通过计算不同尺寸、随机分布散射球的散射以及分析散射球间的高阶散射效应,结果表明:矢量递归T矩阵算法具有很高的计算精度,算法中包含多散射体间的高阶散射效应,因此能够精确计算多散射体总的散射效应.本文所建模型可应用于土壤湿度探测工程中评估地表下掩埋离散随机散射体散射对雷达回波信号产生的影响.     

A high-order gas-kinetic Navier–Stokes flow solver

The foundation for the development of modern compressible flow solver is based on the Riemann solution of the inviscid Euler equations. The high-order schemes are basically related to high-order spatial interpolation or reconstruction. In order to overcome the low-order wave interaction mechanism due to the Riemann solution, the temporal accuracy of the scheme can be improved through the Runge–Kutta method, where the dynamic deficiencies in the first-order Riemann solution is alleviated through the sub-step spatial reconstruction in the Runge–Kutta process. The close coupling between the spatial and temporal evolution in the original nonlinear governing equations seems weakened due to its spatial and temporal decoupling. Many recently developed high-order methods require a Navier–Stokes flux function under piece-wise discontinuous high-order initial reconstruction. However, the piece-wise discontinuous initial data and the hyperbolic-parabolic nature of the Navier–Stokes equations seem inconsistent mathematically, such as the divergence of the viscous and heat conducting terms due to initial discontinuity. In this paper, based on the Boltzmann equation, we are going to present a time-dependent flux function from a high-order discontinuous reconstruction. The theoretical basis for such an approach is due to the fact that the Boltzmann equation has no specific requirement on the smoothness of the initial data and the kinetic equation has the mechanism to construct a dissipative wave structure starting from an initially discontinuous flow condition on a time scale being larger than the particle collision time. The current high-order flux evaluation method is an extension of the second-order gas-kinetic BGK scheme for the Navier–Stokes equations (BGK-NS). The novelty for the easy extension from a second-order to a higher order is due to the simple particle transport and collision mechanism on the microscopic level. This paper will present a hierarchy to construct such a high-order method. The necessity to couple spatial and temporal evolution nonlinearly in the flux evaluation can be clearly observed through the numerical performance of the scheme for the viscous flow computations.     

3-D High-Order Edge-Element Analysis of 3-D Discontinuity Problems in Guiding Wave Structures

A 54-parameter high-order 3-D edge-element construction based on full magnetic vector variational formulas is discussed. This approach eliminates the spurious solutions efficiently with high accuracy. The analyses of the scattering characteristics of some 3-D discontinuity problems in guided wave structures verify the effective and reliability of the present method. The comparison between the numerical results obtained with 54 and 12-parameter edge element method shows that the present approach is more accurate and efficient.