基于单位分解积分的伽辽金无网格方法研究

数值积分是伽辽金无网格方法实施的一个重要环节,提出了一种适合于伽辽金无网格方法的单位分解积分技术．该积分技术建立在有限覆盖和单位分解基础之上,不需要对积分区域进行分解,具有较高的积分精度．并以无单元伽辽金方法为例,详细说明了基于单位分解积分的伽辽金无网格方法的实现过程．这样,在近似函数建立和数值积分过程中都不需要进行网格划分,从而形成一种“真正的”无网格方法．     

薄板的局部Petrov-Galerkin方法

利用薄板控制微分方程的等效积分对称弱形式和对变量(挠度)采用移动最小二乘近似函数进行插值,研究了薄板弯曲问题的无网格局部Petrov-Galerkin方法．这是一种真正的无网格方法,它不需要任何有限元或边界元网格,不管这种网格是用于能量积分还是进行插值的目的．所有的积分都在规则形状的子域及其边界上进行,并用罚因子法施加本质边界条件．数值例子表明,无网格局部Petrov-Galerkin法不但能够求解二阶微分方程的边值问题,而且求解四阶微分方程的边值问题也很有效,也具有收敛快、稳定性好、对挠度和内力都具有精度高的特点．     

《数值计算与计算机应用》2021年第42卷第4期摘要

潘佳佳,李会元,二阶椭圆问题的弱迦辽金四边形谱元方法[J].数值计算与计算机应用,2021,42(4):303-322.摘要:本文对二阶椭圆方程特征值问题的弱伽辽金谱元方法开展相关数值研究.与弱有限元方法类似,弱伽辽金谱元方法的逼近函数空间包括各个单元上的独立内部分量、并辅以各单元边界分量作为单元与单元间的联系.本文聚焦任意凸四边形网格剖分下的弱伽辽金四边形谱元方法,弱逼近函数中的各内部分量与边界分量分别由参考正方形单元与参考单元边界上的正交多项式通过双线性变换来构造;而弱梯度逼近空间则由参考正方形上的正交多项式通过Piola变换构造.在此基础上,本文提出了二阶椭圆方程特征值问题的弱伽辽金四边形谱元方法逼近格式和实现算法,并通过对离散弱梯度核空间的系统研究。     

Hermite三次样条多小波自然边界元法

针对应用自然边界元方法解上半平面的Laplace方程的Neumann边值问题时存在奇异积分的困难,本文提出了Hermite三次样条多小波自然边界元法.Hermite三次样条多小波具有较短的紧支集、很好的稳定性和显式表达式,而且它们在不同层上的导数还是相互正交的.因此,本文将它与自然边界元法相结合,利用小波伽辽金法离散自然边界积分方程,使自然边界积分方程中的强奇异积分化为弱奇异积分,从而降低了问题的复杂性.文中给出的算例表明了该方法的可行性和有效性.     

Pasternak地基上简支板振动问题的准格林函数方法

提出一种新的数值方法——准格林函数方法．以Pasternak地基上简支多边形薄板的振动问题为例,详细阐明了准格林函数方法的思想．即利用问题的基本解和边界方程构造一个准格林函数,这个函数满足了问题的齐次边界条件,采用格林公式将Pasternak地基上薄板自由振动问题的振型控制微分方程化为两个耦合的第二类Fredholm积分方程．边界方程有多种选择,在选定一种边界方程的基础上,可以通过建立一个新的边界方程来表示问题的边界,以克服积分核的奇异性,最后由积分方程的离散化方程组有非平凡解的条件,求得固有频率．数值方法表明,该方法具有较高的精度．     

新的三维力学GELD正演和反演算法

在本文中 ,我们提出了新的整体积分和局部微分GILD的力学正演和反演方法 .我们建立了弹性和塑性力学的体积分微分方程 .我们证明了这个体积分方程和伽辽金虚功原理等价 .新的GILD方法是基于这个体积分微分方程 .GL方法是进一步的发展 ,GL是一种整体场和局部场相互作用的全新方法 .在这个方法中 ,仅仅需要解 3× 3或者 6 × 6的局部小矩阵 .特别是 ,用GL方法求解无限域的偏微分方程时 ,不需要任何人工边界 ,不需要任何吸收边界条件和不需要任何边界积分方程 .新的三维力学GILD正演和反演算法已被应用研究奈米材料的力学性质的模拟计算 .我们获得非常好的奈米材料的力学变形的超拉力的力学性质 .我们提出了新的奈米地球物理新概念和发现了GILD数值量子     

单裂纹问题的虚边界无网格伽辽金法分析

使用含裂纹复变基本解,虚边界无网格伽辽金法被进一步推广应用于弹性材料的单裂纹问题求解.为了清晰地说明单裂纹问题的虚边界元法实现过程,单裂纹问题的虚边界元法示意图、复变坐标平面下含裂纹问题的复变位移和复变面力基本解示意图被展示.含裂纹复变基本解,因自动满足裂纹处边界条件,故使用虚边界无网格伽辽金法计算单裂纹问题,无需在裂纹处布置节点或单元.给出含裂纹复变基本解中的Φ'(x)的详细表达式、裂纹左右裂尖应力强度因子的虚边界无网格离散公式,方便了其他学者使用本方法计算裂纹问题.数值计算两端受拉长方形钢板中心含有裂纹的应力强度因子的算例,计算结果证明了本方法的精确性与稳定性.     

定常Kuramoto-Sivashinsky方程关于非线性伽辽金方法的一个注记

本文讨论了定常Ｋ－Ｓ方程关于伽辽金方法和非线性伽辽金方法的收敛性和最大模估计;对相同模数而言,两者的误差阶完全一致,数值结果表明非线性伽辽金方法同样成功地计算出了Ｋ－Ｓ方程的分歧解,并且在计算时间方面非线性伽辽金方法比伽辽金方法要少得多。     

Winkler地基上固支薄板自由振动问题的准Green函数方法

将准Green函数方法应用于求解Winkler地基上固支薄板的自由振动问题.即利用问题的基本解和边界方程构造一个准Green函数,这个函数满足了问题的齐次边界条件.采用Green公式,将Winkler地基上固支薄板自由振动问题的振型控制微分方程化为第二类Fredholm积分方程.通过边界方程的适当选择,积分方程核的奇异性被克服了.数值算例表明,该方法具有较高的精度,是一种有效的数学方法.     

改进的无奇异局部边界积分方程方法

在局部边界积分方程方法中,当源节点位于分析域的整体边界上时,局部边界积分将出现奇异积分问题,这些奇异积分需要做特别的处理．为此,提出了对域内节点采用局部积分方程,而对边界节点直接采用移动最小二乘近似函数引入边界条件来解决奇异积分问题,这同时也解决了对积分边界进行插值引入近似误差的问题．作为应用和数值实验,对Laplace方程和Helmholtz方程问题进行了分析,取得了很好的数值结果．进而,在Helmholtz方程求解中,采用了含波解信息的修正基函数来代替单项式基函数进行近似．数值结果显示,这样处理是简单高效的,在高波数声传播问题的求解中非常具有前景．     

An improved boundary element method for the charge density of a thin electrified plate in ℝ3

A weakly singular integral equation of the first kind on a plane surface piece Γ is solved approximately via the Galerkin method. The determination of the solution of this integral equation (with the single-layer potential) is a classical problem in physics, since its solution represents the charge density of a thin, electrified plate Γ loaded with some given potential. Using piecewise constant or piecewise bilinear boundary elements we derive asymptotic estimates for the Galerkin error in the energy norm and analyse the effect of graded meshes. Estimates in lower order Sobolev norms are obtained via the Aubin–Nitsche trick. We describe in detail the numerical implementation of the Galerkin method with both piecewise-constant and piecewise-linear boundary elements. Numerical experiments show experimental rates of convergence that confirm our theoretical, asymptotic results.     

Meshless Galerkin algorithms for boundary integral equations with moving least square approximations

In this paper, we first give error estimates for the moving least square (MLS) approximation in the Hk norm in two dimensions when nodes and weight functions satisfy certain conditions. This two-dimensional error results can be applied to the surface of a three-dimensional domain. Then combining boundary integral equations (BIEs) and the MLS approximation, a meshless Galerkin algorithm, the Galerkin boundary node method (GBNM), is presented. The optimal asymptotic error estimates of the GBNM for three-dimensional BIEs are derived. Finally, taking the Dirichlet problem of Laplace equation as an example, we set up a framework for error estimates of the GBNM for boundary value problems in three dimensions.     

An a posteriori error estimate for a first-kind integral equation

In this paper we present a new a posteriori error estimate for the boundary element method applied to an integral equation of the first kind. The estimate is local and sharp for quasi-uniform meshes and so improves earlier work of ours. The mesh-dependence of the constants is analyzed and shown to be weaker than expected from our previous work. Besides the Galerkin boundary element method, the collocation method and the qualocation method are considered. A numerical example is given involving an adaptive feedback algorithm.

     

A Symmetric Galerkin Boundary Element Method for Linear Poroelasticity

In linear poroelasticity so far only collocation boundary element methods have been available. However, in some applications, e.g., when coupling with finite elements is desired, a symmetric formulation is preferable. Choosing a Galerkin approach which involves the second boundary integral equation, such a formulation is possible. Here, a previously presented integration by part technique for the regularization of the first boundary integral equation is extended to the second boundary integral equation as well. While the weakly singular representation of the double layer operator has been presented before, the emphasis lies here on the so called hyper-singular boundary integral operator. Due to the regularization, this operator can be evaluated numerically and, hence, be used within a numerical scheme for the first time. Different numerical studies will be presented to show the behavior of the established symmetric Galerkin boundary element method, also comparing it with collocation boundary element methods. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analytical formulation of meshless local integral equation method

In this paper, the exact forms of integrals in the meshless local boundary integral equation method are derived and implemented for elastostatic problems. A weak form for a set of governing equations with a unit test function or polynomial test functions is transformed into local integral equations. Each node has its own support domain and is surrounded by a local integral domain with different shapes of boundaries. The meshless approximation based on the radial basis function (RBF) is employed for the implementation of displacements. A completed set of closed forms of the local boundary integrals are obtained. As there are no numerical integrations to be carried out the computational time is significantly reduced. Three examples are presented to demonstrate the application of this approach in solid mechanics.     

Modelling wave propagation phenomena in time domain by using the symmetric Galerkin BEM and the Convolution Quadrature Method

The boundary element method (BEM) is known to be well suited to treat wave propagation phenomena. Here, a 3-d elastodynamic body is under consideration. For a numerical treatment of the underlying boundary integral equations appropriate time and spatial discretizations have to be introduced. The time discretization is done via the Convolution Quadrature Method (CQM) as proposed by Lubich while a symmetric Galerkin scheme is applied in space. Unfortunately, the symmetric Galerkin boundary element method (SGBEM) requires the use of the second boundary integral equation involving a hypersingular kernel. Therefore, an appropriate regularization of the hypersingular bilinear form has to be established. Finally, a numerical study show that the presented approach yields good convergence rates as well as good numerical stability properties. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Adaptive meshless Galerkin boundary node methods for hypersingular integral equations

Adaptive refinement techniques are developed in this paper for the meshless Galerkin boundary node method for hypersingular boundary integral equations. Two types of error estimators are derived. One is a perturbation error estimator that is formulated based on the difference between numerical solutions obtained using two consecutive nodal arrangements. The other is a projection error estimator that is formulated based on the difference between the numerical solution itself and its projection. These error estimators are proven to have an upper and a lower bound by the constant multiples of the exact error in the energy norm. A localization scheme is presented to accomodate the non-local property of hypersingular integral operators for the needed computable local error indicators. The convergence of the adaptive meshless techniques is verified theoretically. To confirm the theoretical results and to show the efficiency of the adaptive techniques, numerical examples in 2D and 3D with high singularities are provided.     

A meshless local Galerkin method for the numerical solution of Hammerstein integral equations based on the moving least squares technique

In this paper, a computational scheme is proposed to estimate the solution of one- and two-dimensional Fredholm-Hammerstein integral equations of the second kind. The method approximates the solution using the discrete Galerkin method based on the moving least squares (MLS) approach as a locally weighted least squares polynomial fitting. The discrete Galerkin technique for integral equations results from the numerical integration of all integrals in the system corresponding to the Galerkin method. Since the proposed method is constructed on a set of scattered points, it does not require any background meshes and so we can call it as the meshless local discrete Galerkin method. The implication of the scheme for solving two-dimensional integral equations is independent of the geometry of the domain. The new method is simple, efficient and more flexible for most classes of nonlinear integral equations. The error analysis of the method is provided. The convergence accuracy of the new technique is tested over several Hammerstein integral equations and obtained results confirm the theoretical error estimates.     

Product Gauss quadrature rules vs. cubature rules in the meshless local Petrov–Galerkin method

A crucial point in the implementation of meshless methods such as the meshless local Petrov–Galerkin (MLPG) method is the evaluation of the domain integrals arising over circles in the discrete local weak form of the governing partial differential equation. In this paper we make a comparison between the product Gauss numerical quadrature rules, which are very popular in the MLPG literature, with cubature formulas specifically constructed for the approximation of an integral over the unit disk, but not yet applied in the MLPG method, namely the spherical, the circularly symmetrical and the symmetric cubature formulas. The same accuracy obtained with 64×64 points in the product Gauss rules may be obtained with symmetric quadrature formulas with very few points.