In this paper, we investigate the numerical solution of the two-dimensional fractional Laplacian wave equations. After splitting out the Riesz fractional derivatives from the fractional Laplacian, we treat the Riesz fractional derivatives with an implicit scheme while solving the rest part explicitly. Thanks to the tensor structure of the Riesz fractional derivatives, a splitting alternative direction implicit (S-ADI) scheme is proposed by incorporating an ADI remainder. Then the Gohberg-Semencul formula, combined with fast Fourier transform, is proposed to solve the derived Toeplitz linear systems at each time integration. Theoretically, we demonstrate that the S-ADI scheme is unconditionally stable and possesses second-order accuracy. Finally, numerical experiments are performed to demonstrate the accuracy and efficiency of the S-ADI scheme.

28 pages, 27 figures