各向异性介质波导不连续问题的半解析谱单元法分析

将谱单元法与精细积分法相结合求解各向异性介质的波导不连续问题.从矢量波动方程的单变量变分形式出发,采用基于Gauss-Lobatto-Legendre多项式零点作为插值结点的谱单元,对含有各向异性介质波导的横截面进行离散,然后将问题导入哈密顿体系利用精细积分法进行求解.由于采用了谱单元法,在单元网格数较少时,可获得高准确度的计算结果;又由于利用了精细积分法,结构的纵向长度可以任意设定,克服了当人工边界设置在离介质块较远处时,计算量不断增加的缺点.研究表明半解析谱单元法可有效地应用于各向异性介质的波导不连续问题,在提高准确度的同时可大量节省计算时间.     

标量波传播问题的双互易精细积分法

为避免使用计算多种特征频率下的声场响应，采用双互易方法将边界积分方程中时间二次导数项的域积分转化为边界积分.首先，将计算场点配置在边界上并考虑边界条件，可以获得由内部节点上声压量线性表示的边界节点上的物理量；其次，将计算场点配置于域内离散节点上，将所得边界积分方程组中关于边界物理量用内部节点的声压量线性表示，获得关于声压量的二阶常微分方程组；第三，引入声压变化速度作为未知量，将二阶常微分方程组转化为一阶常微分方程组；最后，采用精细积分法精确求解常微分方程组.数值算例验证了双互易精细积分法的正确性和稳定性.     

基于不连续因子的沸水堆芯三维瞬态数值模拟

基于不连续因子校正的粗网格有限差分法是实现堆芯瞬态三维数值模拟的高效方法之一, 粗网节块的界面不连续因子与边界反照率的计算方法决定了实时数值模拟过程中的精度。在计算不连续因子的过程中, 省去了细网节块计算与粗网均匀化过程, 直接在粗网格划分情况下, 基于节块展开法和非线性迭代策略, 推导了粗网格界面不连续因子比率与边界反照率的计算公式, 并编制了相应的计算程序。沸水堆典型算例的三维瞬态模拟证实该方法可在空间域和时间域两方面, 使静态、瞬态精度均达到与先进节块法相等同的程度, 并且计算效率优于先进节块法, 为核电站全范围模拟机三维堆芯的实时仿真模型开发提供了一种切实可行的选择。     

一种有限元-光滑粒子流体动力学耦合算法

介绍一种新的有限元-光滑粒子流体动力学耦合算法.通过将SPH法中的粒子固连到单元边界上实现有限元计算与SPH计算的耦合,粒子可以连接到单元边界上的任意点,且一条单元边界可以与多个粒子连接,耦合界面处粒子分布与单元分布无相互限制.应用该算法模拟应力波传播和高速冲击问题.结果表明,算法在应力波传播计算中具有良好的精度,在高速冲击问题中计算精度与SPH法相近的,而计算效率大大提高.     

基于Burton-Miller方程的轴对称结构声学边界元方法

提出了一种求取轴对称结构任意边界条件下声辐射特性的边界元方法。采用Burton和Miller改进型公式将高阶奇异项转化为弱奇异项之和,保证声辐射参数的唯一性,且计算简单精确。将结构表面声压与振速按照旋转轴角度进行Fourier级数展开,利用级数的正交性建立各项待定系数的求解公式;然后转化格林函数的法向偏导为切向偏导,方便直接计算各项积分,并将面积分公式表示为沿结构边界的线积分和沿旋转角度的积分;进一步采用二次等参单元离散结构边界线,建立声压与振速的关系矩阵,从而确定结构声辐射参数。以脉动球源和横向振动球源为例计算,与解析解和传统边界元法结果作对比,说明该方法的有效精确性。     

三维静电场线性插值边界元中的解析积分方法

提出求解三维静电场的三角形线性插值边界元解析积分方法.针对含1/R和1/R²的积分项,将单元形状函数分解为常数项、含x的线性项和含y的线性项,从而将边界单元积分简化为6个基本积分组合,并导出其解析计算公式,避免了因形状函数改变而导致的重复计算.该方法不仅可以准确计算远离奇异情况下的边界元积分,而且可以准确计算一阶和二阶接近奇异积分以及一阶奇异积分.计算结果表明,在接近奇异积分和奇异积分比较突出的问题中,当数值积分方法不能给出正确结果时,用同样的边界元网格,解析积分方法可以给出正确的结果,提高了三维静电场线性插值边界元法的计算精度.     

狭缝式裂变探测系统的中子灵敏度解析计算方法

 利用基于解析公式的理论方法来计算狭缝式裂变探测系统的中子灵敏度。采用积分和泰勒展开法推导了铀裂变材料在中子作用下产生的碎片数目及PIN探测器接收碎片的概率,并得到该系统的中子灵敏度理论计算公式。在D-T中子源上对该探测系统的中子灵敏度进行了实验标定,并将实验结果与解析公式计算得到的理论值进行比较。结果表明：理论计算值与实验结果保持一致,两者之间的偏差小于10%。     

接近式光刻衍射光场的有效波前作用规律分析

掩模图形不同区域的波前因衍射光的相干叠加作用而对光刻胶上的场点产生不同的影响。为研究这一作用规律,采用波前分割的方法对掩模图形上的波前进行区域划分,并利用光场相干叠加相互抵消的性质,最终得到了掩模图形上对场点光场影响最大的波前区域。理论分析表明,对于图形内部场点而言,这个区域范围为场点附近约一个半波带的大小;而对靠近图形边缘的场点,其范围与边缘形状有关,一般要稍大于一个半波带。利用该方法可以显著提高仿真效率,对于分析由衍射造成的光刻误差来源和大小具有一定的理论指导意义。该理论计算结果得到了实验验证。     

一类可变禁区的不连续系统的加周期分岔

研究了一类可变禁区不连续系统的加周期分岔行为, 发现由可变禁区导致不同类型的加周期分岔. 研究表明, 系统的迭代轨道和禁区的上下两个边界均可发生边界碰撞, 从而产生加周期分岔. 基于边界碰撞分岔理论, 定义基本的迭代单元, 解析推导出了相应的分岔曲线, 在全参数空间中给出了不同加周期所出现的范围. 与数值模拟结果比较, 理论分析结果与数值结果高度一致.     

基于方向图变换的快速不连续相位展开

几何突变导致的相位不连续一直是相位展开中具有挑战的问题。针对这一问题,提出基于方向与变换的快速不连续相位展开算法。在提出的算法中,利用图像的结构张量估计包裹相位图的方向,将方向图进行变换和差分计算得到的可靠的权重系数图作为加权最小二乘法的权重,并使用预处理共轭梯度法迭代求解。该算法可以快速地找出不连续位置,并在不连续分割后的区域进行单独的相位展开,具有很好的识别和展开效果。详细地描述了算法的原理和实现步骤,并对算法进行仿真和实验数据验证。实验结果表明:其相位展开的均方根误差为0.36,证明该算法能够快速、准确地对不连续的包裹相位进行展开。     

Singular meshless method using double layer potentials for exterior acoustics

Time-harmonic exterior acoustic problems are solved by using a singular meshless method in this paper. It is well known that the source points cannot be located on the real boundary, when the method of fundamental solutions (MFS) is used due to the singularity of the adopted kernel functions. Hence, if the source points are right on the boundary the diagonal terms of the influence matrices cannot be derived. Herein we present an approach to obtain the diagonal terms of the influence matrices of the MFS for the numerical treatment of exterior acoustics. By using the regularization technique to regularize the singularity and hypersingularity of the proposed kernel functions, the source points can be located on the real boundary and therefore the diagonal terms of influence matrices are determined. We also maintain the prominent features of the MFS, that it is free from mesh, singularity, and numerical integration. The normal derivative of the fundamental solution of the Helmholtz equation is composed of a two-point function, which is one of the radial basis functions. The solution of the problem is expressed in terms of a double-layer potential representation on the physical boundary based on the potential theory. The solutions of three selected examples are used to compare with the results of the exact solution, conventional MFS, boundary element method, and Dirichlet-to-Neumann finite element method. Good numerical performance is demonstrated by close agreement with other solutions.     

A boundary integral equation method using auxiliary interior surface approach for acoustic radiation and scattering in two dimensions

This paper presents an effective solution method for predicting acoustic radiation and scattering fields in two dimensions. The difficulty of the fictitious characteristic frequency is overcome by incorporating an auxiliary interior surface that satisfies certain boundary condition into the body surface. This process gives rise to a set of uniquely solvable boundary integral equations. Distributing monopoles with unknown strengths over the body and interior surfaces yields the simple source formulation. The modified boundary integral equations are further transformed to ordinary ones that contain nonsingular kernels only. This implementation allows direct application of standard quadrature formulas over the entire integration domain; that is, the collocation points are exactly the positions at which the integration points are located. Selecting the interior surface is an easy task. Moreover, only a few corresponding interior nodal points are sufficient for the computation. Numerical calculations consist of the acoustic radiation and scattering by acoustically hard elliptic and rectangular cylinders. Comparisons with analytical solutions are made. Numerical results demonstrate the efficiency and accuracy of the current solution method.     

Acoustic scattering by multiple elliptical cylinders using collocation multipole method

This paper presents the collocation multipole method for the acoustic scattering induced by multiple elliptical cylinders subjected to an incident plane sound wave. To satisfy the Helmholtz equation in the elliptical coordinate system, the scattered acoustic field is formulated in terms of angular and radial Mathieu functions which also satisfy the radiation condition at infinity. The sound-soft or sound-hard boundary condition is satisfied by uniformly collocating points on the boundaries. For the sound-hard or Neumann conditions, the normal derivative of the acoustic pressure is determined by using the appropriate directional derivative without requiring the addition theorem of Mathieu functions. By truncating the multipole expansion, a finite linear algebraic system is derived and the scattered field can then be determined according to the given incident acoustic wave. Once the total field is calculated as the sum of the incident field and the scattered field, the near field acoustic pressure along the scatterers and the far field scattering pattern can be determined. For the acoustic scattering of one elliptical cylinder, the proposed results match well with the analytical solutions. The proposed scattered fields induced by two and three elliptical–cylindrical scatterers are critically compared with those provided by the boundary element method to validate the present method. Finally, the effects of the convexity of an elliptical scatterer, the separation between scatterers and the incident wave number and angle on the acoustic scattering are investigated.     

A new method for true and spurious eigensolutions of arbitrary cavities using the combined Helmholtz exterior integral equation formulation method

Integral equation methods have been widely used to solve interior eigenproblems and exterior acoustic problems (radiation and scattering). It was recently found that the real-part boundary element method (BEM) for the interior problem results in spurious eigensolutions if the singular (UT) or the hypersingular (LM) equation is used alone. The real-part BEM results in spurious solutions for interior problems in a similar way that the singular integral equation (UT method) results in fictitious solutions for the exterior problem. To solve this problem, a Combined Helmholtz Exterior integral Equation Formulation method (CHEEF) is proposed. Based on the CHEEF method, the spurious solutions can be filtered out if additional constraints from the exterior points are chosen carefully. Finally, two examples for the eigensolutions of circular and rectangular cavities are considered. The optimum numbers and proper positions for selecting the points in the exterior domain are analytically studied. Also, numerical experiments were designed to verify the analytical results. It is worth pointing out that the nodal line of radiation mode of a circle can be rotated due to symmetry, while the nodal line of the rectangular is on a fixed position.     

Application of the Wiener–Hopf method for describing the propagation of sound in cylindrical and rectangular channels with an impedance jump in the presence of a flow

Exact solutions to problems of the propagation of acoustic modes in lined channels with an impedance jump in the presence of a uniform flow are constructed. Two problems that can be solved by the Wiener- Hopf method—the propagation of acoustic modes in an infinite cylindrical channel with a transverse impedance jump and the propagation of acoustic modes in a rectangular channel with an impedance jump on one of its walls—are considered. On the channel walls, the Ingard–Myers boundary conditions are imposed and, as an additional boundary condition in the vicinity of the junction of the linings, the condition expressing the finiteness of the acoustic energy. Analytical expressions for the amplitudes of the transmitted and reflected fields are obtained.     

A fast direct solver for the integral equations of scattering theory on planar curves with corners

We describe an approach to the numerical solution of the integral equations of scattering theory on planar curves with corners. It is rather comprehensive in that it applies to a wide variety of boundary value problems; here, we treat the Neumann and Dirichlet problems as well as the boundary value problem arising from acoustic scattering at the interface of two fluids. It achieves high accuracy, is applicable to large-scale problems and, perhaps most importantly, does not require asymptotic estimates for solutions. Instead, the singularities of solutions are resolved numerically. The approach is efficient, however, only in the low- and mid-frequency regimes. Once the scatterer becomes more than several hundred wavelengths in size, the performance of the algorithm of this paper deteriorates significantly. We illustrate our method with several numerical experiments, including the solution of a Neumann problem for the Helmholtz equation given on a domain with nearly 10000 corner points.     

半空间球面波函数叠加法重构声源直接辐射声场

提出了基于半空间球面波函数叠加的声场重构方法,以重构含有限声阻抗边界半空间中声源直接辐射的声场。在半空间中多极子声源声压场的解析解的基础上,构造出以边界声阻抗为参量的半空间球面波函数的正交基;通过求逆获得半空间总声压解的基函数系数,同时也获得声源直接辐射声场即自由空间中的基函数系数,进而重构出声源直接辐射的声场。在边界声阻抗已知和边界声阻抗未知两种条件下,对该方法进行了仿真验证和参数分析,并在全消声室内进行了实验验证。结果表明,所提方法能重构出半空间中典型声源即球形声源和平面声源的直接辐射声场;该方法在边界声阻抗已知时的重构精度与稳定性高于在边界声阻抗未知时的情形。      

轴对称体声振耦合的边界子波谱与有限元耦合方法

探讨了子波在Helmholtz积分方程及声振耦合中的应用,在建立了求解轴对称Helmholtz积分方程的子波谱方法的基础上,构造了轴对称子波谱与轴对称有限元的耦合方法,该方法可以处理轴对称问题的任意边界条件.进行了声振耦合问题的模态分析.