A wall-function approach to incorporating Knudsen-layer effects in gas micro flow simulations

For gas flows in microfluidic configurations, the Knudsen layer close to the wall can comprise a substantial part of the entire flow field and has a major effect on quantities such as the mass flow rate through micro devices. The Knudsen layer itself is characterized by a highly nonlinear relationship between the viscous stress and the strain rate of the gas, so even if the Navier-Stokes equations can be used to describe the core gas flow they are certainly inappropriate for the Knudsen layer itself. In this paper we propose a "wall-function" model for the stress/strain rate relations in the Knudsen layer. The constitutive structure of the Knudsen layer has been derived from results from kinetic theory for isothermal shear flow over a planar surface. We investigate the ability of this simplified model to predict Knudsen-layer effects in a variety of configurations. We further propose a semi-empirical Knudsen-number correction to this wall function, based on high-accuracy DSMC results, to extend the predictive capabilities of the model to greater degrees of rarefaction

Solving the Boltzmann equation deterministically by the fast spectral method : application to gas microflows

Based on the fast spectral approximation to the Boltzmann collision operator, we present an accurate and efficient deterministic numerical method for solving the Boltzmann equation. First, the linearised Boltzmann equation is solved for Poiseuille and thermal creep flows, where the influence of different molecular models on the mass and heat flow rates is assessed, and the Onsager-Casimir relation at the microscopic level for large Knudsen numbers is demonstrated. Recent experimental measurements of mass flow rates along a rectangular tube with large aspect ratio are compared with numerical results for the linearised Boltzmann equation. Then, a number of two-dimensional micro flows in the transition and free molecular flow regimes are simulated using the nonlinear Boltzmann equation. The influence of the molecular model is discussed, as well as the applicability of the linearised Boltzmann equation. For thermally driven flows in the free molecular regime, it is found that the magnitudes of the flow velocity are inversely proportional to the Knudsen number. The streamline patterns of thermal creep flow inside a closed rectangular channel are analysed in detail: when the Knudsen number is smaller than a critical value, the flow pattern can be predicted based on a linear superposition of the velocity profiles of linearised Poiseuille and thermal creep flows between parallel plates. For large Knudsen numbers, the flow pattern can be determined using the linearised Poiseuille and thermal creep velocity profiles at the critical Knudsen number. The critical Knudsen number is found to be related to the aspect ratio of the rectangular channel

Modelling thermal flow in a transition regime using a lattice Boltzmann approach

Lattice Boltzmann models are already able to capture important rarefied flow phenomena, such as velocity-slip and temperature jump, provided the effects of the Knudsen layer are minimal. However, both conventional hydrodynamics, as exemplified by the Navier-Stokes-Fourier equations, and the lattice Boltzmann method fail to predict the nonlinear velocity and temperature variations in the Knudsen layer that have been observed in kinetic theory. In the present paper, we propose an extension to the lattice Boltzmann method that will enable the simulation of thermal flows in the transition regime where Knudsen layer effects are significant. A correction function is introduced that accounts for the reduction in the mean free path near a wall. This new approach is compared with direct simulation Monte Carlo data for Fourier flow and good qualitative agreement is obtained for Knudsen numbers up to 1.58

Knudsen Diffusion in Silicon Nanochannels

Measurements on helium and argon gas flow through an array of parallel, linear channels of 12 nm diameter and 200 micrometer length in a single crystalline silicon membrane reveal a Knudsen diffusion type transport from 10^2 to 10^7 in Knudsen number Kn. The classic scaling prediction for the transport diffusion coefficient on temperature and mass of diffusing species,D_He ~ sqrt(T), is confirmed over a T range from 40 K to 300 K for He and for the ratio of D_He/D_Ar ~ sqrt(m_Ar/m_He). Deviations of the channels from a cylindrical form, resolved with transmission electron microscopy down to subnanometer scales, quantitatively account for a reduced diffusivity as compared to Knudsen diffusion in ideal tubular channels. The membrane permeation experiments are described over 10 orders of magnitude in Kn, encompassing the transition flow regime, by the unified flow model of Beskok and Karniadakis.Comment: 4 pages, 3 figure

Asymptotic Theory of Rayleigh Problem in Rarefied Gas

The asymptotic theory of Rayleigh shear flow for large values of time is developed on the basis of the linearized Boltzmann-Krook equation. Asymptotic equations for mean velocity outside the Knudsen layer are obtained by employing the Hilbert expansion. Slip boundary conditions are derived from the analysis of the Knudsen layer adjacent to the wall. A solution of the asymptotic equation is obtained under the slip boundary condition and zero initial condition. Discussions are also made of the flow induced by a slowly oscillating flat plate

Knudsen number, ideal hydrodynamic limit for elliptic flow and QGP viscosity in s\sqrt{s}=62 and 200 GeV Cu+Cu/Au+Au collisions

Taking into account of entropy generation during evolution of a viscous fluid, we have estimated inverse Knudsen number, ideal hydrodynamic limit for elliptic flow and QGP viscosity to entropy ratio in $\sqrt{s}$=62 and 200 GeV Cu+Cu/Au+Au collisions. Viscosity to entropy ratio is estimated as $\eta/s=0.17\pm 0.10\pm 0.20$, the first error is statistical, the second one is systematic. In a central Au+Au collision, inverse Knudsen number is $\approx 2.80\pm 1.63$, which presumably small for complete equilibration. In peripheral collisions it is even less. Ideal hydrodynamic limit for elliptic flow is $\sim$40% more than the experimental flow in a central collision.Comment: 4 pages, 2 figures, 2 tables. Final version to be published in Phys. Rev.