腔光子-自旋波量子耦合系统中各向异性奇异点的实验研究

通过耦合三维微波腔中光子和腔内钇铁石榴石单晶小球中的自旋波量子形成腔-自旋波量子的耦合系统,并通过精确调节系统参数在该实验系统中观测到各向异性奇异点.奇异点对应于非厄米系统中一种特殊状态,在奇异点处,耦合系统的本征值和本征矢均简并,并且往往伴随着非平庸的物理性质.以往大量研究主要集中在各向同性奇异点的范畴,它的特征是在系统参数空间中沿着不同参数坐标趋近该奇异点时具有相同的函数关系.在这篇文章中,主要介绍实验上在腔光子-自旋波量子耦合系统中通过调节系统的耦合强度和腔的耗散衰减系数两条趋近奇异点的路径而实现了各向异性奇异点,具体分别对应于在趋近奇异点时,本征值的虚部的变化与耦合强度和腔的衰减系数的变化会有线性和平方根不同的行为.各向异性奇异点的实现有助于基于腔光子-自旋波量子耦合系统的量子信息处理和精密探测器件的进一步研究.     

极坐标系中的B样条有限元法在求解尖角奇异边值问题中的应用

基于尖角状角点附近的场分布特性，提出了应用极坐标系中的Ｂ样条有限元法求解尖角奇异边值问题。该方法具有精度高和易于实现等优点。通过对Ｍｏｔｚ问题、裂梁问题以及扇形波导本征值问题等尖角奇异边值问题的计算分析，验证方法的有效性。     

混沌背景中信号参数估计的新方法

怎样提取混沌背景中信号的参数具有重要意义.在重构的相空间中,叠加有其他信号的混沌信号时间序列重构的点集会偏离混沌吸引子所在光滑流形,依据这一性质并综合利用混沌背景中信号本身的特性,提出一种参数估计的新方法：最小相对奇异值(MRSV)法.该方法先建立逆滤波器,由其输出重构相空间,然后改变其参数,使输出信号在嵌入空间中作局部奇异值分解的相对奇异值最小,来实现参数估计.AR模型参数和正弦信号频率估计的仿真结果验证了该方法的有效性. 关键词： 混沌 参数估计 最小相对奇异值(MRSV) 逆滤波器     

三维位势问题边界元法中几乎奇异积分的正则化

采用一种半解析正则化算法,计算了三维位势问题边界元法中近边界点的几乎强奇异和几乎超奇异面积分.该算法适用于三角形线性等参元.对高次单元将其细分为几个三节点三角形单元即可应用该算法.由于几乎奇异性,与内点邻近的单元上的积分,采用半解析正则化积分算法计算;而远处单元的积分仍保持常规高斯积分.对三维热传导算例,计算了近边界点的温度和热流.数值结果证明了该算法的有效性和精确性.     

动力边界元法中的强奇异积分问题

本文讨论了动力边界元法中的奇异积分问题,对其中的强奇异积分提出了一个有效的计算方法.该方法从合非零初始态的边界积分方程出发,利用动力方程的特解间接地确定了主系数(即所谓强奇异积分),从而避免了直接计算强奇异积分的困难.根据该方法编制了计算程序,并给出了一个简单算例。     

基于局域均值滤波的Stoilov改进算法

Stoilov算法中相移量的计算依赖于采集到的变形条纹图的光强,存在对光强的减法、除法和开方等运算,使相位计算时在某些位置出现分子分母为零或开方为复数等奇异现象,导致相位展开无意义或出现超大误差。从而使三维面形重构出现畸变、失真,甚至无法进行三维重构。大量实验说明这些奇异点都是孤立的点,本文提出一种基于局域均值滤波的方法对Stoilov算法进行改进,首先对奇异点进行标记,然后依次用其周围8个点中的非奇异点的平均值代替该奇异点的值。实验结果表明,该方法在考虑传感系统带来的误差的情况下,有效抑制了奇异现象,提高了三维测量精度。     

两量子比特系统中相互作用对高阶奇异点的影响

近年来,与环境耦合的非厄米开放系统成为人们研究的热点.非厄米体系中的奇异点会发生本征值和本征态的聚合,是区分厄米体系的重要性质之一.在具有宇称-时间反演对称性的体系中,奇异点通常伴随着对称性的自发破缺,存在很多值得探究的新奇物理现象.以往的研究多关注无相互作用系统中的二阶奇异点,对具有相互作用的多粒子系统,及其中可能出现的高阶奇异点讨论较少,特别是相关的实验工作尚未见报道.本文研究了具有宇称-时间反演对称性的两量子比特体系,证明了该体系中存在三阶奇异点,并且量子比特间的伊辛型相互作用能够诱导体系在三阶奇异点附近出现能量的高阶响应,可通过测量特定量子态占据数随时间的演化拟合体系本征值的方法来验证.其次通过探究该体系本征态的性质,展示了奇异点的态聚合特征,并提出了利用长时间演化后稳态的密度矩阵验证态聚合的方法.此外,还将理论的两量子比特哈密顿量映射到两离子实验系统中,基于¹⁷¹Yb~+囚禁离子系统设计了实现和调控奇异点,进而验证三阶响应的实验方案.这一方案具有极高的可行性,并有望对利用非厄米系统实现精密测量和高灵敏度量子传感器提供新的思路.     

奇异点快速检测在牛奶成分近红外光谱测量中的应用

近红外光谱作为一种依靠模型对物化性质进行分析的技术,对光谱数据的准确性进行快速准确的判断是得到可靠分析结果的前提。但是光谱数据中奇异点的存在会在很大程度上影响多变量校正模型的准确性,从而影响模型的预测效果。文章综合利用半数重采样法(Resampling by Half-Mean,RHM)和最小半球体积法(Smallest Half-Volume,SHV)成功剔除了被测量的牛奶成分近红外光谱中的奇异点,其效果远优于传统的奇异点剔除方法,并且该方法具有简单快速、计算量小、数值稳定等特点,非常适用于在线分析和其他类型的光谱数据中奇异点的检测。     

热应力边界元法中几乎超奇异积分的计算

通过分部积分变换将热弹性力学应力边界积分方程中的超奇异积分转化为强奇异积分,然后与另一个强奇异积分求和,得到仅含几乎强奇异的热应力自然边界积分方程.再对其中的几乎强奇异积分施以正则化,消除了热弹性力学边界元法中的几乎奇异积分,可以准确计算出热弹性力学问题中近边界内点的热应力.算例证明了方法的有效性.     

FGMs稳态热传导分析的重心Lagrange插值配点法

提出一种求解二维功能梯度材料（FGMs）稳态热传导问题的重心Lagrange插值配点法.基于Chebyshev节点构造二维重心Lagrange插值函数及其偏导数,然后基于配点法将其直接代入FGMs热传导问题的控制方程和边界条件,得到系统离散方程.重心Lagrange插值配点法是一种真正的无网格方法,很好地融合了重心Lagrange插值和配点格式的优势,具有高效、稳定、高精度和易于数值实现的优点.采用重心Lagrange插值配点法分别对指数型、二次型和三角型FGMs热传导问题进行数值模拟.结果表明：该方法具有较高的计算效率和计算精度,对材料梯度参数的变化不敏感.可以进一步拓展到FGMs瞬态问题和FGMs的热力耦合分析.     

2-D time-harmonic BEM for solids of general anisotropy with application to eigenvalue problems

We present the direct formulation of the two-dimensional boundary element method (BEM) for time-harmonic dynamic problems in solids of general anisotropy. We split the fundamental solution, obtained by Radon transform, into static singular and dynamics regular parts. We evaluate the boundary integrals for the static singular part analytically and those for the dynamic regular part numerically over the unit circle.We apply the developed BEM to eigenvalue analysis. We determine eigenvalues of full non-symmetric complex-valued matrices, depending non-linearly on the frequency, by first reducing them to the generalized linear eigenvalue problem and then applying the QZ algorithm. We test the performance of the QZ algorithm thoroughly in comparison with the FEM solution. The proposed BEM is not only a strong candidate to replace the FEM for industrial eigenvalue problems, but it is also applicable to a wider class of two-dimensional time-harmonic problems.     

边界元奇异与近奇异数值积分方法及其应用于大规模声学问题

提出了综合处理Burton-Miller方法所导致的奇异积分与近奇异积分问题的数值求积方法,以此改进了基于常量元素的常规边界元和低频快速多极边界元方法。对于奇异积分问题,利用Hadamard有限积分方法进行解决;对于近奇异积分问题,则采用极坐标变换法和PART方法(Projection and Angular&;Radial Transformation)进行克服。与解析解和LMS Virtual.Lab商业软件的结果比较验证了方法的正确性,并对比分析了奇异积分与近奇异积分对计算精度的影响。采用低频快速多极子方法以加速常规边界元法的计算效率,计算分析了计算复杂度,并成功实现了34万自由度大规模问题的计算。结果表明,近奇异积分问题主要由超奇异核函数引起,对计算精度的影响不容忽略;快速多极边界元法的精度与常规边界元法一致,但计算复杂度要远低于后者。     

大规模声学边界元法的GPU并行计算

提出一种大规模声学边界元法的高效率、高精度GPU并行计算方法.基于Burton-Miller边界积分方程,推导适于GPU的并行计算格式并实现了传统边界元法的GPU加速算法.为提高原型算法的效率,研究GPU数据缓存优化方法.由于GPU的双精度浮点运算能力较低,为了降低数值误差,研究基于单精度浮点运算实现的doublesingle精度算法.数值算例表明,改进的算法实现了最高89.8%的GPU使用效率,且数值精度与直接使用双精度数相当,而计算时间仅为其1/28,显存消耗也仅为其一半.该方法可在普通PC机(8GB内存,NVIDIA Ge Force 660 Ti显卡)上快速完成自由度超过300万的大规模声学边界元分析,计算速度和内存消耗均优于快速边界元法.     

Numerical quadrature for singular and near-singular integrals of boundary element method and its applications in large-scale acoustic problems

The numerical quadrature methods for dealing with the problems of singular and near-singular integrals caused by Burton-Miller method are proposed,by which the conventional and fast multipole BEMs(boundary element methods) for 3D acoustic problems based on constant elements are improved.To solve the problem of singular integrals,a Hadamard finite-part integral method is presented,which is a simplified combination of the methods proposed by Kirkup and Wolf.The problem of near-singular integrals is overcome by the simple method of polar transformation and the more complex method of PART(Projection and Angular Radial Transformation).The effectiveness of these methods for solving the singular and near-singular problems is validated through comparing with the results computed by the analytical method and/or the commercial software LMS Virtual.Lab.In addition,the influence of the near-singular integral problem on the computational precisions is analyzed by computing the errors relative to the exact solution.The computational complexities of the conventional and fast multipole BEM are analyzed and compared through numerical computations.A large-scale acoustic scattering problem,whose degree of freedoms is about 340,000,is implemented successfully.The results show that,the near singularity is primarily introduced by the hyper-singular kernel,and has great influences on the precision of the solution.The precision of fast multipole BEM is the same as conventional BEM,but the computational complexities are much lower.     

Multilevel fast multipole algorithm for acoustic wave scattering by truncated ground with trenches

The multilevel fast multipole algorithm (MLFMA) is extended to solve for acoustic wave scattering by very large objects with three-dimensional arbitrary shapes. Although the fast multipole method as the prototype of MLFMA was introduced to acoustics early, it has not been used to study acoustic problems with millions of unknowns. In this work, the MLFMA is applied to analyze the acoustic behavior for very large truncated ground with many trenches in order to investigate the approach for mitigating gun blast noise at proving grounds. The implementation of the MLFMA is based on the Nystrom method to create matrix equations for the acoustic boundary integral equation. As the Nystrom method has a simpler mechanism in the generation of far-interaction terms, which MLFMA acts on, the resulting scheme is more efficient than those based on the method of moments and the boundary element method (BEM). For near-interaction terms, the singular or near-singular integrals are evaluated using a robust technique, which differs from that in BEM. Due to the enhanced efficiency, the MLFMA can rapidly solve acoustic wave scattering problems with more than two million unknowns on workstations without involving parallel algorithms. Numerical examples are used to demonstrate the performance of the MLFMA with report of consumed CPU time and memory usage.     

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

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

Use of nonsingular boundary integral formulation for reducing errors due to near-field measurements in the boundary element method based near-field acoustic holography

In the conformal near-field acoustic holography (NAH) using the boundary element method (BEM), the transfer matrix relating the vibro-acoustic properties of source and field depends solely on the geometrical condition of the problem. This kind of NAH is known to be very powerful in dealing with the sources having irregular shaped boundaries. When the vibro-acoustic source field is reconstructed by using this conformal NAH, one tends to position the sensors as close as possible to the source surface in order to get rich information on the nonpropagating wave components. The conventional acoustic BEM based on the Kirchhoff-Helmholtz integral equation has the singularity problem in the close near field of the source surface. This problem stems from the singular kernel of the Green function of the boundary integral equation (BIE) and the singularity can influence the reconstruction accuracy greatly. In this paper, the nonsingular BIE is introduced to the NAH calculation and the holographic BIE is reformulated. The effectiveness of nonsingular BEM has been investigated for the reduction of reconstruction error. Through interior and exterior examples, it is shown that the resolution of predicted field pressure could be improved in the close near field by employing the nonsingular BIE. Because the BEM-based NAH inevitably requires the field pressure measured in the close proximity to the source surface, the present approach is recommended for improving the resolution of the reconstructed source field.     

声学Burton-Miller方程边界元法GPU并行计算

基于GPU,对声学Burton-Miller积分方程的边界元解法进行并行计算.提出并行计算格式和程序实现方法,以及Burton-Miller方程中各类奇异(包括强奇异、超奇异)积分的GPU计算和局部修正方法.典型算例结果表明,在特征频率处可获得正确的解,具有较高精度,可在普通个人计算机上快速完成自由度超过2×10⁵的声学边界元分析.为计算声学及相关工程领域的中、大规模声场分析问题提供一种快速、高效、简便的数值计算工具.     

Burton–Miller-type singular boundary method for acoustic radiation and scattering

This paper proposes the singular boundary method (SBM) in conjunction with Burton and Miller?s formulation for acoustic radiation and scattering. The SBM is a strong-form collocation boundary discretization technique using the singular fundamental solutions, which is mathematically simple, easy-to-program, meshless and introduces the concept of source intensity factors (SIFs) to eliminate the singularities of the fundamental solutions. Therefore, it avoids singular numerical integrals in the boundary element method (BEM) and circumvents the troublesome placement of the fictitious boundary in the method of fundamental solutions (MFS). In the present method, we derive the SIFs of exterior Helmholtz equation by means of the SIFs of exterior Laplace equation owing to the same order of singularities between the Laplace and Helmholtz fundamental solutions. In conjunction with the Burton–Miller formulation, the SBM enhances the quality of the solution, particularly in the vicinity of the corresponding interior eigenfrequencies. Numerical illustrations demonstrate efficiency and accuracy of the present scheme on some benchmark examples under 2D and 3D unbounded domains in comparison with the analytical solutions, the boundary element solutions and Dirichlet-to-Neumann finite element solutions.