波动切比雪夫谱元模拟的时间积分方法研究

由于在波动切比雪夫谱元模拟中使用隐式时间积分方法存在计算效率较低、不易施加边界条件以及波动输入不便的问题,而中心差分法是能够平衡精度和计算效率的较优选择,并以一维波动模型探讨相应的波动输入方法和时域积分稳定条件。通过输入Ricker子波以及复合正弦波的数值算例证实了方法的有效性,并接着分析不同计算参数对模拟精度的影响。结果表明:空间上一个最短波长尺度内的谱元节点数、单元阶次分别取不低于5和4时能够达到比较理想的精度,而时间步长变化对精度的影响不大。