Japan Geoscience Union Meeting 2019

Presentation information

[J] Oral

S (Solid Earth Sciences ) » S-TT Technology & Techniques

[S-TT47] Creating future of solid Earth science with high performance computing (HPC)

Sun. May 26, 2019 10:45 AM - 12:15 PM 304 (3F)

convener:Takane Hori(R&D Center for Earthquake and Tsunami, Japan Agency for Marine-Earth Science and Technology), Yuji Yagi(Graduate School of Life and Environmental Sciences, University of Tsukuba, Tsukuba), Katsuhiko Shiomi(National Research Institute for Earth Science and Disaster Resilience), Chairperson:Takane Hori

11:15 AM - 11:30 AM

[STT47-03] Optimal spectral element method for elastic wave computation

*Kei Hasegawa1, Nobuaki Fuji2, Kensuke Konishi3 (1.Earthquake Research Institute, University of Tokyo, 2.Institut de Physique du Globe de Paris, 3.Institute of Earth Sciences, Academia Sinica)

Keywords:Elastic wave, Numerical analysis, Finite-Element Method, Spectral Element Method, Error optimization, Computational seismology

We introduce new numerical integration operators which compose the mass and stiffness matrices of a modified spectral element method for simulation of elastic wave propagation. While these operators use the same quadrature nodes as does the original spectral element method, they are designed in order that their lower-order contributions to the numerical dispersion error cancel each other. As a result, the modified spectral element method yields two extra-orders of accuracy, and is comparable to the original method of one order higher. The theoretical results are confirmed by numerical dispersion analysis and examples of computation of waveforms using our operators. Replacing the ordinary operators by those proposed in this study could be a non-expensive solution to improve the accuracy.

*A significant portion of this presentation was published in Computer Methods in Applied Mechanics and Engineering (doi: 10.1016/j.cma.2018.07.025).