高精度特征解及其外推求解位势方程

根据位势理论,基本边界特征值问题可转化为具有对数奇性的边界积分方程.利用机械求积方法求解特征值和特征向量,以及利用这些特征解求解Laplace方程.特征解和Laplace方程的解具有高精度和低的计算复杂度.利用Anselone聚紧和渐近紧理论,证明了方法的收敛性和稳定性.此外,还给出了误差的奇数阶渐近展开.利用h3-Richardson外推,不仅误差近似的精度阶大为提高,而且,得到的后验误差估计可以构造自适应算法.具体的数值例子说明了算法的有效性.     

平面双调和问题的第一类边界积分方程的高精度求积方法与外推

借助位势理论,平面双调和方程的Dirichlet问题被转化为第一类边界积分方程组,本文使用新型的反常积分的求积公式构造出解造解此类边界积分方程的机械求积方法,证明了该方法具有O（h^3)阶精度和误差的h^3幂渐近展开,故借助Richardson外推还能提高精度阶。     

解带有Hilbert核的奇异积分方程的高精度组合算法

该文给出了用求积法解带Hilber核的奇异积分方程的高精度组合算法.把网格点分成互不相交的子集合, 在子集合上并行求解离散方程组, 再利用组合算法求得全局网格点的逼近.如果积分方程的系数属于B_δ, 则求积法的精度可达O(e^-nδ). 此外, 使用组合算法不仅能得到更高的精度阶, 而且能够得到后验误差估计. 数值算例的结果表明组合算法是极其有效的.     

粘弹性方程类Wilson元的整体超收敛和外推

本文借助双线性元积分恒等式技巧,对粘弹性方程的类Wilson元解进行了高精度分析.通过证明类 Wilson元的非协调误差在矩形网格下可以达到O(h3)这一独特性质及利用插值后处理技术给出了H1模意义下O(h2)阶的超逼近和整体超收敛结果.进而通过构造合适的外推格式,得到具有更高阶O(h3)精度的数值逼近解.     

一类分数阶Langevin方程block-by-block算法的数值分析

分数阶Langevin方程有重要的科学意义和工程应用价值,基于经典block-by-block算法,求解了一类含有Caputo导数的分数阶Langevin方程的数值解.Block-by-block算法通过引入二次Lagrange基函数插值,构造出逐块收敛的非线性方程组,通过在每一块耦合求得分数阶Langevin方程的数值解.在0<α<1条件下,应用随机Taylor展开证明block-by-block算法是3+α阶收敛的,数值试验表明在不同α和时间步长h取值下,block-by-block算法具有稳定性和收敛性,克服了现有方法求解分数阶Langevin方程速度慢精度低的缺点,表明block-by-block算法求解分数阶Langevin方程是高效的.     

外部非定常Navier-Stokes方程的非线性Galerkin算法

提出了求解外部非定常Navier-stokes方程的有限元边界元耦合的非线性Galerkin算法,证明了相应变分问题的正则性和数值解的收敛速度。收敛性分析表明如果选取粗网格尺度H是细网格尺度h的开平方数量级,则该算法提供了与古典Galerkin算法同阶的收敛速度。然而非线性Galerkin算法仅仅需要在粗网格解非线性问题,在细网格上解线性问题。因此,该算法可以节省计算工作量。     

非线性湿气迁移方程的超收敛分析

本文用混合有限元方法研究一般的非线性湿气迁移方程.利用双线性元Q₁₁和零阶Raviart-Thomas元(Q₁₀×Q₀₁)证明方程的超收敛性.利用这两个单元插值算子的性质和平均值技巧,得到了方程半离散格式的O(h²)阶超收敛结果.对于方程线性化的Crank-Nicolson(C-N)全离散格式,得到了具有O(h²+τ²)阶的超收敛结果,这里h是空间剖分参数,τ是时间步长.该方法说明如果线性化问题有超收敛性,那么对应的非线性问题有同样的超收敛性.最后,给出数值算例,证实了理论分析的正确性和方法的有效性.     

Adomian分解法求解非线性分数阶Fredholm积分微分方程

为了求解一类非线性分数阶Fredholm积分微分方程的数值解,本文将Adomian分解法(Adomian Decomposition Method,ADM)引入到非线性分数阶Fredholm微积分方程的求解中.将ADM多项式与分数阶积分定义有效结合,得到Adomian级数解.通过收敛性分析证明所得的级数解收敛于精确解,给出最大绝对截断误差.并结合实例,证明方法的有效性和实用性.     

二维分数阶发展型方程的正式的二阶BDF交替方向隐式紧致差分格式

该文将研究二维分数阶发展型方程的正式的二阶向后微分公式(BDF)的交替方向隐式(ADI)紧致差分格式.在时间方向上用二阶向后微分公式离散一阶时间导数,积分项用二阶卷积求积公式近似,在空间方向上用四阶精度的紧致差分离散二阶空间导数得到全离散紧致差分格式.基于与卷积求积相对应的实二次型的非负性,利用能量方法研究了差分格式的稳定性和收敛性,理论结果表明紧致差分格式的收敛阶为O(k~(a+1)+h_1~4+h_2~4),其中k为时间步长,h_1和h_2分别是空间x和y方向的步长.最后,数值算例验证了理论分析的正确性.     

非线性Sine-Gordon方程Hermite型有限元新的超收敛分析及外推

在半离散格式下讨论了一类非线性Sine-Gordon方程的Hermite型矩形元逼近.利用该元的高精度分析和对时间t的导数转移技巧,得到了H1模意义下O(h2)阶的最优误差估计和O(h3)阶的超逼近性.进一步地,通过运用插值后处理方法,给出了超收敛结果.与此同时,借助于构造一个新的外推格式,导出了与线性情形相同的O(h4)阶外推解.     

Boundary value methods for Volterra integral and integro-differential equations

New and effective quadrature rules generated by boundary value methods are introduced. We employ the introduced quadrature rules to construct quadrature methods for the second kind Volterra integral equations and Volterra integro-differential equations. These methods are shown to be effective and possess excellent convergence properties. The nonlinear multigrid method is applied to solve the discrete systems derived from the introduced numerical scheme. Numerical simulations are presented and confirm the efficiency and accuracy of the methods.     

解平面弹性力学问题的边界积分方程的机械求积方法及其外推

１．引 言 考虑平面弹性力学内或外位移边值问题和内或外应力边值问题这里Ω是平面有界开集,Ωｃ是闭包Ω的补集,Γ是Ω或Ωｃ的边界,ｕ＝（ｕ１,ｕ２）是位移,ｎ＝（ｎ１,ｎ２）是Γ的外法向单位向量,δｉｊ＝（ｕｉ,ｊ＋ｕｊ,ｉ）／２是应变张量,λ和μ是Ｌａｍｅ常数,并且按张量计算规则：重复下标蕴含对该下标从１到２的求和． 使用直接边界元方法（１．１）与（１．２）皆可被转换为边界积分方程组这里αｉｊ（ｙ）是取决于ｙ∈Γ的常数,当ｙ是Γ的光滑点时,;式中是ｋｅｌｖｉｎ基本解,有以下表达式［５,７］这里ｒ＝…     

Extrapolation algorithms for solving nonlinear boundary integral equations by mechanical quadrature methods

We study the numerical solution procedure for two-dimensional Laplace’s equation subjecting to non-linear boundary conditions. Based on the potential theory, the problem can be converted into a nonlinear boundary integral equations. Mechanical quadrature methods are presented for solving the equations, which possess high accuracy order O(h ³) and low computing complexities. Moreover, the algorithms of the mechanical quadrature methods are simple without any integration computation. Harnessing the asymptotical compact theory and Stepleman theorem, an asymptotic expansion of the errors with odd powers is shown. Based on the asymptotic expansion, the h ³ −Richardson extrapolation algorithms are used and the accuracy order is improved to O(h ⁵). The efficiency of the algorithms is illustrated by numerical examples.     

The Mechanical Quadrature Methods and Their Extrapolation for Solving BIE of Steklov Eigenvalue Problems

By means of the potential theory Steklov eigenvalue problems are transformed into general eigenvalue problems of boundary integral equations (BIE) with the logarithmic singularity.Using the quadrature rules, the paper presents quadrature methods for BIE of Steklov eigenvalue problem, which possess high accuracies O(h^3) and low computing complexities. Moreover, an asymptotic expansion of the errors with odd powers is shown. Using h^3-Richardson extrapolation, we can not only improve the accuracy order of approximations, but also derive a posterior estimate as adaptive algorithms. The efficiency of the algorithm is illustrated by some examples.     

Saturation-free numerical scheme for computing the flow past a lattice of airfoils and the determination of separation points in a viscous fluid

A numerical method for computing the potential flow past a lattice of airfoils is described. The problem is reduced to a linear integrodifferential equation on the lattice contour, which is then approximated by a linear system of equations with the help of specially derived quadrature formulas. The quadrature formulas exhibit exponential convergence in the number of points on an airfoil and have a simple analytical form. Due to its fast convergence and high accuracy, the method can be used to directly optimize the airfoils as based on any given integral characteristics. The shear stress distribution and the separation points are determined from the velocity distribution at the airfoil boundary calculated by solving the boundary layer equations. The method proposed is free of laborious grid generation procedures and does not involve difficulties associated with numerical viscosity at high Reynolds numbers.     

Numerical solution of nonlinear two-dimensional Fredholm integral equations of the second kind using Gauss product quadrature rules

The Gauss product quadrature rules and collocation method are applied to reduce the second-kind nonlinear two-dimensional Fredholm integral equations (FIE) to a nonlinear system of equations. The convergence of the proposed numerical method is proved under certain conditions on the kernel of the integral equation. An iterative method for approximating the solution of the obtained nonlinear system is provided and its convergence is proved. Also, some numerical examples are presented to show the efficiency and accuracy of the proposed method.     

Block boundary value methods for solving Volterra integral and integro-differential equations

Reducible quadrature rules generated by boundary value methods are considered in block version and applied to solve the second kind Volterra integral equations and Volterra integro-differential equations. These extended block boundary value methods are shown to possess both excellent stability properties and high accuracy for Volterra-type equations. Numerical experiments are presented and the efficiency, accuracy and stability of the schemes are confirmed.     

Dirac delta methods for Helmholtz transmission problems

In this paper we use a boundary integral method with single layer potentials to solve a class of Helmholtz transmission problems in the plane. We propose and analyze a novel and very simple quadrature method to solve numerically the equivalent system of integral equations which provides an approximation of the solution of the original problem with linear convergence (quadratic in some special cases). Furthermore, we also investigate a modified quadrature approximation based on the ideas of qualocation methods. This new scheme is again extremely simple to implement and has order three in weak norms.      

Hybrid Galerkin boundary elements: theory and implementation

Summary. In this paper we present a new quadrature method for computing Galerkin stiffness matrices arising from the discretisation of 3D boundary integral equations using continuous piecewise linear boundary elements. This rule takes as points some subset of the nodes of the mesh and can be used for computing non-singular Galerkin integrals corresponding to pairs of basis functions with non-intersecting supports. When this new rule is combined with standard methods for the singular Galerkin integrals we obtain a “hybrid” Galerkin method which has the same stability and asymptotic convergence properties as the true Galerkin method but a complexity more akin to that of a collocation or Nystr?m method. The method can be applied to a wide range of singular and weakly-singular first- and second-kind equations, including many for which the classical Nystr?m method is not even defined. The results apply to equations on piecewise-smooth Lipschitz boundaries, and to non-quasiuniform (but shape-regular) meshes. A by-product of the analysis is a stability theory for quadrature rules of precision 1 and 2 based on arbitrary points in the plane. Numerical experiments demonstrate that the new method realises the performance expected from the theory. Received January 22, 1998 / Revised version received May 26, 1999 / Published online April 20, 2000 –? Springer-Verlag 2000     

Integration of singular Galerkin-type boundary element integrals for 3D elasticity problems

Summary. A Galerkin approximation of both strongly and hypersingular boundary integral equation (BIE) is considered for the solution of a mixed boundary value problem in 3D elasticity leading to a symmetric system of linear equations. The evaluation of Cauchy principal values (v. p.) and finite parts (p. f.) of double integrals is one of the most difficult parts within the implementation of such boundary element methods (BEMs). A new integration method, which is strictly derived for the cases of coincident elements as well as edge-adjacent and vertex-adjacent elements, leads to explicitly given regular integrand functions which can be integrated by the standard Gauss-Legendre and Gauss-Jacobi quadrature rules. Problems of a wide range of integral kernels on curved surfaces can be treated by this integration method. We give estimates of the quadrature errors of the singular four-dimensional integrals. Received June 25, 1995 / Revised version received January 29, 1996