提出了波动时域数值模拟的一种新的显式方法,并通过一维非规则网格节点递推公式的建立说明该方法.为了阐明应用该法构建高精度、稳定递推公式的可行性,详细论述了一维规则网格标量波动数值模拟的精度和稳定性;提出了构建时空精度皆为2M阶(M为正整数)的稳定递推公式的技术途径,并以构建二阶(M=1)和四阶(M=2)公式为例予以说明.最后,通过算例详细说明了理论结果.
A new explicit method for numerical simulation of wave motion in time domain is first presented and illustrated by constructing the recursion formulas at nodal points of an irregular grid for 1-D wave equation. To demonstrate the possibility to develop stable and high-order formulas with the method, accuracy and stability of the numerical simulation are discussed in detail for the 1-D scalar wave equation in a uniform grid, and an approach is then proposed to construct the stable formulas with 2M~order accuracy both in time and space. The theoretical result is illustrated by constructing the formulas of second order (M = 1) and fourth order (M = 2), and is validated by a Series of numerical tests. The implication of this study is finally disscussed for further works.