Weidong Lei¹, Hai Zhou^1,*, Hongjun Li², Rui Chen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.126, No.2, pp. 577-597, 2021, DOI:10.32604/cmes.2021.011857 - 21 January 2021

Abstract In order to simulate the propagation process of subway vibration of parallel tunnels in semi-infinite rocks or soils, time domain boundary element method (TD-BEM) formulation for analyzing the dynamic response of twin-parallel circular tunnels in an elastic semi-infinite medium is developed in this paper. The time domain boundary integral equations of displacement and stress for the elastodynamic problem are presented based on Betti’s reciprocal work theorem, ignoring contributions from initial conditions and body forces. In the process of establishing time domain boundary integral equations, some virtual boundaries are constructed between finite boundaries and the free… More >

Liu Liqi¹, Wang Haitao^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.109-110, No.4, pp. 303-324, 2015, DOI:10.3970/cmes.2015.109.303

Abstract 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 More >

Mikael A. Langthjem¹, Masami Nakano²

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.4, pp. 227-241, 2013, DOI:10.3970/cmes.2013.096.227

Abstract This paper is concerned with a mathematical model of a simple axisymmetric silencer-like model, consisting of a hole-tone feedback system equipped with a tailpipe. The unstable shear layer is modeled via a discrete vortex method, based on axisymmetric vortex rings. The aeroacoustic model is based on the Powell- Howe theory of vortex sound. Boundary integrals are discretized via the boundary element method; but the tailpipe is represented by the exact (one-dimensional) solution. It is demonstrated though numerical examples that this numerical model can display lock-in of the self-sustained flow oscillations to the resonant acoustic oscillations. More >

D. Soares Jr.¹, M. P. Vinagre²

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.1, pp. 1-8, 2008, DOI:10.3970/cmes.2008.025.001

Abstract This work presents an alternative procedure to compute time-domain electromagnetic fields. The Boundary Element Method is here adopted to numerically analyze wave propagation problems, computing just a so-called primary field (either the electric or the magnetic field can be selected as primary field; the complementary field is here named secondary field). The secondary field is obtained following Maxwell's equations, i.e., considering space derivatives of the primary field (computed by the Complex Variable Method) and time integration procedures. This methodology is more efficient and flexible since fewer systems of equations must be solved at each time-step. More >

Ch. Zhang², A. Savaidis³

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.5, pp. 603-618, 2003, DOI:10.3970/cmes.2003.004.603

Abstract A novel non-hypersingular time-domain traction BEM is presented for three-dimensional (3-D) transient elastodynamic crack analysis. The initial-boundary value problem is formulated as a set of non-hypersingular time-domain traction boundary integral equations (BIEs). To solve the time-domain traction BIEs, a time-stepping scheme based on the convolution quadrature formula of Lubich (1988a,b; 1994) for temporal discretization and a collocation method for spatial discretization is adopted. Numerical examples are given for an unbounded solid with a penny-shaped crack under a tensile and shear impact loading. A comparison of the present time-domain BEM with the conventional one shows that More >