The Effect of Creep Flow on Two-Dimensional Isoflux Microchannels

Author(s):  
Jennifer van Rij ◽  
Todd Harman ◽  
Tim Ameel

Micro channel convective heat transfer and friction loss characteristics are numerically evaluated for gaseous, two-dimensional, steady state, laminar, constant wall heat flux flows. The effects of Knudsen number, accommodation coefficients, second order slip boundary conditions, creep flow, and thermal/hydrodynamic developing flow are considered. These effects are compared through the Poisuelle number and Nusselt number. Numerical values for the Poisuelle and Nusselt numbers are obtained using a continuum based three-dimensional, unsteady, compressible computational fluid dynamics algorithm that has been modified with slip boundary conditions. To verify the numerical results, analytic solutions for the hydrodynamically and thermally fully developed Poisuelle and Nusselt numbers have been derived. The fully developed analytic Poisuelle and Nusselt numbers are given as a function of Knudsen number, the first and second order velocity slip and temperature jump coefficients, the Brinkman number, and the ratio of the thermal creep velocity to the mean velocity. Excellent agreement between the numerical and analytical data is demonstrated. Second order slip terms and creep velocity are shown to have significant effects on the Poisuelle and Nusselt numbers.

2017 ◽  
Vol 821 ◽  
pp. 31-58 ◽  
Author(s):  
Pierre-Yves Passaggia ◽  
Alberto Scotti ◽  
Brian White

The linear instability mechanisms of horizontal convection in a rectangular cavity forced by a horizontal buoyancy gradient along its surface are investigated using local and global stability analyses for a Prandtl number equal to unity. The results show that the stability of the base flow, a steady circulation characterized by a narrow descending plume and a broad upwelling region, depends on the Rayleigh number, $Ra$. For free-slip boundary conditions at a critical value of $Ra\approx 2\times 10^{7}$, the steady base flow becomes unstable to three-dimensional perturbations, characterized by counter-rotating vortices originating within the plume region. A Wentzel–Kramers–Brillouin (WKB) method applied along closed streamlines demonstrates that this instability is of a Rayleigh–Taylor type and can be used to accurately reconstruct the global instability mode. In the case of no-slip boundary conditions, the base flow also becomes unstable to a self-sustained two-dimensional instability whose critical Rayleigh number is $Ra\approx 1.7\times 10^{8}$. Beyond this critical $Ra$, two-dimensional equilibrium stationary states of the Navier–Stokes equations are computed using the selective frequency damping method. The two-dimensional onset of instability is shown to be characterized by a family of modes also originating within the plume. A local spatio-temporal stability analysis shows that the flow becomes absolutely unstable at the origin of the plume. Taken together, these results illustrate the mechanisms behind the onset of turbulence that has been observed in horizontal convection.


2021 ◽  
Vol 65 (1) ◽  
pp. 1-23
Author(s):  
Ranis Ibragimov ◽  
◽  
Vesselin Vatchev ◽  

We examine the viscous effects of slip boundary conditions for the model describing two-dimensional Navier-Stokes flows in a plane diffuser. It is shown that the velocity profile is related to a half period shifted Weierstrass function with two parameters. This allows to approximate the explicit solution by a Taylor series expansion with two new micro- parameters, that can be measured in physical experiments. It is shown that the assumption for no-slip boundary conditions is stable in the sense that a small perturbation of the boundary values result in a small perturbation in the solutions.


2019 ◽  
Vol 29 (5) ◽  
pp. 1786-1814 ◽  
Author(s):  
A.A. Avramenko ◽  
N.P. Dmitrenko ◽  
I.V. Shevchuk ◽  
A.I. Tyrinov ◽  
V.I. Shevchuk

Purpose The paper aims to consider heat transfer in incompressible flow in a rotating flat microchannel with allowance for boundary slip conditions of the first and second order. The novelty of the paper encompasses analytical and numerical solutions of the problem, with the latter based on the lattice Boltzmann method (LBM). The analytical solution of the problem includes relations for the velocity and temperature profiles and for the Nusselt number depending on the rotation rate of the microchannel and slip velocity. It was demonstrated that the velocity profiles at high rotation rates transform from parabolic to M-shaped with a minimum at the channel axis. The temperature profiles tend to become uniform (i.e. almost constant). An increase in the channel rotation rate contributes to the increase in the Nusselt number. An increase in the Prandtl number causes a similar effect. The trend caused by the effect of the second-order slip boundary conditions depends on the closure hypothesis. It is shown that heat transfer in a flat microchannel can be successfully modeled using the LBM methodology, which takes into account the second-order boundary conditions. Design/methodology/approach The paper is based on the comparisons of an analytical solution and a numerical solution, which employs the lattice Boltzmann method. Both mathematical approaches used the first-order and second-order slip boundary conditions. The results obtained using both methods agree well with each other. Findings The analytical solution of the problem includes relations for the velocity and temperature profiles and for the Nusselt number depending on the rotation rate of the microchannel and slip velocity. It was demonstrated that the velocity profiles at high rotation rates transform from parabolic to M-shaped with a minimum at the channel axis. The temperature profiles tend to become uniform (i.e. almost constant). The increase in the channel rotation rate contributes to the increase in the Nusselt number. An increase in the Prandtl number causes the similar effect. The trend caused by the effect of the second-order slip boundary conditions depends on the closure hypothesis. It is shown that heat transfer in a flat microchannel can be successfully modeled using the LBM methodology, which considers the second-order boundary conditions. Originality/value The novelty of the paper encompasses analytical and numerical solutions of the problem, whereas the latter are based on the LBM.


Author(s):  
Jennifer van Rij ◽  
Tim Ameel ◽  
Todd Harman

Microchannel convective heat transfer characteristics in the slip flow regime are numerically evaluated for two-dimensional, steady state, laminar, constant wall heat flux and constant wall temperature flows. The effects of Knudsen number, accommodation coefficients, viscous dissipation, pressure work, second-order slip boundary conditions, axial conduction, and thermally/hydrodynamically developing flow are considered. The effects of these parameters on microchannel convective heat transfer are compared through the Nusselt number. Numerical values for the Nusselt number are obtained using a continuum based three-dimensional, unsteady, compressible computational fluid dynamics algorithm that has been modified with slip boundary conditions. Numerical results are verified using analytic solutions for thermally and hydrodynamically fully developed flows. The resulting analytical and numerical Nusselt numbers are given as a function of Knudsen number, the first- and second-order velocity slip and temperature jump coefficients, the Peclet number, and the Brinkman number. Excellent agreement between numerical and analytical data is demonstrated. Viscous dissipation, pressure work, second-order slip terms, and axial conduction are all shown to have significant effects on Nusselt numbers in the slip flow regime.


2014 ◽  
Vol 136 (5) ◽  
Author(s):  
Leila Issa ◽  
Issam Lakkis

We present reduced-order models of unsteady low-Mach-number ideal gas flows in two-dimensional rectangular microchannels subject to first-order slip-boundary conditions. The pressure and density are related by a polytropic process, allowing for isothermal or isentropic flow assumptions. The Navier–Stokes equations are simplified using low-Mach-number expansions of the pressure and velocity fields. Up to first order, this approximation results in a system that is subject to no-slip condition at the solid boundary. The second-order system satisfies the slip-boundary conditions. The resulting equations and the subsequent pressure-flow-rate relationships enable modeling the flow using analog circuit components. The accuracy of the proposed models is investigated for steady and unsteady flows in a two-dimensional channel for different values of Mach and Knudsen numbers.


Author(s):  
M. Darbandi ◽  
Y. Daghighi

A new finite-volume-based finite-element method using the quadratic elements is developed in the present study, to analyze the flow in micro and nano sizes with higher-order slip boundary conditions. The method is applied to gaseous flow in micro and nanoscale-channels. The developed method is carried out over a wide range of Knudsen numbers, which cover not only the continuum slip flow regime with 0≤Kn≤0.1 but also it entire the range of transient regime with 0.1<Kn≤10. To make the present computational model capable of simulating micro and nano sizes with the help of the Navier-Stokes equations, the modified high-order slip boundary conditions are applied which need utilizing the advantages of general quadratic second order elements in the computational domain. In other words, this paper introduces a new developed method, which is applied on higher-order elements, and employing reliable boundary conditions that all of these issues are used for the first time in the Micro/Nano study as well. The results reveal excellent agreement with those represented by analytical, DSMC, and Boltzmann calculations. The proposed method (using finite-volume-element strategy which benefits from the advantages of general quadratic second-order elements) is proved to be an efficient, practical, and accurate tool, which robustly extends the capability of our primitive large scale Navier-Stokes solver to micro and nano-scale flow predictions in slip and transient regimes. It can be regarded as a super alternative to classical molecular dynamics-based methods.


2009 ◽  
Vol 131 (10) ◽  
Author(s):  
Xiaohong Yan ◽  
Qiuwang Wang

In this paper, first, the Navier–Stokes equations for incompressible fully developed flow in microchannels and microtubes with the first-order and second-order slip boundary conditions are analytically solved. Then, the compressible Navier–Stokes equations are numerically solved with slip boundary conditions. The numerical methodology is based on the control volume scheme. Numerical results reveal that the compressibility effect increases the velocity gradient near the wall and the friction factor. On the other hand, the increment of velocity gradient near the wall leads to a much larger slip velocity than that for incompressible flow with the same value of Knudsen number and results in a corresponding decrement of friction factor. General correlations for the Poiseuille number (fRe), the Knudsen number (Kn), and the Mach number (Ma) containing the first-order and second-order slip coefficients are proposed. Correlations are validated with available experimental and numerical results.


Sign in / Sign up

Export Citation Format

Share Document