Numerical simulation of the large-scale high frequency sound field is a computational challenging task. To solve the difficulty that the traditional boundary collocation methods are not easy to be applied to large-scale problems thanks to the resulting large-scale fully-populated matrix, a dual-level singular boundary method is proposed in this study. By introducing a dual level structure,the fully-populated matrix is transformed to a large-scale locally supported sparse matrix on fine mesh. The bottleneck ...