A fast Fourier spectral method for wave kinetic equation

The central object in wave turbulence theory is the wave kinetic equation (WKE), which is an evolution equation for wave action density and acts as the wave analog of the Boltzmann kinetic equations for particle interactions. Despite recent exciting progress in the theoretical aspects of the WKE, numerical developments have lagged behind. In this paper, we introduce a fast Fourier spectral method for solving the WKE. The key idea lies in reformulating the high-dimensional nonlinear wave kinetic operator as a spherical integral, analogous to classical Boltzmann collision operator. The conservation of mass and momentum leads to a double convolution structure in Fourier space, which can be efficiently handled using the fast Fourier transform. We demonstrate the accuracy and efficiency of the proposed method through several numerical tests in both 2D and 3D, revealing some interesting and unique features of this equation.