High-order kinetic relaxation schemes as high-accuracy Poisson solvers

We present a new approach to find accurate solutions to the Poisson equation, as obtained from the steady-state limit of a diffusion equation with strong source terms. For this purpose, we start from Boltzmann's kinetic theory and investigate the influence of higher-order terms on the resulting macroscopic equations. By performing an appropriate expansion of the equilibrium distribution, we provide a method to remove the unnecessary terms up to a desired order and show that it is possible to find, with high level of accuracy, the steady-state solution of the diffusion equation for sizeable Knudsen numbers. In order to test our kinetic approach, we discretize the Boltzmann equation and solve the Poisson equation, spending up to six order of magnitude less computational time for a given precision than standard lattice Boltzmann methods (LBMs).