Discrete boundary element methods on general meshes in 3D

Abstract. This paper is concerned with the stability and convergence of fully discrete Galerkin methods for boundary integral equations on bounded piecewise smooth surfaces in . Our theory covers equations with very general operators, provided the associated weak form is bounded and elliptic on , for some . In contrast to other studies on this topic, we do not assume our meshes to be quasiuniform, and therefore the analysis admits locally refined meshes. To achieve such generality, standard inverse estimates for the quasiuniform case are replaced by appropriate generalised estimates which hold even in the locally refined case. Since the approximation of singular integrals on or near the diagonal of the Galerkin matrix has been well-analysed previously, this paper deals only with errors in the integration of the nearly singular and smooth Galerkin integrals which comprise the dominant part of the matrix. Our results show how accurate the quadrature rules must be in order that the resulting discrete Galerkin method enjoys the same stability properties and convergence rates as the true Galerkin method. Although this study considers only continuous piecewise linear basis functions on triangles, our approach is not restricted in principle to this case. As an example, the theory is applied here to conventional “triangle-based” quadrature rules which are commonly used in practice. A subsequent paper [14] introduces a new and much more efficient “node-based” approach and analyses it using the results of the present paper. Received December 10, 1997 / Revised version received May 26, 1999 / Published online April 20, 2000 –? Springer-Verlag 2000     

Hypersingular kernel integration in 3D Galerkin boundary element method

We consider Cauchy singular and Hypersingular boundary integral equations associated with 3D potential problems defined on polygonal domains, whose solutions are approximated with a Galerkin boundary element method, related to a given triangulation of the boundary. In particular, for constant and linear shape functions, the most frequently used basis functions, we give explicit results of the analytical inner integrations and suggest suitable quadrature schemes to evaluate the outer integrals required to form the Galerkin matrix elements. These numerical indications are given after an analysis of the singularities arising in the whole integration process, which is valid also for shape functions of higher degrees.     

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     

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.     

A primal mixed formulation for exterior transmission problems in ${\bf R}^2$

Summary.   We combine a primal mixed finite element approach with a Dirichlet-to-Neumann mapping (arising from the boundary integral equation method) to study the weak solvability and Galerkin approximations of a class of linear exterior transmission problems in potential theory. Our results are mainly based on the Babuska-Brezzi theory for variational problems with constraints. We establish the uniqueness of solution for the continuous and discrete formulations, and show that finite element subspac es of Lagrange type satisfy the discrete compatibility conditions. In addition, we provide the error analysis, including polygonal approximations of the domain, and prove strong convergence of the Galerkin solutions. Moreover, under additional regularity assumptions on the solution of the continuous formulation, we obtain the asymptotic rate of convergence O(h). Received August 25, 1998 / Revised version received March 8, 2000 / Published online October 16, 2000     

Discrete Galerkin Methods for Fredholm Integral Equations of the Second Kind

The discrete Galerkin and discrete iterated Galerkin methodsarise when the integrals required in the Galerkin and iteratedGalerkin methods are calculated using numerical integration.In this paper, prolongation and restriction operators are usedto give an error analysis for these two discrete Galerkin methods.From this analysis, we can then give conditions on the quadratureerrors, under which the two discrete Galerkin solutions havethe same order of convergence as their exact counterparts.     

Analysis of fully discrete approximations to parabolic problems with rough or distribution-valued initial data

Summary. We consider fully discrete approximations to a parabolic initial-boundary value problem with rough or distribution-valued initial data in two space dimensions. For discretization in time and space, we apply single step methods and the standard Galerkin method with piecewise linear test functions, respectively. For spatial discretization of the initial condition, we are however forced to use more involved constructions. Our main result is stability and error estimates of the discrete solutions. Received October 21, 1999 / Revised version received May 3, 2001 / Published online December 18, 2001     

Variational time discretizations of higher order and higher regularity

We consider a family of variational time discretizations that are generalizations of discontinuous Galerkin (dG) and continuous Galerkin–Petrov (cGP) methods. The family is characterized by two parameters. One describes the polynomial ansatz order while the other one is associated with the global smoothness that is ensured by higher order collocation conditions at both ends of the subintervals. Applied to Dahlquist’s stability problem, the presented methods provide the same stability properties as dG or cGP methods. Provided that suitable quadrature rules of Hermite type are used to evaluate the integrals in the variational conditions, the variational time discretization methods are connected to special collocation methods. For this case, we present error estimates, numerical experiments, and a computationally cheap postprocessing that allows to increase both the accuracy and the global smoothness by one order.

     

A novel asymptotic extraction technique for the efficient evaluation of a class of double Sommerfeld integrals

When the Galerkin method of moments (MoM) is used in the numerical analysis of microwave microstrip networks, elements of a network impedance matrix are required. Two dimensional (2D) Sommerfeld (spectral) integrals must be evaluated in such cases. When semi-infinite circuit feed lines are used, certain specific types of Sommerfeld integrals arise whose integrands oscillate rapidly and present serious numerical difficulties for their evaluation. In this paper, a novel, accurate, and much more efficient asymptotic extraction technique (AET) is developed that involves extracting an inner asymptotic angular integral, which is then evaluated analytically before evaluating the outer integral numerically.     

Stochastic discrete Hamiltonian variational integrators

Variational integrators are derived for structure-preserving simulation of stochastic Hamiltonian systems with a certain type of multiplicative noise arising in geometric mechanics. The derivation is based on a stochastic discrete Hamiltonian which approximates a type-II stochastic generating function for the stochastic flow of the Hamiltonian system. The generating function is obtained by introducing an appropriate stochastic action functional and its corresponding variational principle. Our approach permits to recast in a unified framework a number of integrators previously studied in the literature, and presents a general methodology to derive new structure-preserving numerical schemes. The resulting integrators are symplectic; they preserve integrals of motion related to Lie group symmetries; and they include stochastic symplectic Runge–Kutta methods as a special case. Several new low-stage stochastic symplectic methods of mean-square order 1.0 derived using this approach are presented and tested numerically to demonstrate their superior long-time numerical stability and energy behavior compared to nonsymplectic methods.     

A Galerkin method for retarded boundary integral equations with smooth and compactly supported temporal basis functions

We consider retarded boundary integral formulations of the three-dimensional wave equation in unbounded domains. Our goal is to apply a Galerkin method in space and time in order to solve these problems numerically. In this approach the computation of the system matrix entries is the major bottleneck. We will propose new types of finite-dimensional spaces for the time discretization. They allow variable time-stepping, variable order of approximation and simplify the quadrature problem arising in the generation of the system matrix substantially. The reason is that the basis functions of these spaces are globally smooth and compactly supported. In order to perform numerical tests concerning our new basis functions we consider the special case that the boundary of the scattering problem is the unit sphere. In this case explicit solutions of the problem are available which will serve as reference solutions for the numerical experiments.     

Full asymptotics of spline Petrov–Galerkin methods for some periodic pseudodifferential equations

In this paper we study the existence of a formal series expansion of the error of spline Petrov–Galerkin methods applied to a class of periodic pseudodifferential equations. From this expansion we derive some new superconvergence results as well as alternative proofs of already known weak norm optimal convergence results. As part of the analysis the approximation of integrals of smooth functions multiplied by splines by rectangular rules is analyzed in detail. Finally, some numerical experiments are given to illustrate the applicability of Richardson extrapolation as a means of accelerating the convergence of the methods.     

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.     

Nonconforming Galerkin methods for the Helmholtz equation

Nonconforming Galerkin methods for a Helmholtz‐like problem arising in seismology are discussed both for standard simplicial linear elements and for several new rectangular elements related to bilinear or trilinear elements. Optimal order error estimates in a broken energy norm are derived for all elements and in L² for some of the elements when proper quadrature rules are applied to the absorbing boundary condition. Domain decomposition iterative procedures are introduced for the nonconforming methods, and their convergence at a predictable rate is established. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 475–494, 2001     

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.     

A Petrov–Galerkin RBF method for diffusion equation on the unit sphere

This paper concerns a numerical solution for the diffusion equation on the unit sphere. The given method is based on the spherical basis function approximation and the Petrov–Galerkin test discretization. The method is meshless because spherical triangulation is not required neither for approximation nor for numerical integration. This feature is achieved through the spherical basis function approximation and the use of local weak forms instead of a global variational formulation. The local Petrov–Galerkin formulation allows to compute the integrals on small independent spherical caps without any dependence on a connected background mesh. Experimental results show the accuracy and the efficiency of the new method.     

Semi-analytic integration of hypersingular Galerkin BIEs for three-dimensional potential problems

An accurate and efficient semi-analytic integration technique is developed for three-dimensional hypersingular boundary integral equations of potential theory. Investigated in the context of a Galerkin approach, surface integrals are defined as limits to the boundary and linear surface elements are employed to approximate the geometry and field variables on the boundary. In the inner integration procedure, all singular and non-singular integrals over a triangular boundary element are expressed exactly as analytic formulae over the edges of the integration triangle. In the outer integration scheme, closed-form expressions are obtained for the coincident case, wherein the divergent terms are identified explicitly and are shown to cancel with corresponding terms from the edge-adjacent case. The remaining surface integrals, containing only weak singularities, are carried out successfully by use of standard numerical cubatures. Sample problems are included to illustrate the performance and validity of the proposed algorithm.     

Quadrature Formulas for Refinable Functions and Wavelets II: Error Analysis

This paper is concerned with the construction and the analysis of Gauss quadrature formulas for computing integrals of (smooth) functions against refinable functions and wavelets. The main goal of this paper is to develop rigorous error estimates for these formulas. For the univariate setting, we derive asymptotic error bounds for a huge class of weight functions including spline functions. We also discuss multivariate quadrature rules and present error estimates for specific nonseparable refinable functions, i.e., for some special box splines.     

Stability and convergence of optimum spectral non‐linear Galerkin methods

Our objective in this article is to present some numerical schemes for the approximation of the 2‐D Navier–Stokes equations with periodic boundary conditions, and to study the stability and convergence of the schemes. Spatial discretization can be performed by either the spectral Galerkin method or the optimum spectral non‐linear Galerkin method; time discretization is done by the Euler scheme and a two‐step scheme. Our results show that under the same convergence rate the optimum spectral non‐linear Galerkin method is superior to the usual Galerkin methods. Finally, numerical example is provided and supports our results. Copyright © 2001 John Wiley & Sons, Ltd.     

椭圆PDE-约束优化问题的一个预条件子

针对由Galerkin有限元离散椭圆PDE-约束优化问题产生的具有特殊结构的3×3块线性鞍点系统,提出了一个预条件子并给出了预处理矩阵特征值及特征向量的具体表达形式.数值结果表明了该预条件子能够有效地加速Krylov子空间方法的收敛速率,同时也验证了理论结果.