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

将一种通用算法应用于平面位势问题边界元法中近边界点几乎奇异积分的正则化。对线性单元,位势问题近边界点的几乎强和超奇异积分可归纳为两种形式。通过分部积分,将引起奇异的积分元素变换到积分号之外,从而对这两种积分分别给出了无奇异的正则化计算公式。除了线性元,二次元也应用于该算法。与近边界点临近的二次单元划分为两段线性单元,该算法仍然适用。算例证明了方法的有效性和精确性。对曲线边界问题,联合二次元和线性元可提高计算结果精确度。     

三维快速多极边界元法分析地埋管群传热问题北大核心CSCD

基于三节点三角形线性单元,为克服单元跨叶子积分难题,将三维位势问题快速多极边界元法与几乎奇异积分的半解析算法相结合,实现了三维边界元法中几乎奇异积分的准确计算,该方法适用于U型地埋管薄体结构的换热分析.在制冷、制热两种工况下研究了U型地埋管壁厚对换热量的影响,并进一步分析了管群间的热相互作用.计算结果显示,当管壁导热系数一定时,管壁越厚,对管内流体和土壤之间的换热影响越大.当钻孔间距一定时,管群中埋管数量越多,热干扰现象越强烈,提高管群换热量的主要措施是降低管群间热干扰.因准确计算了几乎奇异积分,三维快速多极边界元法可以有效计算薄体和厚体耦合的三维热传导问题.该文方法和分析结果可为地埋管换热器系统的工程应用提供参考.     

边界元法中奇异积分计算的极坐标变换法

在边界元法中,奇异积分的处理是一个极为引人注目的问题.本文提出了一种在单元状态作极坐标变换的新的处理方法,它能显式地消除奇异积分的奇异性,使之成为常规积分,因而易于在边界元法中使用高次单元.计算实例表明,本文所提出的方法是有效的、方便的.     

边界元近奇异积分计算的迭代Sinh-Sigmoidal组合式变换法

精确有效地消除积分的近奇异性是三维边界元法在工程应用中的首要问题.当源点与三角形积分单元间的距离无限趋近于零时,会出现近奇异积分问题,积分单元的形状和投影点的位置都是影响近奇异积分计算精度的重要因素.现有的非线性变换法大多只关注径向上积分的近奇异性,而忽略了角度方向和积分单元形状的影响,在投影点接近三角形积分单元边界的情况下,无法获得令人满意的计算精度,并且对子三角形积分单元的形状非常敏感.因此提出了一种改进的基于自适应分块技术和不同坐标变换的迭代sinh sigmoidal组合式变换法,分别消除径向和角度方向积分的近奇异性,在确保计算精度的同时,大大减小了计算规模.数值算例验证了该方法的有效性.     

用于解析函数复分析的共轭边界元法

由2个共轭的实调和函数构建1个复解析函数,其复分析在应用数学和力学领域具有重要的作用.提出了一个加权残数方程组,证明了该方程组为2个共轭函数的域内控制方程、边界条件和边界上Cauchy Riemann(柯西-黎曼)条件的近似解,等效为复解析函数的逼近方程.在离散空间中,由该加权残数方程分别推导出2个位势问题的直接边界积分方程和1个表示Cauchy-Riemann条件的有限差分方程,随后解决了弱奇异线性方程组的求解难题,并提出用Cauchy积分公式求内点值的方法,从而建立了一种用于复分析的常单元共轭边界元法.最后,用3个算例证明了所提出方法适用于域内或域外的幂函数、指数函数或对数函数形式的解析函数,而且其误差与2维位势问题是同等量级的.     

二维Helmholtz方程的插值型边界无单元法

针对二维Helmholtz方程的内外边值问题,提出了插值型边界无单元法(interpolating boundary element-free method).在间接位势理论的基础上,利用Laplace方程基本解的特性,建立了求解Helmholtz方程Neumann边值内外问题的正则化形式,有效消除了强奇异积分的计算.其次通过引入全局距离展开成局部距离的幂级数,详细推导了距离函数的导数和法向导数差值的极限表达式.最后给出了4个插值型边界无单元法的数值算例,表明了该方法可取得较高的可行性和有效性.     

角度-距离复合变换法消除边界积分方程近奇异性

针对薄型结构边界单元分析过程中出现的近奇异积分问题,研究了采用一种角度变换和距离变换相结合的方法,节省了计算量,提高了计算精度.研究发现,当积分单元上与配置点距离最近的点落在积分单元的边沿或者顶点附近时,经过基础变换后的积分在两个方向都表现出奇异性,因此,对两个方向同时使用一维非线性变换能够切实消除近奇异性.数值算例验证了复合变换对近奇异积分计算精度的提升效果.     

复超球面上奇异积分方程的正则化

<正> 在一维奇异积分方程论中,复变函数的 Cauchy 型积分起着十分重要的作用;但在研究高维奇异积分方程时,利用多复变数 Cauchy 型积分作为工具者,至今尚少(参看[2]—[6]).本文是用复超球的 Cauchy 型积分边界性质,处理复超球面上含 Cauchy 核、B 核与 h 核的奇异积分方程的正则化问题.     

圆内平面弹性问题的边界积分公式

根据双解析函数可以得到单位圆内平面弹性问题应力函数的边界积分公式,但式中包含强奇异积分,不能用于直接计算．将边界上的应力函数展开为Fourier级数,再利用广义函数论中的几个公式进行卷积计算,可以得到不含强奇异积分核的边界积分公式,通过边界的应力函数值和法向导数的积分,直接得到圆内应力函数值,并给出几个算例,表明该结果用于求解单位圆内平面弹性问题十分方便．     

反演二维瞬态热传导问题随温度变化的导热系数

基于边界元法反演二维瞬态热传导问题随温度变化的导热系数.采用Kirchhoff变换将非线性的控制方程转变为线性方程.边界元法用于构建二维瞬态热传导问题的数值分析模型.将反演参数作为优化变量,测点温度计算值与测量值之间的残差平方和作为优化目标函数.引入复变量求导法求解目标函数的梯度矩阵,梯度正则化法用于优化目标函数获得反演结果.探讨时间步长、测点数量和随机偏差对反演结果的影响.减小步长、增加测点数量收敛速度加快.降低了随机偏差,计算结果更精确.算例证明了算法的有效性与稳定性.     

A stable boundary elements method for magnetohydrodynamic channel flows at high Hartmann numbers

The article is devoted to extension of boundary element method (BEM) for solving coupled equations in velocity and induced magnetic field for time dependent magnetohydrodynamic (MHD) flows through a rectangular pipe. The BEM is equipped with finite difference approach to solve MHD problem at high Hartmann numbers up to 10⁶. In fact, the finite difference approach is used to approximate partial derivatives of unknown functions at boundary points respect to outward normal vector. It yields a numerical method with no singular boundary integrals. Besides, a new approach is suggested in this article where transforms 2D singular BEM's integrals to 1D nonsingular ones. The new approach reduces computational cost, significantly. Note that the stability of the numerical scheme is proved mathematically when computational domain is discretized uniformly and Hartmann number is 40 times bigger than length of boundary elements. Numerical examples show behavior of velocity and induced magnetic field across the sections.     

Spherical Singular Integrals, Monogenic Kernels and Wavelets on the Three-Dimensional Sphere

The construction of wavelets relies on translations and dilations which are perfectly given in . On the sphere translations can be considered as rotations but it is difficult to say what are dilations. For the 2-dimensional sphere there exist two different approaches to obtain wavelets which are worth to be considered. The first concept goes back to W. Freeden and collaborators who define wavelets by means of kernels of spherical singular integrals. The other concept developed by J.P. Antoine and P. Vandergheynst is a purely group theoretical approach and defines dilations as dilations in the tangent plane. Surprisingly both concepts coincides for zonal functions. We will define singular integrals and kernels of singular integrals on the three dimensional sphere which are also approximate identities. In particular the Cauchy kernel in Clifford analysis is a kernel of a singular integral, the singular Cauchy integral, and an approximate identity. Furthermore, we will define wavelets on the 3-dimensional sphere by means of kernels of singular integrals. This paper is dedicated to the memory of our friend and colleague Jarolim Bureš Received: October, 2007. Accepted: February, 2008.     

Numerical integration of logarithmic and nearly logarithmic singularity in BEMs

A new semi-analytical integration scheme is proposed for evaluation of logarithmically singular and/or nearly singular integrals occurring in 2D BEM formulations. Extensive numerical experiments are performed to study the accuracy by the proposed and other schemes of numerical integration. The accuracy of the numerical integration is almost exact even for curvilinear elements and the accuracy of numerical computation is determined by the accuracy of the approximation of the boundary density and geometry.     

Numerical evaluation of various levels of singular integrals, arising in BEM and its application in hydrofoil analysis

The sound implementation of the boundary element method (BEM) is highly dependent on an accurate numerical integration of singular integrals. In this paper, a set of various types of singular domain integrals with three-dimensional boundary element discretization is evaluated based on a transformation integration technique. In the BEM, the integration domain (body surface) needs to be discretized into small elements. For each element, the integral I(x^p, x) is calculated on the domain dS. Several types of integrals IB_α and IC_α are numerically and analytically computed and compared with the relative error. The method is extended to evaluate singular integrals which arise in the solution of the three-dimensional Laplace’s equation. An example of the elliptic hydrofoil is performed to study the physical accuracy. The results obtained using both numerical and analytical methods are shown in good agreement with the experimental data.     

二重奇异积分近似计算的一种优化算法

本文对于矩形区域上某一内点为奇点的奇异积分的近似计算给出了优化中心数值算法，它在迭代计算过程中避免了函数值的重复计算．采用外推法减少迭代次数．     

Singularity preserving Galerkin method for Hammerstein equations with logarithmic kernel

In a recent paper [3], Cao and Xu established the Galerkin method for weakly singular Fredholm integral equations that preserves the singularity of the solution. Their Galerkin method provides a numerical solution that is a linear combination of a certain class of basis functions which includes elements that reflect the singularity of the solution. The purpose of this paper is to extend the result of Cao and Xu and to establish the singularity preserving Galerkin method for Hammerstein equations with logarithmic kernel. The iterated singularity preserving Galerkin method is also discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

Exact Singularity Subtraction from Boundary Integral Equations in Modeling Vesicles and Red Blood Cells

The study of vesicles, capsules and red blood cells (RBCs) under flow is a field of active research, belonging to the general problematic of fluid/structure interactions. Here, we are interested in modeling vesicles, capsules and RBCs using a boundary integral formulation, and focus on exact singularity subtractions of the kernel of the integral equations in 3D. In order to increase the precision of singular and near-singular integration, we propose here a refinement procedure in the vicinity of the pole of the Green-Oseen kernel. The refinement is performed homogeneously everywhere on the source surface in order to reuse the additional quadrature nodes when calculating boundary integrals in multiple target points. We also introduce a multi-level look-up algorithm in order to select the additional quadrature nodes in vicinity of the pole of the Green-Oseen kernel. The expected convergence rate of the proposed algorithm is of order$\mathcal{O}(1/N^2)$ while the computational complexity is of order$\mathcal{O}$($N^2$ln$N$), where $N$ is the number of degrees of freedom used for surface discretization. Several numerical tests are presented to demonstrate the convergence and the efficiency of the method.     

Singularity cancellation method for time-domain boundary element formulation of elastodynamics: A direct approach

In the implementation of time-domain boundary element method for elasto-dynamic problems, there are two types of singularities: the wave front singularity arising when the product of wave velocity and time is equal to the distance between the source point and the field point, and the spatial singularity arising when the source point coincides with the field point. In this paper, the singularity of the first type in the integrand is eliminated by an analytical integration over time, Cauchy principal value and Hadamard finite part integral. Four types of singularities with different orders appear in the integrand after analytical time integration. In order to accurately calculate the integral, in which the integrand is piecewise continuous, the integral domain is subdivided into several patches based on the relation between the product of wave velocity and time and the distance. In singular patches, the integrands are separated into a regular part and a singular part. The regular part can be computed by traditional numerical integration method such as Gaussian integration, while the singular part can be analytically integrated. Using the proposed method, the spatial singular integrals can be calculated directly. Numerical examples using various kinds of elements are presented to verify the proposed method.