A discrete Galerkin method for a hyperslngular boundary integral equation

Received on 14 August 1995. Revised on 20 August 1996. Consider solving the interior Neumann problem with D a simply-connected planar region and S=D a smooth curve.A double-layer potential is used to represent the solution,and it leads to the problem of solving a hypersingular integralequation. This integral equation is reformulated as a Cauchysingular integral equation. A discrete Galerkin method withtrigonometric polynomials is then given for its solution. Anerror analysis is given, and numerical examples complete thepaper.     

A boundary element Galerkin method for a hypersingular integral equation on open surfaces

A hypersingular boundary integral equation of the first kind on an open surface piece Γ is solved approximately using the Galerkin method. As boundary elements on rectangles we use continuous, piecewise bilinear functions which vanish on the boundary of Γ. We show how to compensate for the effect of the edge and corner singularities of the true solution of the integral equation by using an appropriately graded mesh and obtain the same convergence rate as for the case of a smooth solution. We also derive asymptotic error estimates in lower-order Sobolev norms via the Aubin–Nitsche trick. Numerical experiments for the Galerkin method with piecewise linear functions on triangles demonstrate the effect of graded meshes and show experimental rates of convergence which underline the theoretical results.     

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.     

Error estimates for a discontinuous Galerkin method for elliptic problems

In this paper, we consider an elliptic problem with the homogeneous Dirichlet boundary condition and introduce discontinuous Galerkin approximations of the problem. Optimal error estimates of discontinuous Galerkin approximations are obtained.     

Error estimates for a Galerkin numerical scheme applied to a generalized BBM equation

Herein we obtain error estimates for a new Galerkin spectral scheme to approximate the solutions of a dispersive-type equation which is a model to describe the propagation of a wave on the surface of a channel with a flat bottom.     

Energy estimates for Galerkin semidiscretizations of time domain boundary integral equations

In this paper we present a battery of results related to how Galerkin semidiscretization in space affects some formulations of wave scattering and propagation problems when retarded boundary integral equations are used.     

Error estimates for discontinuous Galerkin method for nonlinear parabolic equations

We consider the nonlinear parabolic partial differential equations. We construct a discontinuous Galerkin approximation using a penalty term and obtain an optimal L^∞(L²) error estimate.     

Error estimates for a projection-difference method for a linear differential-operator equation

We study a projection-difference method for approximately solving the Cauchy problem u′(t) + A(t)u(t) + K(t)u(t) = h(t), u(0) = 0 for a linear differential-operator equation in a Hilbert space, where A(t) is a self-adjoint operator and K(t) is an operator subordinate to A(t). Time discretization is based on a three-level difference scheme, and space discretization is carried out by the Galerkin method. Under certain smoothness conditions on the function h(t), we obtain estimates for the convergence rate of the approximate solutions to the exact solution.     

Error estimates for a Galerkin method for a class of model equations for long waves

The Galerkin method, together with a second order time discretization, is applied to the periodic initial value problem for $$\frac{\partial }{{\partial t}}(u - (a(x)u_x )_x ) + (f(x,u))_x = 0$$ . Heref(x, ·) may be highly nonlinear, but a certain cancellation effect is assumed for∫f(x, u) _x u. Optimal order error estimates inL ₂,H ₁, andL _∞ are derived for a general class of piecewise polynomial spaces.     

Collocation method for the natural boundary integral equation

In this work, we construct the Legendre wavelet and apply it to investigate the numerical solution of the natural boundary integral equation of the Laplace equation in the upper half-plane by the collocation method. In our algorithm the coefficient matrix of the linear algebraic system is sparse when the order of the matrix is large. Two test examples show that our algorithm yields very accurate results at less computational cost.     

Error estimates for the Galerkin method as applied to time-dependent equations

A projection method is studied as applied to the Cauchy problem for an operator-differential equation with a non-self-adjoint operator. The operator is assumed to be sufficiently smooth. The linear spans of eigenelements of a self-adjoint operator are used as projection subspaces. New asymptotic estimates for the convergence rate of approximate solutions and their derivatives are obtained. The method is applied to initial-boundary value problems for parabolic equations.     

A fast numerical method for a natural boundary integral equation for the Helmholtz equation

A Neumann boundary value problem of the Helmholtz equation in the exterior circular domain is reduced into an equivalent natural boundary integral equation. Using our trigonometric wavelets and the Galerkin method, the obtained stiffness matrix is symmetrical and circulant, which lead us to a fast numerical method based on fast Fourier transform. Furthermore, we do not need to compute the entries of the stiffness matrix. Especially, our method is also efficient when the wave number k in the Helmholtz equation is very large.     

Error estimates for the third order explicit Runge-Kutta discontinuous Galerkin method for a linear hyperbolic equation in one-dimension with discontinuous initial data

In this paper we present an error estimate for the explicit Runge-Kutta discontinuous Galerkin method to solve a linear hyperbolic equation in one dimension with discontinuous but piecewise smooth initial data. The discontinuous finite element space is made up of piecewise polynomials of arbitrary degree $k\ge 1$ , and time is advanced by the third order explicit total variation diminishing Runge-Kutta method under the standard CFL temporal-spatial condition. The $L^2(\mathbb R \backslash \mathcal R _T)$ -norm error at the final time $T$ is optimal in both space and time, where $\mathcal R _T$ is the pollution region due to the initial discontinuity with the width $\mathcal O (\sqrt{T\beta }h^{1/2}\log (1/h))$ . Here $h$ is the maximum cell length and $\beta $ is the flowing speed. These results are independent of the time step and hold also for the semi-discrete discontinuous Galerkin method.     

A boundary element method for the solution of the transient diffusion equation in two dimensions

The main advantage of using potential or boundary element methods for generating numerical solutions to linear elliptic partial differential equations is that the dimension of the problem is reduced by one. This feature not only reduces the number of algebraic equations to be solved, but, in two and three space dimensions, it also eliminates the problems of mesh generation. Both of these factors contribute to reducing the computational cost. A potential method is described which provides numerical solutions of the transient diffusion equation for two-dimensional problems. The solutions are generally as accurate as conventional techniques whilst retaining the advantage of a reduction in dimensionality throughout the time domain. The numerical study, however, also indicates that the method may have limitations, particularly when there are discontinuities at the boundaries.     

A multiscale Galerkin method for the hypersingular integral equation reduced by the harmonic equation

The aim of this paper is to investigate the numerical solution of the hypersingular integral equation reduced by the harmonic equation. First, we transform the hypersingular integral equation into 2π-periodic hypersingular integral equation with the map x=cot(θ/2). Second, we initiate the study of the multiscale Galerkin method for the 2π-periodic hypersingular integral equation. The trigonometric wavelets are used as trial functions. Consequently, the 2j+1 × 2j+1 stiffness matrix Kj can be partitioned j×j block matrices. Furthermore, these block matrices are zeros except main diagonal block matrices. These main diagonal block matrices are symmetrical and circulant matrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform and the inverse fast Fourier transform instead of the inverse matrix. Finally, we provide several numerical examples to demonstrate our method has good accuracy even though the exact solutions are multi-peak and almost singular.     

Error estimates of a Galerkin method for some nonlinear Sobolev equations in one space dimension

Summary A semidiscrete Galerkin finite element method is defined and analyzed for nonlinear evolution equations of Sobolev type in a single space variable. Optimal orderL ^p error estimates are derived for 2p. And it is shown that the rates of convergence of the approximate solution and its derivative are one order better than the optimal order at certain spatial Jacobi and Gauss points, respectively. Also the standard nodal superconvergence results are established. Futher, it is considered that an a posteriori procedure provides superconvergent approximations at the knots for the spatial derivatives of the exact solution.     

A direct boundary integral equation method for transmission problems

A system of integral equations for the field and its normal derivative on the boundary in acoustic or potential scattering by a penetrable homogeneous object in arbitrary dimensions is presented. The system contains the operators of the single and double layer potentials, of the normal derivative of the single layer, and of the normal derivative of the double layer potential. It defines a strongly elliptic system of pseudodifferential operators. It is shown by the method of Mellin transformation that a corresponding property, namely a Gårding's inequality in the energy norm, holds also in the case of a polygonal boundary of a plane domain. This yields asymptotic quasioptimal error estimates in Sobolev spaces for the corresponding Galerkin approximation using finite elements on the boundary only.     

A direct boundary integral equation method for transmission problems

A system of integral equations for the field and its normal derivative on the boundary in acoustic or potential scattering by a penetrable homogeneous object in arbitrary dimensions is presented. The system contains the operators of the single and double layer potentials, of the normal derivative of the single layer, and of the normal derivative of the double layer potential. It defines a strongly elliptic system of pseudodifferential operators. It is shown by the method of Mellin transformation that a corresponding property, namely a Gårding's inequality in the energy norm, holds also in the case of a polygonal boundary of a plane domain. This yields asymptotic quasioptimal error estimates in Sobolev spaces for the corresponding Galerkin approximation using finite elements on the boundary only.