Quadrature formulas for periodic functions and their application to the boundary element method

Two-dimensional and axisymmetric boundary value problems for the Laplace equation in a domain bounded by a closed smooth contour are considered. The problems are reduced to integral equations with a periodic singular kernel, where the period is equal to the length of the contour. Taking into account the periodicity property, high-order accurate quadrature formulas are applied to the integral operator. As a result, the integral equations are reduced to a system of linear algebraic equations. This substantially simplifies the numerical schemes for solving boundary value problems and considerably improves the accuracy of approximation of the integral operator. The boundaries are specified by analytic functions, and the remainder of the quadrature formulas decreases faster than any power of the integration step size. The examples include the two-dimensional potential inviscid circulation flow past a single blade or a grid of blades; the axisymmetric flow past a torus; and free-surface flow problems, such as wave breakdown, standing waves, and the development of Rayleigh-Taylor instability.     

Modelling two-dimensional flow past arbitrary cylindrical bodies using boundary element formulations

This paper presents an efficient method of solving Queen's linearized equations for steady plane flow of an incompressible, viscous Newtonian fluid past a cylindrical body of arbitrary cross-section. The numerical solution technique is the well known direct boundary element method. Use of a fundamental solution of Oseen's equations, the ‘Oseenlet’, allows the problem to be reduced to boundary integrals and numerical solution then only requires boundary discretization. The formulation and solution method are validated by computing the net forces acting on a single circular cylinder, two equal but separated circular cylinders and a single elliptic cylinder, and comparing these with other published results. A boundary element representation of the full Navier-Stokes equations is also used to evaluate the drag acting on a single circular cylinder by matching with the numerical Oseen solution in the far field.     

Approximation of boundary values in the boundary element method

Various techniques may be applied to the approximation of the unknown boundary functions involved in the boundary element method (BEM). Several techniques have been examined numerically to find the most efficient. Techniques considered were: Lagrangian polynomials of the zeroth, first and second orders; spline functions; and the novel weighted minimization technique used successfully in the finite difference method (FDM) for arbitrarily irregular meshes. All these approaches have been used in the BEM for the numerical analysis of plates with various boundary conditions.Both coarse and fine grids on the boundary have been assumed. Maximal errors of the deflections of each plate and the bending moments have been found and the effective computer CPU times determined.Analysis of the results showed that, for the same computer time, the greatest accuracy was obtained by the weighted FDM approach. In the case of the Lagrange approximation, higher order polynomials have proved more efficient. The spline technique yielded more accurate results, but with a higher CPU time.Two discretization approaches have been investigated: the least-squares technique and the collocation method. Despite the fact that the simultaneous algebraic equations obtained were not symmetric, the collocation approach has been confirmed as clearly superior to the least-squares technique, because of the amount of computation time used.     

An extrapolation method for a class of boundary integral equations

Boundary value problems of the third kind are converted into boundary integral equations of the second kind with periodic logarithmic kernels by using Green's formulas. For solving the induced boundary integral equations, a Nyström scheme and its extrapolation method are derived for periodic Fredholm integral equations of the second kind with logarithmic singularity. Asymptotic expansions for the approximate solutions obtained by the Nyström scheme are developed to analyze the extrapolation method. Some computational aspects of the methods are considered, and two numerical examples are given to illustrate the acceleration of convergence.

     

An implementation of stress discontinuity in the boundary element method and application to gear teeth

The boundary element method as it applies to three-dimensional elasto-static problems is implemented using isoparametric quadratic elements. A scheme for resolving the problems of traction vector discontinuity is presented and the required additional equations are derived. Example problems considered, including the stress analysis of a three-dimensional gear tooth, demonstrate that high accuracy may be achieved using a relatively small number of elements if continuity of displacements and discontinuity of tractions are properly implemented.     

A convolution quadrature Galerkin boundary element method for the exterior Neumann problem of the wave equation

The numerical solution of the Neumann problem of the wave equation on unbounded three‐dimensional domains is calculated using the convolution quadrature method for the time discretization and a Galerkin boundary element method for the spatial discretization. The mathematical analysis that has been built up for the Dirichlet problem is extended and developed for the Neumann problem, which is important for many modelling applications. Numerical examples are then presented for one of these applications, modelling transient acoustic radiation. Copyright © 2009 John Wiley & Sons, Ltd.     

The boundary integral equation method applied to shallow membrane shells of positive Gaussian curvature

The boundary integral equation method (BIEM) is developed for the analysis of shallow membrane shells with positive Gaussian curvatures. Shells with constant thickness and constant curvatures are considered. In the infinite domain, fundamental solutions are obtained which correspond to generalized concentrated tangential forces in the x and y coordinate directions. The Betti-Maxwell reciprocal theorem and Green's second identity are used to obtain the boundary integral equations of the solution presented.This approach, which is applied for the first time in membrane shell theory, seems to be a powerful alternative to domain type methods. Shells with various boundary conditions, loadings and arbitrary plan forms can be considered. It is also possible to add the effects of thermal fields and openings in the shells.The potential of the method is demonstrated by means of a worked example.     

Nonlocal overdetermined boundary value problem for stationary Navier-Stokes equations

A stationary system of Stokes and Navier-Stokes equations describing the flow of a homogeneous incompressible fluid in a bounded domain is considered. The vector of the flow velocity and a finite number of nonlocal conditions are defined at a part of the domain boundary. It is proved that, in the linear case, the problem has at least one stable solution. In the nonlinear case, the local solvability of the problem is proved.     

An optimal asymptotic-numerical method for convection dominated systems having exponential boundary layers

The convection dominated diffusion problems are studied. Higher order accurate numerical methods are presented for problems in one and two dimensions. The underlying technique utilizes a superposition of given problem into two independent problems. The first one is the reduced problem that refers to the outer or smooth solution. Stretching transformation is used to obtain the second problem for inner layer solution. The method considered for outer or degenerate problems are based on higher order Runge–Kutta methods and upwind finite differences. However, inner problem is solved analytically or asymptotically. The schemes presented are proved to be consistent and stable. Possible extensions to delay differential equations and to nonlinear problems are outlined. Numerical results for several test examples are illustrated and a comparative analysis is presented. It is observed that the method presented is highly accurate and easy to implement. Moreover, the numerical results obtained are not only comparable with the exact solution but also in agreement with the theoretical estimates.     

A new boundary element method for the solution of plane steady-state thermoelastic fracture mechanics problems

An efficient numerical technique is developed for plane, homogeneous, isotropic, steady-state thermoelasticity problems involving arbitrary internal smooth and/or kinkedcracks. The thermal stress intensity factors and relative crack surface displacements due to steady-state temperature distributions are determined and compared to available solutions obtained by other methods. In these analyses the thermal boundary conditions across the crack surface are assumed to be insulated. The present approach involves coupling the direct boundary integral equations to newly developed crack integral equations.     

The boundary integral method for the linearized rotating Navier-Stokes equations in exterior domain

In this paper, we apply the boundary integral method to the linearized rotating Navier-Stokes equations in exterior domain. Introducing some open ball which decomposes the exterior domain into a finite domain and an infinite domain, we obtain a coupled problem by the linearized rotating Navier-Stokes equations in finite domain and a boundary integral equation without using the artificial boundary condition. For the coupled problem, we show the existence and uniqueness of solution. Finally, we study the finite element approximation for the coupled problem and obtain the error estimate between the solution of the coupled problem and its approximation solution.     

Remarks concerning the boundary element method in potential theory

The BEM has become a well-known tool in the numerical treatment of potential and elesticity problems. In this paper some modifications concerning the treatment of singularities and of Neumann problems are presented, which lead to improved accuracy or reduced computing time.     

Development of boundary element method to dynamic problems for porous media

Biot's consolidation theory is extended to a general class of viscoelastic bodies defined by Riemann-Stieltjes integral convolutions. From a new reciprocity theorem, proved for the governing equations including the inertia terms, the basic integral representations of the displacement fields and pore pressure are obtained. It is shown that, in the absence of internal inputs, a formulation of the dynamic problem in terms of the boundary unknown fields only is possible.     

Numerical solution of thermal convection problems using the multidomain boundary element method

The multidomain dual reciprocity method (MD‐DRM) has been effectively applied to the solution of two‐dimensional thermal convection problems where the momentum and energy equations govern the motion of a viscous fluid. In the proposed boundary integral method the domain integrals are transformed into equivalent boundary integrals by the dual reciprocity approach applied in a subdomain basis. On each subregion or domain element the integral representation formulas for the velocity and temperature are applied and discretised using linear continuous boundary elements, and the equations from adjacent subregions are matched by additional continuity conditions. Some examples showing the accuracy, the efficiency and flexibility of the proposed method are presented. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 469–489, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10016     

The solution of the non-homogeneous Helmholtz equation by means of the boundary element method

One of the most important advantages of the Boundary Element Method (BEM) is that no internal discretization of the domain is required. This advantage, however, is generally lost when source terms are present in the governing differential equation. It is shown here that for the non-homogeneous Helmholtz equation with a harmonic source term, it is possible to transform the volume integral into a surface integral thus retaining this feature. The transformation is achieved using the Green formula. The technique is applied to solve numerically a test problem with known simple analytical solution.     

Energetic boundary element method analysis of wave propagation in 2D multilayered media

The analysis of scalar wave propagation in 2D zonewise homogeneous media with vanishing initial and mixed boundary conditions is carried out. The problem is formulated in terms of time‐dependent boundary integral equations, and then it is set in a weak form, based on a natural energy identity satisfied by the differential problem solution. Several numerical results have been obtained by means of the related energetic Galerkin boundary element method showing accuracy and stability of the method. Copyright © 2012 John Wiley & Sons, Ltd.