On the effect of numerical integration in the Galerkin boundary element method

Summary. We discuss the effect of cubature errors when using the Galerkin method for approximating the solution of Fredholm integral equations in three dimensions. The accuracy of the cubature method has to be chosen such that the error resulting from this further discretization does not increase the asymptotic discretization error. We will show that the asymptotic accuracy is not influenced provided that polynomials of a certain degree are integrated exactly by the cubature method. This is done by applying the Bramble-Hilbert Lemma to the boundary element method. Received May 24, 1995     

Sobolev方程的时间间断Galerkin有限元方法

引入Sobolev方程的等价积分方程,构造Sobolev方程的新的时间间断Galerkin有限元格式.该格式不仅保持有限元解在时间剖分点处的间断特性,而且避免了传统时空有限元格式中跳跃项的出现,从而降低了格式理论分析和数值模拟的复杂性.证明了Sobolev方程的时间间断而空间连续的时空有限元解的稳定性、存在唯一性、L2...     

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.     

Time discretisation of monotone nonlinear evolution problems by the discontinuous Galerkin method

A class of discontinuous Galerkin methods is studied for the time discretisation of the initial-value problem for a nonlinear first-order evolution equation that is governed by a monotone, coercive, and hemicontinuous operator. The numerical solution is shown to converge towards the weak solution of the original problem. Furthermore, well-posedness of the time-discrete problem as well as a priori error estimates for sufficiently smooth exact solutions are studied.     

A comparison of numerical integration rules for the meshless local Petrov–Galerkin method

The meshless local Petrov–Galerkin (MLPG) method is a mesh-free procedure for solving partial differential equations. However, the benefit in avoiding the mesh construction and refinement is counterbalanced by the use of complicated non polynomial shape functions with subsequent difficulties, and a potentially large cost, when implementing numerical integration schemes. In this paper we describe and compare some numerical quadrature rules with the aim at preserving the MLPG solution accuracy and at the same time reducing its computational cost.     

A note on the convergence of product integration and Galerkin method for weakly singular integral equations

We consider weakly singular integral equations of Fredholm-type whose kernels satisfy certain algebraic estimates with their derivatives. In particular, we establish optimal convergence order estimates for product integration and Galerkin method applied on suitable grading mesh for the solution of such equations. Some superconvergence results are also derived.     

The re-iterated Galerkin method

The iterated Galerkin method introduced by Sloan is modifiedto allow the improvement in the Galerkin approximation whichiteration produces to be repeated. It is shown in principle,and verified in illustration, that the repeated improvementsmay be achieved under weaker conditions than those requiredby Sloan. The sequence of approximations obtained is appliedto variational principles and in the case when a maximum principleis available it is shown that each step in the process improvesthe estimate of the maximum value of the associated functional.Applications to integral equations are given.     

Interior estimates for the wavelet Galerkin method

Summary. In this paper we derive an interior estimate for the Galerkin method with wavelet-type basis. Such an estimate follows from interior Galerkin equations which are common to a class of methods used in the solution of elliptic boundary value problems. We show that the error in an interior domain can be estimated with the best order of accuracy possible, provided the solution is sufficiently regular in a slightly larger domain, and that an estimate of the same order exists for the error in a weaker norm (measuring the effects from outside the domain ). Examples of the application of such an estimate are given for different problems. Received May 17, 1995 / Revised version received April 26, 1996     

A multilevel discontinuous Galerkin method

A variable V-cycle preconditioner for an interior penalty finite element discretization for elliptic problems is presented. An analysis under a mild regularity assumption shows that the preconditioner is uniform. The interior penalty method is then combined with a discontinuous Galerkin scheme to arrive at a discretization scheme for an advection-diffusion problem, for which an error estimate is proved. A multigrid algorithm for this method is presented, and numerical experiments indicating its robustness with respect to diffusion coefficient are reported.Received: 5 June 2001, Revised: 12 December 2001, Published online: 4 April 2003Mathematics Subject Classification (1991): 65F10, 65N55, 65N30This research was supported in part by Institute for Mathematics and its Applications, Supercomputing Institute of University of Minnesota, and Deutsche Forschungsgemeinschaft.     

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.     

Galerkin method for the linear-radiator integral equation

The article examines the general Galerkin method scheme for the singular integral equation in the problem of radiation from a finite-thickness linear radiator. Translated from Obratnye Zadachi Estestvoznaniya, Published by Moscow University, Moscow, 1997, pp. 150–158.     

Justification of the Galerkin method for hypersingular equations

The paper presents a theoretical study of hypersingular equations of the general form for problems of electromagnetic-wave diffraction on open surfaces of revolution. Justification of the Galerkin is given. The method is based on the separation of the principal term and its analytic inversion. The inverse of the principal operator is completely continuous. On the basis of this result, the equivalence of the initial equation to a Fredholm integral equation of the second kind is proven. An example of numerical solution with the use of Chebyshev polynomials of the second kind is considered.     

Some details of the Galerkin finite element method

The incorporation of the Galerkin technique in the finite element method has removed the constraint of finding a variational formulation for many problems of mathematical physics. The method has been successfully applied to many areas and has received wide acceptance. However, in the process of transplanting the concept from the Galerkin method for the entire domain to the Galerkin finite element method, some formal details have been overlooked or glossed over in the literature. This paper considers some of these details, including a possible reason for integration by parts and the contribution of interelement discontinuity terms.     

Nonlinear Galerkin method and two-step method for the Navier-Stokes equations

This article represents a new nonlinear Galerkin scheme for the Navier-Stokes equations. This scheme consists of a nonlinear Galerkin finite element method and a two-step difference method. Moreover, we also provide a Galerkin scheme. By convergence analysis, two numerical schemes have the same second-order convergence accuracy for the spatial discretization and time discretization if H is chosen such that H = O(h^2/3). However, the nonlinear Galerkin scheme is simpler than the Galerkin scheme, namely, this scheme can save a large amount of computational time. © 1996 John Wiley & Sons, Inc.     

Observations on the numerical stability of the Galerkin method

A simple factorization of the finite-dimensional Galerkin operators motivates a study of the numerical stability of a Galerkin procedure on the basis of its “potential stability” and the “conditioning” of its coordinate functions. Conditions sufficient for stability and conditions leading to instability are thereby identified. Numerical examples of stability and instability occurring in the application of the Galerkin method to boundary-integral equations arising in simple scattering problems are provided and discussed within this framework. Numerical instabilities reported by other authors are examined and explained from the same point of view. This revised version was published online in June 2006 with corrections to the Cover Date.     

A discontinuous Galerkin method for the Cahn-Hilliard equation

The Cahn-Hilliard equation is modeled to describe the dynamics of phase separation in glass and polymer systems. A priori error estimates for the Cahn-Hilliard equation have been studied by the authors. In order to control accuracy of approximate solutions, a posteriori error estimation of the Cahn-Hilliard equation is obtained by discontinuous Galerkin method.     

A non-linear version of the amir-lindenstrauss method

A non-linear version of the Amir-Lindenstrauss method of projections for weakly compactly generated Banach spaces is proved, that implies immediately the Benyamini-Rudin-Wage result on continuous images of Eberlein compact spaces.     

The Solution of Time Dependent Problems by Galerkin Methods

Evolutionary problems involving parabolic and hyperbolic equationsin any number of space dimensions are formulated in such a waythat Galerkin methods, using a local basis, can be used to findapproximate solutions. These methods are particularly suitedto problems involving general regions and derivative boundaryconditions.