三维波动方程的高精度显式格式

摘要

本文基于波速有限物理概念，通过波动方程的解析解构造高精度显式格式的一种显式方法推广到三维波动的时域数值模拟。运用此显式方法构造规则和非规则网格均适用的系列高精度显式格式基本形式，并针对均匀立方体划分。规则网格，详细论述了时空精度皆为2M阶(M为正整数)的稳定显式格式构造途径，最后以构建二阶显式格式(M=1)为例说明整个构建过程，并通过数值试验检验了该显式格式的稳定性和精度。