An enhancement of the fast time-domain boundary element method for the three-dimensional wave equation

被引:11
作者
Takahashi, Toru [1 ]
Tanigawa, Masaki [2 ]
Miyazawa, Naoya [1 ]
机构
[1] Nagoya Univ, Dept Mech Syst Engn, Chikusa Ku, Nagoya, Aichi 4648603, Japan
[2] Shimizu Corp, Koto Ku, Inst Technol, 3-4-17 Etchujima, Tokyo 1358530, Japan
关键词
Boundary element method; Fast multipole method; Wave equation; Time domain; Interpolation; Parameter optimisation; 2-DIMENSIONAL SCATTERING PHENOMENA; TRANSIENT ANALYSIS; ALGORITHM; SOLVERS;
D O I
10.1016/j.cpc.2021.108229
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Our objective is to stabilise and accelerate the time-domain boundary element method (TDBEM) for the three-dimensional wave equation. To overcome the potential time instability, we considered using the Burton-Miller-type boundary integral equation (BMBIE) instead of the ordinary boundary integral equation (OBIE), which consists of the single- and double-layer potentials. In addition, we introduced a smooth temporal basis, i.e. the B-spline temporal basis of order d, whereas d = 1 was used together with the OBIE in a previous study[1]. Corresponding to these new techniques, we generalised the interpolation-based fast multipole method that was developed in [1]. In particular, we constructed the multipole-to-local formula (M2L) so that even for d >= 2 we can maintain the computational complexity of the entire algorithm, i.e. O(N-s(1+delta) N-t), where N-s and N-t denote the number of boundary elements and the number of time steps, respectively, and delta is theoretically estimated as 1/3 or 1/2. The numerical examples indicated that the BMBIE is indispensable for solving the homogeneous Dirichlet problem, but the order d cannot exceed 1 owing to the doubtful cancellation of significant digits when calculating the corresponding layer potentials. In regard to the homogeneous Neumann problem, the previous TDBEM based on the OBIE with d = 1 can be unstable, whereas it was found that the BMBIE with d = 2 can be stable and accurate. The present study will enhance the usefulness of the TDBEM for 3D scalar wave (C) 2021 Elsevier B.V. All rights reserved.
引用
收藏
页数:25
相关论文
共 33 条
[1]   An energy approach to space-time Galerkin BEM for wave propagation problems [J].
Aimi, A. ;
Diligenti, M. ;
Guardasoni, C. ;
Mazzieri, I. ;
Panizzi, S. .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2009, 80 (09) :1196-1240
[2]  
Asakura J., 2009, JSIAM Letters, V1, P52, DOI [10.14495/jsiaml.1.52, DOI 10.14495/JSIAML.1.52]
[3]   A fast hybrid field-circuit simulator for transient analysis of microwave circuits [J].
Aygün, K ;
Fischer, BC ;
Meng, J ;
Shanker, B ;
Michielssen, E .
IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, 2004, 52 (02) :573-583
[4]  
Bowman J.J., 1987, ELECTROMAGNETIC ACOU
[5]   APPLICATION OF INTEGRAL EQUATION METHODS TO NUMERICAL SOLUTION OF SOME EXTERIOR BOUNDARY-VALUE PROBLEMS [J].
BURTON, AJ ;
MILLER, GF .
PROCEEDINGS OF THE ROYAL SOCIETY OF LONDON SERIES A-MATHEMATICAL AND PHYSICAL SCIENCES, 1971, 323 (1553) :201-&
[6]  
Chew W. C., 2001, Fast and Efficient Algorithms in Computational Electromagnetics
[7]  
Chiyoda S., 2019, P C COMPUT ENG SCI, V24
[8]   Analysis of transient wave scattering from rigid bodies using a Burton-Miller approach [J].
Ergin, AA ;
Shanker, B ;
Michielssen, E .
JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA, 1999, 106 (05) :2396-2404
[9]   The plane-wave time-domain algorithm for the past analysis of transient wave phenomena [J].
Ergin, AA ;
Shanker, B ;
Michielssen, E .
IEEE ANTENNAS AND PROPAGATION MAGAZINE, 1999, 41 (04) :39-52
[10]   Fast transient analysis of acoustic wave scattering from rigid bodies using a two-level plane wave time domain algorithm [J].
Ergin, AA ;
Shanker, B ;
Michielssen, E .
JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA, 1999, 106 (05) :2405-2416