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

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

子波变换在声辐射和声散射数值解中的应用

通过将边界变量用于波展开,获得了求解二维及三维轴对称声辐射和声散射的边界积分方程的子波谱方法,既可求解Ｄｉｒｉｃｈｌｅｔ、Ｎｅｕｍａｎｎ问题,也可求解混合过值问题;它能处理任意边界条件的油对称体．用三维元方法解决了三维轴对称问题边界于波谱方法的奇异积分。给出了二维问题奇异积分的近似积分公式。给出了子波谱方法的系数计算方法,它与传统的边界元系数计算方法相似,易于计算机程序实现,能处理复杂的边界几何形状。该方法的优点是可以获得稀疏的系数矩阵。算例表明：该方法收敛较快,精度高。     

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

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

声学时域边界元法的时间卷积计算方法研究

核函数中保留Dirac函数的原型,形成关于时间的卷积积分,是声学时域边界元法中一种稳定、有效的时间数值积分计算方法 (CQ-BEM)。然而,传统CQ-BEM中卷积积分系数的获取有计算量大、耗时长,且对不同单元需要重新计算的问题,极大地降低了CQ-BEM法计算时域声场的效率。针对传统CQ-BEM积分系数计算效率低的问题,本文利用多项式展开定理给出了待求函数泰勒系数的解析表达与数值计算方法,建立了不同单元间待求系数的转换理论,可以在一次循环迭代内完成不同单元的积分系数的计算,大幅降低了计算量,提高了CQ-BEM方法的声场计算效率。脉动球源数值算例结果表明,在相同要求下,本文方法计算时间较传统方法减少50%以上,相对误差小5个数量级以上,且计算时间随单元数的增长率仅为传统方法的2.34%。因此,本文提出的系数计算方法能够有效提高CQ-BEM方法的时域声场计算效率,拓展了CQ-BEM在大型机电设备时域声场模拟的计算规模。     

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

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

变系数瞬态热传导问题边界元格式的特征正交分解降阶方法

提出了一种基于边界元法求解变系数瞬态热传导问题的特征正交分解(POD)降阶方法,重组并推导出变系数瞬态热传导问题适合降阶的边界元离散积分方程,建立了变系数瞬态热传导问题边界元格式的POD降阶模型,并用常数边界条件下建立的瞬态热传导问题的POD降阶模态,对光滑时变边界条件瞬态热传导问题进行降阶分析.首先,对一个变系数瞬态热传导问题,建立其边界域积分方程,并将域积分转换成边界积分;其次,离散并重组积分方程,获得可用于降阶分析的矩阵形式的时间微分方程组;最后,用POD模态矩阵对该时间微分方程组进行降阶处理,建立降阶模型并对其求解.数值算例验证了本文方法的正确性和有效性.研究表明:1)常数边界条件下建立的低阶POD模态矩阵,能够用来准确预测复杂光滑时变边界条件下的温度场结果;2)低阶模型的建立,解决了边界元法中采用时间差分推进技术求解大型时间微分方程组时求解速度慢、算法稳定性差的问题.     

有限元/边界积分方法在海面及其上方弹体目标电磁散射中的应用

本文将有限元/边界积分方法(FE/BIM)结合区域分解方法引入到粗糙海面及其上方目标 的电磁散射问题的研究中. 由于积分边界可以以任意形状设置在距模型表面任意远的距离处, 故本文采用共形人工边界结合区域分解建模方法截断模型的开放计算区域以减少求解未知量, 在截断区域内部采用有限元方法求解, 而计算区域的边界条件通过边界积分方程方法得到. 通过与矩量法获得的数值计算结果进行比较, 证明了该混合算法及模型处理方法的正确性, 进而研究了海面上方弹体目标的电磁散射特性, 并讨论了其双站散射系数随电磁波入射角度、目标高度、海面风速以及弹体尺寸的电磁散射特性变化情况. 本文结果可用于反演复杂背景下的目标信息及目标探测等领域. 关键词： 电磁散射 粗糙海面 目标 有限元/边界积分方法     

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

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

聚焦换能器自由场互易校准的冲击响应积分法衍射系数计算

在聚焦换能器的会聚球面波自由场互易校准法中，衍射损失修正系数是获得互易常数的一个关键参数。本文利用冲击响应函数积分法系统地计算了超声治疗头、机械扫描B超探头以及多普勒诊断探头的衍射损失系数，得出了一些典型参数下的计算曲线，可供不同类型的聚焦换能器会聚球面波互易校准时引用。与传统的瑞利积分法分析比较显示，本文中所采用的冲击响应函数积分法显著地提高了运算速度，同时基本保持了与瑞利积分算法计算结果的一致性。     

三维声学多层快速多极子边界元及其应用

快速多极子边界元算法可以加速矩阵和向量乘法运算, 将传统边界元算法的计算量和内存占用量分别降为O(N log²N)和O(N), 适用于大型声学模型模拟计算. 本文发展了一种基于Burton-Miller方程的三维多层声学快速多极子边界元算法. 将新的自适应树状算法应用到对角形式的快速多极子边界元算法, 并使用最新提出的解析式源点矩计算公式, 进一步提高了快速多极子边界元的计算效率. 绝对软球体在内部共振频率处的散射声场计算, 验证了所发展算法在共振频率处求解的正确性. 与Bapat所提供的程序在多脉动球体辐射声场计算精度的比较, 验证了算法及程序在大型模型声学计算中的准确性, 同时显示了其求解的高效性. 最后, 将该算法用于车内声场及水下声学探测的分析计算.     

白光和红绿蓝光LED可调光源研制

为了提高LED光源色温和亮度的调节精度和准确度,结合色温由低向高变化时光色所呈现的渐变特点,提出了一种低色温白光LED灯珠、高色温白光LED灯珠加红绿蓝光LED灯珠补偿式调光的方法.将色温分成三个部分进行调节,每个部分选用不同的LED灯珠组合来进行调光.实验结果表明:不同组合情况下的LED光源的初始输出色温相对于目标值的偏差范围在1％以内;亮度可以在保证色温不变的情况下独立进行调节,初始输出值与目标值的偏差范围在1％以内;经过微调之后可以达到目标值;达到了色温和亮度独立调节的要求;光源发光稳定,不会因为长时间工作而影响调节精度.     

Muonic and pionic L- and M-series in carbon and oxygen and the pionic 2p level shift and width of oxygen

Muonic and pionic X-rays of the L- and M-series in C and O have been measured with a Si(Li) detector in the energy range between 7 keV and 60 keV. The target consisted of mylar (C₅H₄O₂). Energies and intensities of 21 transitions have been determined. The strong interaction shift of the pionic 2p level in O was measured and found to be +4.1 ±2.3 eV. The measured width of this level is 11±6 eV. The measured yields have been compared with cascade calculations.     

Visualization and interferometry for cylindrical and spherical sample tomography and profilometry

A technique combining image processing and laser interferometry for visualizing and detecting the deformation of transparent cylindrical and spherical sample is proposed. This deformation includes geometric deformation such as volume transition in profilometry and physical deformation such as refractive index change in tomography. Phase contour lines are used for quantitative analysis and graphical representation of the deformation. This method allows us to visually detect the spatial variation of the deformation field and to evaluate the test quality such as misalignment of optical system. A theoretical analysis using phase contour map to characterize the deformation field is described in detail. A method using phase contour map to qualify the interferometric test is proposed. Analysis of test examples is carried out. Suggestions on using phase contour line method to ameliorate test system design are finally discussed.     

Nonlocal potentials and their exact and approximate local and velocity-dependent equivalents

We show that good approximations to the exact equivalent local potential (ELP) and damping factor of a nonlocal Perey-Buck potential can be calculated in the partial wave WKB approximation of Horiuchi. The exact ELP and damping factor are obtained by means of a method previously given by one of us. We also confirm that an approximate ELP proposed by Bauhoff et al. is of comparable accuracy as the Horiuchi approximation. Thesel-dependent ELP's exhibit reduced attraction in the interior and provide a test for higher order WKB approximations. We subsequently obtain an equivalent velocity dependent potential (EVDP) which is even exactly wave function equivalent to the original nonlocal potential. This almost local potential, unlike the trivial equivalent local potential, is smooth and well-behaved and is therefore particularly useful in nuclear reactions where the off-shell behaviour of the potential is important.     

Quantum and Classical Disparity and Accord

Discrepancies and accords between quantum (QM) and classical mechanics (CM) related to expectation values and periods are generally found for both the harmonic oscillator (SHO) and a free particle in a box (FPB), which may apply generally. These indicate non-locality is expected throughout QM. The FPB energy states violate the Correspondence Principle. Previously unexpected accords are found and proven that 〈x ²〉 _CM =〈x ²〉 _QM and τ _CM =τ _QMb (beat period i.e. beats between the phases for adjoining energy states) for the SHO for all quantum numbers, n. However, for the FPB the beat periods differ at small n. It is shown that a particle’s velocity in an infinite square well varies, no matter how wide the box, nor how far the particle is from the walls. The quantum free particle variances share an indirect commonality with the Aharonov-Bohm and Aharonov-Casher effects in that there is a quantum action in the absence of a force. The concept of an “Expectation Value over a Partial Well Width” is introduced. This paper raises the question as to whether these inconsistencies are undetectable, or can be empirically ascertained. These inherent variances may need to be fixed, or nature is manifestly more non-classical than expected.     

Spin and Statistics and First Principles

It was shown in the early seventies that, in Local Quantum Theory (that is the most general formulation of Quantum Field Theory, if we leave out only the unknown scenario of Quantum Gravity) the notion of Statistics can be grounded solely on the local observable quantities (without assuming neither the commutation relations nor even the existence of unobservable charged field operators); one finds that only the well known (para)statistics of Bose/Fermi type are allowed by the key principle of local commutativity of observables. In this frame it was possible to formulate and prove the Spin and Statistics Theorem purely on the basis of First Principles.