3D multidomain BEM for solving the Laplace equation

An efficient 3D multidomain boundary element method (BEM) for solving problems governed by the Laplace equation is presented. Integral boundary equations are discretized using mixed boundary elements. The field function is interpolated using a continuous linear function while its derivative in a normal direction is interpolated using a discontinuous constant function over surface boundary elements. Using a multidomain approach, also known as the subdomain technique, sparse system matrices similar to the finite element method (FEM) are obtained. Interface boundary conditions between subdomains leads to an over-determined system matrix, which is solved using a fast iterative linear least square solver. The accuracy and robustness of the developed numerical algorithm is presented on a scalar diffusion problem using simple cube geometry and various types of meshes. Efficiency is demonstrated with potential flow around the complex geometry of a fighter airplane using tetrahedral mesh with over 100,000 subdomains on a personal computer.     

A multidomain boundary element method for two equation turbulence models

The paper deals with the multidomain Boundary Element Method (BEM) for modelling 2D complex turbulent flow using low Reynolds two equation turbulence models. While the BEM is widely accepted for laminar flow this is the first case, where this method is applied for a complex flow problems using k–ε turbulence model. The integral boundary domain equations are discretised using mixed boundary elements and a multidomain method also known as subdomain technique. The resulting system matrix is overdetermined, sparse, block banded and solved using fast iterative linear least squares solver. The simulation of turbulent flow over a backward step is in excellent agreement with the finite volume method using the same turbulent model.     

Analytical integration of the moments in the diagonal form fast multipole boundary element method for 3-D acoustic wave problems

A diagonal form fast multipole boundary element method (BEM) is presented in this paper for solving 3-D acoustic wave problems based on the Burton-Miller boundary integral equation (BIE) formulation. Analytical expressions of the moments in the diagonal fast multipole BEM are derived for constant elements, which are shown to be more accurate, stable and efficient than those using direct numerical integration. Numerical examples show that using the analytical moments can reduce the CPU time by a lot as compared with that using the direct numerical integration. The percentage of CPU time reduction largely depends on the proportion of the time used for moments calculation to the overall solution time. Several examples are studied to investigate the effectiveness and efficiency of the developed diagonal fast multipole BEM as compared with earlier p³ fast multipole method BEM, including a scattering problem of a dolphin modeled with 404,422 boundary elements and a radiation problem of a train wheel track modeled with 257,972 elements. These realistic, large-scale BEM models clearly demonstrate the effectiveness, efficiency and potential of the developed diagonal form fast multipole BEM for solving large-scale acoustic wave problems.     

3D web-splines solution to human eye heat distribution using bioheat equation

A bioheat transfer model of the human eye is constructed using weighted extended b-splines as shape functions for the finite element method. This newly developed computational approach is employed to calculate the steady-state temperature distribution in a normal human eye. Firstly, the human eye is evaluated in two dimensions. The simulation results which are verified using the values reported in the literature, point out to better efficiency in terms of the accuracy level. Consequently, to give a more precise representation of the actual human eye, three dimensional modeling is simulated using these new finite elements in conjunction with linear, quadratic and cubic b-splines. Grid convergence number estimates are derived for both sets of simulations. It is shown in this paper that 3D cubic web-spline model is 0.03 °C closer to the mean value of the previous studies including the experimental studies as compared with 0.17 °C improvement obtained with the standard finite element method. Our findings indicate that weighted extended b-spline solutions improve the computational methods for health care.     

Comparison of boundary and finite element methods for moving-boundary problems governed by a potential

Three formulations of the boundary element method (BEM) and one of the Galerkin finite element method (FEM) are compared according to accuracy and efficiency for the spatial discretization of two-dimensional, moving-boundary problems based on Laplace's equation. The same Euler-predictor, trapezoid-corrector scheme for time integration is used for all four methods. The model problems are on either a bounded or a semi-infinite strip and are formulated so that closed-form solutions are known. Infinite elements are used with both the BEM and FEM techniques for the unbounded domain. For problems with the bounded region, the BEM using the free-space Green's function and piecewise quadratic interpolating functions (QBEM) is more accurate and efficient than the BEM with linear interpolation. However, the FEM with biquadratic basis functions is more efficient for a given accuracy requirement than the QBEM, except when very high accuracy is demanded. For the unbounded domain, the preferred method is the BEM based on a Green's function that satisfies the lateral symmetry conditions and which leads to discretization of the potential only along the moving surface. This last formulation is the only one that reliably satisfies the far-field boundary condition.     

A fast multipole boundary element method for solving two-dimensional thermoelasticity problems

A fast multipole boundary element method (BEM) for solving general uncoupled steady-state thermoelasticity problems in two dimensions is presented in this paper. The fast multipole BEM is developed to handle the thermal term in the thermoelasticity boundary integral equation involving temperature and heat flux distributions on the boundary of the problem domain. Fast multipole expansions, local expansions and related translations for the thermal term are derived using complex variables. Several numerical examples are presented to show the accuracy and effectiveness of the developed fast multipole BEM in calculating the displacement and stress fields for 2-D elastic bodies under various thermal loads, including thin structure domains that are difficult to mesh using the finite element method (FEM). The BEM results using constant elements are found to be accurate compared with the analytical solutions, and the accuracy of the BEM results is found to be comparable to that of the FEM with linear elements. In addition, the BEM offers the ease of use in generating the mesh for a thin structure domain or a domain with complicated geometry, such as a perforated plate with randomly distributed holes for which the FEM fails to provide an adequate mesh. These results clearly demonstrate the potential of the developed fast multipole BEM for solving 2-D thermoelasticity problems.     

Development of circular arc boundary elements method

The boundary element method (BEM) is a popular technique to solve engineering problems. We compare the circular arc elements (CAE) discretization to both linear and quadratic discretizations. The main aim of this paper is to determine analytical expressions for the discretization error in 2D BEM for the Laplace equation using CAE discretization. The results are validated by numerical examples.     

A coupled boundary element-finite difference solution of the elliptic modified mild slope equation

The modified mild slope equation of [5] is solved using a combination of the boundary element method (BEM) and the finite difference method (FDM). The exterior domain of constant depth and infinite horizontal extent is solved by a BEM using linear or quadratic elements. The interior domain with variable depth is solved by a flexible order of accuracy FDM in boundary-fitted curvilinear coordinates. The two solutions are matched along the common boundary of two methods (the BEM boundary) to ensure continuity of value and normal flux. Convergence of the individual methods is shown and the combined solution is tested against several test cases. Results for refraction and diffraction of waves from submerged bottom mounted obstacles compare well with experimental measurements and other computed results from the literature.     

Combined single domain and subdomain BEM for 3D laminar viscous flow

A subdomain boundary element method (BEM) using a continuous quadratic interpolation of function and discontinuous linear interpolation of flux is presented for the solution of the vorticity transport equation and the kinematics equation in 3D. By employing compatibility conditions between subdomains an over-determined system of linear equations is obtained, which is solved in a least squares manner. The method, combined with the single domain BEM, is used to solve laminar viscous flows using the velocity vorticity formulation of Navier–Stokes equations. The versatility and accuracy of the method are proven using the 3D lid driven cavity test case.     

Integral and geometrical means in the analytical evaluation of the BEM integrals for a 3D Laplace equation

By using integral theorems and geometrical interpretations, the analytical formulas for the coefficients occurring in boundary element method (BEM) equations for a 3D Laplace equation were found for arbitrary planar polygonal boundary elements with constant approximation. Closed forms for the gradients of the coefficients were also obtained. In addition, an analytical formula for planar triangular boundary elements with a linear approximation of potential was given. The formulas obtained are appropriate especially in singular and nearly singular cases.     

Dynamic analysis of interfacial crack problems in anisotropic bi-materials by a time-domain BEM

A time-domain boundary element method (BEM) together with the sub-domain technique is applied to study dynamic interfacial crack problems in two-dimensional (2D), piecewise homogeneous, anisotropic and linear elastic bi-materials. The bi-material system is divided into two homogeneous sub-domains along the interface and the traditional displacement boundary integral equations (BIEs) are applied on the boundary of each sub-domain. The present time-domain BEM uses a quadrature formula for the temporal discretization to approximate the convolution integrals and a collocation method for the spatial discretization. Quadratic quarter-point elements are implemented at the tips of the interface cracks. A displacement extrapolation technique is used to determine the complex dynamic stress intensity factors (SIFs). Numerical examples for computing the complex dynamic SIFs are presented and discussed to demonstrate the accuracy and the efficiency of the present time-domain BEM.     

Boundary element analysis for an interface crack between dissimilar elastoplastic materials

The boundary element method (BEM) is presented for elastoplastic analysis of cracks between two dissimilar materials. The boundary integral equations and integral representation of stress rates are written in such a form that all integrals can be evaluated by the regular Gaussian quadrature rule. An advanced multidomain BEM formulation is suggested for the solution of analysed problems where the substantial reduction of stiffness matrix is observed. The elastoplastic behaviour is modelled through the use of an approximation for the plastic component of the stresses. The boundary and the yielding zone are discretized by elements with quadratic approximations. In numerical examples the path independence of the J- and L-integrals for a straight interface crack and a circular arc-shaped interface crack are investigated, respectively. The influence of the different values of Young's modulus on the J-integral, shape and size of plastic zones is treated too.     

Adaptive fast multipole boundary element method for three-dimensional half-space acoustic wave problems

A new adaptive fast multipole boundary element method (BEM) for solving 3-D half-space acoustic wave problems is presented in this paper. The half-space Green's function is employed explicitly in the boundary integral equation (BIE) formulation so that a tree structure of the boundary elements only for the boundaries of the real domain need to be applied, instead of using a tree structure that contains both the real domain and its mirror image. This procedure simplifies the implementation of the adaptive fast multipole BEM and reduces the CPU time and memory storage by about a half for large-scale half-space problems. An improved adaptive fast multipole BEM is presented for the half-space acoustic wave problems, based on the one developed recently for the full-space problems. This new fast multipole BEM is validated using several simple half-space models first, and then applied to model 3-D sound barriers and a large-scale windmill model with five turbines. The largest BEM model with 557470 elements was solved in about an hour on a desktop PC. The accuracy and efficiency of the BEM results clearly show the potential of the adaptive fast multipole BEM for solving large-scale half-space acoustic wave problems that are of practical significance.     

An adaptive fast multipole boundary element method for three-dimensional acoustic wave problems based on the Burton–Miller formulation

The high solution costs and non-uniqueness difficulties in the boundary element method (BEM) based on the conventional boundary integral equation (CBIE) formulation are two main weaknesses in the BEM for solving exterior acoustic wave problems. To tackle these two weaknesses, an adaptive fast multipole boundary element method (FMBEM) based on the Burton–Miller formulation for 3-D acoustics is presented in this paper. In this adaptive FMBEM, the Burton–Miller formulation using a linear combination of the CBIE and hypersingular BIE (HBIE) is applied to overcome the non-uniqueness difficulties. The iterative solver generalized minimal residual (GMRES) and fast multipole method (FMM) are adopted to improve the overall computational efficiency. This adaptive FMBEM for acoustics is an extension of the adaptive FMBEM for 3-D potential problems developed by the authors recently. Several examples on large-scale acoustic radiation and scattering problems are presented in this paper which show that the developed adaptive FMBEM can be several times faster than the non-adaptive FMBEM while maintaining the accuracies of the BEM.     

Numerical p‐version refinement studies for the regularized stress‐BEM

In the development of the boundary element method (BEM) and the finite element method (FEM) researchers have typically selected similar basis functions. That is, both methods typically employ low‐order interpolations such as piece‐wise linear or piece‐wise quadratic and rely on h‐version refinement to increase accuracy as required. In the case of the FEM, the decision to use low‐order elements is made for computational efficiency as an attractive compromise between local modeling accuracy and sparseness of the resulting linear system. However, in many BEM formulations, low‐order elements may be the only practical choice given the complexity of using analytic integration formulae in conjunction with special integral interpretations. Unlike their efficient use in the FEM, fine meshes of low‐order elements in the BEM are highly inefficient from a computational standpoint given the dense nature of BEM systems. Moreover, owing to singularities in the kernel functions, the BEM should be expected to benefit more so than the FEM from very high levels of local accuracy. Through the use of regularized algorithms which only require numerical integration, p‐version refinement in the BEM is easily extended to include any set of basis functions with no significant increase in programming complexity. Numerical results show that by using interpolations as high as 12th and 16th order, one can expect reductions in error by as many as five orders of magnitude over comparable algorithms based on similar system size. For two‐dimensional problems, it is also shown that, for a given level of error, one can expect reductions in system size by an order of magnitude, thus leading to a reduction in computational expense for conventional algorithms by three orders of magnitude. Copyright © 2003 John Wiley & Sons, Ltd.     

Modeling of fatigue crack propagation using dual boundary element method and Gaussian Monte Carlo method

This paper studies the modeling of fatigue crack propagation on a multiple crack site of a finite plate using deterministic and probabilistic methods. Stress intensity factor has been calculated by the combined deterministic approach of the dual boundary element method (DBEM) and the probabilistic approach of the Gaussian Monte Carlo method. The Gaussian Monte Carlo method has been incorporated to simulate the random process of the fatigue crack propagation. A finite plate of aluminum alloy 2024-T3 with a thickness of 1.6 mm and 14 holes is analyzed and the fatigue life of the plate is predicted by following a linear elastic law of fracture mechanics. The results of fatigue life predicted by DBEM-Monte Carlo method are in good agreement with experimental ones. The same approach is also applied to two other engineering applications of a gear tooth and a bracket.     

A fast boundary cloud method for 3D exterior electrostatic analysis

An accelerated boundary cloud method (BCM) for boundary‐only analysis of 3D electrostatic problems is presented here. BCM uses scattered points unlike the classical boundary element method (BEM) which uses boundary elements to discretize the surface of the conductors. BCM combines the weighted least‐squares approach for the construction of approximation functions with a boundary integral formulation for the governing equations. A linear base interpolating polynomial that can vary from cloud to cloud is employed. The boundary integrals are computed by using a cell structure and different schemes have been used to evaluate the weakly singular and non‐singular integrals. A singular value decomposition (SVD) based acceleration technique is employed to solve the dense linear system of equations arising in BCM. The performance of BCM is compared with BEM for several 3D examples. Copyright © 2004 John Wiley & Sons, Ltd.     

On the effective implementation of a boundary element code on graphics processing units using an out-of-core LU algorithm

A collocation boundary element code for solving the three-dimensional Laplace equation, publicly available from http://intetec.org, has been adapted to run on an Nvidia Tesla general-purpose graphics processing unit (GPU). Global matrix assembly and LU factorization of the resulting dense matrix are performed on the GPU. Out-of-core techniques are used to solve problems larger than the available GPU memory. The code achieved about 10 times speedup in matrix assembly over a single CPU core and about 56 Gflops/s in the LU factorization using only 512 Mbytes of GPU memory. Details of the GPU implementation and comparisons with the standard sequential algorithm are included to illustrate the performance of the GPU code.     

Coupled FE–BE method for eigenvalue analysis of elastic structures submerged in an infinite fluid domain

For thin elastic structures submerged in heavy fluid, e.g., water, a strong interaction between the structural domain and the fluid domain occurs and significantly alters the eigenfrequencies. Therefore, the eigenanalysis of the fluid–structure interaction system is necessary. In this paper, a coupled finite element and boundary element (FE–BE) method is developed for the numerical eigenanalysis of the fluid–structure interaction problems. The structure is modeled by the finite element method. The compressibility of the fluid is taken into consideration, and hence the Helmholtz equation is employed as the governing equation and solved by the boundary element method (BEM). The resulting nonlinear eigenvalue problem is converted into a small linear one by applying a contour integral method. Adequate modifications are suggested to improve the efficiency of the contour integral method and avoid missing the eigenvalues of interest. The Burton–Miller formulation is applied to tackle the fictitious eigenfrequency problem of the BEM, and the optimal choice of its coupling parameter is investigated for the coupled FE–BE method. Numerical examples are given and discussed to demonstrate the effectiveness and accuracy of the developed FE–BE method. Copyright © 2016 John Wiley & Sons, Ltd.     

A BEM analysis of fracture mechanics in 2D anisotropic piezoelectric solids

This paper presents a single-domain boundary element method (BEM) analysis of fracture mechanics in 2D anisotropic piezoelectric solids. In this analysis, the extended displacement (elastic displacement and electrical potential) and extended traction (elastic traction and electrical displacement) integral equations are collocated on the outside boundary (no-crack boundary) of the problem and on one side of the crack surface, respectively. The Green's functions for the anisotropic piezoelectric solids in an infinite plane, a half plane, and two joined dissimilar half-planes are also derived using the complex variable function method. The extrapolation of the extended relative crack displacement is employed to calculate the extended `stress intensity factors' (SIFs), i.e., K_I, K_II, K_III and K_IV. For a finite crack in an infinite anisotropic piezoelectric solid, the extended SIFs obtained with the current numerical formulation were found to be very close to the exact solutions. For a central and inclined crack in a finite and anisotropic piezoelectric solid, we found that both the coupled and uncoupled (i.e., the piezoelectric coefficient e_ijk=0) cases predict very similar stress intensity factors K_I and K_II when a uniform tension σ_yy is applied, and very similar electric displacement intensity factor K_IV when a uniform electrical displacement D_y is applied. However, the relative crack displacement and electrical potential along the crack surface are quite different for the coupled and uncoupled cases. Furthermore, for a inclined crack within a finite domain, we found that while a uniform σ_yy (=1 N m⁻²) induces only a very small electrical displacement intensity factor (in the unit of Cm^−3/2), a uniform D_y (=1 C m⁻²) can produce very large stress intensity factors (in the unit of Nm^−3/2).