牛顿科茨公式计算超奇异积分的误差估计

超奇异积分的数值计算是边界元方法中的重要的课题之一,本文得到了牛顿科茨公式计算任意阶超奇异积分误差估计,当误筹函数中的Sκ(p)(Τ)=0时,便得到超收敛现象,并给出了Sκ(p)(Τ)之间的相互关系.相应的数值算例验证了理论分析的正确性.     

圆周上带希尔伯特核的柯西奇异积分方程的数值解(英文)

本文研究了圆周上带希尔伯特核的柯西奇异积分的复合梯型公式.利用连续的分片线性函数逼近被积函数,得到积分公式的误差估计;然后用积分公式构造求解对应奇异积分方程的两种格式;最后给出数值实验验证理论分析结果.     

Legendre小波求解超奇异积分

超奇异积分的数值算法一直是近些年来研究的重要课题. 基于超奇异积分的 Hadamard 有限部分积分定义, 本文给出了利用 Legendre 小波计算超奇异积分的方法. 当奇异点位于区间内时, 由于 Legendre 小波具有很好的正交性、显式表达式以及小波函数的可计算性, 将区间内的奇异点变换到区间端点处, 再利用区间端点处 Hadamard 有限部分积分的定义,进而可以计算 p+1(p∈N⁺) 阶超奇异积分. 文中最后给出的算例表明了该方法的可行性和有效性.     

偏微分方程边值反问题的数值方法研究

本文研究了奇异积分方程在反边值问题中的应用问题.利用圆周上的自然积分方程及其反演公式,把Laplace方程的边值反问题转化为一对超奇异积分方程和弱奇异积分方程的组合,通过选取三角插值近似奇异积分的计算并构造相应的配置格式,并使用Tikhonov正则化方法求解所得到的线性方程组.数值实验表明了该方法的有效性.     

薄体结构温度场的高阶边界元分析

分析了二维问题边界元法3节点二次单元的几何特征,区分和定义了源点相对高阶单元的Ⅰ型和Ⅱ型接近度.针对二维位势问题高阶边界元中奇异积分核,构造出具有相同Ⅱ型几乎奇异性的近似核函数,在几乎奇异积分单元上分离出积分核中主导的奇异函数部分.原积分核扣除其近似核函数后消除几乎奇异性,成为正则积分核函数,并采用常规Gauss数值方法计算该正则积分;对奇异核函数的积分推导出解析公式,从而建立了一种新的边界元法高阶单元几乎奇异积分半解析算法.应用该算法计算了二维薄体结构温度场算例,计算结果表明高阶单元半解析算法能充分发挥边界元法优势,显著提高计算精度.     

开口曲线上的奇异积分近似求积公式

本文利用复插值样函数讨论了开口光滑曲线上的奇异积分在被积函数分别属于 H类和 H* 类时的近似求积公式 ,给出了误差估计和收敛性     

边界元方法中超奇异积分的计算方法 献给林群教授80华诞

超奇异积分的近似计算是边界元方法,特别是自然边界元理论中必须面对的难题之一.经典的数值方法,如Gauss求积公式和Newton-Cotes积分公式等数值方法,都不能直接用于超奇异积分的近似计算.本文将介绍超奇异积分基于不同定义的Gauss积分公式、S型变换公式、Newton-Cotes积分公式和外推法近似计算超奇异积分的思路,重点阐述Newton-Cotes积分公式和基于有限部分积分定义的外推法近似计算超奇异积分的主要结论.     

一类多节对偶积分方程组正则化为超(强)奇异积分方程组求解法

基于Mellin变换法,首先方程组进行Mellin变换,然后,通过引入新的未知函数的Mellin变换代换原来未知函数的Mellin变换,使对偶积分方程组退耦正则化为超(强)奇异积分方程组．将未知函数分解并表示成未知函数和已知幂函数的乘积,幂指数(a_i,v_i)需使超(强)奇异积分方程组中的超(强)奇异积分,在端点(a_i,b_i)有界或可积奇异,求解超(强)奇异积分方程组可以使用有限部分积分式．将未知函数展成任意完备函数系(?)_n*(u)的级数,将超(强)奇异积分方程组,化成线性代数方程组,通过求解级数中的各项系数,由此给出对偶积分方程组的一般性解．并严格证明了对偶积分方程组和由它化成的超(强)奇异积分方程组的等价性,解的存在性和解的表示形式不唯一性．本文给出的理论解和解法,可供求解数学,物理,力学中的混合边值问题应用．     

圆周上超奇异积分计算及其误差估计

许多科学和工程计算问题归化为积分方程。这些积分方程往往是奇异的,有些甚至是超奇异的。特别是自然边界归化无一例外都导致超奇异积分方程。由于通常的数值积分方法对计算超奇异积分都无效,故长期以来边界元研究中总是回避超奇异积分方程。但近年来     

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

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

二维Fredholm方程的Taylor展开式解法

本文讨论了一类二维Fredholm方程的一种近拟解,通过利用二元函数的Taylor展开式,积分方程转化成一个关于未知函数及其相应的偏导数的线性代数方程组.数值例子表明了该方法的有效性.     

On the convergence of the vortex loop method with regularization for the Neumann boundary value problem on a plane screen

We consider a linear integral equation with a supersingular integral treated in the sense of the Hadamard finite value, which arises in the solution of the Neumann boundary value problem for the Laplace equation with the representation of the solution in the form of a doublelayer potential. We consider the case in which the exterior boundary value problem is solved outside a plane surface (a screen). For the integral operator in the above-mentioned equation, we suggest quadrature formulas of the vortex loop method with regularization, which provide its approximation on the entire surface when using an unstructured partition. In the problem in question, the derivative of the unknown density of the double-layer potential, as well as the errors of quadrature formulas, has singularities in a neighborhood of the screen edge. We construct a numerical scheme for the integral equation on the basis of the suggested quadrature formulas and prove an estimate for the norm of the inverse matrix of the resulting system of linear equations and the uniform convergence of the numerical solutions to the exact solution of the supersingular integral equation on the grid.     

FINITE ELEMENT NONLINEAR GALERKIN COUPLING METHOD FOR THE EXTERIOR STEADY NAVIER-STOKES PROBLEM

1.IntroductionNonlinearGalerkinmethodsaremultilevelschemesforthedissipativeevolutionpartialdifferentialequations.Theycorrespondtothesplittingsoftheunknownu:u=y z)wherethecomponentsareofdifferentorderofmagnitudewithrespecttoaparameterrelatedtothespati...     

An approximation method using approximate approximations

The aim of this article is to extend the method of approximate approximations to boundary value problems. This method was introduced by V. Maz'ya in 1991 and has been used until now for the approximation of smooth functions defined on the whole space and for the approximation of volume potentials. In the present article we develop an approximation procedure for the solution of the interior Dirichlet problem for the Laplace equation in two dimensions using approximate approximations. The procedure is based on potential theoretical considerations in connection with a boundary integral equations method and consists of three approximation steps as follows. In the first step, the unknown source density in the potential representation of the solution is replaced by approximate approximations. In the second, the decay behavior of the generating functions is used to gain a suitable approximation for the potential kernel, and in the third, Nyström's method leads to a linear algebraic system for the approximate source density. For every step a convergence analysis is established and corresponding error estimates are given.     

Thin plate spline Galerkin scheme for numerically solving nonlinear weakly singular Fredholm integral equations

The present work proposes a numerical method to obtain an approximate solution of non-linear weakly singular Fredholm integral equations. The discrete Galerkin method in addition to thin-plate splines established on scattered points is utilized to estimate the solution of these integral equations. The thin-plate splines can be regarded as a type of free shape parameter radial basis functions which create an efficient and stable technique to approximate a function. The discrete Galerkin method for the approximate solution of integral equations results from the numerical integration of all integrals in the method. We utilize a special accurate quadrature formula via the non-uniform composite Gauss-Legendre integration rule and employ it to compute the singular integrals appeared in the scheme. Since the approach does not need any background meshes, it can be identified as a meshless method. Error analysis is also given for the method. Illustrative examples are shown clearly the reliability and efficiency of the new scheme and confirm the theoretical error estimates.     

A COLLOCATION METHOD FOR THE CONDUCTIVITY PROBLEM WITH DISCONTINUOUS COEFFICIENT

In this paper, a new collocation BEM for the Robin boundary value problem of the conductivity equation ▽(γ▽u) = 0 is discussed, where the 7 is a piecewise constant function. By the integral representation formula of the solution of the conductivity equation on the boundary and interface, the boundary integral equations are obtained. We discuss the properties of these integral equations and propose a collocation method for solving these boundary integral equations. Both the theoretical analysis and the error analysis are presented and a numerical example is given.     

一道用元素法求立体体积习题的探讨

定积分的主要思想是用近似的方法获得微元的表示,然后用积分得到精确值.合理选取积分元素是运用定积分元素法解决问题的关键.对一道用元素法求立体体积的习题进行了探讨.     

A numerical study of parabolic problems with nonlinear boundary conditions

In this paper, we consider an initial‐boundary value problem for a parabolic equation with nonlinear boundary conditions. The solution to the problem can be expressed as a convolution integral of a Green's function and two unknown functions. We change the problem to a system of two nonlinear Volterra integral equations of convolution type. By using an explicit procedure on the basis of Sinc‐function properties, the resulting integral equations are replaced by a system of nonlinear algebraic equations, whose solution yields an accurate approximate solution to the parabolic problem. Some examples are considered to illustrate the ability of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.