GENERALIZED MULTISCALE FINITE ELEMENT METHOD FOR THE STEADY STATE LINEAR BOLTZMANN EQUATION

The Boltzmann equation, as a model equation in statistical mechanics, is used to describe the statistical behavior of a large number of particles driven by the same physics laws. Depending on the media and the particles to be modeled, the equation has slightly different forms. In this article, we investigate a model Boltzmann equation with highly oscillatory media in the small Knudsen number regime and study the numerical behavior of the generalized multiscale finite element method (GMsFEM) in the fluid regime when high oscillation in the media presents. The GMsFEM is a general approach [E. Chung, Y. Efendiev, and T. Y. Hou, T. Comput. Phys., 320 (2016), pp. 69-95] to numerically treat equations with multiscale structures. The method is divided into the offline and online steps. In the offline step, basis functions are prepared from a snapshot space via a well-designed generalized eigenvalue problem (GEP), and these basis functions are then utilized to patch up for a solution through DG formulation in the online step to incorporate specific boundary and source information. We prove the well-posedness of the method on the Boltzmann equation and show that the GEP formulation provides a set of optimal basis functions that achieve spectral convergence. Such convergence is independent of the oscillation in the media, or the smallness of the Knudsen number, making it one of the few methods that simultaneously achieve numerical homogenization and asymptotic preserving properties across all scales of oscillations and the Knudsen number.