scholarly journals Simulation of 2D Elastic Bodies with Randomly Distributed Circular Inclusions Using the BEM

2007 ◽  
Vol 1 (2) ◽  
Author(s):  
Yao Zhenhan ◽  
Kong Fanzhong ◽  
Zheng Xiaoping

Based on the Rizzo’s direct boundary integral equation formulation for elasticity problems, elastic bodies with randomly distributed circular inclusions are simulated using the boundary element method. The given numerical examples show that the boundary element method is more accurate and more efficient than the finite element method for such type of problems. The presented approach can be successfully applied to estimate the equivalent elastic properties of many composite materials.

1987 ◽  
Vol 109 (1) ◽  
pp. 15-21 ◽  
Author(s):  
A. F. Seybert ◽  
C. Y. R. Cheng

This paper is concerned with the application of the Boundary Element Method (BEM) to interior acoustics problems governed by the reduced wave (Helmholtz) differential equation. The development of an integral equation valid at the boundary of the interior region follows a similar formulation for exterior problems, except for interior problems the Sommerfeld radiation condition is not invoked. The boundary integral equation for interior problems does not suffer from the nonuniqueness difficulty associated with the boundary integral equation formulation for exterior problems. The boundary integral equation, once obtained, is solved for a specific geometry using quadratic isoparametric surface elements. A simplification for axisymmetric cavities and boundary conditions permits the solution to be obtained using line elements on the generator of the cavity. The present formulation includes the case where a node may be placed at a position on the boundary where there is not a unique tangent plane (e.g., at an edge or a corner point). The BEM capability is demonstrated for two types of classical interior axisymmetric problems: the acoustic response of a cavity and the transmission loss of a muffler. For the cavity response comparison data are provided by an analytical solution. For the muffler problem the BEM solution is compared to data obtained by a finite element method analysis.


1983 ◽  
Vol 18 (4) ◽  
pp. 199-205 ◽  
Author(s):  
R T Fenner

The principles of the boundary integral equation (BIE) or boundary element method (BEM) are discussed in a non-mathematical way. The technique is compared with other numerical methods, particularly the finite element method (FEM), on the basis of computational efficiency, and the main advantages and disadvantages of the BIE approach are outlined.


2017 ◽  
Vol 2017 ◽  
pp. 1-12
Author(s):  
Jui-Hsiang Kao

This research develops an Advance-Tracing Boundary Element Method in the time domain to calculate the waves that radiate from an immersed obstacle moving with random acceleration. The moving velocity of the immersed obstacle is multifrequency and is projected along the normal direction of every element on the obstacle. The projected normal velocity of every element is presented by the Fourier series and includes the advance-tracing time, which is equal to a quarter period of the moving velocity. The moving velocity is treated as a known boundary condition. The computing scheme is based on the boundary integral equation in the time domain, and the approach process is carried forward in a loop from the first time step to the last. At each time step, the radiated pressure on each element is updated until obtaining a convergent result. The Advance-Tracing Boundary Element Method is suitable for calculating the radiating problem from an arbitrary obstacle moving with random acceleration in the time domain and can be widely applied to the shape design of an immersed obstacle in order to attain security and confidentiality.


Author(s):  
Yijun Liu ◽  
Milind Bapat

Some recent development of the fast multipole boundary element method (BEM) for modeling acoustic wave problems in both 2-D and 3-D domains are presented in this paper. First, the fast multipole BEM formulation for 2-D acoustic wave problems based on a dual boundary integral equation (BIE) formulation is presented. Second, some improvements on the adaptive fast multipole BEM for 3-D acoustic wave problems based on the earlier work are introduced. The improvements include adaptive tree structures, error estimates for determining the numbers of expansion terms, refined interaction lists, and others in the fast multipole BEM. Examples involving 2-D and 3-D radiation and scattering problems solved by the developed 2-D and 3-D fast multipole BEM codes, respectively, will be presented. The accuracy and efficiency of the fast multipole BEM results clearly demonstrate the potentials of the fast multipole BEM for solving large-scale acoustic wave problems that are of practical significance.


1996 ◽  
Vol 63 (2) ◽  
pp. 278-286 ◽  
Author(s):  
A. Nagarajan ◽  
S. Mukherjee ◽  
E. Lutz

This paper presents a novel variant of the boundary element method, here called the boundary contour method, applied to three-dimensional problems of linear elasticity. In this work, the surface integrals on boundary elements of the usual boundary element method are transformed, through an application of Stokes’ theorem, into line integrals on the bounding contours of these elements. Thus, in this formulation, only line integrals have to be numerically evaluated for three-dimensional elasticity problems—even for curved surface elements of arbitrary shape. Numerical results are presented for some three-dimensional problems, and these are compared against analytical solutions.


2017 ◽  
Vol 140 (1) ◽  
Author(s):  
Sofia Sarraf ◽  
Ezequiel López ◽  
Laura Battaglia ◽  
Gustavo Ríos Rodríguez ◽  
Jorge D'Elía

In the boundary element method (BEM), the Galerkin weighting technique allows to obtain numerical solutions of a boundary integral equation (BIE), giving the Galerkin boundary element method (GBEM). In three-dimensional (3D) spatial domains, the nested double surface integration of GBEM leads to a significantly larger computational time for assembling the linear system than with the standard collocation method. In practice, the computational time is roughly an order of magnitude larger, thus limiting the use of GBEM in 3D engineering problems. The standard approach for reducing the computational time of the linear system assembling is to skip integrations whenever possible. In this work, a modified assembling algorithm for the element matrices in GBEM is proposed for solving integral kernels that depend on the exterior unit normal. This algorithm is based on kernels symmetries at the element level and not on the flow nor in the mesh. It is applied to a BIE that models external creeping flows around 3D closed bodies using second-order kernels, and it is implemented using OpenMP. For these BIEs, the modified algorithm is on average 32% faster than the original one.


2016 ◽  
Vol 19 (4) ◽  
pp. 970-997 ◽  
Author(s):  
Brian Zinser ◽  
Wei Cai

AbstractIn this paper, we will present a high-order, well-conditioned boundary element method (BEM) based on Müller's hypersingular second kind integral equation formulation to accurately compute electrostatic potentials in the presence of inhomogeneity embedded within layered media. We consider two types of inhomogeneities: the first one is a simple model of an ion channel which consists of a finite height cylindrical cavity embedded in a layered electrolytes/membrane environment, and the second one is a Janus particle made of two different semi-spherical dielectric materials. Both types of inhomogeneities have relevant applications in biology and colloidal material, respectively. The proposed BEM gives condition numbers, allowing fast convergence of iterative solvers compared to previous work using first kind of integral equations. We also show that the second order basis converges faster and is more accurate than the first order basis for the BEM.


Sign in / Sign up

Export Citation Format

Share Document