The boundary element method for the solution of Stokes equations in two-dimensional domains

A boundary element method for the solution of Stokes equations governing creeping flow or Stokes flow in the interior of an arbitrary two-dimensional domain is presented. A procedure for introducing pressure data on the boundary of the domain is also included and the integral coefficients of the resulting linear algebraic equations are evaluated analytically. Calculations are performed in a circular domain using a variety of different boundary conditions, including a combination of the fluid velocity and the pressure. Results are presented both on the boundary and inside the solution domain in order to illustrate that the boundary element method developed here provides an efficient technique, in terms of accuracy and convergence, to investigate Stokes flow numerically.     

Modeling of 3D transversely piezoelectric and elastic bimaterials using the boundary element method

 This paper presents a numerical model for three-dimensional transversely isotropic bimaterials based on the boundary element formulation. The point force solutions expressed in a united-form for distinct eigenvalues are studied for transversely isotropic piezoelectricity and pure elasticity. A boundary integral formulation is implemented for the modeling of two-phase materials. In this study, the stress distributions are computed for a near interface flaw. The influences of the shape and location of the flaw on the the stress concentration are examined. The accuracy of the numerical procedures is validated through selected example problems and comparison studies. Received 3 October 2001 / Accepted 9 April 2002     

An iterative boundary element method for Cauchy inverse problems

 We are interested in this paper in recovering an harmonic function from the knowledge of Cauchy data on some part of the boundary. A new inversion method is introduced. It reduces the Cauchy problem resolution to the determination of the resolution of a sequence of well-posed problems. The sequence of these solutions is proved to converge to the Cauchy problem solution. The algorithm is implemented in the framework of boundary elements. Displayed numerical results highlight its accuracy, as well as its robustness to noisy data. Received 6 November 2000     

Numerical analysis of the inviscid incompressible flow in two-dimensional radial-flow pump impellers

The flow through a two-dimensional centrifugal impeller fitted with equiangular blades of arbitrary geometry is investigated using a combination of conformal mapping with a boundary element technique. The blades can be thin or thick of arbitrary cross-section. A theoretical analysis and a numerical procedure are developed to determine the pressure distributions along the blade.     

A meshfree formulation of local radial point interpolation method (LRPIM) for incompressible flow simulation

 The LRPIM method is adopted to simulate the two-dimensional natural convection problems within enclosed domain of different geometries. In this paper, the vorticity-stream function form of N-S equations is taken as the governing equations. It was observed that the obtained results agreed very well with others available in the literatures, and with the same nodal density, the accuracy achieved by the LRPIM method is much higher than that of the finite difference (FD) method. The numerical examples show that the present LRPIM method can successfully deal with incompressible flow problems on randomly distributed nodes. Received: 2 April 2002 / Accepted: 6 January 2003 The authors would like to thank Mr. Y. T. Gu for his contribution to this work.     

The least-squares meshfree method for solving linear elastic problems

 A meshfree method based on the first-order least-squares formulation for linear elasticity is presented. In the authors' previous work, the least-squares meshfree method has been shown to be highly robust to integration errors with the numerical examples of Poisson equation. In the present work, conventional formulation and compatibility-imposed formulation for linear elastic problems are studied on the convergence behavior of the solution and the robustness to the inaccurate integration using simply constructed background cells. In the least-squares formulation, both primal and dual variables can be approximated by the same function space. This can lead to higher rate of convergence for dual variables than Galerkin formulation. In general, the incompressible locking can be alleviated using mixed formulations. However, in meshfree framework these approaches involve an additional use of background grids to implement lower approximation space for dual variables. This difficulty is avoided in the present method and numerical examples show the uniform convergence performance in the incompressible limit. Therefore the present method has little burden of the requirement of background cells for the purposes of integration and alleviating the incompressible locking. Received: 16 December 2001 / Accepted: 4 November 2002     

The traction boundary contour method for linear elasticity

This paper presents a further development of the boundary contour method. The boundary contour method is extended to cover the traction boundary integral equation. A traction boundary contour method is proposed for linear elastostatics. The formulation of traction boundary contour method is regular for points except the ends of the boundary element and corners. The present approach only requires line integrals for three‐dimensional problems and function evaluations at the ends of boundary elements for two‐dimensional cases. The implementation of the traction boundary contour method with quadratic boundary elements is presented for two‐dimensional problems. Numerical results are given for some two‐dimensional examples, and these are compared with analytical solutions. This method is shown to give excellent results for illustrative examples. Copyright © 1999 John Wiley & Sons, Ltd.     

A complex variable boundary collocation method for plane elastic problems

 Lagrange interpolation is extended to the complex plane in this paper. It turns out to be composed of two parts: polynomial and rational interpolations of an analytical function. Based on Lagrange interpolation in the complex plane, a complex variable boundary collocation approach is constructed. The method is truly meshless and singularity free. It features high accuracy obtained by use of a small number of nodes as well as dimensionality advantage, that is, a two-dimensional problem is reduced to a one-dimensional one. The method is applied to two-dimensional problems in the theory of plane elasticity. Numerical examples are in very good agreement with analytical ones. The method is easy to be implemented and capable to be able to give the stress states at any point within the solution domain. Received: 20 August 2002 / Accepted: 31 January 2003     

A Trefftz function approximation in local boundary integral equations

 In the present paper the Trefftz function as a test function is used to derive the local boundary integral equations (LBIE) for linear elasticity. Since Trefftz functions are regular, much less requirements are put on numerical integration than in the conventional boundary integral method. The moving least square (MLS) approximation is applied to the displacement field. Then, the traction vectors on the local boundaries are obtained from the gradients of the approximated displacements by using Hooke's law. Nodal points are randomly spread on the domain of the analysed body. The present method is a truly meshless method, as it does not need a finite element mesh, either for purposes of interpolation of the solution variables, or for the integration of the energy. Two ways are presented to formulate the solution of boundary value problems. In the first one the local boundary integral equations are written in all nodes (interior and boundary nodes). In the second way the LBIE are written only at the interior nodes and at the nodes on the global boundary the prescribed values of displacements and/or tractions are identified with their MLS approximations. Numerical examples for a square patch test and a cantilever beam are presented to illustrate the implementation and performance of the present method. Received 6 November 2000     

An advanced meshless LBIE/RBF method for solving two-dimensional incompressible fluid flows

The present work presents a meshless local boundary integral equation (LBIE) method for the solution of two-dimensional incompressible fluid flow problems governed by the Navier–Stokes equations. The method uses, for its meshless implementation, nodal points spread over the analyzed domain and employs in an efficient way the radial basis functions (RBF) for the interpolation of the interior and boundary variables. The inverse matrix of the RBF is computed only once for every nodal point and the interpolation functions are evaluated by the inner product of the inverse matrix with the weight vector associated to the integration point. This technique leads to a fast and efficient meshless approach, the locality of the method is maintained and the system matrices are banded with small bandwidth. The velocity–vorticity approach of the Navier–Stokes equations is adopted and the LBIEs are derived for the velocity and the vorticity field, resulting in a very stable and accurate implementation. The evaluation of the volume integrals is accomplished via a very efficient and accurate technique by triangularizing the local area of the nodal point to the minimum number of well formed triangles. Numerical examples illustrate the proposed methodology and demonstrate its accuracy.     

Evolutionary optimization in thermoelastic problems using the boundary element method

 The paper is devoted to application of evolutionary algorithms and the boundary element method to shape optimization of structures for various thermomechanical criteria, inverse problems of finding an optimal distribution of temperature on the boundary and identification of unknown boundary. Design variables are specified by Bezier curves. Several numerical examples of evolutionary computation are presented. Received 6 November 2000     

A study on time schemes for DRBEM analysis of elastic impact wave

 The precise integration and differential quadrature methods are two new unconditionally stable numerical schemes to approximate time derivative with more than the second order accuracy. Recent studies showed that compared with the Houbolt and Newmark methods, they produced more accurate solutions with large time step for the problems where response is primarily dominated by low and intermediate frequency modes. This paper aims to investigate these time schemes in the context of the dual reciprocity BEM (DRBEM) formulation of various shock-excited scalar elastic wave problems, where high modes have important affect on traction response. The Houbolt method was widely recommended in many literatures for such DRBEM dynamic formulations. However, this study found that the damped Newmark algorithm was the most efficient and accurate for impact traction analysis in conjunction with the DRBEM. The precise integration and differential quadrature methods are shown inapplicable for such shock-excited problems due to the absence of numerical damping. On the other hand, we also found that to achieve the same order of accuracy, the differential quadrature method required much less computing effort than the precise integration method due to the use of the Bartels–Stewart algorithm solving the resulting Lyapunov matrix analogization equation. Received 6 November 2000     

Meshless analysis of two-dimensional Stokes flows with the Galerkin boundary node method

In this paper, the Galerkin boundary node method (GBNM) is developed for the solution of stationary Stokes problems in two dimensions. The GBNM is a boundary only meshless method that combines a variational form of boundary integral formulations for governing equations with the moving least-squares (MLS) approximations for construction of the trial and test functions. Boundary conditions in this approach are included into the variational form, thus they can be applied directly and easily despite the MLS shape functions lack the property of a delta function. Besides, the GBNM keeps the symmetry and positive definiteness of the variational problems. Convergence analysis results of both the velocity and the pressure are given. Some selected numerical tests are also presented to demonstrate the efficiency of the method.     

Domain decomposition in the symmetric boundary element analysis

 Recent developments in the symmetric boundary element method (SBEM) have shown a clear superiority of this formulation over the collocation method. Its competitiveness has been tested in comparison to the finite element method (FEM) and is manifested in several engineering problems in which internal boundaries are present, i.e. those in which the body shows a jump in the physical characteristics of the material and in which an appropriate study of the response must be used. When we work in the ambit of the SBE formulation, the body is subdivided into macroelements characterized by some relations which link the interface boundary unknowns to the external actions. These relations, valid for each macroelement and characterized by symmetric matricial operators, are similar in type to those obtainable for the FEM. The assembly of the macroelements based on the equilibrium conditions, or on the compatibility conditions, or on both of these conditions leads to three analysis methods: displacement, force, and mixed-value methods, respectively. The use of the fundamental solutions involves the punctual satisfaction of the compatibility and of the equilibrium inside each macroelement and it causes a stringent elastic response close to the actual solution. Some examples make it possible to perform numerical checks in comparison with solutions obtainable in closed form. These checks show that the numerical solutions are floating ones when the macroelement geometry obtained by subdividing the body changes. Received 26 January 2001     

Shape design sensitivity analysis and optimization of two-dimensional periodic thermal problems using BEM

 A general procedure to perform shape design sensitivity analysis for two-dimensional periodic thermal diffusion problems is developed using boundary integral equation formulation. The material derivative concept to describe shape variation is used. The temperature is decomposed into a steady state component and a perturbation component. The adjoint variable method is used by utilizing integral identities for each component. The primal and adjoint systems are solved by boundary element method. The sensitivity results compared with those by finite difference show good accuracy. The shape optimal design problem of a plunger model for the panel of a television bulb, which operates periodically, is solved as an example. Different objectives and amounts of heat flux allowed are studied. Corresponding optimum shapes of the cooling boundary of the plunger are obtained and discussed. Received 15 August 2001 / Accepted 28 February 2002     

Interpolation functions in control volume finite element method

 The main contribution of this paper is the study of interpolation functions in control volume finite element method used in equal order and applied to an incompressible two-dimensional fluid flow. Especially, the exponential interpolation function expressed in the elemental local coordinate system is compared to the classic linear interpolation function expressed in the global coordinate system. A quantitative comparison is achieved by the application of these two schemes to four flows that we know the analytical solutions. These flows are classified in two groups: flows with privileged direction and flows without. The two interpolation functions are applied to a triangular element of the domain then; a direct comparison of the results given by each interpolation function to the exact value is easily realized. The two functions are also compared when used to solve the discretized equations over the entire domain. Stability of the numerical process and accuracy of solutions are compared. Received: 20 October 2002 / Accepted: 2 December 2002     

On probabilistic viscous incompressible flow of some composite fluids

 The perturbation-based stochastic finite element formulation for the viscous incompressible fluid flow with heat transfer is proposed below. Analyzed viscoelastic fluid contains elastic spherical particles with randomly treated radii and total number in the fluid volume; these random variables are defined using their expected values and variances. Starting from these parameters and statistically described fluid viscosity, the probabilistic moments of the effective viscosity for a fluid with these suspensions are derived thanks to the second order perturbation second central probabilistic moment approach. So defined random effective fluid is next studied in the incompressible isobaric Couette flow between parallel plates, where the heat transfer effects are included. Defining boundary velocities and temperatures by their first two moments, the expected values and cross-covariances of the relevant homogenized fluid state functions are calculated. The engineering application of the approach can be modeling of polymers in fluid state (during their processing), reinforced with rubber particles to strengthen the entire composite as well as in general computational simulation of the composite fluids. The stochastic methodology applied to fluids with random solid suspensions can be extended on random flows with stochastic bubbles as well as multiphase coupled flow problems, where partially saturated media are considered [7], too. The main value of the stochastic approach is that it creates the opportunity to determine reliability indices, quite analogously to the research in the domain of solids and structures made of composites. Received 20 October 2001 / Accepted 15 February 2002     

Boundary-integral equation analysis of twisted internally cracked axisymmetric bimaterial elastic solids

 Green's function is obtained for the infinite bimaterial elastic solid, containing an internal circular interface crack, loaded by a unit tangential co-axial circular source. An axisymmetric direct boundary integral equation (BIE) is used for the analysis of a finite bimaterial axisymmetric body containing an internal circular interface crack and a finite homogeneous cracked cylinder, both under torsional loading. Using the proposed technique, no discretization of the crack surface is necessary. Numerical results for both examples as obtained by the proposed method are presented and discussed. Received: 29 October 2001 / Accepted: 29 May 2002     

Identification of material properties and cavities in two-dimensional linear elasticity

 In this paper, the simultaneous identification of the Poisson ratio, the shear modulus and a circular cavity embedded in an isotropic linear elastic material from boundary measurements is investigated. The numerical method proposed is based on the least-squares minimisation of the errors between the measured and the calculated tractions on the outer boundary using the boundary element method (BEM). Received: 26 July 2002 / Accepted: 4 February 2003 L. Marin would like to acknowledge the financial support received from the EPSRC. The very constructive comments made by the referees on the first version of this paper are gratefully acknowledged.