A linearization of Mieussens's discrete velocity model for kinetic equations

A linearization is developed for Micussens's discrete velocity model (see, e.g., [L. Mieussens, Discrete-velocity models and numerical schemes for the Boltzmann-BGK equation in plane and axisymmetric geometries, J. Comput. Phys. 162 (2000) 429-466]) for kinetic equations. The basic idea is to use a linearized expression of the reference distribution function in the kinetic equation, instead of its exact expression, in the numerical scheme. This modified scheme is applied to various kinetic models, which include the BGK model, the ES-BGK model, the BGK model with velocity-dependent collision frequency, and the recently proposed ES-BGK model with velocity-dependent collision frequency. One-dimensional stationary shock waves and stationary planar Couette flow, which are two benchmark problems for rarefied gas flows, are chosen as test examples. Molecules are modeled as Maxwell molecules and hard sphere molecules. It is found that results from the modified scheme are very similar to results from the original Mieussens's numerical scheme for various kinetic equations in almost all tests we did, while, depending on the test case, 20-40 percent of computational time can be saved. The application of the method is not affected by the Knudsen number and molecular models, but is restricted to lower Mach numbers for the BGK (or the ES-BGK) model with velocity-dependent collision frequency. (c) 2006 Elsevier Masson SAS. All rights reserved.