Regularization Techniques Applied to Boundary Element Methods

1994 ◽  
Vol 47 (10) ◽  
pp. 457-499 ◽  
Author(s):  
Masataka Tanaka ◽  
Vladimir Sladek ◽  
Jan Sladek
Keyword(s):  
Boundary Element ◽  
Integral Equations ◽  
Boundary Integral Equations ◽  
Scalar Potential ◽  
Fundamental Solutions ◽  
Boundary Integral ◽  
Boundary Element Methods ◽  
Element Formulation ◽  
Hypersingular Integrals ◽  
Regularization Techniques

This review article deals with the regularization of the boundary element formulations for solution of boundary value problems of continuum mechanics. These formulations may be singular owing to the use of two-point singular fundamental solutions. When the physical interpretation is irrelevant for this topic of computational mechanics, we consider various mechanical problems simultaneously within particular sections selected according to the main topic. In spite of such a structure of the paper, applications of the regularization techniques to many mechanical problems are described. There are distinguished two main groups of regularization techniques according to their application to singular formulations either before or after the discretization. Further subclassification of each group is made with respect to basic principles employed in individual regularization techniques. This paper summarizes the substances of the regularization procedures which are illustrated on the boundary element formulation for a scalar potential field. We discuss the regularizations of both the strongly singular and hypersingular integrals, occurring in the boundary integral equations, as well as those of nearly singular and nearly hypersingular integrals arising when the source point is near the integration element (as compared to its size) but not on this element. The possible dimensional inconsistency (or scale dependence of results) of the regularization after discretization is pointed out. Finally, we discuss the numerical approximations in various boundary element formulations, as well as the implementations of solutions of some problems for which derivative boundary integral equations are required.

Dual Boundary Element Analysis for Time-Dependent Fracture Problems in Creeping Materials

2008 ◽  
Vol 383 ◽  
pp. 109-121 ◽  
Author(s):  
E. Pineda ◽  
M.H. Aliabadi
Keyword(s):  
Boundary Element ◽  
Integral Equations ◽  
Boundary Integral Equations ◽  
Creep Behaviour ◽  
Boundary Integral ◽  
Element Formulation ◽  
Secondary Creep ◽  
Element Analysis ◽  
Power Law Creep ◽  
Element Method

This paper presents the development of a new boundary element formulation for analysis of fracture problems in creeping materials. For the creep crack analysis the Dual Boundary Element Method (DBEM), which contains two independent integral equations, was formulated. The implementation of creep strain in the formulation is achieved through domain integrals in both boundary integral equations. The domain, where the creep phenomena takes place, is discretized into quadratic quadrilateral continuous and discontinuous cells. The creep analysis is applied to metals with secondary creep behaviour. This is con…ned to standard power law creep equations. Constant applied loads are used to demonstrate time e¤ects. Numerical results are compared with solutions obtained from the Finite Element Method (FEM) and others reported in the literature.

A Higher-Order Poly-Region Boundary Element Method for Steady Thermoviscous Fluid Flows

2004 ◽  
Author(s):  
M. . M. Grigoriev ◽  
G. F. Dargush
Keyword(s):  
Boundary Element ◽  
Boundary Integral Equations ◽  
Fluid Flows ◽  
Boundary Elements ◽  
Higher Order ◽  
Boundary Integral ◽  
Iterative Solvers ◽  
Boundary Element Methods ◽  
Element Formulation ◽  
Thermoviscous Fluid

In this presentation, we re-visit the poly-region boundary element methods (BEM) proposed earlier for the steady Navier-Stokes [1] and Boussinesq [2] flows, and develop a novel higher-order BEM formulation for the thermoviscous fluid flows that involves the definition of the domains of kernel influences due to steady Oseenlets. We introduce region-by-region implementation of the steady-state Oseenlets within the poly-region boundary element fequatramework, and perform integration only over the (parts of) higher-order boundary elements and volume cells that are influenced by the kernels. No integration outside the domains of the kernel influences are needed. Owing to the properties of the convective Oseenlets, the kernel influences are very local and propagate upstream. The localization becomes more prominent as the Reynolds number of the flow increases. This improves the conditioning of the global matrix, which in turn, facilitates an efficient use of the iterative solvers for the sparse matrices [3]. Here, we consider quartic boundary elements and bi-quartic volume cells to ensure a high level resolution in space. Similar to the previous developments [4–6], coefficients of the discrete boundary integral equations are evaluated with the sufficient precision using semi-analytic approach to ensure exceptional accuracy of the boundary element formulation. To demonstrate the attractiveness of the poly-region BEM formulation, we consider a numerical example of the well-known Rayleigh-Benard problem governed by the Boussinesq equations.

scholarly journals Quotient-space boundary element methods for scattering at complex screens

2021 ◽  
Author(s):  
Xavier Claeys ◽  
Lorenzo Giacomel ◽  
Ralf Hiptmair ◽  
Carolina Urzúa-Torres
Keyword(s):  
Boundary Element ◽  
Integral Equations ◽  
Quotient Space ◽  
Boundary Integral Equations ◽  
Boundary Integral ◽  
Boundary Element Methods ◽  
Link Type ◽  
Trace Spaces ◽  
Space Boundary ◽  
Single Trace

AbstractA complex screen is an arrangement of panels that may not be even locally orientable because of junction lines. A comprehensive trace space framework for first-kind variational boundary integral equations on complex screens has been established in Claeys and Hiptmair (Integr Equ Oper Theory 77:167–197, 2013. https://doi.org/10.1007/s00020-013-2085-x) for the Helmholtz equation, and in Claeys and Hiptmair (Integr Equ Oper Theory 84:33–68, 2016. https://doi.org/10.1007/s00020-015-2242-5) for Maxwell’s equations in frequency domain. The gist is a quotient space perspective that allows to make sense of jumps of traces as factor spaces of multi-trace spaces modulo single-trace spaces without relying on orientation. This paves the way for formulating first-kind boundary integral equations in weak form posed on energy trace spaces. In this article we extend that idea to the Galerkin boundary element (BE) discretization of first-kind boundary integral equations. Instead of trying to approximate jumps directly, the new quotient space boundary element method employs a Galerkin BE approach in multi-trace boundary element spaces. This spawns discrete boundary integral equations with large null spaces comprised of single-trace functions. Yet, since the right-hand-sides of the linear systems of equations are consistent, Krylov subspace iterative solvers like GMRES are not affected by the presence of a kernel and still converge to a solution. This is strikingly confirmed by numerical tests.

scholarly journals Stress Analysis of Three-Dimensional Media Containing Localized Zone by FEM-SGBEM Coupling

2011 ◽  
Vol 2011 ◽  
pp. 1-27
Author(s):  
Jaroon Rungamornrat ◽  
Sakravee Sripirom
Keyword(s):  
Stress Analysis ◽  
Boundary Element ◽  
Integral Equations ◽  
Boundary Integral Equations ◽  
Numerical Technique ◽  
Computational Cost ◽  
Three Dimensional ◽  
Boundary Integral ◽  
Boundary Element Methods ◽  
Element Method

This paper presents an efficient numerical technique for stress analysis of three-dimensional infinite media containing cracks and localized complex regions. To enhance the computational efficiency of the boundary element methods generally found inefficient to treat nonlinearities and non-homogeneous data present within a domain and the finite element method (FEM) potentially demanding substantial computational cost in the modeling of an unbounded medium containing cracks, a coupling procedure exploiting positive features of both the FEM and a symmetric Galerkin boundary element method (SGBEM) is proposed. The former is utilized to model a finite, small part of the domain containing a complex region whereas the latter is employed to treat the remaining unbounded part possibly containing cracks. Use of boundary integral equations to form the key governing equation for the unbounded region offers essential benefits including the reduction of the spatial dimension and the corresponding discretization effort without the domain truncation. In addition, all involved boundary integral equations contain only weakly singular kernels thus allowing continuous interpolation functions to be utilized in the approximation and also easing the numerical integration. Nonlinearities and other complex behaviors within the localized regions are efficiently modeled by utilizing vast features of the FEM. A selected set of results is then reported to demonstrate the accuracy and capability of the technique.

Numerical Implementation

1994 ◽  
Author(s):  
Nhan Phan-Thien ◽  
Sangtae Kim
Keyword(s):  
Boundary Element ◽  
Numerical Integration ◽  
Integral Equations ◽  
Large Scale ◽  
Boundary Integral Equations ◽  
Numerical Technique ◽  
Boundary Integral ◽  
Practical Implementation ◽  
Computing Environment ◽  
The Individual

Analytical solutions to a set of boundary integral equations are rare, even with simple geometries and boundary conditions. To make any reasonable progress, a numerical technique must be used. There are basically four issues that must be discussed in any numerical scheme dealing with integral equations. The first and most basic one is how numerical integration can be effected, together with an effective way of dealing with singular kernels of the type encountered in elastostatics. Numerical integration is usually termed numerical quadrature, meaning mathematical formulae for numerical integration. The second issue is the boundary discretization: when integration over the whole boundary is replaced by a sum of the integrations over the individual patches on the boundary. Each patch would be a finite element, or in our case, a boundary element on the surface. Obviously a high-order integration scheme can be devised for the whole domain, thus eliminating the need for boundary discretization. Such a scheme would be problem dependent and therefore would not be very useful to us. The third issue has to do with the fact that we are constrained by the very nature of the numerical approximation process to search for solutions within a certain subspace of L2, say the space of piecewise constant functions in which the unknowns are considered to be constant over a boundary element. It is the order of this subspace, together with the order and the nature of the interpolation of the geometry, that gives rise to the names of various boundary element schemes. Finally, one is faced with the task of solving a set of linear algebraic equations, which is usually dense (the system matrix is fully populated) and potentially ill-conditioned. A direct solver such as Gauss elimination may be very efficient for small- to medium-sized problems but will become stuck in a large-scale simulation, where the only feasible solution strategy is an iterative method. In fact, iterative solution strategies lead naturally to a parallel algorithm under a suitable parallel computing environment. This chapter will review various issues involved in the practical implementation of the CDL-BIEM on a serial computer and on a distributed computing environment.

Sign in / Sign up

Export Citation Format

Share Document