Numerical Models of Gas Flow and Heat Transfer in Microscale Channels: Capturing Rarefaction Behaviour Using a Constitutive Scaling Approach

Author(s):  
Lynne O’Hare ◽  
Jason M. Reese

In this paper we discuss the physics of rarefied gas flows at the micro-scale, including discontinuities of momentum and energy at solid boundaries, and the Knudsen layer (a region within one to two mean free paths of any solid surface where non-equilibrium flow features are dominant.) We describe how scaling the constitutive relations of the Navier-Stokes-Fourier equation set can capture key rarefaction behaviour observed in gas microsystems, and how a new implementation of this approach in fully compressible, non-isothermal CFD facilitates the analysis of “real-world” engineering problems. Details of our implementation are given, as are the results of a compressible Couette flow case study, successfully validated against available data sources. We also discuss the relative merits of two published constitutive scaling models, comparing their micro-flow predictions for half-space problems, and contrasting their individual means of application. Some practical implications of using constitutive-relation scaling are explained, and some advantages of the technique compared to alternative methods are outlined. To conclude, we examine some limitations of the method, and outline avenues of research that could potentially broaden the scope of what is a flexible and efficient approach to gas microsystem design using CFD.

Author(s):  
Angelos G. Klothakis ◽  
Georgios N. Lygidakis ◽  
Ioannis K. Nikolos

During the past decade considerable efforts have been exerted for the simulation of rarefied gas flows in a wide range of applications, like the flow over suborbital vehicles, in microelectromechanical systems, etc. Such flows appear to be significantly different from those at the continuum regime, making the Navier-Stokes equations to fail without further amendment. In this study an in-house academic CFD solver, named Galatea, is modified appropriately to account for rarefied gases. The no-slip condition on solid walls is no longer valid, hence, velocity slip and temperature jump boundary conditions are applied instead. Additionally, a second-order accurate slip model has been incorporated, namely, this of Beskok and Karniadakis, increasing the accuracy in the same area but avoiding simultaneously the numerical difficulties, entailed by the computation of the second derivative of slip velocity when complex geometries and unstructured grids are coupled. The proposed solver is validated against rarefied laminar flow over a suborbital shuttle, designed by the Azim’UTBM team. The obtained results are compared with those extracted with the parallel open-source kernel SPARTA, which is based on the DSMC method. A satisfactory agreement is reported between the two methodologies, demonstrating the potential of the modified solver to simulate effectively such flows.


Author(s):  
Duncan A. Lockerby ◽  
Jason M. Reese ◽  
Henning Struchtrup

Switching criteria for hybrid hydrodynamic/molecular gas flow solvers are developed, and are demonstrated to be more appropriate than conventional criteria for identifying thermodynamic non-equilibrium. For switching from a molecular/kinetic solver to a hydrodynamic (continuum-fluid) solver, the criterion is based on the difference between the hydrodynamic near-equilibrium fluxes (i.e. the Navier–Stokes stress and Fourier heat flux) and the actual values of stress and heat flux as computed from the molecular solver. For switching from hydrodynamics to molecular/kinetic, a similar criterion is used but the values of stress and heat flux are approximated through higher order constitutive relations; in this case, we use the R13 equations. The efficacy of our proposed switching criteria is tested within an illustrative hybrid kinetic/Navier–Stokes solver. For the test cases investigated, the results from the hybrid procedure compare very well with the full kinetic solution, and are obtained at a fraction of the computational cost.


2013 ◽  
Vol 13 (5) ◽  
pp. 1330-1356 ◽  
Author(s):  
G. H. Tang ◽  
G. X. Zhai ◽  
W. Q. Tao ◽  
X. J. Gu ◽  
D. R. Emerson

AbstractGases in microfluidic structures or devices are often in a non-equilibrium state. The conventional thermodynamic models for fluids and heat transfer break down and the Navier-Stokes-Fourier equations are no longer accurate or valid. In this paper, the extended thermodynamic approach is employed to study the rarefied gas flow in microstructures, including the heat transfer between a parallel channel andpressure-driven Poiseuille flows through a parallel microchannel andcircular microtube. The gas flow characteristics are studied and it is shown that the heat transfer in the non-equilibrium state no longer obeys the Fourier gradient transport law. In addition, the bimodal distribution of streamwise and spanwise velocity and temperature through a long circular microtube is captured for the first time.


2020 ◽  
Vol 34 (14n16) ◽  
pp. 2040076
Author(s):  
Zhen-Yu Yuan ◽  
Zhong-Zheng Jiang ◽  
Wen-Wen Zhao ◽  
Wei-Fang Chen

This paper is focused on the gas properties over a cylinder from continuum to rarefied regimes based on the non-Newtonian constitutive model. This new constitutive model is first derived from Eu’s nonequilibrium ensemble method, which is intended for accurate description of nonequilibrium flows. Some assumptions and simplifications are made during the establishing progress of the new constitutive model by both Eu and Myong. To verify its accuracy, temperature contours and skin frictions around the cylinder are simulated by this new model. The inflow Mach number is equal to 10 and the Knudsen number ranges from 0.002 to 0.05. All simulation results are compared with Navier–Stokes (NS) and the direct simulation Monte Carlo (DSMC) methods in detail. The comparisons of friction around the surface show that the non-Newtonian constitutive models are better than the linear constitutive relations of NS equations for the prediction of nonequilibrium flow and much more close to DSMC simulation results.


2015 ◽  
Author(s):  
K. Farber ◽  
P. Farber ◽  
J. Gräbel ◽  
S. Krick ◽  
J. Reitz ◽  
...  

Author(s):  
Deepak Nabapure ◽  
Ram Chandra Murthy

Abstract The present study investigates the flow behavior of the rarefied gas over a wall-mounted cube. The problem is studied for different cube heights (h) of 9mm and 18mm in the slip and transition regimes. The Direct Simulation Monte Carlo (DSMC) method is employed to evaluate the properties such as velocity, pressure and temperature fields. The Reynolds number (Re) ranges from 403 to 807, and the Knudsen number (Kn) is in the range from 0.05 to 0.103. A typical shock wave is formed in front of the cube. The recirculation length of the vortices normalized with respect to the respective cube heights for Kn = 0.05 and Kn = 0.103 are about 1.11 and 1.95 respectively. Similarly, the center of the vortices is located at about 3.33 and 6.11 times the respective cube heights upstream, for Kn = 0.05 and Kn = 0.103. The local temperature and pressure variations observed upstream of the cube are two orders higher in magnitude and are primarily attributed to strong compressibility effects. The present study paves the way for benchmarking, and forms a basis for understanding the rarefied gas flows over complex geometries.


Author(s):  
M. Hossein Gorji ◽  
Stephan Küchlin ◽  
Patrick Jenny

In this work, we present a hybrid algorithm based on the Fokker-Planck (FP) kinetic model and direct simulation Monte Carlo (DSMC) for studies of rarefied gas flows. A particle based FP solution algorithm for rarefied gas flow simulations has recently been devised by the authors. The motivation behind the FP approximation is purely computational, i.e. due to the fact that the resulting random processes are continuous in time the computational cost of the corresponding time integration becomes independent of the Knudsen number. However, the method faces limitations for flows with very high Knudsen numbers (larger than approximately 5). In the method presented here, the FP approach is coupled with DSMC in order to gain from the efficiency of the FP model and from the accuracy of DSMC at small and large cell based Knudsen numbers, respectively.


Author(s):  
Mehti Koklu ◽  
Nurhak Erbas ◽  
Oktay Baysal

Effectiveness of two-dimensional synthetic jet is studied using numerical simulations. A Navier-Stokes (NS) solver for moving and deforming meshes has been modified to investigate numerically the diaphragm-driven flow in and out of two synthetic jet cavity geometries. Compressible flow simulations are required for rarefied gas flows to accurately predict the micro flow field. The solver is modified to accommodate slip wall boundary condition proposed in literature for micro scale flow problems. The piezoelectric-driven diaphragm of the cavity is modeled in a realistic manner as a moving boundary to accurately compute the flow inside the jet cavity. The primary focus of the proposed paper will be on the analysis of the design space determined by the geometric and flow-type design variables that identify the effectiveness of the synthetic jet by means of the orifice jet velocity and local jet momentum rate. The design variables are the membrane oscillation frequency (f), membrane oscillation amplitude (A), orifice width (d), and membrane width (W). The present computations for jet discharging into quiescent medium reveal that these variables have determining effects on the flow control parameters, which are the jet exit velocity, local momentum rate, as well as vortex shedding from the orifice.


Sign in / Sign up

Export Citation Format

Share Document