On the modelling of isothermal gas flows at the microscale

This paper makes two new propositions regarding the modelling of rarefied (non-equilibrium) isothermal gas flows at the microscale. The first is a new test case for benchmarking high-order, or extended, hydrodynamic models for these flows. This standing time-varying shear-wave problem does not require boundary conditions to be specified at a solid surface, so is useful for assessing whether fluid models can capture rarefaction effects in the bulk flow. We assess a number of different proposed extended hydrodynamic models, and we find the R13 equations perform the best in this case. Our second proposition is a simple technique for introducing non-equilibrium effects caused by the presence of solid surfaces into the computational fluid dynamics framework. By combining a new model for slip boundary conditions with a near-wall scaling of the Navier--Stokes constitutive relations, we obtain a model that is much more accurate at higher Knudsen numbers than the conventional second-order slip model. We show that this provides good results for combined Couette/Poiseuille flow, and that the model can predict the stress/strain-rate inversion that is evident from molecular simulations. The model's generality to non-planar geometries is demonstrated by examining low-speed flow around a micro-sphere. It shows a marked improvement over conventional predictions of the drag on the sphere, although there are some questions regarding its stability at the highest Knudsen numbers

The usefulness of higher-order constitutive relations for describing the Knudsen layer

The Knudsen layer is an important rarefaction phenomenon in gas flows in and around microdevices. Its accurate and efficient modeling is of critical importance in the design of such systems and in predicting their performance. In this paper we investigate the potential that higher-order continuum equations may have to model the Knudsen layer, and compare their predictions to high-accuracy DSMC (direct simulation Monte Carlo) data, as well as a standard result from kinetic theory. We find that, for a benchmark case, the most common higher-order continuum equation sets (Grad's 13 moment, Burnett, and super-Burnett equations) cannot capture the Knudsen layer. Variants of these equation families have, however, been proposed and some of them can qualitatively describe the Knudsen layer structure. To make quantitative comparisons, we obtain additional boundary conditions (needed for unique solutions to the higher-order equations) from kinetic theory. However, we find the quantitative agreement with kinetic theory and DSMC data is only slight

Capturing the Knudsen layer in continuum-fluid models of non-equilibrium gas flows

In hypersonic aerodynamics and microflow device design, the momentum and energy fluxes to solid surfaces are often of critical importance. However, these depend on the characteristics of the Knudsen layer - the region of local non-equilibrium existing up to one or two molecular mean free paths from the wall in any gas flow near a surface. While the Knudsen layer has been investigated extensively using kinetic theory, the ability to capture it within a continuum-fluid formulation (in conjunction with slip boundary conditions) suitable for current computational fluid dynamics toolboxes would offer distinct and practical computational advantages