Higher-order triangular spectral element method with optimized cubature points for seismic wavefield modeling
文献类型：期刊论文
| 作者 | Liu, Youshan1; Teng, Jiwen1; Xu, Tao1,2; Badal, Jose3 |
| 刊名 | JOURNAL OF COMPUTATIONAL PHYSICS
 |
| 出版日期 | 2017-05-01 |
| 卷号 | 336页码:458-480 |
| 关键词 | Triangular spectral element method Dispersion relation Optimized cubature points Positive integration weights Courant-Friedrichs-Lewy numbers |
| ISSN号 | 0021-9991 |
| DOI | 10.1016/j.jcp.2017.01.069 |
| 英文摘要 | The mass-lumped method avoids the cost of inverting the mass matrix and simultaneously maintains spatial accuracy by adopting additional interior integration points, known as cubature points. To date, such points are only known analytically in tensor domains, such as quadrilateral or hexahedral elements. Thus, the diagonal-mass-matrix spectral element method (SEM) in non-tensor domains always relies on numerically computed interpolation points or quadrature points. However, only the cubature points for degrees 1 to 6 are known, which is the reason that we have developed a p-norm-based optimization algorithm to obtain higher-order cubature points. In this way, we obtain and tabulate new cubature points with all positive integration weights for degrees 7 to 9. The dispersion analysis illustrates that the dispersion relation determined from the new optimized cubature points is comparable to that of the mass and stiffness matrices obtained by exact integration. Simultaneously, the Lebesgue constant for the new optimized cubature points indicates its surprisingly good interpolation properties. As a result, such points provide both good interpolation properties and integration accuracy. The Courant-FriedrichsLewy (CFL) numbers are tabulated for the conventional Fekete-based triangular spectral element (TSEM), the TSEM with exact integration, and the optimized cubature-based TSEM (OTSEM). A complementary study demonstrates the spectral convergence of the OTSEM. A numerical example conducted on a half-space model demonstrates that the OTSEM improves the accuracy by approximately one order of magnitude compared to the conventional Fekete-based TSEM. In particular, the accuracy of the 7th-order OTSEM is even higher than that of the 14th-order Fekete-based TSEM. Furthermore, the OTSEM produces a result that can compete in accuracy with the quadrilateral SEM (QSEM). The high accuracy of the OTSEM is also tested with a non-flat topography model. In terms of computational efficiency, the OTSEM is more efficient than the Fekete-based TSEM, although it is slightly costlier than the QSEM when a comparable numerical accuracy is required. (C) 2017 Elsevier Inc. All rights reserved. |
| WOS关键词 | LUMPED FINITE-ELEMENTS ; GAUSS-LOBATTO INTEGRATION ; CONJUGATE-GRADIENT METHOD ; DISPERSION ANALYSIS ; POLYNOMIAL INTERPOLATION ; DIFFERENCE SCHEMES ; SYMPLECTIC SCHEME ; QUADRATURE POINTS ; COMPUTING FEKETE ; PROPAGATION |
| 资助项目 | National Key Research and Development Program of China[2016YFC0600101] ; National Key Research and Development Program of China[2016YFC0600201] ; National Natural Science Foundation of China[41604076] ; National Natural Science Foundation of China[41674102] ; National Natural Science Foundation of China[41674095] ; National Natural Science Foundation of China[41522401] ; National Natural Science Foundation of China[41374062] ; National Natural Science Foundation of China[41404073] ; China Earthquake Administration[201408023] ; China Postdoctoral Science Foundation[2016M600128] |
| WOS研究方向 | Computer Science ; Physics |
| 语种 | 英语 |
| WOS记录号 | WOS:000397362800023 |
| 出版者 | ACADEMIC PRESS INC ELSEVIER SCIENCE |
| 资助机构 | National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; China Earthquake Administration ; China Earthquake Administration ; China Earthquake Administration ; China Earthquake Administration ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; China Earthquake Administration ; China Earthquake Administration ; China Earthquake Administration ; China Earthquake Administration ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; China Earthquake Administration ; China Earthquake Administration ; China Earthquake Administration ; China Earthquake Administration ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Key Research and Development Program of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; National Natural Science Foundation of China ; China Earthquake Administration ; China Earthquake Administration ; China Earthquake Administration ; China Earthquake Administration ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation ; China Postdoctoral Science Foundation |
| 源URL | [http://ir.iggcas.ac.cn/handle/132A11/93370]  |
| 专题 | 地质与地球物理研究所_岩石圈演化国家重点实验室 |
| 通讯作者 | Liu, Youshan |
| 作者单位 | 1.Chinese Acad Sci, Inst Geol & Geophys, State Key Lab Lithospher Evolut, Beijing 100029, Peoples R China 2.CAS Ctr Excellence Tibetan Plateau Earth Sci, Beijing 100101, Peoples R China 3.Univ Zaragoza, Sci B, Phys Earth, Pedro Cerbuna 12, E-50009 Zaragoza, Spain |
| 推荐引用方式 GB/T 7714 | Liu, Youshan,Teng, Jiwen,Xu, Tao,et al. Higher-order triangular spectral element method with optimized cubature points for seismic wavefield modeling[J]. JOURNAL OF COMPUTATIONAL PHYSICS,2017,336:458-480. |
| APA | Liu, Youshan,Teng, Jiwen,Xu, Tao,&Badal, Jose.(2017).Higher-order triangular spectral element method with optimized cubature points for seismic wavefield modeling.JOURNAL OF COMPUTATIONAL PHYSICS,336,458-480. |
| MLA | Liu, Youshan,et al."Higher-order triangular spectral element method with optimized cubature points for seismic wavefield modeling".JOURNAL OF COMPUTATIONAL PHYSICS 336(2017):458-480. |
入库方式： OAI收割
来源：地质与地球物理研究所
浏览0
下载0
收藏0
其他版本
除非特别说明，本系统中所有内容都受版权保护，并保留所有权利。