Evaluating high order discontinuous Galerkin discretization of the Boltzmann collision integral in O(N^2) operations using the discrete Fourier transform

From MaRDI portal
Publication:2323402

DOI10.3934/KRM.2019027zbMATH Open1421.76201arXiv1801.05892OpenAlexW2946764832MaRDI QIDQ2323402FDOQ2323402

Jeffrey Limbacher, Alexander M. Alekseenko

Publication date: 2 September 2019

Published in: Kinetic and Related Models (Search for Journal in Brave)

Abstract: We present a numerical algorithm for evaluating the Boltzmann collision operator with O(N2) operations based on high order discontinuous Galerkin discretizations in the velocity variable. To formulate the approach, Galerkin projection of the collision operator is written in the form of a bilinear circular convolution. An application of the discrete Fourier transform allows to rewrite the six fold convolution sum as a three fold weighted convolution sum in the frequency space. The new algorithm is implemented and tested in the spatially homogeneous case, and results in a considerable improvement in speed as compared to the direct evaluation. Simultaneous and separate evaluations of the gain and loss terms of the collision operator were considered. Less numerical error was observed in the conserved quantities with simultaneous evaluation.


Full work available at URL: https://arxiv.org/abs/1801.05892




Recommendations





Cited In (5)





This page was built for publication: Evaluating high order discontinuous Galerkin discretization of the Boltzmann collision integral in \(\mathcal{O}(N^2)\) operations using the discrete Fourier transform

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2323402)