Resumen
The finite element method is a popular numerical method in engineering applications. However, there is not enough research about the finite element method in underwater sound propagation. The finite element method can achieve high accuracy and great universality. We aim to develop a three-dimensional finite element model focusing on underwater sound propagation. As the foundation of this research, we put forward a finite element model in the Cartesian coordinate system for a sound field in a two-dimensional environment. We firstly introduce the details of the implementation of the finite element model, as well as different methods to deal with boundary conditions and a comparison of these methods. Then, we use four-node quadrilateral elements to discretize the physical domain, and apply the perfectly matched layer approach to deal with the infinite region. After that, we apply the model to underwater sound propagation problems including the wedge-shaped waveguide benchmark problem and the problem where the bathymetry consists of a sloping region and a flat region. The results by the presented finite element model are in excellent agreement with analytical and benchmark numerical solutions, implying that the presented finite element model is able to solve complex two-dimensional underwater sound propagation problems accurately. In the end, we compare the finite element model with the popular normal mode model KRAKEN by calculating sound fields in Pekeris waveguides, and find that the finite element model has better universality than KRAKEN.