Abstract:In this paper, we propose an algorithm that imposes macroscopic physical constraints with Lagrange multiplier approach in implementing the Maxwell boundary condition within the framework of the discrete velocity method. For the simulation of rarefied gas flows in the presence of solid walls with complex geometry, the distribution function in the reflection region of the wall surface needs to be constructed in the discrete velocity space, to fulfill the specular reflection in the Maxwell boundary condition. The construction process should not consist of interpolation only, but include certain macroscopic physical constraints at the wall surface, so as to correctly account for gas-surface interaction on a macroscopic level. We demonstrate that for the specular reflection, keeping the symmetry of the first three moments of the distribution function between the incident and reflected region is sufficient for maintaining the conservation of mass, momentum, and energy at the wall surface. Furthermore, to strictly satisfy macroscopic physical constraints, a Lagrange multiplier method is introduced into the construction of the distribution function to correct the pure interpolation solution. In addition, the construction process requires the inversion of a large and sparse matrix (of dimension N × N, where N is the number of points in the velocity space). To improve the computational efficiency, the matrix inversion is converted into that of a much smaller matrix, i.e. (D+2) × (D+2) in the D-dimensional physical space. A series of numerical experiments are conducted to examine the performance of the proposed algorithm under different flow conditions. We demonstrate that the results obtained by the proposed algorithm are more accurate than the pure interpolation solution, comparing with the benchmark data. Moreover, after the validation of our results with previous studies, we find that the method significantly enhances the conservation of total mass and energy, especially for flows in an enclosed domain.

A boundary matching micro/macro decomposition for kinetic equations

Micro-macro decomposition based asymptotic-preserving numerical schemes and numerical moments conservation for collisional nonlinear kinetic equations

A Micro-Macro Decomposition-Based Asymptotic-Preserving Scheme for the Multispecies Boltzmann Equation

Kinetic/Fluid Micro-Macro Numerical Scheme for a Two Component Gas Mixture

WELL-BALANCED AND ASYMPTOTIC PRESERVING SCHEMES FOR KINETIC MODELS

Boundary Conditions for Kinetic Theory Based Models I: Lattice Boltzmann Models

A class of asymptotic-preserving schemes for kinetic equations and related problems with stiff sources

Accurate macroscale modelling of spatial dynamics in multiple dimensions

Boundary Scheme for a Discrete Kinetic Approximation of the Navier–Stokes Equations

Two-Dimensional System of Moment Equations and Macroscopic Boundary Conditions Depending on the Velocity of Movement and the Surface Temperature of a Body Moving in Fluid

A new local and explicit kinetic method for linear and non-linear convection-diffusion problems with finite kinetic speeds: II. Multi-dimensional case

Maxwell boundary condition for discrete velocity methods: macroscopic physical constraints and Lagrange multiplier-based implementation

Microscopic determination of macroscopic boundary conditions in Newtonian liquids

A Numerical Scheme for the Quantum Boltzmann Equation Efficient in the Fluid Regime

Hydrodynamic limits for kinetic equations preserving mass, momentum and energy: a spectral and unified approach in the presence of a spectral gap

A Multiscale Kinetic-Fluid Solver with Dynamic Localization of Kinetic Effects

An asymptotic preserving scheme for front propagation in a kinetic reaction-transport equation

An Asymptotic Preserving Scheme for the ES-BGK Model of the Boltzmann Equation

Kinetic equations for systems with long-range interactions: a unified description

High Order Asymptotic Preserving DG-IMEX Schemes for Discrete-Velocity Kinetic Equations in a Diffusive Scaling

A Smooth Transition Model Between Kinetic and Diffusion Equations