Numerical Solution of Range-Dependent Acoustic Propagation

The direct global matrix approach can be applied to modeling of range-dependent sound propagation in order to achieve numerically stable and accurate solutions. By solving the global system directly, this method features high efficiency as well as accuracy by avoiding error accumulation. It is an important issue to solve linear systems numerically in the direct global matrix approach, especially for the large-scale problems. An efficient and memory-saving algorithm is developed for solving the global system, in which the global coefficient matrix is treated as a block pentadiagonal matrix. As a result, this numerical model has the ability to solve large-scale problems on regular computers. Numerical examples are also presented to demonstrate the accuracy and efficiency of this method.