谱元法与透射边界的配合使用及其稳定性研究.pdfVIP

第32 卷第11 期 Vol.32 No.11 工 程 力 学 2015 年 11 月 Nov. 2015 ENGINEERING MECHANICS 40 文章编号：1000-4750(2015)11-0040-11 谱元法与透射边界的配合使用及其稳定性研究 1 1 2 戴志军 ，李小军 ，侯春林 (1. 中国地震局地球物理研究所，北京 100080 ；2. 环境保护部核与辐射安全中心，北京 100082) 摘 要：在无限域波动模拟中引入透射边界条件时，目前多将边界上的透射公式与内域的有限元法结合使用，其 计算精度由有限元方法决定，而谱元法因结合有限元和频谱法的优势则比有限元空间域积分具有更高的计算精 度。该文基于谱元法非等距网格划分特性，研究了内域的谱元法与边界上的透射公式结合的理论方法，给出了相 应的透射公式使用方法，并基于建立的谱元法波动数值模型探讨了透射公式的稳定性问题。研究表明：空间域插 值系数需控制在一个合理范围内，空间域插值方法相对于时间域插值方法更为稳定，高频失稳出现可能性相对较 小；Gamma 算子的使用可提高模拟的精度，采用Gamma 算子后对于高阶透射公式仍可出现低频漂移现象，可结 合降阶消漂的方式实现稳定精度高的透射边界应用。 关键词：谱元法；有限元法；透射公式；透射边界；数值稳定性 中图分类号：P315.9 文献标志码：A doi: 10.6052/j.issn.1000-4750.2014.03.0196 A COMBINATION USAGE OF TRANSMITTING FORMULA AND SPECTRAL ELEMENT METHOD AND THE STUDY OF ITS STABILITY 1 1 2 DAI Zhi-jun , LI Xiao-jun , HOU Chun-lin (1. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China; 2. Nuclear and Radiation Safety Centre, Ministry of Environmental Protection of the People’s Republic of China, Beijing 100082, China) Abstract: For the simulation of the problem on infinite domain, finite element method (FEM) is usually combined with transmitting formulas, and the precision of the simulation is mainly determined by the FEM methods. Spectral element method (SEM) takes advantages of both spectral method and FEM, and achieves higher precision for integration in space domain compared with FEM. This paper investigates the application of transmitt