Videos

High order conservative semilagrangian methods for kinetic equations

Presenter
January 13, 2022
Abstract
In this talk, we overview a class of high order conservative semi-Lagrangian schemes for BGK model. The schemes are constructed by coupling the conservative non-oscillatory reconstruction technique with a conservative treatment of the collision term, obtained by either a discrete Maxwellian or by an L2-minimization technique. Due to the semi-Lagrangian nature, the time step is not restricted by a CFL-type condition, while the implicit treatment of the relaxation term based on Runge-Kutta or BDF time discretization enables us to avoid the stiffness problem coming from a small Knudsen number. Because of L-stability and exact conservation, the resulting scheme is asymptotic preserving for the underlying fluid dynamic limit. Several test cases confirm the accuracy and robustness of the methods, and the AP property of the schemes. The method has been extended to the treatment of inert gas mixtures, and applied to compare different models in various regimes. In general, such approaches use fixed velocity grids, and one must secure a sufficient number of grid points in phase space to resolve the structure of the distribution function. When dealing with high Mach number problems, where large variation of mean velocity and temperature are present in the domain under consideration, the computational cost and memory allocation requirements become prohibitively large. Local velocity grid methods have been developed to overcome such difficulty in the context of Eulerian based schemes. In this talk, we introduce a velocity adaption technique for the semi-Lagrangian scheme applied to the BGK model. The velocity grids will be set locally in time and space. We apply a weighted minimization approach to impose global conservation, generalizing the L2-minimization technique introduced in. The efficiency of the proposed scheme is illustrated numerically. An additional application of conservative SL schemes concerns the numerical simulation of Vlasov-type equations. Here conservation of the scheme will provide some advantage over standard non conservative schemes for long time computation. The research is conducted with the following collaborators: S. Boscarino, S. Y. Cho, M. Groppi, S.-B. Yun, JM Qiu, and T. Xiong.