Towards an ultra efficient kinetic scheme Part II: The high order case
arXiv:1212.0204 · doi:10.1016/j.jcp.2013.07.017
Abstract
In a recent paper we presented a new ultra efficient numerical method for solving kinetic equations of the Boltzmann type (G. Dimarco, R. Loubere, Towards an ultra efficient kinetic scheme. Part I: basics on the 689 BGK equation, J. Comp. Phys., (2013), http://dx.doi.org/10.1016/j.jcp.2012.10.058). The key idea, on which the method relies, is to solve the collision part on a grid and then to solve exactly the transport part by following the characteristics backward in time. On the contrary to classical semi-Lagrangian methods one does not need to reconstruct the distribution function at each time step. This allows to tremendously reduce the computational cost and to perform efficient numerical simulations of kinetic equations up to the six dimensional case without parallelization. However, the main drawback of the method developed was the loss of spatial accuracy close to the fluid limit. In the present work, we modify the scheme in such a way that it is able to preserve the high order spatial accuracy for extremely rarefied and fluid regimes. In particular, in the fluid limit, the method automatically degenerates into a high order method for the compressible Euler equations. Numerical examples are presented which validate the method, show the higher accuracy with respect to the previous approach and measure its efficiency with respect to well known schemes (Direct Simulation Monte Carlo, Finite Volume, MUSCL, WENO).
References in corpus (6)
- Locally refined discrete velocity grids for stationary rarefied flow simulations
- Local discrete velocity grids for deterministic rarefied flow simulations
- Arbitrarily high order Convected Scheme solution of the Vlasov-Poisson system
- A high order cell-centered semi-Lagrangian scheme for multi-dimensional kinetic simulations of neutral gas flows
- A hybrid OpenMP and MPI implementation of a conservative spectral method for the Boltzmann equation
- High Performance Computing with a Conservative Spectral Boltzmann Solver
Cited by in corpus (5)
- Approximation of the Boltzmann Collision Operator Based on Hermite Spectral Method
- Conservative semi-Lagrangian schemes for kinetic equations Part II: Applications
- Addressing the gas kinetics Boltzmann equation with branching-path statistics
- Implicit Discontinuous Galerkin Method for the Boltzmann Equation
- High order finite volume schemes with IMEX time stepping for the Boltzmann model on unstructured meshes