Forward and inverse problem for cardiac magnetic field and electric potential using two boundary element methods

This paper discusses the forward and inverse problem for cardiac magnetic fields and electric potentials.A torso-heart model established by boundary element method(BEM) is used for studying the distributions of cardiac magnetic fields and electric potentials.Because node-to-node and triangle-to-triangle BEM can lead to discrepant field distributions,their properties and influences are compared.Then based on constructed torso-heart model and supposed current source functional model-current dipole array,the magnetic and electric imaging by optimal constrained linear inverse method are applied at the same time.Through figure and reconstructing parameter comparison,though the magnetic current dipole array imaging possesses better reconstructing effect,however node-to-node BEM and triangleto-triangle BEM make little difference to magnetic and electric imaging.     

Simulating regular wave propagation over a submerged bar by boundary element method model and Boussinesq equation model

Numerical models based on the boundary element method and Boussinesq equation are used to simulate the wave transform over a submerged bar for regular waves.In the boundary-element-method model the linear element is used,and the integrals are computed by analytical formulas.The Boussinesq-equation model is the well-known FUNWAVE from the University of Delaware.We compare the numerical free surface displacements with the laboratory data on both gentle slope and steep slope,and find that both the models simulate the wave transform well.We further compute the agreement indexes between the numerical result and laboratory data,and the results support that the boundary-element-method model has a stable good performance,which is due to the fact that its governing equation has no restriction on nonlinearity and dispersion as compared with Boussinesq equation.     

快速多极边界元法用于扩散光学断层成像研究

在求解扩散光学断层成像中的正向问题时, 目前普遍采用有限元法, 但是随着实际模型规模的增大, 有限元法的计算量问题日益显著, 而边界元法则由于可以降低计算维度使计算量减少而备受关注. 本文以均匀的高散射介质为模型, 研究了将快速多极边界元法用于扩散光学断层成像的正向问题. 快速多极边界元法利用核函数的多极展开, 将常规边界元法中系数矩阵和迭代矢量的乘积项等价为相应四叉树结构的一次递归, 再结合广义最小残量法进行迭代求解. 将计算结果和蒙特卡罗法的模拟结果进行了比较, 表明利用快速多极边界元法的模拟结果和蒙特卡罗法的结果有很好的一致性. 研究结果验证了快速多极边界元法可以用于扩散光学断层成像, 为其大规模和实时成像带来可观的前景. 关键词： 扩散光学断层成像 边界元法 快速多极边界元法     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势, 具有良好的工程应用前景.     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势,     

用有限电流元法识别HL-2A 装置的等离子体边界

在电流丝模型对HL-2A装置等离子体边界重建取得满意结果[3]、[5]的基础上，本文用有限电流元模型对HL-2A等离子体边界进行了重建研究。计算和实验结果表明，在通常情况下，有限电流元模型比固定电流丝模型的重建误差稍小，前者误差小于3mm，后者误差小于6mm。当部分有限电流元分布在等离子体边界之外时，用有限电流元模型仍然可以成功重建边界。有限电流元的位置在一定范围内变化时，重建的误差都很小。用普通奔腾4 2.4GHz PC机计算一组等离子体放电边界的时间不超过1ms。有限电流元法能准确而快速地识别等离子体的偏滤器位形，这对于HL-2A装置的实时位形控制是基本和重要的。     

回旋管双阳极磁控注入枪的边界元法计算

 根据边界元法建立物理与数学模型,编制了计算磁控注入枪的程序。使用该程序设计和模拟一个工作在35 GHz,70 kV,10 A基波回旋行波管放大器的双阳极磁控注入枪,获得较好的模拟结果。并对电子注在不同工作电压和磁场下的质量进行了分析计算。计算表明,边界元法在分析回旋器件电子光学方面是一种非常有效的方法。     

计算波导加载谐振腔谐振频率的新方法

基于边界元法,提出了一种新的简单有效的计算速调管输出系统波导加载腔的谐振频率、场分布的方法。计算了一种矩形波导加载矩形谐振腔的常见结构。利用主模TE10波激励,在波导的端口处采用推导出的Robin-type边界条件,通过扫频的方式,观测选定场点的值,极值对应系统的谐振频率。数值结果表明:波导加载谐振腔的谐振频率要低于同尺寸谐振腔的频率,谐振腔内的电场分布也受到耦合孔的影响。利用边界元法计算波导加载腔这样的半开放系统,具有实现简单、计算精确度高、数据输入少、计算时间短等优点。     

回旋管双阳极磁控注入枪的边界元法计算

根据边界元法建立物理与数学模型,编制了计算磁控注入枪的程序。使用该程序设计和模拟一个工作在35 GHz,70 kV,10 A基波回旋行波管放大器的双阳极磁控注入枪,获得较好的模拟结果。并对电子注在不同工作电压和磁场下的质量进行了分析计算。计算表明,边界元法在分析回旋器件电子光学方面是一种非常有效的方法。     

基于边界元方法的超导核磁共振成像设备高阶轴向匀场线圈优化算法

在磁共振成像设备中,为了消除目标区域内的高阶谐波磁场分量,传统方法采用无源匀场,但该方法匀场精度较低,针对性较差,适用于全局匀场,而有源匀场则可以通过优化线圈分布来产生所需要的特定的磁场分布.但是,由于匀场线圈线型的复杂度会随着线圈阶数的增加而增加,难以满足设计需要,因此本文提出了一种用于磁共振成像超导匀场线圈系统的多变量非线性优化设计方法.该方法基于边界元方法,将匀场线圈所产生的磁场与目标磁场之间的偏差作为目标函数,线匝间距、线圈半径等作为约束条件,通过非线性优化算法,得到满足设计要求的线圈分布.通过一个中心磁场为0.5 T的开放式双平面磁共振成像超导轴向匀场线圈的设计案例,说明本方法具有计算效率高、灵活性好的特点.     

基于边界元法的近平板圆孔气泡动力学行为研究

研究了带有圆孔的平板附近气泡动力学特性. 基于不可压缩势流理论, 建立了平板圆形破口附近气泡运动数值模型, 并针对气泡初始位置距离破口很近而导致计算结果发散的数值缺陷, 采用气泡壁和壁面融合的方法, 将流场分离为两个半无限域问题进行求解, 实现了在不同无量纲参数范围内的数值模拟, 数值结果与实验结果符合良好. 通过对圆孔附近气泡运动特性的研究发现, 圆孔对气泡的影响基本与壁面相反, 在膨胀阶段对气泡产生腔吸作用, 收缩阶段产生排斥, 在特定的工况下会产生对射流现象. 最后分析了气泡壁与壁面融合, 流场分离后的气泡动态特性以及各工况参数对其影响规律. 关键词： 气泡 边界元 射流 圆孔     

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

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

Fast multipole accelerated boundary element method for the Helmholtz equation in acoustic scattering problems

We apply the fast multipole method (FMM) accelerated boundary element method (BEM) for the three-dimensional (3D) Helmholtz equation, and as a result, large-scale acoustic scattering problems involving 400000 elements are solved efficiently. This is an extension of the fast multipole BEM for two-dimensional (2D) acoustic problems developed by authors recently. Some new improvements are obtained. In this new technique, the improved Burton-Miller formulation is employed to over-come non-uniqueness difficultie...     

A new method for orbit determination: Unit vector method

In this paper, a method of orbit determination is presented according to the principle of unit vector method (UVM). The model and arithmetic are improved and it not only suits initial orbit determination with short arc data, it also suits orbit improvement with data longer. It is also suitable for orbit of any eccentricity and any inclination. It omits most partial derivatives of all the elements which must be calculated in classical differential orbit improvement (DOI), so, it is more efficient than DOI, and the accuracy of orbit determination and convergence of algorithm are also improved appreciably.     

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

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

用边界元法分析非均匀介质中的传输线

 从静电场边值问题的积分解出发,推导出用边界元法求解分区均匀介质填充传输线问题的矩阵表达式,给出传输线电容参数的计算公式,介绍用边界元法求解分区均匀介质填充传输线问题的基本原理和求解过程。对两类传输线的计算结果表明：用边界元法求解分区均匀介质填充传输线问题,不仅具有较高的计算精度,而且可以很方便地应用于各类复杂截面分区均匀介质填充传输线问题的工程设计与计算,边界元法是求解分区均匀介质填充传输线问题的一种有效方法。     

一种新型光滑粒子动力学固壁边界施加模型

由于Lagrange粒子法的本质, 固壁边界条件的施加一直是光滑粒子动力学方法的难点之一. 本文从固壁边界的物理原理出发, 应用多层虚粒子表征固壁边界, 提出了一种新型固壁边界施加模型. 将虚粒子看作流体的扩展, 计算中虚粒子密度保持不变, 压力、速度等参数通过对流体粒子的插值获得, 虚粒子有条件的参与控制方程的计算, 对流体的密度/压力产生影响, 通过压力梯度隐式地表征壁面与流体之间的作用强度并对流体粒子施加沿壁面法线方向的斥力作用, 防止流体粒子对壁面的穿透. 数值算例测试结果表明, 与现有固壁边界施加方法相比, 本文方法更加符合流体与固壁边界作用的物理原理, 可以简单、有效地施加固壁边界条件, 方便地应用于具有复杂几何边界的问题, 获得稳定的流场形态、规则的粒子秩序及良好的速度、压力等参量的分布.     

Deriving correction functions to model the efficiency of noise barriers with complex shapes using boundary element simulations

The boundary element method (BEM) is a commonly used method to compute the insertion loss of noise barriers having arbitrary cross-sections. For large scale three-dimensional problems, however, the BEM is not feasible. On the other hand, standardized calculation methods for noise mapping are efficient, but shapes other than the straight barrier cannot be properly calculated. Attempts to merge these two approaches by using BEM to derive correction functions based on geometrical quantities such as source and target angle as well as the path length elongation between source and receiver caused by the barrier were usually focused on a small set of barrier types, dimensions, absorptive configurations, source or receiver positions. The main objective of this study is to investigate which functions based on the most common geometrical parameters are well suited for approximating the efficiency of different types of barriers, dimensions and absorptive configurations. To achieve this, numerous combinations of 7 different barrier types, different heights and widths as well as 3 different absorptive configurations were simulated using the 2D BEM for 8 different source positions. The octave-band-wise efficiency, i.e. the frequency-dependent gain in insertion loss compared to an equally high, fully reflective straight barrier was used as a basis for the correction functions. Linear as well as polynomial models were compared yielding a polynomial of third degree in the source and fourth degree in the target angle as the best model. Effects on the error using uniform sampling in the target angle instead of a uniform receiver grid as a basis for the correction functions are also investigated. Furthermore, wide-band efficiencies based on standardized traffic emission spectra are calculated showing small errors compared to single-band errors, in particular in the high-frequency range. A linear interpolation scheme is suggested to deal with barriers having dimensions not simulated in this work.     

矢势的边值关系A1t=A2t是怎样得出来的

指出某些文献在推导矢势边值关系A1t=A2t时的不妥之处.示明得出这一边值关系的正确作法.论述了en·(×A2-×A1)=0与A1t=A2t两者之间的关系.     

Position calculation of the frequency band gap in the finite photonic crystal with multiple alternations using boundary element method

In this paper, the wave transmission from finite photonic crystals with multiple alternations is investigated using boundary element method (BEM). Since that, in these structures the alternation is not in all directions of space; the investigations of the frequency band gap with desired accuracy are not practical by analytical methods. Also, the frequency dispersion of dielectric rods is an effective parameter in photonic crystals, which this effect in our calculations has been considered. Due to the high capabilities of the BEM, the transmitted wave spectrum in the photonic crystal is calculated by changing the geometrical and optical parameters of the photonic crystal and applying more alternation in its structure and the position and width of the frequency band gap is investigated. Then, it is assumed that the photonic crystal with an arbitrary angle is rotated around the axis which is perpendicular on the crystal cross section and then, it is irradiated with a plan wave. The band gap of the photonic crystals with the desired structure, desired rotation angle and multiple alternations have been solved. Very low information volume, high speed and accuracy during the calculation and useable for any desired structures are the characteristics of this method.