N. Simões^1,2, A. Tadeu²

CMES-Computer Modeling in Engineering & Sciences, Vol.9, No.3, pp. 221-234, 2005, DOI:10.3970/cmes.2005.009.221

Abstract The use of the Boundary Element Method (BEM) to formulate the 3D transient heat transfer through cylindrical structures with irregular cross-sections, bounded by a homogeneous elastic medium, is described in this paper. In this formulation, both the conduction and the convection phenomena are modeled. This system can be subjected to heat emitted by either point or line sources located somewhere in the media. The solution is first obtained in the frequency domain for a wide range of frequencies and axial wavenumbers. Time domain responses are later calculated by means of (fast) inverse Fourier transforms into More >

J. Purbolaksono¹, M. H. Aliabadi^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.8, No.1, pp. 73-90, 2005, DOI:10.3970/cmes.2005.008.073

Abstract This paper presents the dual boundary integral equations for the buckling analysis of the shear deformable cracked plates. The domain integrals which appear in this formulation are transferred to boundary integrals using the dual reciprocity method. The plate buckling displacement and hypersingular traction integral equations are presented as a standard eigenvalue problem, which would allow direct evaluation of the critical load factor and buckling modes. Several examples with different geometries and boundary conditions are presented to demonstrate the accuracy of the proposed formulation. More >

Y.C. Shiah¹, T.L. Guao¹, C.L. Tan²

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.3, pp. 321-338, 2005, DOI:10.3970/cmes.2005.007.321

Abstract It is well known in elastic stress analysis using the boundary element method (BEM) that an additional volume integral appears in the basic form of the boundary integral equation if thermal effects are considered. In order to restore this general numerical tool as a truly boundary solution technique, it is perhaps most desirable to transform this volume integral exactly into boundary ones. For general 2D anisotropic thermo-elastostatics without heat sources, this was only achieved very recently. The presence of concentrated heat sources in the domain, however, leads to singularities at these points that pose additional More >

Zai You Yan^1,2, Fang Sen Cui², Kin Chew Hung²

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.1, pp. 97-106, 2005, DOI:10.3970/cmes.2005.007.097

Abstract Taking the normal derivative of solid angles on the surface into account, a modified Burton and Miller's formulation is derived. From which, a more reasonable expression of the hypersingular operator is obtained. To overcome the hypersingular integral, the regularization scheme developed recently is employed. Plane acoustic wave scattering from a rigid sphere is computed to show the correctness of the modified formulation with the regularization scheme. In the computation, eight-nodded isoparametric element is applied. More >

Haitao Wang¹, Zhenhan Yao¹

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.1, pp. 85-96, 2005, DOI:10.3970/cmes.2005.007.085

Abstract This paper addresses a new boundary element method (BEM) for the numerical analysis of mechanical properties in 3D particle-reinforced composites. The BEM is accelerated by a new version fast multipole method (FMM) in order to perform large scale simulation of a representative volume element (RVE) containing up to several hundred randomly distributed elastic spherical particles on only one personal computer. The maximum number of degrees of freedom (DOF) reaches more than 300,000. Efficiency of the developed new version fast multipole BEM code is evaluated compared with other conventional solutions for BEM. The effects of micro-structural More >

B. Amaziane¹, A. Naji², D. Ouazar³, A. H.-D. Cheng⁴

CMC-Computers, Materials & Continua, Vol.2, No.2, pp. 85-96, 2005, DOI:10.3970/cmc.2005.002.085

Abstract In this paper the optimization of groundwater pumping in coastal aquifers under the threat of saltwater intrusion is investigated. The aquifer is inhomogeneous and contains several hydraulic conductivities zones. The aquifer data such as the hydraulic conductivities are uncertain, but with their expected mean and standard deviation values given. A stochastic boundary element method based on the perturbation technique is employed as the simulation tool. The stochastic optimization is handled by the chance-constrained programming. Genetic algorithm is selected as the optimization tool. Numerical examples of deterministic and stochastic problems are provided to demonstrate the feasibility More >

P. A. C. Porto¹, A. B. Jorge¹, G. O. Ribeiro²

CMES-Computer Modeling in Engineering & Sciences, Vol.10, No.1, pp. 65-78, 2005, DOI:10.3970/cmes.2005.010.065

Abstract This work deals with a numerical solution technique for the self-regular gradient form of Green's identity, the flux boundary integral equation (flux-BIE). The required C^1,α inter-element continuity conditions for the potential derivatives are imposed in the boundary element method (BEM) code through a non-symmetric variational formulation. In spite of using Lagrangian C⁰ elements, accurate numerical results were obtained when applied to heat transfer problems with singular or quasi-singular conditions, like boundary points and interior points which may be arbitrarily close to the boundary. The numerical examples proposed show that the developed algorithm based on the self-regular More >

Eduardo Divo¹, Alain J. Kassab²

CMC-Computers, Materials & Continua, Vol.2, No.4, pp. 277-288, 2005, DOI:10.3970/cmc.2005.002.277

Abstract A Dual Reciprocity Boundary Element Method is formulated to solve non-linear heat conduction problems. The approach is based on using the Kirchhoff transform along with lagging of the effective non-linear thermal diffusivity. A posteriori error estimate is used to provide effective estimates of the temporal and spatial error. A numerical example is used to demonstrate the approach. More >

Yoshihiro Ochiai¹, Vladimir Sladek²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.6, pp. 525-536, 2004, DOI:10.3970/cmes.2004.006.525

Abstract The conventional boundary element method (BEM) uses internal cells for the domain integralsCwhen solving nonlinear problems or problems with domain effects. This paper is concerned with conversion of the domain integral into boundary ones and some non-integral terms in a three-dimensional BIEM, which does not require the use of internal cells. This method uses arbitrary internal points instead of internal cells. The method is based on a three-dimensional interpolation method in this paper by using a polyharmonic function with volume distribution. In view of this interpolation method, the three-dimensional numerical integration is replaced by boundary More >

S. Callsen¹, O. von Estorff¹, O. Zaleski²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.5, pp. 421-430, 2004, DOI:10.3970/cmes.2004.006.421

Abstract In standard boundary element formulations, singular integrals need to be solved as soon as the considered sources coincide with the collocation points at the boundary. Using a desingularized boundary element approach, the sources are distributed on a surface outside the acoustic domain which means that they are never located at the boundary. Consequently, all the resulting kernels are nonsingular which reduces the complexity of the numerical treatment of the boundary integral equations considerably. In the current contribution a desingularized formulation is given for both, the direct and the indirect boundary element method used to solve More >