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.

Maximum Principle in Linear Finite Element Approximations of Anisotropic Diffusion–convection–reaction Problems

A study on the conditioning of finite element equations with arbitrary anisotropic meshes via a density function approach

Finite Element Methods Respecting the Discrete Maximum Principle for Convection-Diffusion Equations

Adaptive least-squares methods for convection-dominated diffusion-reaction problems

Third order Maximum-Principle-Satisfying Direct discontinuous Galerkin methods for time dependent convection diffusion equations on unstructured triangle mesh

A symmetric and coercive finite volume scheme preserving the discrete maximum principle for anisotropic diffusion equations on star-shaped polygonal meshes

Maximum principle preserving nonlocal diffusion model with Dirichlet boundary condition

The Finite Volume Scheme Preserving Maximum Principle for Two-Dimensional Time-Fractional Fokker–Planck Equations on Distorted Meshes

A study on anisotropic mesh adaptation for finite element approximation of eigenvalue problems with anisotropic diffusion operators

Weak discrete maximum principle of isoparametric finite element methods in curvilinear polyhedra

Maximum Principles for P1-Conforming Finite Element Approximations of Quasi-Linear Second Order Elliptic Equations

Discrete strong extremum principles for finite element solutions of advection‐diffusion problems with nonlinear corrections

A monotone $Q^1$ finite element method for anisotropic elliptic equations

Optimal H1 spatial convergence of a fully discrete finite element method for the time-fractional Allen-Cahn equation

Crouzeix–Raviart and Raviart–Thomas finite-element error analysis on anisotropic meshes violating the maximum-angle condition

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

Nonlocal diffusion model with maximum principle

Anisotropic mesh adaptation for 3D anisotropic diffusion problems with application to fractured reservoir simulation

Maximum principle preserving space and time flux limiting for Diagonally Implicit Runge-Kutta discretizations of scalar convection-diffusion equations

High order maximum principle preserving discontinuous Galerkin method for convection-diffusion equations

An enhanced flux continuity three-dimensional finite element method for heterogeneous and anisotropic diffusion problems on general meshes