|Source||CMES: Computer Modeling in Engineering & Sciences, Vol. 109, No. 4, pp. 303-324, 2015|
|Download||Full length paper in PDF format. Size = 1,334,833 bytes|
|Keywords||Boundary element method, Radial integration method, Time domain, Non-homogeneous problems, Wave propagation.|
This paper presents a three-dimensional (3-D) boundary element method (BEM) scheme based on the Radial Integration Method (RIM) for wave propagation analysis of continuously non-homogeneous problems. The Kelvin fundamental solutions are adopted to derive the boundary-domain integral equation (BDIE). The RIM proposed by Gao (Engineering Analysis with Boundary Elements 2002; 26(10):905-916) is implemented to treat the domain integrals in the BDIE so that only boundary discretization is required. After boundary discretization, a set of second-order ordinary differential equations with respect to time variable are derived, which are solved using the Wilson-q method. Main advantages of the proposed method are that 1) it can treat wave propagations in non-homogeneous domains with only boundary mesh required, and that 2) coefficient matrices arising from the BEM are evaluated and stored only once so that solving large-scale problems with huge time steps is possible. In the numerical examples, the present method is tested in terms of accuracy, capacity to treat non-homogeneous problems and large-scale potentials.