It is well established that non-equilibrium flows cannot properly be described by traditional hydrodynamics, namely, the Navier-Stokes-Fourier (NSF) equations. Such flows occur, for example, in micro-electro-mechanical systems (MEMS), and ultra vacuum systems, where the dimensions of the devices are comparable to the mean free path of a gas molecule. Therefore, the study of non-equilibrium effects in gas flows is extremely important. The general interest of the present study is to explore boundary value problems for moderately rarefied gas flows, with an emphasis on numerical solutions of the regularized 13--moment equations (R13). Boundary conditions for the moment equations are derived based on either phenomenological principles or on microscopic gas-surface scattering models, e.g., Maxwell's accommodation model and the isotropic scattering model. Using asymptotic analysis, several non-linear terms in the R13 equations are transformed into algebraic terms. The reduced equations allow us to obtain numerical solutions for multidimensional boundary value problems, with the same set of boundary conditions for the linearized and fully non-linear equations. Some basic flow configurations are employed to investigate steady and unsteady rarefaction effects in rarefied gas flows, namely, planar and cylindrical Couette flow, stationary heat transfer between two plates, unsteady and oscillatory Couette flow. A comparison with the corresponding results obtained previously by the DSMC method is performed. The influence of rarefaction effects in the lid driven cavity problem is investigated. Solutions obtained from several macroscopic models, in particular the classical NSF equations with jump and slip boundary conditions, and the R13--moment equations are compared. The R13 results compare well with those obtained from more costly solvers for rarefied gas dynamics, such as the Direct Simulation Monte Carlo (DSMC) method. Flow and heat transfer in a bottom heated square cavity in a moderately rarefied gas are investigated using the R13 and NSF equations. The results obtained are compared with those from the DSMC method with emphasis on understanding thermal flow characteristics from the slip flow to the early transition regime. The R13 theory gives satisfying results including flow patterns in fair agreement with DSMC in the transition regime, which the conventional Navier-Stokes-Fourier equations are not able to capture. / Graduate / 0548 / anirudh@uvic.ca

The fast-paced growth in microelectromechanical systems (MEMS), microfluidic fabrication, porous media applications, biomedical assemblies, space propulsion, and vacuum technology demands accurate and practical transport equations for rarefied gas flows. It is well-known that in rarefied situations, due to strong deviations from the continuum regime, traditional fluid models such as Navier-Stokes-Fourier (NSF) fail. The shortcoming of continuum models is rooted in nonequilibrium behavior of gas particles in miniaturized and/or low-pressure devices, where the Knudsen number (Kn) is sufficiently large. Since kinetic solutions are computationally very expensive, there has been a great desire to develop macroscopic transport equations for dilute gas flows, and as a result, several sets of extended equations are proposed for gas flow in nonequilibrium states. However, applications of many of these extended equations are limited due to their instabilities and/or the absence of suitable boundary conditions. In this work, we concentrate on regularized 13-moment (R13) equations, which are a set of macroscopic transport equations for flows in the transition regime, i.e., Kn≤1. The R13 system provides a stable set of equations in Super-Burnett order, with a great potential to be a powerful CFD tool for rarefied flow simulations at moderate Knudsen numbers. The goal of this research is to implement the R13 equations for problems of practical interest in arbitrary geometries. This is done by transformation of the R13 equations and boundary conditions into general curvilinear coordinate systems. Next steps include adaptation of the transformed equations in order to solve some of the popular test cases, i.e., shear-driven, force-driven, and temperature-driven flows in both planar and curved flow passages. It is shown that inexpensive analytical solutions of the R13 equations for the considered problems are comparable to expensive numerical solutions of the Boltzmann equation. The new results present a wide range of linear and nonlinear rarefaction effects which alter the classical flow patterns both in the bulk and near boundary regions. Among these, multiple Knudsen boundary layers (mechanocaloric heat flows) and their influence on mass and energy transfer must be highlighted. Furthermore, the phenomenon of temperature dip and Knudsen paradox in Poiseuille flow; Onsager's reciprocity relation, two-way flow pattern, and thermomolecular pressure difference in simultaneous Poiseuille and transpiration flows are described theoretically. Through comparisons it is shown that for Knudsen numbers up to 0.5 the compact R13 solutions exhibit a good agreement with expensive solutions of the Boltzmann equation.