Domain decomposition based on the spectral element method for frequency-domain computational elastodynamics  

作　　者：Linlin SHI Mingwei ZHUANG Yuanguo ZHOU Na LIU Qinghuo LIU 

机构地区：[1]Institute of Electromagnetics and Acoustics and Key Laboratory of Electromagnetic Wave Science and Detection Technology,Xiamen University,Xiamen 361005,China [2]The Fifth Electronics Research Institute of Ministry of Industry and Information Technology,Guangzhou 510610,China [3]College of Communication and Information Engineering,Xi’an University of Science and Technology,Xi’an 710054,China [4]Department of Electrical and Computer Engineering,Duke University,Durham,NC 27708,USA

出　　处：《Science China Earth Sciences》2021年第3期388-403,共16页中国科学（地球科学英文版）

基　　金：the National Natural Science Foundation of China(Grant Nos.41390453,62001408);the National Key R&D Program of the Ministry of Science and Technology of China(Grant No.2018YFC0603503);the Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund(Grant No.U1501501)。

摘　　要：We propose a domain decomposition method based on the spectral element method(DDM-SEM)for elastic wave computation in frequency domain.It combines the high accuracy of the spectral element method and the high degree of parallelism of a domain decomposition technique,which makes this method suitable for accurate and efficient simulations of large scale problems in elastodynamics.In the DDM-SEM,the original large-scale problem is divided into a number of well designed subdomains.We use the spectral element method independently for each subdomain,and the neighboring subdomains are connected by a frequency-domain version of Riemann transmission condition(RTC)for elastic waves.For the proposed method,we can employ the non-conforming meshes and different interpolation orders in different subdomains to maximize the efficiency.By separating the internal and boundary unknowns of each subdomain,an efficient and naturally parallelizable block LDU direct solver is developed to solve the final system matrix.Numerical experiments verify its accuracy and efficiency,and show that the proposed DDM-SEM can be a promising numerical tool for accurately and effectively solving large and multi-scale problems of elastic waves.It is potentially valuable for the frequency domain seismic inversion where multiple source illuminations are required.

关 键 词：Elastic waves Domain decomposition Spectral element method Frequency domain Numerical simulations 

分 类 号：P618.13[天文地球—矿床学] P631.4[天文地球—地质学]