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

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

计算物理第23卷（2006）总目录

第1期计算中子增殖率的时间强迫碰撞方法……………………………王瑞宏邓力许海燕裴鹿成(1)Navier-Stokes方程的线性化微分求积法(英)……………………………………………李晨吴雄华(18)统一坐标系下多介质流体力学计算的一种方法…………………………………………………贾鹏彦(19)SRLW方程的多辛Fourier拟谱格式及其守恒律…………………孔令华曾文平刘儒勋孔令健(25)一类正交增强的阶谱六面体矢量单元构造……………………………………………李江海孙秦(32)SALE速度重映算法的改进……………………………………………………熊俊周海…     

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

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

弹性波传播模拟的Chebyshev谱元法

通过在每一个单元中采用谱展开近似,Chebyshev谱元法兼具了有限元处理边界及复杂结构的灵活性和谱方法的快速收敛特性,为弹性波传播的数值模拟提供了一种有效工具。从加权余量原理出发,详细阐述了Chebyshev谱元法用于弹性波传播模拟的基本理论及相应数学公式.给出了使用Chebyshev正交多项式展开得到的,存在等参变换时有关单元质量矩阵和单元刚度矩阵的精确积分公式。同时应用逐元技术极大地减少了内存和计算需求.最后,两个数值算例被用于验证这种谱元方法的高精度和强适应性。     

多块对接网格技术在电磁场散射问题中的应用

采用多块对接网格技术和时域有限体积法(FVTD)研究了典型目标和多体的电磁场散射问题.控制方程采用三维一般曲线坐标系下的时变麦克斯韦方程组.时间方向采用四步Runge-Kutta方法,空间离散采用基于通量雅可比矩阵特征结构的通矢量差分分裂,依赖变量采用MUSCL插值.时间方向的计算精度为2阶,空间方向的计算精度可达3阶.典型算例的雷达散射截面(RCS)的计算结果与理论解吻合很好.对于多体问题计算与文献结果相当一致,说明该算法具有对复杂拓扑结构外形(包括多体问题)进行数值模拟的能力.     

小波元方法在声子晶体周期结构中的应用

将双正交小波系统和谱元法的思想结合得到一般有界区域中的双正交小波元,将小波元的边界适应性推广到一阶微分的情形,通过匹配得到严格满足边界条件的小波基函数;基于小波元发展一种一维声子晶体能带计算方法.该方法利用声子晶体本身的结构特点,兼顾小波在数值分析中的优势和边界条件的满足,与周期小波法相比,具有更高的计算精度和计算效率.     

体心和面心立方晶体中高对称方向格波的计算

本文在单式体心和面心立方晶体的一些具有高对称性质的方向上，利用偏振矢量的正交完备性，找出偏振矢量的形式，采用矩阵运算的方法，解本征值方程，求出晶格波的频率ω（ｋ），所得结果与用其它方法计算的结果相同．     

非均匀散射层矢量辐射传输(VRT)方程高阶散射解的迭代法

将散射介质层在z轴方向划分成薄层，用薄层的一阶散射强度、Fourier变换和迭代方法求解散射介质整层的矢量辐射传输(VRT)方程的高阶散射解.该方法将一阶散射与高阶散射迭代结合起来，计算公式简明，可计算高阶迭代解，计算时间少.计算结果与一层均匀散射介质的VRT方程一阶Mueller矩阵解、半空间均匀散射介质二阶Mueller矩阵解、以及离散坐标-特征值特征矢量法的VRT热辐射的数值解作了全面的比较.提出并讨论了非均匀散射层主动与被动VRT方程的高阶解.本计算程序可以通用于非球形粒子多层结构及非均匀介质的散射和热辐射计算. 关键词： VRT方程 分层 迭代解     

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

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

非协调数值流形方法的稳定性和收敛性分析

在流形元的基础上，提出了非协调数值流形方法，非协调数值流形方法的优点是在不增加广义节点自由度的前提下，大大提高数值流形方法的计算精度和计算效率.利用内部自由度静力凝聚处理，推导了消除内参后的单元应变矩阵和单元刚度矩阵.在Hilbert空间内，从最小势能原理出发对非协调数值流形方法的稳定性和收敛性进行了分析和讨论，得到了保证非协调流形元解唯一存在和收敛的基本条件，完善了非协调数值流形方法的理论基础.数值试验表明，新单元构造过程简单，有较高的精度，从而证明了本方法的可行性.     

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

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

一种二阶混合有限体元格式的GAMG预条件子

针对一种含跳系数椭圆问题的二阶混合有限体元格式,讨论求解相应离散系统PGMRES法的预条件子构造问题.通过严格的理论分析,建立分层基下该二阶混合有限体元刚度矩阵和二次有限元刚度矩阵的谱等价关系,并利用关于二次有限元刚度矩阵的一种基于分层思想的GAMG预条件子,为二阶混合有限体元刚度矩阵设计一种高效GAMG预条件子.数值结果验证理论分析的正确性和新预条件子的高效性与稳定性.     

Coupling transfer matrix method to finite element method for analyzing the acoustics of complex hollow body networks

This paper exposes a procedure to couple multiport transfer matrices to finite elements for analyzing the acoustics of automotive hollow body networks with a minimum of memory requirements and computational time. Generally, hollow body networks are made up from a series of elongated fluid partitions similar to ducts or waveguides. These fluid partitions generally contain complex elements: junctions, noise control elements, and cavities. The location and type of these elements in the network, mainly the noise control elements (e.g., sealing parts), may impact the noise inside a car. In the proposed hybrid method, the elongated fluid partitions are modeled with fluid finite elements. All complexities are modeled with two-port or multiport transfer matrices. The coupling of these matrices to finite elements is naturally done at the weak integral formulation stage of the acoustical problem. The coupling does not add any degrees of freedom to, nor modify, the original finite element matrix system. Consequently, changing locations and types of noise control elements in the hollow body network is fast and does not require rebuilding the finite element system. This enables optimizing the acoustics of a complex network on a desktop computer. The hybrid method is compared to experimental results on a tee-shaped hollow body networks. Good correlations are obtained.     

Vector parameterization of the Lorentz group transformations and polar decomposition of Mueller matrices

It is shown that the representation of the coherence matrix (the polarization density matrix) of beams of electromagnetic waves as a biquaternion corresponding to the four-vector of a pseudo-Euclidean space whose components are the intensity and the Stokes parameters provides a possibility of introducing the group transformations of these quantities isomorphic to SO(3.1) group. These transformations are a subset of the set of Mueller polarization matrices which, generally speaking, form a semigroup. The reduction of the semigroup of Mueller matrices to the group of transformations opens the possibility to use the vector parameterization of SO(3.1) group for interpretation of the polar decomposition of Mueller matrices. In particular, in this approach, the elements of the Mueller matrices corresponding to phase elements and polarizers turn out to be most simply and naturally related to their eigenpolarizations.     

Estimation of the loss factor of viscoelastic laminated panels from finite element analysis

A wave finite element (WFE) method is applied for predicting wave dispersion, wave attenuation and dissipation in viscoelastic laminated panels. The method involves postprocessing (using periodic structure theory) of element matrices of a small segment of the structure, which is modelled using a stack of three-dimensional finite elements meshed through the cross-section. Each layer can be discretised using either one solid element or more solid elements in order to more accurately represent interlaminar stress and strain. The finite element model of the segment of the structure is typically very small, resulting in very small computation cost. Formulations for the evaluation of the global loss factor using the WFE approach are given. In particular a formulation to calculate the average loss factor in the general case of an anisotropic component is proposed. Numerical examples are then shown. These concern the evaluation of the dispersion curves and of the global loss factor for damped laminated panels of different constructions.     

Improved eigenvalues for combined dynamical systems using a modified finite element discretization scheme

New approaches are presented to discretize an arbitrarily supported linear structure carrying various lumped attachments. Specifically, the exact eigendata, i.e., the exact natural frequencies and mode shapes, of the linear structure without the lumped attachments are first used to modify its finite element mass and stiffness matrix so that the eigensolutions of the discretized system coincide with the exact modes of vibration. This is achieved by identifying a set of minimum changes in the finite element system matrices and enforcing certain constraint conditions. Once the updated matrices for the linear structure are found, the finite element assembling technique is then used to include the lumped attachments by adding their parameters to the appropriate elements in the modified mass and stiffness matrices. Numerical experiments show that for the same number of elements, the proposed scheme returns higher natural frequencies that are substantially more accurate than those given by the finite element model. Alternatively, the proposed discretization scheme allows one to efficiently and accurately determine the higher natural frequencies of a combined system without increasing the number of elements in the finite element model.     

Efficient gradient estimation using finite differencing and likelihood ratios for kinetic Monte Carlo simulations

While many optimization and control methods for stochastic processes require gradient information from the process of interest, obtaining gradient information from experiments is prohibitively expensive and time-consuming. As a result, such information is often obtained from stochastic process simulations. Computing gradients efficiently and accurately from stochastic simulations is challenging, especially for simulations involving computationally expensive models with significant inherent noise. In this work, we analyze and characterize the applicability of two gradient estimation methods for kinetic Monte Carlo simulations: finite differencing and likelihood ratio. We developed a systematic method for choosing an optimal perturbation size for finite differencing and discuss, for both methods, important implementation issues such as scaling with respect to the number of elements in the gradient vector. Through a series of numerical experiments, the methods were compared across different time and size regimes to characterize the precision and accuracy associated with each method. We determined that the likelihood ratio method is appropriate for estimating gradients at short (transient) times or for systems with small population sizes, whereas finite differencing is better-suited for gradient estimation at long times (steady state) or for systems with large population sizes.     

On an Integrated Transfer Matrix method for multiply connected mufflers

The commercial automotive mufflers are generally of a complicated shape with multiply connected parts and complex acoustic elements. The analysis of such complex mufflers has always been a great challenge. In this paper, an Integrated Transfer Matrix method has been developed to analyze complex mufflers. Integrated transfer matrix relates the state variables across the entire cross-section of the muffler shell, as one moves along the axis of the muffler, and can be partitioned appropriately in order to relate the state variables of different tubes constituting the cross-section. The paper presents a generalized one-dimensional (1-D) approach, using the transfer matrices of simple acoustic elements, which are available from the literature. The present approach is robust and flexible owing to its capability to construct an overall matrix of the muffler with the transfer matrices of individual acoustic elements and boundary conditions, which can then be used to evaluate the transmission loss, insertion loss, etc. Results from the present approach have been validated through comparisons with the available experimental and three-dimensional finite element method (FEM) based results. The results show good agreement with both measurements and FEM analysis up to the cut-off frequency.     

随机介质背散射二维Mueller矩阵的数值模拟与分析

推导了随机介质背散射Mueller矩阵的直接计算公式,并运用矢量Monte Carlo方法进行了数值模拟.结果表明随机介质背散射二维Mueller矩阵方位关系随散射系数的减小而增强,而与微粒大小关系不大;Mueller矩阵元素绝对值的空间分布随径向呈近似指数规律衰减,矩阵元素的方位变化具有周期性.对称系统的二维Mueller矩阵的花样图中仅有7幅独立,其余9幅可通过对称、旋转变换得到.