Reducing the Dispersion Error in Acoustic Wave Modeling with Modified Integration Rules

Bin Yue,Murthy N. Guddati
DOI: https://doi.org/10.1121/1.4785055
2004-01-01
Abstract:Numerical analysis of transient acoustic wave propagation is often performed using finite element or finite difference methods along with central difference time stepping. Such analysis results in wave velocities that are different from exact ones. For conventional algorithms, this dispersion error is second order with respect to discretization parameters. In this paper, modified integration rules are combined with a modified central difference scheme to reduce this error to fourth order. Essentially, the locations of integration points used for evaluating the system matrices are selected such that the wave velocity error is minimized. It was found that the optimal integration points are square roots of 2/3 for the stiffness matrix, which eliminates the anisotropy in the second-order error. The optimal integration points for mass matrix, which depend on the wave velocity, the mesh size, and the time-step size, are also identified to completely eliminate the second-order error. The resulting dispersion error is thus reduced to fourth order. Numerical experiments illustrate that the proposed method has superior performance over existing methods, not only for structured square meshes, but also for unstructured meshes.
What problem does this paper attempt to address?