TY - JOUR
T1 - A pseudospectral penalty scheme for 2D isotropic elastic wave computations
AU - Feng, Ko An
AU - Teng, Chun Hao
AU - Chen, Min Hung
N1 - Funding Information:
Acknowledgements The authors would like to thank the support from the National Science Council of Taiwan (grant number: NSC 95-2120-M-001-003). Chun-Hao Teng would like to express his gratitude to Dr. Yuh-Lin Wang and Dr. Juen-Kai Wang in the Institute of Atomic and Molecular Science, Academia Sinica, Taiwan, for their great help. The authors also would like to thank an anonymous reviewer for useful comments that helped to improve both the contents and the style of the paper. The high-performance computing support from Computer and Network Center at National Cheng Kung University is highly acknowledged.
PY - 2007/12
Y1 - 2007/12
N2 - In this paper, we present a pseudospectral scheme for solving 2D elastic wave equations. We start by analyzing boundary operators leading to the well-posedness of the problem. In addition, equivalent characteristic boundary conditions of common physical boundary conditions are discussed. These theoretical results are further employed to construct a Legendre pseudospectral penalty scheme based on a tensor product formulation for approximating waves on a general curvilinear quadrilateral domain. A stability analysis of the scheme is conducted for the case where a straight-sided quadrilateral element is used. The analysis shows that, by properly setting the penalty parameters, the scheme is stable at the semi-discrete level. Numerical experiments for testing the performance of the scheme are conducted, and the expected p- and h-convergence patterns are observed. Moreover, the numerical computations also show that the scheme is time stable, which makes the scheme suitable for long time simulations.
AB - In this paper, we present a pseudospectral scheme for solving 2D elastic wave equations. We start by analyzing boundary operators leading to the well-posedness of the problem. In addition, equivalent characteristic boundary conditions of common physical boundary conditions are discussed. These theoretical results are further employed to construct a Legendre pseudospectral penalty scheme based on a tensor product formulation for approximating waves on a general curvilinear quadrilateral domain. A stability analysis of the scheme is conducted for the case where a straight-sided quadrilateral element is used. The analysis shows that, by properly setting the penalty parameters, the scheme is stable at the semi-discrete level. Numerical experiments for testing the performance of the scheme are conducted, and the expected p- and h-convergence patterns are observed. Moreover, the numerical computations also show that the scheme is time stable, which makes the scheme suitable for long time simulations.
UR - http://www.scopus.com/inward/record.url?scp=35548948031&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=35548948031&partnerID=8YFLogxK
U2 - 10.1007/s10915-007-9154-8
DO - 10.1007/s10915-007-9154-8
M3 - Article
AN - SCOPUS:35548948031
VL - 33
SP - 313
EP - 348
JO - Journal of Scientific Computing
JF - Journal of Scientific Computing
SN - 0885-7474
IS - 3
ER -