A FAST SPECTRAL METHOD FOR THE BOLTZMANN COLLISION OPERATOR WITH GENERAL COLLISION KERNELS

作者：Gamba Irene M; Haack Jeffrey R; Hauck Cory D; Hu Jingwei

来源：SIAM Journal on Scientific Computing, 2017, 39(4): B658-B674.

DOI：10.1137/16M1096001

• Boltzmann collision integral
• spectral method
• convolution
• fast Fourier transform
• Lebedev quadrature

We propose a simple fast spectral method for the Boltzmann collision operator with general collision kernels. In contrast to the direct spectral method [L. Pareschi and G. Russo, SIAM J. Numer. Anal., 37 (2000), pp. 1217-1245; I. M. Gamba and S. H. Tharkabhushanam, J. Comput. Phys., 228 (2009), pp. 2012-2036], which requires O(N-6) memory to store precomputed weights and has O(N-6) numerical complexity, the new method has complexity O(MN4 log N), where N is the number of discretization points in each of the three velocity dimensions and M is the total number of discretization points on the sphere and M << N-2. Furthermore, it requires no precomputation for the variable hard sphere model and only O(MN4) memory to store precomputed functions for more general collision kernels. Although a faster spectral method is available [C. Mouhot and L. Pareschi, Math. Comp., 75 (2006), pp. 1833-1852] (with complexity O(MN3 log N)), it works only for hard sphere molecules, thus limiting its use for practical problems. Our new method, on the other hand, can apply to arbitrary collision kernels. A series of numerical tests is performed to illustrate the efficiency and accuracy of the proposed method.

• 出版日期2017